我的代码计算与牛顿-Rhapson方法的非线性方程系统,和我试图显示为“误差矢量”(H1)的函数“iterration”(I)的曲线图。
也许是因为我是一个标量和H1是矢量?
如何解决制造阴谋像我想要的代码?
Syms.X Yž
U =(X * Y)-z ^ 3-2;
V =(X * Y * Z)-x ^ 2 + Y ^ 2-4;
W = EXP(X)-exp(Y)+ COS(Z)-6-;
XO = 1;哟= 1; ZO = 1;
设k = 10;
我= 0;
一个雅可比=([U; V; W],[X,Y,Z]);
尽管我
J =潜艇(A,[X,Y,Z],[XO,YO,ZO]);
j_i = j的^ -1;
UO =潜艇(U,[X,Y,Z],[XO,YO,ZO]);
VO =潜艇(V,[X,Y,Z],[XO,YO,ZO]);
WO =潜艇(W,[X,Y,Z],[XO,YO,ZO]);
溶胶= [XO;哟; ZO] -j_i * [UO; VO; WO];
溶胶=双(溶胶);
X1 =溶胶(1,1);
H1 = j_i * [UO; VO; WO];
XO =溶胶(1,1);
哟=溶胶(2,1);
ZO =溶胶(3,1);
我= i + 1;
图(I,H1,'r');
抓住在
Xlabel(“误差矢量”)
ylabel('Iterration')
结尾
X = XO;
Y =哟;
Z = ZO;
x,y,z