[Δομημένος Πρ.] Εργασία E

[Thmmy.gr portal]

  Thmmy.gr portal
   Forum
   Downloads
   Ενεργ. Λογαριασμού
   Επικοινωνία
  
  Χρήσιμα links
   Σελίδα τμήματος
   Βιβλιοθήκη Τμήματος
   Elearning
   Φοιτητικά fora
   Πρόγραμμα Λέσχης
   Πρακτική Άσκηση
   Ηλεκτρονική Εξυπηρέτηση Φοιτητών
   Διανομή Συγγραμμάτων
   Ψηφιακό Καταθετήριο Διπλωματικών
   Πληροφορίες Καθηγητών
   Instagram @thmmy.gr
   mTHMMY
  
  Φοιτητικές Ομάδες
   ACM
   Aristurtle
   ART
   ASAT
   BEAM
   BEST Thessaloniki
   EESTEC LC Thessaloniki
   EΜΒ Auth
   IAESTE Thessaloniki
   IEEE φοιτητικό παράρτημα ΑΠΘ
   SpaceDot
   VROOM
   Panther
  

THMMY.gr Web

Helios_MultiScribbles2SMF Default ThemeSMFone_Blue

[Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads]

Είσαι πρωτοετής;... Καλώς ήρθες! Μπορείς να βρεις πληροφορίες εδώ. Βοήθεια για τους καινούργιους μέσω χάρτη.

Κατεβάστε εδώ το Android Application για εύκολη πρόσβαση στο forum.

Ανεβάζετε τα θέματα των εξετάσεων στον τομέα Downloads με προσοχή στα ονόματα των αρχείων!

Για ανανέωση (ή προσθήκη νέου) avatar, πρέπει η μεγαλύτερη διάσταση της εικόνας να είναι 110 pixels.

[alexpaxi]

paidia poia einai h sunarthsh pou theloume na broume tis rizes ths?emeis epilegoume poia tha xrhsimopoihsoume?

[Exomag]

paidia poia einai h sunarthsh pou theloume na broume tis rizes ths?emeis epilegoume poia tha xrhsimopoihsoume?

Ναι, όποια θες εσύ...

ΥΣ: Να βλέπετε και τα προηγούμενα posts

[alexpaxi]

Τα ειδα αλλα δεν εβγαλα ακρη..thanks anyway..

[jthois]

Παιδιά, εκτός από την f(x), θα δώσουμε και την f '(x);
Γιατί αλλιώς θα πρέπει να την υπολογίσουμε.

[thanospr]

Γινεται να το κανουμε με define η οχι?

[Exomag]

Γινεται να το κανουμε με define η οχι?

Αυτό που θα πρέπει να κάνεις είναι να ορίσεις δύο (μια για την f(x) και μια για την παράγωγο της) δικιές σου συναρτήσεις που να επιστρέφουν τιμές...
πχ

Code:

float f(float x) {
 return x*x;
}

float df_dx(float x) {
 return 2*x;
}

[manos]

