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

Fit mit fminsearch an Daten nach exponentialverteilung
 

Gregor

Gast

Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 26.01.2010, 11:49     Titel: Fit mit fminsearch an Daten nach exponentialverteilung
  Antworten mit Zitat      
Ich möchte gern an Daten eine Exponentielle Kurve anfitten. Habe dazu einen Code von http://www.cae.tntech.edu/help/programs/matlab/curve_fitting für mein Problem adaptiert. Funktioniert aber leider nicht, ich verstehe nicht warum.

Hier der Code im ersten m-file:
Code:

mydata=[%0       62
           %15       61
           30       62
           45   76
      60   96
   75   125
   90   143
   105   153
        120   166
   135   177
   150   181
   165   187
   180   190
   195   194
   210   191
   225   186];


xdata=mydata(:,1);
   ydata=mydata(:,2);
   w_start=[60 196 105 100];  %startparameter
   options=optimset('MaxFunEvals',10e9);
   w_guess=fminsearch('myfitfunc',w_start,options);
   w1_guess=w_guess(1); %schwellwert
   w2_guess=w_guess(2); %sättigungswert
   w3_guess=w_guess(3); %t_0 (Beginn Enhancement)
   w4_guess=w_guess(4); %Tau
   ystart=w_start(2)-w_start(1)*exp(-(xdata-w_start(3))/w_start(4));
   yfit=w1_guess-w1_guess*exp(-(xdata-w3_guess)/w4_guess);
   plot(xdata,ydata,'ro',xdata,ystart,'g',xdata,yfit,'b');
   legend('Original Data','Original Fit','Final Fit',4);
   title(sprintf('Curve Fit : yfit=%.3f+%.3f*exp(-(x-%.3f)/%.3f)', ...
    w2_guess,w1_guess,w3_guess,w4_guess));
   xlabel('x')
   ylabel('y')
 

im zweiten file:
Code:

function err=myfitfunc(w)
  global xdata ydata;
  w1=w(1);
  w2=w(2);
  w3=w(3);
  w4=w(4);
    yfit=w2-w1*exp(-(xdata-w3)/w4);
  err=norm(ydata-yfit);

 


Eigentlich sollten die Parameter minimiert werden. Aber das Resultat verschiebt immer nur die Fitfunktion entlang der y-achse, unabhängig von den gewählten Startparametern. Vielleicht hat jemand eine Idee.

Danke!

Gregor


Nacken

Gast

Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 10.07.2013, 14:21     Titel:
  Antworten mit Zitat      
Hi,

die Verschiebung kommt durch einen Schreibfehler bei yfit
und das eigentliche Fitten funktioniert bei dir nicht, das mit den globalen Variablen so nicht funktioniert. Sie müssen bereits im Hauptteil als global definiert werden. Dann klappts!

Grüße
 
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.