| 查看: 1766 | 回复: 2 | ||
[求助]
请各位大神看看MATLAB运算报错原因
|
|
function hhh %初值（0.000016 0.000002 0.0125 11252 359 0.9 0.595 0.593 0.937 0.003 2.2125 300） DDR=1.6; TS=333; ms=0.01454; us=1069.6; Ps=99763; As=0.000015588; A0=DDR^2*As; MS=0.00018; x2s=0.99965; Cvs=(1-0.810031)*1.6+0.810031*2.16; eq1=sym('MS-MV')； eq2=sym('x2s*MS-x20*M0')； eq6=sym('1/P20'); eq12=sym('0.5*TS*(ms*Cvs+m0*Cv0)+ms*0.5*us^2-m0*0.5*u0^2+(Ps*As*us-P0*A0*u0)-mv*(719.7+1.28*40))/(0.5*(ms*Cvs+m0*Cv0)-mv*2.16)'); [M0,MV,x20,T0]=solve(eq1,eq2,eq6,eq12,'M0,MV,x20,T0'); M0=subs(M0),MV=subs(MV),x20=subs(x20),T0=subs(T0) fsolve(@myfuns,x0) function eqs=myfuns(x) M0=x(1); MV=x(2); m0=x(3); P0=x(4); u0=x(5); x20=x(6); Psat=x(7); P20=x(8); V2=x(9); V1=x(10); Cv0=x(11); T0=x(12); eq1=MS-MV; eq2=x2s*MS-x20*M0; eq3=M0*((1-x20)*65175+x20*58.1); eq4=m0*u0/A0*10; eq5=m0/A0/0.876; eq6=1/P20; eq7=(EXP(16.6513-2940.49/(T0-35.93)))/760; eq8=Psat*V2*EXP(V1*(1-73.53/31147-V1*V2*0.11)); eq9=x20*1.31/(x20*1.31+(1-x20)*0.79); eq10=1-V2; eq11=(1-V2)*1.6+V2*2.16; eq12=(0.5*TS*(ms*Cvs+m0*Cv0)+ms*0.5*us^2-m0*0.5*u0^2+(Ps*As*us-P0*A0*u0)-mv*(719.7+1.28*40))/(0.5*(ms*Cvs+m0*Cv0)-mv*2.16); eqs=[eq1;eq2;eq3;eq4;eq5;eq6;eq7;eq8;eq9;eq10;eq11;eq12]; >> fsolve(@myfuns,[0.000016 0.000002 0.0125 11252 359 0.9 0.595 0.593 0.937 0.003 2.2125 300]) Error using feval Undefined function 'myfuns' for input arguments of type 'double'. Error in fsolve (line 243) fuser = feval(funfcn{3},x,varargin{:}); Caused by: Failure in initial user-supplied objective function evaluation. FSOLVE cannot continue. 请各位大神看看原因，谢谢了。 |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有230人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有185人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab出错。。。求大神解答。。。
已经有3人回复
- 求大神帮我改一下下面的程序！matlab迭代法的 不知道错在哪里！
已经有7人回复
- MATLAB关于lsqcurvefit程序几个错误，求大神指导一下
已经有6人回复
- linux下运行python报错，不知原因
已经有12人回复
libralibra
至尊木虫 (著名写手)
骠骑将军
- 计算强帖: 3
- 应助: 817 (博后)
- 金币: 12914.1
- 红花: 64
- 帖子: 2238
- 在线: 287.3小时
- 虫号: 696514
- 注册: 2009-02-05
- 专业: 计算机软件
2楼2013-03-19 21:08:23
3楼2013-03-20 09:55:58
30