求个非线性方程组,六个参数,急.用MATLAB 求不出来哈,代码如下>> syms qn dn dr pn pr cn cr sn sr hn hr t e a o A T u1 k>> syms Fn Fr fn fr>> syms zn zr>>o=0.8;>>a=0.1;>>cn=0.4;>>cr=0.2;>>k=2;>>T=1;>>hn=0.1;>>sn=0.1;>>hr=0.05;>>sr=0.
来源：学生作业帮助网 编辑：六六作业网 时间：2024/12/18 18:12:14
求个非线性方程组,六个参数,急.用MATLAB 求不出来哈,代码如下>> syms qn dn dr pn pr cn cr sn sr hn hr t e a o A T u1 k>> syms Fn Fr fn fr>> syms zn zr>>o=0.8;>>a=0.1;>>cn=0.4;>>cr=0.2;>>k=2;>>T=1;>>hn=0.1;>>sn=0.1;>>hr=0.05;>>sr=0.
求个非线性方程组,六个参数,急.
用MATLAB 求不出来哈,代码如下
>> syms qn dn dr pn pr cn cr sn sr hn hr t e a o A T u1 k>> syms Fn Fr fn fr>> syms zn zr
>>o=0.8;
>>a=0.1;
>>cn=0.4;
>>cr=0.2;
>>k=2;
>>T=1;
>>hn=0.1;
>>sn=0.1;
>>hr=0.05;
>>sr=0.05;
>>A=0.1;
>> dn=1-(pn-pr)/(1-o)>> dr=(o*pn-pr)/(o-o^2)>> Fn=20*(qn-dn)>> Fr=20*(qn*e*(a*t+1)-dr)>> zn=qn-dn>> zr=qn*e*(a*t+1)-dr>> eq1='-(pn+hn-cn-sn)*Fn+pn+hn-cn-e*(a*t+1)*(pr+hr+A*(T-t)-cn-sr)*Fr+e*(a*t+1)*(pr+hr+A*(T-t)-cn)+u1*(1-e*(a*t+1))=0'>> eq2='qn+(hn-hr)/(1-o)+(1+qn-(2*pn-pr-sn-cn+hn)/(1-o))*Fn+Fr*(pr+hr+A*(T-t)-cn-sr)/(1-o)+10*zn^2=0'>> eq3='e*qn*(a*t+1)+(hr-o*hn)/(o-o^2)+Fn*(pn+hn-cn-sn)/(1-o)+Fr*(o*pn-2*pr+cn+sr-A*(T-t)-hr-o*(1-o)*e*qn*(a*t+1))/(o-o^2)+10*0.05^2-10*zr^2=0'>> eq4='a*e*qn*(hr+pr+A*(T-2*t-1/a)-cn)+a*e*qn*(sr+cn-hr-A*(T-2*t-1/a)-pr)*Fr-A*(20*dr*zr+10*zr^2)-u1*qn*e*a=0'>> eq5='qn*(a*t+1)*(pr+hr+A*(T-t)-cn)-qn*(a*t+1)*(pr+hr+A*(T-t)-cn-sr)*Fr-2*k*e-u1*a*t*qn=0'>> eq6='qn*(1-e*(a*t+1))=0'>> [qn,pn,pr,t,e,u1]=solve(eq1,eq2,eq3,eq4,eq5,eq6)
算不出来,有多少钱全打赏了.
显示错误是In solve at 140
if isempty(R) | findstr(R,'fsolve') warning('symbolic:solve:warnmsg3','Explicit solution could not be found.'); varargout = cell(1,nargout); varargout{1} = sym([]);
求个非线性方程组,六个参数,急.用MATLAB 求不出来哈,代码如下>> syms qn dn dr pn pr cn cr sn sr hn hr t e a o A T u1 k>> syms Fn Fr fn fr>> syms zn zr>>o=0.8;>>a=0.1;>>cn=0.4;>>cr=0.2;>>k=2;>>T=1;>>hn=0.1;>>sn=0.1;>>hr=0.05;>>sr=0.
我这边完全可以求解啊,把你的语句做成了一个文件,就是附件里面的那个 .m 文件.
把结果做成了一个word文档,也见附件.