[求助]关于非线性代数方程组求解的问题

无水1324

请各位高手指教！
要解一个非线性的方程组，其中有一个参数变化时解的曲线。
我自己编了一个程序运行时出错，程序如下：
function q=myfun(p,flag,wh)
qa2=p(1);qm2=p(2);
q(1)=-qa2+2.*wh^2/(100.*(.5000-.3183*(2.-2.*qa2)*(1.-.2146*(1.-1.*qm2)^2/qa2^2)/qm2-1.*wh^2)^2+wh^2)^(1/2);
q(2)=-qm2+.1000/(1.+.5000*qa2*(.6366+.3635*(1.-1.*qm2)^2/qa2^2-1.*(1.+qm2)/qa2)/qm2);
% 定义方程组
wh=0.001:0.01:3;
x=[];
for i=1:length(wh)
x=fsolve('myfun',[1,1],[],wh(i));
y(i)=x(1);yy(i)=x(2);
end
plot(wh,y)
%%%%%%解、画图
请高手帮我看一下,谢谢！

无水1324

fsolve('')中是不是不能像解微分方程一样用如下的调用方法：
fsolve('myfun',[1,1],[],wh(i));
到底错在哪里？

lovelyvickey

使用fsolve函数和求解非线性方程组有以下几种：
1、X=fsolve(fun,x0)，x0为初始矩阵求解fun
2、X=fsolve(fun,x0,options),以options为选择参数的输入变量，详细看optimset的帮助信息。
3、X=fsolve(fun,x0,options，P1,P2）将问题定性参数P1,P2等直接赋值给函数fun(X,P1,P2,...)
4、[X,fval]=fsolve(fun,P1,P2,P3....),返回原方程在X处的值。