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

Polynom Fit

 

Detox
Forum-Fortgeschrittener

Forum-Fortgeschrittener


Beiträge: 68
Anmeldedatum: 25.10.14
Wohnort: Berlin
Version: ---
     Beitrag Verfasst am: 09.06.2017, 13:42     Titel: Polynom Fit
  Antworten mit Zitat      
Hi Leute,

ich bin gerade dabei King's Law für die Hot-Wire Anemometry zu implementieren bzw. es geht darum die Koeffizienten dafür herauszufinden. Meine Idee war mit fittype da ran zu gehen. Das Problem was ich aber habe ist,dass ich meinen Exponenten glaube ich nicht definieren kann.

Die Formel sieht wie folgt aus E^2 = A + B*u^n wobei n in meinem Beispiel gegeben ist und die Vektoren E und u sind ebenfalls gegeben.

Code:

start = [3;0.5;0.45];  
up = [inf;inf;inf];
low = [-inf;-inf;-inf];
options = fitoptions('Method','NonlinearLeastSquares','Upper',up,'Lower',low,'Display','iter','TolFun',10^-10,'TolX',10^-10,'StartPoint',start);

f = fittype('sqrt(A+B*x^n)');

res1 = fit(u,E,f,options)
 


Nun die Frage, ob das nicht irgendwie eleganter lösbar ist das Problem und ich habe noch nicht so ganz verstanden was der StartPoint genau macht bzw. ob dieser wichtig ist.

Besten Dank für eure Hilfe! Very Happy
Private Nachricht senden Benutzer-Profile anzeigen


Harald
Forum-Meister

Forum-Meister


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

mit fitlm (das Modell ist ja linear in A und B) sollte das deutlich einfacher gehen. Unter der Annahme von Spaltenvektoren:
Code:
x = u.^n;
y = E.^2;
mdl = fitlm(x, y)


Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
Detox
Themenstarter

Forum-Fortgeschrittener

Forum-Fortgeschrittener


Beiträge: 68
Anmeldedatum: 25.10.14
Wohnort: Berlin
Version: ---
     Beitrag Verfasst am: 10.06.2017, 14:24     Titel:
  Antworten mit Zitat      
Besten Dank Harald!

Werde ich gleich mal austesten Smile

Edit: Mit Polyfit müssten A und B bekannt sein, oder?
Private Nachricht senden Benutzer-Profile anzeigen
 
Harald
Forum-Meister

Forum-Meister


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

wenn A und B bekannt wäre, was würde dann noch gesucht?
Du kannst ebensogut
Code:
coeff = polyfit(x, y, 1)
verwenden, um A und B zu bestimmen.

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