Μήπως η αναδρομική σχέση για την froot θα έπρεπε να είναι xi+1 = xi-(xi-xi-1)f(xi)/(f(xi)-f(xi-1)) και όχι               
xi+1 = xi-(xi-xi-1)f(xi)/(f(xi-f(xi-1))?

[thanospr]

Αμα το κανω μεσα στις nroot και froot πειραζει?

[c0ndemn3d]

Αμα το κανω μεσα στις nroot και froot πειραζει?

Δεν νομίζω να υπάρχει πρόβλημα... Αλλά ο τρόπος που λέει ο exomag είναι ο προτεινόμενος από Κορτέσιο και είναι και γενικά καλύτερος από άποψη οργάνωσης του κώδικα. Αλλά όπως και να το κάνεις πιστεύω δεν θα ενοχλήσει κανέναν.

[Exomag]

Αμα το κανω μεσα στις nroot και froot πειραζει?

Σκέψου μετά, για παράδειγμα, να θες να αλλάξεις την συνάρτηση f(x). Θα πρέπει να ψάξεις όλο τον κώδικα για να βρεις που (μέσα στις nroot και froot) έχει γράψει τη σχέση της f(x). Ενώ με τον τρόπο που εξήγησα (δήλωση ως ξεχωριστή συνάρτηση) είναι πιο εύκολο/γρήγορο...

Μήπως η αναδρομική σχέση για την froot θα έπρεπε να είναι xi+1 = xi-(xi-xi-1)f(xi)/(f(xi)-f(xi-1)) και όχι               
xi+1 = xi-(xi-xi-1)f(xi)/(f(xi-f(xi-1))?

Ναι, είναι τυπογραφικό...

[thanospr]

Μα εγω ειπα να το κανω με define για να μπορω να αλλαξω μονο τους συντελεστες.Για παραδειγμα 5χ να μπορω να το κανω 3χ.Εσεις αυτο το λετε για να μπορω να αλλαξω την συναρτηση και να βαλω οποια θελω μετα?

[Exomag]

Μα εγω ειπα να το κανω με define για να μπορω να αλλαξω μονο τους συντελεστες.Για παραδειγμα 5χ να μπορω να το κανω 3χ.Εσεις αυτο το λετε για να μπορω να αλλαξω την συναρτηση και να βαλω οποια θελω μετα?

Ναι, εγώ αυτό που λέω είναι για να αλλάξεις εύκολα όλη την συνάρτηση. Βεβαίως και μπορείς να κάνεις define τους συντελεστές που θα χρησιμοποιήσεις στη συνάρτηση για να μπορείς να τους αλλάζεις ευκολότερα...

[thanospr]

Οκ.Ευχαριστω.

[c0ndemn3d]

Παιδιά, καλό θα ήταν να δείτε και αυτό(ο κώδικας παρακάτω): http://alexander.ee.auth.gr:8083/eTHMMY/archive/102/downloadFile/4471/pr22.cpp

Code:

/* Επίλυση της εξίσωσης f(x)=0 με τη μέθοδο Regula-falsi */

#include <stdio.h>
#include <process.h>

double f(double x)
{
  return 2*x*x-18;
}

double root(double (*p)(double),double x1,double x0,double e,long n)
{
  long i;
  double x;
  for(i=0; i<n; i++){
    if(p(x1)==p(x0)){
      printf("f(x at step i) = f(x at step i-1)\n");
      exit(1);
    }
    x = x1-((x1-x0)*p(x1)/(p(x1)-p(x0)));
    if((x1<x?x-x1:x1-x)<e) break;
    x0=x1;
    x1=x;
  }
  if(i>=n)
    printf("Δε συγκλίνει μετά από %ld επαναλήψεις\n",n);
  return x;
}

void main()
{ 
  double x1,x0,r,e=0.00000000001;
  long n=10000;
  scanf("%lf %lf",&x0,&x1);
  r=root(f,x1,x0,e,n);
  printf("root = %lf\n",r);
}

Είναι από τα παραδείγματα του Κορτέσιου στο ethmmy.

[Exomag]

Παιδιά, καλό θα ήταν να δείτε και αυτό(ο κώδικας παρακάτω): http://alexander.ee.auth.gr:8083/eTHMMY/archive/102/downloadFile/4471/pr22.cpp

Code:

/* Επίλυση της εξίσωσης f(x)=0 με τη μέθοδο Regula-falsi */

#include <stdio.h>
#include <process.h>

double f(double x)
{
  return 2*x*x-18;
}

double root(double (*p)(double),double x1,double x0,double e,long n)
{
  long i;
  double x;
  for(i=0; i<n; i++){
    if(p(x1)==p(x0)){
      printf("f(x at step i) = f(x at step i-1)\n");
      exit(1);
    }
    x = x1-((x1-x0)*p(x1)/(p(x1)-p(x0)));
    if((x1<x?x-x1:x1-x)<e) break;
    x0=x1;
    x1=x;
  }
  if(i>=n)
    printf("Δε συγκλίνει μετά από %ld επαναλήψεις\n",n);
  return x;
}

void main()
{ 
  double x1,x0,r,e=0.00000000001;
  long n=10000;
  scanf("%lf %lf",&x0,&x1);
  r=root(f,x1,x0,e,n);
  printf("root = %lf\n",r);
}

Είναι από τα παραδείγματα του Κορτέσιου στο ethmmy.

Απλά εδώ υλοποιεί την διαδικασία χωρίς αναδρομή, ενώ η εκφώνηση της εργασίας Ε ζητάει ρητά αναδρομή...