fit fittet überhaupt nicht

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

pospiech

Forum-Century


	Beiträge: 130

	Anmeldedatum: 03.09.07

	Wohnort: ---

	Version: ---

Verfasst am: 04.03.2010, 15:05 Titel: fit fittet überhaupt nicht

Ich habe eine fit-Funktion wie folgt definiert:

Code:

function [Fit Error] = cgFit(XAXIS, DATA, STARTVALUES)

%Filter für NaN Werte
ok = isfinite(XAXIS) & isfinite(DATA);
%Fit Funktion definieren
f_type = fittype('rectpuls(0.5*(sin(2*pi*((x-xo)+c*(x-xo)^2)/ps * pi/180)+1))',...
'dependent',{'y'},'independent',{'x'},...
'coefficients',{'c', 'xo', 'ps'});
%Fit ausführen
[Fit Error] = fit(XAXIS(ok),DATA(ok),f_type,'Startpoint',STARTVALUES);

Funktion ohne Link?

das ganze generiert eine Rechteckfunktion mit variabler Periodenlänge.

Wenn ich das ganze zum fitten aufrufe, mit:

Code:

[Fit Error] = cgFit(XAXIS', DATA, STARTVALUES)

Funktion ohne Link?

bekomme ich aber nur exakt die Startwerte zurück.

Warum versucht Matlab gar nicht erst zu fitten?

Matthias

pospiech

Themenstarter

Forum-Century


	Beiträge: 130

	Anmeldedatum: 03.09.07

	Wohnort: ---

	Version: ---

Verfasst am: 04.03.2010, 15:45 Titel:

Ich kann den fit sogar umschreiben zu:

Code:

[Fit Error] = cgFit(XAXIS', DATA, STARTVALUES);

function [Fit Error] = cgFit(XAXIS, DATA, STARTVALUES)

options=optimset('Display','iter');
options.TolX = 1e-6;
options.TolFun = 1e-6;
Estimates=fminsearch(@cgfit2,STARTVALUES,options,XAXIS,DATA);

Fit = Estimates;
Error = cgfit2(Estimates,XAXIS,DATA);

function sse=cgfit2(params,Input,Actual_Output)
c=params(1);
xo=params(2);
ps=params(3);

Fitted_Curve = zeros([1 length(Input)]);
for xi=1:length(Input)
x = Input(xi);
Fitted_Curve(xi) = pi*rectpuls(0.5*(-sin(2*pi*((x-xo)+c*(x-xo)^2)/ps * pi/180)+1)) - pi;
end
Error_Vector=Fitted_Curve - Actual_Output';
% When curvefitting, a typical quantity to
% minimize is the sum of squares error
sse=sum(Error_Vector.^2);

Funktion ohne Link?

und ich erhalte nur exakt die Fiteingabedaten zurück.

Anhand der Ausgabe sehe ich nur das der Fit praktisch nichts macht:

Zitat:

Iteration Func-count min f(x) Procedure
0 1 5245.89
1 4 5225.04 initial simplex
2 6 5142.36 reflect
3 7 5142.36 reflect
4 12 5142.36 shrink
5 13 5142.36 reflect
6 15 5135.08 reflect
7 17 5135.08 contract inside
8 18 5135.08 reflect
9 20 5135.08 contract inside
10 22 5135.08 contract inside
11 24 5135.08 contract inside
...

woran kann ich noch drehen?

EDIT: Ich erhalte nicht wirklich die Eingabeparameter zurück, aber diese sind nur minimal variiert worden, so dass überhaupt keine Verbesserung stattgefunden hat.

Wie kann ich die Variation vergrößern?


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.