WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Newtown Verfahren und Iterationen plotten

 

MaxChuckNorris

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 18.05.2022, 08:35     Titel: Newtown Verfahren und Iterationen plotten
  Antworten mit Zitat      
Hi, ich habe die Aufgabenstellung:

"Schreiben Sie ein Programm, dass das Newtown verfahren bis auf einen Fehler
Epsilon > 0 berechnet und die Iterationen plottet.

Ermitteln Sie dann den Plot für die Approximation der Lösung von e^sin(x/2)=x

im Intervall von [0,5] mit dem Startwert x_0=0 und den Approximationsfehler Epsilon=10^-8. "

Ein Kommilitone hat diesen Code bisher gefunden/erstellt:

Code:
clear all;
 
prompt={'Funktion','Ableitung','Startwert','Iteration','Min','Max'}
defaultanswer={'exp(sin(x/2))-x; ','(cos(x/2)*exp(sin(x/2)))/2 - 1;','0','50','0','0.5' };
name='Newton';
numlines=1;


Abfrage=inputdlg(prompt,name,numlines,defaultanswer)



x(1)=[str2num(Abfrage{3})]
n=[str2num(Abfrage{4})]
minx=[str2num(Abfrage{5})]
maxx=[str2num(Abfrage{6})]


fx=eval(['@(x) ' (Abfrage{1})])
f1x=eval(['@(x) ' (Abfrage{2})])


eps=1e-8;

i=1;
while abs(fx(x(i)))>eps
 x(i+1)=x(i)-(fx(x(i))/f1x(x(i)))
 f(i+1)=fx(x(i))
 if i==n
     break
 end
 i=i+1
end

y=fx(minx:0.1:maxx)
plot(x,f,'b*-')
hold on
plot(minx:0.1:maxx,y,'r-')
hold off
 


Wir verstehen leider nicht, was damit gemeint ist, dass wir die Iterationen plotten sollen und leider sieht das, was wir geplottet haben, sehr merkwürdig aus, sind uns auch nicht sicher, ob wir die Aufgabenstellung richtig verstanden haben.

Also haben auch nicht verstanden, was mit Itertaionen plotten gemeint ist.


MaxChuckNorris

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 18.05.2022, 08:38     Titel:
  Antworten mit Zitat      
Im Code oben hat sich ein Fehler eingeschlichen
defaultanswer={'exp(sin(x/2))-x; ','(cos(x/2)*exp(sin(x/2)))/2 - 1;','0','50','0','0.5' };

hier sollte natürlich am Ende 5 stehen, also:
defaultanswer={'exp(sin(x/2))-x; ','(cos(x/2)*exp(sin(x/2)))/2 - 1;','0','50','0','5' }
 
MaxChuckNorris

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 18.05.2022, 08:51     Titel:
  Antworten mit Zitat      
Ein anderes Problem ist zudem, ich habe 1e^-8, bekomme jedocch immer nur 5 Nachkommastellen raus?
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 24.448
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 18.05.2022, 10:24     Titel:
  Antworten mit Zitat      
Hallo,

der Code plottet doch schon die Iterationen (= Schritte des Algorithmus) in blau.

Die zweite Zeile in der while-Schleife sollte m.E. eher sein:
Code:
f(i+1)=fx(x(i+1)); % x(i+1) statt x(i)

vor der while-Schleife
Code:
f(1) = fx(x(1));

Dann sieht man, wie sich das Verfahren an der Funktion "entlang tastet".

Im Editor werden mehrere Code Analyzer-Warnungen in orange angezeigt. Diese sollte man beheben. Zudem ist eval nicht empfehlenswert. Stattdessen könnte man hier str2func verwenden.

Die 5 Nachkommastellen sind eine Frage der Anzeige und haben mit dem gesetzten eps nichts zu tun. Die Anzeige ändern kannst du mit format .

Grüße,
Harald
_________________

1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
Private Nachricht senden Benutzer-Profile anzeigen
 
MaxChuckNorris

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 18.05.2022, 17:43     Titel:
  Antworten mit Zitat      
Danke, wenn ich jedoch die Verbesserungen bezüglich der Whileschleife vollziehe, startet das Programm nicht mehr? Und ich wollte eigentlich, dass man so lange iteriert, bis der x wert 10^-8 entspricht. Hier geht das aber durch den y-wert oder?

Habe keine Ahnung von matlab, es ist eigentlich so, dass ich die Aufgaben mit einem Kommilitonen mache, ich übernehme Java und C Aufgaben, er die Matlab Aufgaben (Er ist nur leider krank zurzeit). Verstehe deshalb bei dem Rest nicht genau, wie ich das implementiere
 
