matlab 最速下降法求解clc;clear;disp('最速下降法')syms x1 x2 lam;x=[x1 x2];t=0; X0=[1;1]; y=(x1^2+x2-11)^2+(x1+x2^2-7)^2;e=0.01;for n=1:10 a=subs(diff(y,x1),x,X0); b=subs(diff(y,x2),x,X0); yj=[a;b]; z=subs(y,x,X0-lam*y

来源:学生作业帮助网 编辑:作业帮 时间:2024/05/04 16:29:19
matlab 最速下降法求解clc;clear;disp('最速下降法')syms x1 x2 lam;x=[x1 x2];t=0;         X0=[1;1]; y=(x1^2+x2-11)^2+(x1+x2^2-7)^2;e=0.01;for n=1:10    a=subs(diff(y,x1),x,X0);    b=subs(diff(y,x2),x,X0);    yj=[a;b];     z=subs(y,x,X0-lam*y

matlab 最速下降法求解clc;clear;disp('最速下降法')syms x1 x2 lam;x=[x1 x2];t=0; X0=[1;1]; y=(x1^2+x2-11)^2+(x1+x2^2-7)^2;e=0.01;for n=1:10 a=subs(diff(y,x1),x,X0); b=subs(diff(y,x2),x,X0); yj=[a;b]; z=subs(y,x,X0-lam*y
matlab 最速下降法求解
clc;
clear;
disp('最速下降法')
syms x1 x2 lam;
x=[x1 x2];
t=0;
X0=[1;1];
y=(x1^2+x2-11)^2+(x1+x2^2-7)^2;
e=0.01;
for n=1:10
a=subs(diff(y,x1),x,X0);
b=subs(diff(y,x2),x,X0);
yj=[a;b];
z=subs(y,x,X0-lam*yj);
z1=diff(z,lam);
lambda=solve(z1,'lam');
lambda=double(lambda);
a=subs(diff(y,x1),x,X0);
b=subs(diff(y,x2),x,X0);
fan=norm(yj);
if fan>e
X1=X0-lambda(1)*yj;
X0=X1;
t=t+1;
else
break
end
end
X0
y=(subs(y,x,X0));
y
t
这是我编的谁能帮我看看有啥问题
明白最速下降法的帮我看看 得不到结果 谢谢了

matlab 最速下降法求解clc;clear;disp('最速下降法')syms x1 x2 lam;x=[x1 x2];t=0; X0=[1;1]; y=(x1^2+x2-11)^2+(x1+x2^2-7)^2;e=0.01;for n=1:10 a=subs(diff(y,x1),x,X0); b=subs(diff(y,x2),x,X0); yj=[a;b]; z=subs(y,x,X0-lam*y
你写的程序我已经运行过了,没有任何语法错误