Summe Abstandsquadrate als Funktion von Parametern

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

matlabel

Forum-Fortgeschrittener


	Beiträge: 83

	Anmeldedatum: 28.01.16

	Wohnort: ---

	Version: 2024a

Verfasst am: 06.08.2024, 16:00 Titel: Summe Abstandsquadrate als Funktion von Parametern

Hallo,
ich habe Messwertpaare (x,y). Zusätzlich habe ich eine Funktion f(a,b,x). a und b sollen im Sinne der kleinsten Abstandsquadrate an diese Messwerte angepasst werden.Für a bzw. b gebe ich einen Wertebereich a_min...a_max bzw. b_min...b_max vor.
Es geht mir aber nicht direkt um die Anpassung, sondern - für Lernzwecke - um eine Grafik, in der die Summe der Abstandsquadrate [f(a,b,x-y]² als Funktion von a und b dargestellt wird.

Code:

fStr='@(a,b,x)a*x+b*sqrt(x)'; %Beispiel für eine Funktion
x=[1,2,3,4];
y=[8,8.5,9,9.2]; %völlig fiktive Werte
f=str2func(fStr);
a_min=1; a_max=2;
b_min=2; b_max=2.8;
Steps=50;
da=(a_max - a_min)/Steps;
a=a_min :da: a_max;
db=(b_max -b_min )/Steps;
b=b_min : db : b_max;
[A,B]=meshgrid(a,b);

% Wie baue ich hier den Zugriff auf a und b ein bzw. wie durchlaufe ich das Meshgrid?

resStr='@(a,b,x,y)(f(a,b,x)-y).^2';
resfunc=str2func(resStr);
res=resfunc(a,b,x,y);

% Wie baue ich zu jedem A,B den jeweiligen Wert von res in eine Funktion G(A,B) so ein, dass ich ich G plotten kann?

Funktion ohne Link?

Harald

Forum-Meister


	Beiträge: 24.495

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 07.08.2024, 17:38 Titel:

Hallo,

zunächst würde ich str2func nur verwenden, wenn es nötig ist - beispielsweise weil die Eingabe der Funktion aus einer App erfolgt. Das zweite str2func erscheint aber unnötig:
resfunc = @(a,b,x,y)(f(a,b,x)-y).^2;

Eine komplette Vektorisierung halte ich hier für sehr schwierig. Am einfachsten dürfte eine doppelte for-Schleife sein:

Code:

for ka = 1:length(a)
for kb = 1:length(b)
res(ka, kb) = sum(resfunc(a(ka), b(kb), x, y));
end
end

Funktion ohne Link?

Statt

Code:

da=(a_max - a_min)/Steps;
a=a_min :da: a_max;

Funktion ohne Link?

wäre übrigens die einfachere Möglichkeit

Code:

a = linspace(a_min, a_max, Steps + 1);

Funktion ohne Link?

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 ;)

matlabel

Themenstarter

Forum-Fortgeschrittener


	Beiträge: 83

	Anmeldedatum: 28.01.16

	Wohnort: ---

	Version: 2024a

Verfasst am: 10.08.2024, 11:29 Titel: Danke Harald, so klappt's


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.