Fittingwerte Plotten

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

hackhuebi

Forum-Anfänger


	Beiträge: 24

	Anmeldedatum: 28.04.15

	Wohnort: ---

	Version: ---

Verfasst am: 17.11.2015, 20:06 Titel: Fittingwerte Plotten

Hallo Leute,
ich habe folgenden Code:

Code:

xtest_reduced=zeros(n_red,10);
for loop2=1:10
for loop=1:n_red
if loop==n_red
xxxx=0;
end

ibase=(loop-1)*n_dist+halfwidth_weight;

i1=ibase-halfwidth_weight;
if i1<1 i1=1; end
i2=ibase+halfwidth_weight;
if i2>nt i2=nt; end
wi1=i1-ibase+halfwidth_weight+1;
wi2=i2-ibase+halfwidth_weight+1;
int1=trapz(weightfun(wi1:wi2));
xtest_reduced(loop,loop2)=sum(xtest(i1:i2,loop2).*weightfun(wi1:wi2))/int1;
end
end

figure('NumberTitle','off','Name','Messwerte+Reduzierung','units','normalized','outerposition',[0 0 1 1]);
hold on;
plot(ttest,xtest(:,1),'.','MarkerSize',1);
plot(ttest,xtest(:,2),'.','MarkerSize',1);
plot(ttest_reduced,xtest_reduced(:,1),'r.-');
plot(ttest_reduced,xtest_reduced(:,2),'c.-');

titlestr=sprintf('Channel %d',channel);
title(titlestr);

linfit=fittype(@(a,b,x) a*x+b);
index_slide_step=1;
deltal=5;
jdet_start_idx=17000/main_deltat; % exclude the first 1.7e4 s from jump detection
jdet_start_idx=deltal+1;

n_max=uint32(n_red-2*deltal);
tmid=zeros(n_max,1);
step_detect=zeros(n_max,1);
ni=2*deltal+1;
stride0=10;

Rolling Eyes

for loop=1:10
detect_range=jdet_start_idx-deltal:stride0:n_red;

s2=fitoptions('Method','NonlinearLeastSquares','Startpoint',[28 0.0016 0.23 11 1.2e4 2.54 0.26 1.95e4],...
'Lower',[0 ,0 ,0 ,1e-4,1 ,1e-3 , 1e-3, 0],...
'Upper',[Inf ,1 ,1 ,1e4 ,1e7,1e2 ,1e2, 1e9],'TolFun',1e-8,'TolX',1e-8,...
'MaxFunEvals',12000,'MaxIter',2000,'Algorithm','Trust-Region');
f2= fittype('a./(1+b*x.^c)+d*((e./(x+h)).^(f+g)-(e./(x+h)).^f)','options',s2);
Rolling Eyes

[fit2,gof2,out2]=fit(ttest_reduced(detect_range),xtest_reduced(detect_range,loop),f2);
if out2.firstorderopt>1
warning('Fit is probably not converged!\nFirst order optimality: %.1f',out2.firstorderopt);
end
end
xtest_reduced(detect_range,1)

Rolling Eyes

plot(ttest_reduced,fit2(ttest_reduced),'g-','Linewidth',2);

Funktion ohne Link?

xtest_reduced hat 10 Spalten, da es 10 Messkanäle enthält. Soweit so gut. Ich möchte jetzt aber die Kurven in den Plots für mehrere Kanäle gleichzeitig darstellen. Vorher war das Skript so geschrieben, dass nur für einen Kanal, der zuvor angegeben wurde, alle Kurven geplottet wurden.

Jetzt habe ich aber ein Problem. Ich bin mir nicht sicher ob das mit den Loops für das Fitten der einzelnen Kurven so möglich ist. Sollte aber laut meinem Verständnis passen. Nur weiß ich dann nicht wie ich im Plot die Y-Werte, also xtest_reduced, der einzelnen Spalten aufrufen soll. Habe vor die fraglichen Zeilen mal einen Smiley gesetzt. Vllt kann mir ja wer von euch helfen. Vielen Dank


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.