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

Messwerte approximieren

 

therion
Forum-Anfänger

Forum-Anfänger


Beiträge: 22
Anmeldedatum: 14.01.13
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 03.03.2015, 15:19     Titel: Messwerte approximieren
  Antworten mit Zitat      
Ich versuche Messwerte (im Anhang) zu approximieren, leider schient mein Ansatz bzw. die Anfangswerte nicht gut zu sein -- hat jemand vielleicht eine gute Idee für einen Ansatz?

Code:

Data =csvread('werte.m');
uout0 = Data(:,1);
t0 = Data(:,2)

options=optimset('TolFun',1e-13,'TolX',1e-13,'MaxFunEvals',Inf,'MaxIter',Inf);
costfun = @(c,t) (c(1) ./ (1 + c(2)* exp(c(3)*t)))
 
c = lsqcurvefit(costfun,[.5 .5 .5],t0,uout0,[],[],options)
 
figure(2)
plot(t0,uout0,'O','MarkerSize',3);
hold on
plot(t0,costfun(c,t0),'LineWidth',1.5)
 


und

Code:

Data =csvread('werte.m');
uout0 = Data(:,1);
t0 = Data(:,2)

options=optimset('TolFun',1e-13,'TolX',1e-13,'MaxFunEvals',Inf,'MaxIter',Inf);
costfun = @(c,t) (c(1) ./ (1 + c(2)* exp(c(3)*t)))

c = lsqnonlin(@(c) (uout - (c(1) ./ (1 + c(2)* exp(c(3)*t)))),[1 1 1],[],[],options)

plot(t0,uout0,'O','MarkerSize',3);
hold on
plot(t0,costfun(c,t0),'LineWidth',1.5)
 


werte.m
 Beschreibung:

Download
 Dateiname:  werte.m
 Dateigröße:  623 Bytes
 Heruntergeladen:  305 mal
Capture.PNG
 Beschreibung:

Download
 Dateiname:  Capture.PNG
 Dateigröße:  6.68 KB
 Heruntergeladen:  367 mal
Private Nachricht senden Benutzer-Profile anzeigen


Winkow
Moderator

Moderator



Beiträge: 3.842
Anmeldedatum: 04.11.11
Wohnort: Dresden
Version: R2014a 2015a
     Beitrag Verfasst am: 03.03.2015, 15:47     Titel:
  Antworten mit Zitat      
deine x werte sind sehr klein.
Code:
Data =csvread('werte.m');
uout0 = Data(:,1);
t0 = Data(:,2)

options=optimset('TolFun',1e-13,'TolX',1e-13,'MaxFunEvals',Inf,'MaxIter',Inf);
costfun = @(c,t) (c(1) ./ (1 + c(2)* exp(10^8*c(3)*t)))
 
c = lsqcurvefit(costfun,[.5 .5 .5],t0,uout0,[],[],options)
 
figure(2)
plot(t0,uout0,'O','MarkerSize',3);
hold on
plot(t0,costfun(c,t0),'LineWidth',1.5)
 

funktioniert da schon weitaus besser.
wobei das mir aber ehr nach tanh aussiht als nach einer exponentialfunktion
Code:
Data =csvread('werte.m');
uout0 = Data(:,1);
t0 = Data(:,2)

options=optimset('TolFun',1e-13,'TolX',1e-13,'MaxFunEvals',Inf,'MaxIter',Inf);
costfun = @(c,t) (c(1) ./ (1 + c(2)* tanh(10^8*c(3)*t)))
 
c = lsqcurvefit(costfun,[.5 .5 .5],t0,uout0,[],[],options)
 
figure(2)
plot(t0,uout0,'O','MarkerSize',3);
hold on
plot(t0,costfun(c,t0),'LineWidth',1.5)
 

liefert da bessere ergebnisse finde ich.
wenn man an sowas rumspielen will kann ich cftool empfehlen. da kann man schneller änderungen vornehmen.
_________________

richtig Fragen
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 - 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.