Parameter Gaussche Glockenkurve - fminsearch - Mein MATLAB Forum

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Parameter Gaussche Glockenkurve - fminsearch

Maria:)

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 04.05.2019, 17:06 Titel: Parameter Gaussche Glockenkurve - fminsearch

Hallo zusammen,

vor mir liegt folgende Aufgabenstellung. Ich soll ein Programm schreiben, welches Daten einliest und anschließend eine Guass'sche Glockenkurve anfittet.

Es sollen vier Parameter mit Hilfe der Funktion fminsearch bestimmt werden.

Mein Code sieht aktuell wie folgt aus:

Code:

clear
clear all

%Daten einlesen
[file,~]=uigetfile('*.dat','Datei auswaehlen');

[FID]=fopen(file,'r');
input=textscan(FID,'%s','Headerlines',3);
[nanz,~]=size(input{1,1});

x=zeros(nanz,1); %Zeit in ms
p=zeros(nanz,1); %Leistung in mW

%Daten sortieren
for n=1:nanz
string=char(input{1,1}(n));
k=strfind(string,',');
x1=string(1:k-1);
p1=string(k+1:end);

x1=str2double(x1);
p1=str2double(p1);
x(n)=x1;
p(n)=p1;
end

%Eingelesene Daten plotten
figure (1)
subplot(2,2,1)
plot(x,p);

title('Ursprüngliche Daten');
xlabel('Zeit [ms]');
ylabel('Leistung [mW]');

a1=1; a2=1;
c1=1;
t1=1;
a0=[a1,c1,t1,a2];

[a0neu]= fminsearch(@gauss_min,a0,[],x,p);
f=a0neu(1).*exp(-(x-a0neu(3)).^2/a0neu(2)^2)+a0neu(4);
subplot(2,2,2)

plot(x,f,'r')

function q = gauss_min(a0,x,p)
f=a0(1).*exp(-(x-a0(3)).^2/a0(2)^2)+a0(4);
q=sum((p-f).^2);
end

Funktion ohne Link?

Ich bekomme keine richtige Fehlermeldung, sonder nur folgenden "Hinweis":
Exiting: Maximum number of function evaluations has been exceeded
- increase MaxFunEvals option.
Current function value: 5620749.629258

Die zweite Kurve, die ich plotte - sieht der ersten überhaupt nicht ähnlich...
Ich bin bezüglich Matlab noch komplett unerfahren und würde mich wirklich freuen, wenn mir jemand helfen kann!

Liebe Grüße
Maria Smile

Harald

Forum-Meister


	Beiträge: 24.448

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 04.05.2019, 20:34 Titel:

Hallo,

bessere Startwerte können helfen.

Meine Empfehlung wäre, statt fminsearch eine Funktion für Kleinste-Quadrate-Optimierung, z.B. lsqcurvefit zu verwenden.

Zudem ist die verwendete Syntax von fminsearch veraltet und nicht mehr dokumentiert, sollte daher nicht mehr verwendet werden.

Grüße,
Harald
_________________

1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)

Maria:)

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 05.05.2019, 11:27 Titel: Parameter Gaussche Glockenkurve - fminsearch

Hallo Harald,

an das hab ich schon gedacht, allerdings wird in der Aufgabenstellung explizit von fminsearch gesprochen - muss ich also verwenden:/

Ich hab jetzt allerdings meine Startwerte angepasst und siehe da - es funktioniert Very Happy

Vielen lieben Dank!

Liebe Grüße
Maria:)


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