杭州交通信息网舞台节目表演:matlab 求解
来源：百度文库 编辑：杭州交通信息网 时间：2024/04/30 18:18:10
已知某非线性方程组如下：
ff(1)=(3-5*x(1))*x(1)+1-2*x(2)＝0
for k=2:9
ff(k)=(3-5*x(k))*x(k)+1-x(k-1)-2*x(k+1)＝0
end
ff(10)=(3-5*x(10))*x(10)+1-x(9)＝0
试求该方程组的解。
ff(1)=(3-5*x(1))*x(1)+1-2*x(2)＝0
for k=2:9
ff(k)=(3-5*x(k))*x(k)+1-x(k-1)-2*x(k+1)＝0
end
ff(10)=(3-5*x(10))*x(10)+1-x(9)＝0
试求该方程组的解。
function y=sol()
x0=ones(10,1);
[x,fval]=fsolve(@myfun,x0);
y=x;
function y=myfun(x)
y=[(3-5*x(1))*x(1)+1-2*x(2)];
for k=2:9
y =[y;(3-5*x(k))*x(k)+1-x(k-1)-2*x(k+1)];
end
y=[y;(3-5*x(10))*x(10)+1-x(9)];
结果：
ans =
0.8810
-0.1195
-0.1527
0.2847
0.8509
-0.1930
-0.3155
-0.1282
0.4237
0.7527
function y=sol()
x0=ones(10,1);
[x,fval]=fsolve(@myfun,x0);
y=x;
function y=myfun(x)
y=[(3-5*x(1))*x(1)+1-2*x(2)];
for k=2:9
y =[y;(3-5*x(k))*x(k)+1-x(k-1)-2*x(k+1)];
end
y=[y;(3-5*x(10))*x(10)+1-x(9)];
结果：
ans =
0.8810
-0.1195
-0.1527
0.2847
我怎么算出来是
ans =
0.7916
0.1168
0.2646
0.7067
0.1560
0.3281
0.7036
0.1207
0.3092
0.7599