Hilfe benötigt zu non-linear fitting

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

TeChierys

Forum-Anfänger


	Beiträge: 37

	Anmeldedatum: 14.11.13

	Wohnort: ---

	Version: ---

Verfasst am: 05.12.2015, 23:49 Titel: Hilfe benötigt zu non-linear fitting

Hallo zusammen,

ich sitze an meiner Masterarbeit und bin immer noch ein Noob in Matlab. Jede Hilfe ist erwünscht!

Ich suche mittels fzero-Funktion eine Zahlenvektor vola mit gegebenen V0, H, b, und Zahlenvektor mty und q. Einige von den habe ich aus einfachheitsgründen in Code reindefiniert.

Jedoch läuft das ganze nicht so richtig. Ich bekomme aus fzero einfach keinen Nullstelle für die Gleichung delta (unten). Ich bitte um Fehlersuche. Vermutung nach könnte ich jetzt 3 Fehlerarten vorstellen:
1. Vektorrechnung gescheitert. Ich weiß die Indizies hier sind schlampig gesetzt, z.b. der timediff ist ein [1*8] Vektor und vola ein [8*1].
2. quadrat auf Vola macht fzero schwach? Vielleicht soll ich ein lineare Variable die ersetzen?
3. Fehler in Schreibweise in fzero-Funktion

Ich bitte um euren kritischen Blick.

Code:

function [ vola] = calcvola2( V0,H,b)
q=[1,0.98,0.96,0.95,0.93,0.91,0.87,0.8];
mty=[0.5,1,2,3,4,5,7,10];

%%%%%%%%%Voreingabe%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
msize=length(mty);
int=zeros(msize+1,1);
vola=nan(msize,1);
tol = 1e-6;

volamax=1;
volarange=[0.01, 0.99];
%%%%%%%%%Diskretisierung Zeitabstand%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

timediff(1)=0.5;
for i=2:msize
timediff(i)=mty(i)-mty(i-1);
end

%Startwert
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
options = optimset('fzero');
options = optimset(options, 'TolX', tol(1), 'Display', 'off');

%Termstruktur
for t=1:msize
if isnan(int(t))
break
else
try
[vola(t),fval,exitflag]=fzero(@objfcn,volarange,options,...
V0,H,b,timediff(t),int(t),q(t))

catch
vola(t)=nan;
end

int(t+1)=int(t)+vola(t)^2*timediff(t);

end
end

vola = reshape(vola, msize, 1);

function delta=objfcn(vola,V0,H,b,timediff,int,q)
logterm=log(V0/H);
variation=(2*b-1)/2;
freeterm=(H/V0)^(2*b-1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Kontrolle:
if isnan(sqrt(int+vola^2*timediff))
break
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
dfirst=(logterm+variation*(int+vola^2*timediff))/sqrt(int+vola^2*timediff);
dsec=dfirst-2*logterm/sqrt(int+vola^2*timediff);

modellvalue=normcdf(dfirst)-freeterm*normcdf(dsec);

delta=q-modellvalue;

Funktion ohne Link?

Danke!

Harald

Forum-Meister


	Beiträge: 24.501

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 06.12.2015, 15:48 Titel:

Hallo,

für den Anfang hilft es sehr, wenn du sagst, was "läuft nicht so richtig" bedeutet.
Bekommst du eine Fehlermeldung? Welche?
Ein unerwartetes Ergebnis? Inwiefern?

Grüße,
Harald

TeChierys

Themenstarter

Forum-Anfänger


	Beiträge: 37

	Anmeldedatum: 14.11.13

	Wohnort: ---

	Version: ---

Verfasst am: 10.12.2015, 22:29 Titel:

Hallo Harald,

sorry für die verspätete Antwort. Weil ich selbst keiner Ahnung hatte, was schief ging. Anscheind liegt das ganze Problem an Startwertfestlegung. Der fzero.m habe ich immer noch nicht durchgeblickt, aber er fordert die Voraussetzung von Zwischenwertsatz, also dass an der Intervallgrenzen (volarange) eine negative UND eine positive Funktionswert existiert. Somit funktioniert die fzero nicht immer für meine Auswertungsfunktion delta...

Inzwischen habe ich das Problem aufgehoben. In dem ich entweder H groß, oder b klein wähle. Numerisch ist fzero halt instabil. Was könnte man anstatt fzero nehmen?

Harald

Forum-Meister


	Beiträge: 24.501

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 11.12.2015, 05:37 Titel:

Hallo,

Zitat:

Weil ich selbst keiner Ahnung hatte, was schief ging.

... aber irgendwoher musst du ja wissen, dass etwas schief geht.

Zitat:

Anscheind liegt das ganze Problem an Startwertfestlegung.

Welches Problem?

Zitat:

Was könnte man anstatt fzero nehmen?

z.B. fsolve .

Grüße,
Harald

TeChierys

Themenstarter

Forum-Anfänger


	Beiträge: 37

	Anmeldedatum: 14.11.13

	Wohnort: ---

	Version: ---

Verfasst am: 11.12.2015, 11:06 Titel:

Harald hat Folgendes geschrieben:

Hallo,

Zitat:

Weil ich selbst keiner Ahnung hatte, was schief ging.

... aber irgendwoher musst du ja wissen, dass etwas schief geht.

fzero untersucht zunächst zu angegebenen Intervallgrenzen

Code:

[0.01, 0.99]

Funktion ohne Link?

ein Vorzeichendifferenz. Als Fehlermeldung immer Vorzeichenfehler.

Zitat:

Anscheind liegt das ganze Problem an Startwertfestlegung.

Welches Problem?

Das Problem war, der fzero-Code lief bis o.g. Fehlermeldung und geht nicht in Nullstellesuche ein. Durch Startwertanpassung ist es nun möglich. Ich soll H über 50% von V0 wählen oder b unter 1 setzen.

Zitat:

Was könnte man anstatt fzero nehmen?

z.B. fsolve .

Danke für den Tipp, das werde ich nun probieren.
Grüße,
Harald

Harald

Forum-Meister


	Beiträge: 24.501

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 11.12.2015, 15:12 Titel:

Hallo,

... und es ist eben wichtig, dass du die Fehlermeldung nennst - am besten kopieren und nicht übersetzen.

Grüße,
Harald

TeChierys

Themenstarter

Forum-Anfänger


	Beiträge: 37

	Anmeldedatum: 14.11.13

	Wohnort: ---

	Version: ---

Verfasst am: 11.12.2015, 18:17 Titel:

Danke, werde ich beim nächsten Mal machen. Die Fehlermeldung hier hieß

"Error using ==> fzero at 283
The function values at the interval endpoints must differ in sign."

Weshalb ich zunächst keine Ahnung von Fehler hatte: ich habe die Anweisung try...catch dringelassen, die versichert für einen Output (nan...) in jeder Hinsicht .... Werde ich zunächst Finger davon lassen.

[EDITED, Jan, Bitte kein Top-Quoting der gesamten vorherigen Nachricht - Danke]


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 - 2025 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.