#求助:matlab中fsolve指令该如何用
先谈一些我认识的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 S0F(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*nx+2,1)=w(nx)-2;