%Programm: Levenberg-Marquardt-Algorithmus

%Modell: (t-a)^2+exp(b*(t^2+y^2)-5=0

%Fummelparameter: a,b
%Gesucht Beste Parameter: a,b
clc
clear all
format Longg
%----------------Stützpunkte----------------
t=[2 3 4];
y=[0 2 0];
%-------------------------------------------
%----------------Startvektor----------------
a=4;
b=0;
x0=[a b];

 

m=length(t);
beta0=0.2;
beta1=0.8;
mu0=1
%-------------------------------------------
%-------------Zugriff auf Symbolik----------
syms a b real
%-------------------------------------------
%--------------Erzeugung-Startwertematrix---
for i=1:m
f(i)=(t(i)-a)^2+exp(b*(t(i)^2+y(i)^2))-5; % <- Modellfunktion
end
f=f';

Parameter=[a b];
%-------------------------------------------
%--------------Erzeugung-Jacobimatrix-------
v=[a b];
j=jacobian(f,v);
%-------------------------------------------
%--------------Matrix-Substitution----------
DF=subs(j,Parameter,x0);
F=subs(f,Parameter,x0);
%-------------------------------------------

maxit=10;
tol=norm(F)/mu0;