GUI Funktion aus Edit plotten und ableiten, Fehler irgendwo!

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Lelin

Forum-Newbie


	Beiträge: 1

	Anmeldedatum: 19.08.15

	Wohnort: ---

	Version: 2015

Verfasst am: 22.08.2015, 16:18 Titel: GUI Funktion aus Edit plotten und ableiten, Fehler irgendwo!

Hallo, ich bin ein Matlab-Neuling und habe die Aufgabe das Newtonverfahren mit GUI zu programmieren. Dabei wird der Graph geplottet und die Tangenten mittels Schleife bis zu einer gewissen Genauigkeit eingezeichnet. Es ist auch soweit fast fertig.

Mein Problem ist, dass alles nur funktioniert, wenn die Funktion keinen Summanden.
Also bei x^2 ist alles schön und bei x^2-1 funktioniert nix mehr.
Desweitern kann ich keine Ableitungen bilden mit diesem f1=diff(f);
Habe dies zwar umgangen, aber mit wäre mir lieber.

Meine Vermutung ist, dass ich etwas falsch gemacht habe, da wo die Funktion als solche definiert wird. Aber ohne Erfahrung hab ich da keine Ahnung was das sein könnte.

Bitte schaut es euch mal an

Code:

function buttonplot_Callback(hObject, eventdata, handles)

% Alte Werte löschen
cla(handles.axes1);

%Bereich festlegen
xmin=str2double(get(handles.exmin,'string'));
xmax=str2double(get(handles.exmax,'string'));

ymin=str2num(get(handles.eymin,'string'));
ymax=str2num(get(handles.eymax,'string'));
x=xmin:0.1:xmax;
set(handles.axes1, 'XLimMode', 'manual');
set(handles.axes1, 'YLimMode', 'manual');
Xlim=[xmin xmax];
Ylim=[ymin ymax];
set(gca,'XLim',Xlim, 'YLim',Ylim);

%Funktion holen
[b]get(handles.efunc,'String');
f = eval(strcat('@(x) ',get(handles.efunc,'String'))); [/b]

%Startwert x aus der Edit holen
xT=str2double(get(handles.estart,'String'));
a=str2double(get(handles.eab,'String'));

%Diagrammeinstellungen
cla;
grid on;
axis equal;
hold on;

xlabel('x')
ylabel('y')

if get(handles.checkbox2,'value') == 0.0
fplot(f, Xlim);
end

x0=0;
yT = f(xT); % Punkt
zaehler=0;
NST=fsolve(f,x0);
abstand=abs(NST-xT);

tableData = cell(3);

%--------------------------------------------------------------------
if yT~=0 % prüfen ob Startpunkt schon Nullstelle
%while zaehler<15 % Abbruch damit nur endlich viele Durchläufe
while abstand>a || zaehler<15; %Betrag(f(x)=0-tangente(x)=0)

% Tangente
yT=f(xT);
x1 = xT - 1e-10;
x2 = xT + 1e-10;
y1 = f(x1);
y2 = f(x2);
m = (y2 - y1) / (x2 - x1);

n = yT - m * xT;
% Tangente als Funktion definieren
t = @(x) m * x + n;
fplot(t, Xlim, 'color','c')

xT=fsolve(t, 0); %neues xT ist Nullstelle von Tangente
abstand=abs(NST-xT);

tableData(zaehler + 1, 1) = num2cell(zaehler + 1);
tableData(zaehler + 1, 2) = num2cell(xT);
tableData(zaehler + 1, 3) = cellstr(strcat(num2str(m),'*x + ',num2str(n)));

zaehler=zaehler+1;

end

set(handles.table1,'data',tableData);

set(gca,'XLim',Xlim, 'YLim',Ylim);

else
errordlg('Nullstelle ist Startwert','File error')
end

%______________________________________________________________________

Funktion ohne Link?

Ich wär sehr dankbar für Hinweise. Kann auch wenn nötig das ganze Programm posten oder Fragen beantworten (Manches könnte ein wenig doppelt gemoppelt sein, da ich so ohne Erfahrung rumprobiert habe und dann einfach alles gelassen habe, was einmal funktioniert Smile

also bitte nicht wundern.

Liebe Grüße

[EDITED, Jan, Bitte Code-Umgebung verwenden - Danke!]

Harald

Forum-Meister


	Beiträge: 24.495

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 22.08.2015, 19:15 Titel:

Hallo,

reproduzierbarer Code ist immer einfacher als für sich nicht lauffähige Teile.
Zudem bitte immer posten, was nicht funktioniert: gibt es eine Fehlermeldung? Ein unerwartetes Verhalten? Welches?
Statt eval würde ich für den Zweck der Umwandlung str2func empfehlen.

Differenzieren kannst du entweder symbolisch oder indem du Gradienten über Sekantensteigungen diff(y) ./ diff(x) annäherst.

Grüße,
Harald


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 |

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.