MaxChuckNorris

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 18.05.2022, 17:46     Titel:
  Antworten mit Zitat      
Wobei sorry, habe es nun hinbekommen, also die Ersetzungen der while-schleife, aber ist es hier so, dass der y-wert gegen 10^-8 iteriert wird oder der x wert?
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 24.448
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 18.05.2022, 17:54     Titel:
  Antworten mit Zitat      
Hallo,

Zitat:
wenn ich jedoch die Verbesserungen bezüglich der Whileschleife vollziehe, startet das Programm nicht mehr?

Bei mir schon. Mit deinen Verbesserungen ist also anscheinend etwas schief gegangen. Was, kann ich dir nur sagen, wenn du deinen aktuellen Code teilst.

Zitat:
Und ich wollte eigentlich, dass man so lange iteriert, bis der x wert 10^-8 entspricht. Hier geht das aber durch den y-wert oder?

Ja, und so läuft nunmal das Newton-Verfahren. Alternativ kannst du iterieren, bis die Änderung des x-Werts < 10^-8 ist. Die Umsetzung ist dann in MATLAB auch nicht viel anders als in C oder Java.

Grüße,
Harald
_________________

1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
Private Nachricht senden Benutzer-Profile anzeigen
 
MaxChuckNorris

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 18.05.2022, 21:38     Titel:
  Antworten mit Zitat      
Danke, konnte das auch ergänzen deine Verbesserung, nun sieht es jedoch so aus, als hätte man 3 Nullstellen und nich tnur eine?
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 24.448
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 18.05.2022, 22:41     Titel:
  Antworten mit Zitat      
Hallo,

[quote]nun sieht es jedoch so aus, als hätte man 3 Nullstellen und nich tnur eine?/quote]
Für mich nicht. Die rote Kurve kreuzt die Null-Linie einmal, und die Iterationen in blau nähern sich der einen Nullstelle an.
Wo siehst du 3 Nullstellen?

Grüße,
Harald
_________________

1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
Private Nachricht senden Benutzer-Profile anzeigen
 
chuckynorisi

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 19.05.2022, 02:50     Titel:
  Antworten mit Zitat      
Danke, ja passt doch! Aber wie funktioniert das mit dem x-Wert, damit ich 10^-6 Stellen anzeige, habs leider nicht hinbekommen
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 24.448
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 19.05.2022, 09:25     Titel:
  Antworten mit Zitat      
Hallo,

zum einen kannst du die Anzeige wie gesagt über format anpassen.
Wenn dir die Anpassungsmöglichkeiten von format nicht reichen, kannst du fprintf verwenden, z.B.
Code:


Grüße,
Harald
_________________

1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
Private Nachricht senden Benutzer-Profile anzeigen
 
MacChuckNoris

Gast


Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 19.05.2022, 10:41     Titel:
  Antworten mit Zitat      
Danke dir, wir sollen nicht den Code abgeben, nur die Ausgabe. Ich fände es jedoch etwas viel, wenn ich alles, was rauskommt, rauskopiere. Hättwst Du vielleicht eine Empfehlung? Was würdest Du alles von der Ausagbe abgeben, bzw. wie viel und was würdest Du nicht abgeben? Wir sollen die Ergebnise abgeben, also die Iterationen, aber wäre das da nicht zu viel?
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 24.448
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2017b
     Beitrag Verfasst am: 19.05.2022, 10:56     Titel:
  Antworten mit Zitat      
Hallo,

die Iterationen bestehen an sich nur aus den blauen Datenpunkten, also würde ich x und f kopieren:
Code:
fprintf("%.6f %.6f \n", [x; f])


Grüße,
Harald
_________________

1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
Private Nachricht senden Benutzer-Profile anzeigen
 
Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

Du kannst Beiträge in dieses Forum schreiben.
Du kannst auf Beiträge in diesem Forum antworten.
Du kannst deine Beiträge in diesem Forum nicht bearbeiten.
Du kannst deine Beiträge in diesem Forum nicht löschen.
Du kannst an Umfragen in diesem Forum nicht mitmachen.
Du kannst Dateien in diesem Forum posten
Du kannst Dateien in diesem Forum herunterladen
.





 Impressum  | Nutzungsbedingungen  | Datenschutz | FAQ | goMatlab RSS Button RSS

Hosted by:


Copyright © 2007 - 2024 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.