24小时热门版块排行榜    

查看: 1750  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 9/450 2026-08-11 22:12 by 医学老男孩
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
信息提示
请填处理意见