disser1.pdf (2006 KB) - Ernst-Moritz-Arndt-Universität Greifswald

C. MATLAB-Programme
end; ps=[ps1;ps2;ps3;ps4;ps5],
f1=@(x1,x2,x3) ps(1,1)*x1^2+ps(2,1)*x1+ps(3,1) ...
+ps(4,1)*x1*x2+ps(5,1)*x1*x3;
f2=@(x1,x2,x3) ps(1,2)*x2^2+ps(2,2)*x2+ps(3,2) ...
+ps(4,2)*x1*x2+ps(5,2)*x2*x3;
f3=@(x1,x2,x3) ps(1,3)*x3^2+ps(2,3)*x3+ps(3,3) ...
+ps(4,3)*x1*x3+ps(5,3)*x2*x3;
f1x1=@(x1,x2,x3) ps(1,1)*2*x1+ps(2,1) ...
+ps(4,1)*x2+ps(5,1)*x3;
f1x2=@(x1) ps(4,1)*x1;
f1x3=@(x1) ps(5,1)*x1;
f2x1=@(x2) ps(4,2)*x2;
f2x2=@(x1,x2,x3) ps(1,2)*2*x2+ps(2,2) ...
+ps(4,2)*x1+ps(5,2)*x3;
f2x3=@(x2) ps(5,2)*x2;
f3x1=@(x3) ps(4,3)*x3;
f3x2=@(x3) ps(5,3)*x3;
f3x3=@(x1,x2,x3) ps(1,3)*2*x3+ps(2,3)...
+ps(4,3)*x1+ps(5,3)*x2;
x=[1;1;1];
i=0;
epsilon=1e-6
while norm([f1(x(1),x(2),x(3));f2(x(1),x(2),x(3)); ...
f3(x(1),x(2),x(3))]) > epsilon
sol=[f1x1(x(1),x(2),x(3)),f1x2(x(1)),f1x3(x(1)); ...
f2x1(x(2)),f2x2(x(1),x(2),x(3)),f2x3(x(2));...
f3x1(x(3)),f3x2(x(3)),f3x3(x(1),x(2),x(3))] \ ...
-[f1(x(1),x(2),x(3));f2(x(1),x(2),x(3)); ...
f3(x(1),x(2),x(3))];
x(1)=x(1)+sol(1);
x(2)=x(2)+sol(2);
x(3)=x(3)+sol(3);
disp([’Iteration: ’ num2str(i)]);
disp([’Lösung : ’ num2str(x’,’%1.7f\t’)]);
i=i+1;
end;
C.2. Programm 2. Berechnung der Adjungierten.
clear all;
format long;
e=[0.4, 0.6, 0.6];
146