作业帮 > 综合 > 作业

求个非线性方程组,六个参数,急.

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/10/01 10:43:28
求个非线性方程组,六个参数,急.
用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([]);
我这边完全可以求解啊,把你的语句做成了一个文件,就是附件里面的那个  .m 文件.
把结果做成了一个word文档,也见附件.
再问: ��磬���ϱ��Ѿ�����ֵ�ˣ�������˵Ӧ���������������ľ�����ֵ��ɡ� ��word�ļ���Ľ�Ӧ����û����ģ����һ���������ļ�������ȥ��ʾ����һ������ �ܲ��ܰ�ֵ�����ң�֮��Ὣ���ͷ��ϡ� �ղ�������һ�£����ֲ���ֵ�Ļ����ǿ�����������Ǹ�����Ǹ�ֵ���㣬��û�н���ˣ�����Ϊ��ֵ�����޽��