24小时热门版块排行榜

>论坛更新日志 (4617)
>虫友互识 (537)
>文献求助 (439)
>休闲灌水 (401)
>导师招生 (399)
>考博 (136)
>硕博家园 (133)
>博后之家 (104)
>基金申请 (101)
>招聘信息布告栏 (100)
>论文投稿 (95)
>考研 (88)
>教师之家 (65)
>找工作 (50)
>公派出国 (48)
>绿色求助(高悬赏) (46)

返回列表

waixxi

金虫 (初入文坛)

应助: 0 (幼儿园)
金币: 953.5
帖子: 16
在线: 14.9小时
虫号: 576940
注册: 2008-06-21
专业: 纤维

[求助] 请各位大神看看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.

请各位大神看看原因，谢谢了。

回复此楼