24小时热门版块排行榜    

查看: 1754  |  回复: 2

waixxi

金虫 (初入文坛)

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

请各位大神看看原因,谢谢了。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

感谢参与,应助指数 +1
ben_ladeng: 专家考核 2013-03-20 16:55:24
Undefined function 'myfuns' for input arguments of type 'double'.

这不是很明确吗?没有定义myfuns函数
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
2楼2013-03-19 21:08:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sun198607

银虫 (初入文坛)

myfuns函数的M文件错误
3楼2013-03-20 09:55:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 waixxi 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 重要来源:本周末出结果 +6 瞬息宇宙 2026-08-12 6/300 2026-08-13 08:51 by zhengxianzb
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 关于Filecode分析方法 +8 majunge000 2026-08-10 9/450 2026-08-12 22:17 by 虫友是什么虫
[基金申请] 不应该看fileCode +5 且听虎啸 2026-08-12 6/300 2026-08-12 16:26 by Tide man
[基金申请] 长年满屏的广告,版主太不责任了。基金也等的急 +4 gltch 2026-08-12 4/200 2026-08-12 16:08 by eulota
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 10/500 2026-08-12 09:12 by 布布和一二
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
信息提示
请填处理意见