Geradenapproximation

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Lisa

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 02.06.2013, 14:30 Titel: Geradenapproximation

Hallo,
ich habe so ein fertiges Programm bekommen wo ich Geraden annähern soll für verschiedene Spektren, allerdings versteh ich nicht genau wie ich da die Koordinaten ändern kann. Es nähert mir immer die gleiche Gerade an.

Code:

function [cut,parameter,Approximationsfehler] = piecewiseLinearApproximation(x,y)
penalty = 1;
step = 10;

%%
y = y(:);
Nsamples = length(y);

startLength = round(0.1*Nsamples); % length for the first and last part

NfirstPart = Nsamples-2*startLength+1; % number how often the length of the first part may be varied
NmiddlePartMax = Nsamples-2*startLength+1;
Napprox = (ceil(NfirstPart/step)+1)/2 * ceil(NmiddlePartMax/step); % bei step > 1 weniger....

results = zeros(Napprox,3+3*2+2);

i = 0;
for iFirst = 1: step: NfirstPart

ypart = y(1:startLength+iFirst-1);
xpart = x(1:startLength+iFirst-1);

cut1 = startLength+iFirst-1;

[ai1,~,~,e1] = curvefit(1,xpart, ypart);

%NmiddlePart = Nsamples-3*startLength+1 - (iFirst+1);
NmiddlePart = NmiddlePartMax - (iFirst + 1);

for iMiddle = 1: step : NmiddlePart

ypart = y(startLength+iFirst : startLength+iFirst+startLength);
xpart = x(startLength+iFirst : startLength+iFirst+startLength);

cut2 = startLength+iFirst+startLength;

[ai2,~,~,e2] = curvefit(1,xpart, ypart);

% last part

ypart = y(startLength+iFirst+startLength+1 : Nsamples);
xpart = x(startLength+iFirst+startLength+1 : Nsamples);

[ai3,~,~,e3] = curvefit(1,xpart, ypart);

i = i+1;
results(i,:) = [e1 e2 e3 ai1' ai2' ai3' cut1 cut2]; %[e1 e2 e3 ai1(1) ai1(2) ai2(1) ai2(2) ai3(1) ai3(3) cut1 cut2]
end
end

results(i+1:end,: ) = [];

Ehelp = [results(:,10) results(:,11)-results(:,10) Nsamples-results(:,11)]; % [cut1 cut2-cut1 Nsamples-cut2];
results(:,1) = results(:,1)*penalty;

E = sum(results(:,1:3)./Ehelp,2);
[~,iBest] = min(E);

Approximationsfehler = results(:,1:3); % führe die Variable nachträglich nochmal nach draußen....
% Variablennamen überarbeiten! - wurde aber eigentlich nur für BMT2013 benötigt

% figure; plot(x,y); hold on
% for i = 1:3
% plotFittingCurve([x(1) x(end)],results(iBest,3+(i-1)*2+1:3+(i-1)*2+2)');
% end

cut = [results(iBest,10:11)];
% convert cut to x
dx = x(2)-x(1);
cut = cut*dx;
parameter = reshape(results(iBest,4:9),2,3)';

Funktion ohne Link?

Jan S

Moderator


	Beiträge: 11.057

	Anmeldedatum: 08.07.10

	Wohnort: Heidelberg

	Version: 2009a, 2016b

Verfasst am: 03.06.2013, 01:11 Titel: Re: Geradenapproximation

Hallo Lisa,

Ok. Und was ist nun Deine Frage?
Das Programm enthält keine Dokumentation. Deshalb ist das Programm unbenutzbar. Es ist zwar möglich, dass Matlab-Experten herum-Raten, was das Programm wie machen soll. Aber da sinnvolle Kommentare im Code fehlen, lässt sich das Programm weder erweitern, noch ändern und erst recht nicht debuggen.

Deshalb lautet mein wohl deprimierender Vorschlag: Löschen und neu schreiben.

Gruß, Jan


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.