| 查看: 1692 | 回复: 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. 请各位大神看看原因，谢谢了。 |
回复此楼» 猜你喜欢
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有171人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有12人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有13人回复
-
桂林理工大学物理学专业招收调剂
已经有13人回复
-
VASP 的一组 GPU / CPU 基准测试记录
已经有0人回复
-
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢
已经有7人回复
-
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名
已经有0人回复
-
【原创讨论】从电子约束到物质编辑：一套可迭代的环形磁场科技树
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
5