#求助:matlab中fsolve指令该如何用
先谈一些我认识的fsolve指令
近期在使用fsolve指令求解一个非线性代数方程式,总共有2002个未知数,2002个方程,我有自己去了解一些关于fsolve的用法,首先他主要有三种算法,trust-region-dogleg(信任区间狗腿法,默认方法)、trust-region(信任区间法)、levenberg-marquardt(简称LM),其他的帮助会有在其他博主的文章中找到,再次不赘述
讲述一下自己的问题
我将自己的函数放在b1的matlab中,
function F=b1(w)
global dt a nx y0 g S0
F(1:nx,1)=1/dt*a*w(1:nx)+1/2*w(nx+1:2*nx)-y0(1:nx)-w(2*nx+1)*ones(nx,1);
F(nx+1:2*nx,1)=1/dt*a*w(nx+1:2*nx)+1/2*w(nx+1:2*nx).^2./w(1:nx)+g/4*w(1:nx).^2...
-1/2*a*g*w(1:nx).*S0-y0(nx+1:2*nx)-w(2*nx+2)*ones(nx,1);
F(2*nx+1,1)=w(nx+1)-4.42;
F(2