| 5 | 1/1 | 返回列表 |
| 查看: 1681 | 回復(fù): 2 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
waixxi金蟲 (初入文壇)
|
[求助]
請(qǐng)各位大神看看MATLAB運(yùn)算報(bào)錯(cuò)原因
|
|
|
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. 請(qǐng)各位大神看看原因,謝謝了。 |
銀蟲 (初入文壇)
至尊木蟲 (著名寫手)
驃騎將軍

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 083000學(xué)碩274求調(diào)劑 +5 | Li李魚 2026-03-26 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 320分,材料與化工專業(yè),求調(diào)劑 +5 | 一定上岸aaa 2026-03-27 | 7/350 |
|
|
[考研] 295材料工程專碩求調(diào)劑 +5 | 1428151015 2026-03-27 | 5/250 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 315分求調(diào)劑 +6 | 26考研上岸版26 2026-03-26 | 6/300 |
|
|
[考研] 0703化學(xué) +6 | 妮妮ninicgb 2026-03-27 | 6/300 |
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +8 | 路癡小琪 2026-03-27 | 8/400 |
|
|
[考研] 328求調(diào)劑 +6 | 嗯滴的基本都 2026-03-27 | 6/300 |
|
|
[考博] 26申博 +3 | 加油沖啊! 2026-03-26 | 3/150 |
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
[考研] 343求調(diào)劑 +4 | 贈(zèng)我一本書 2026-03-23 | 4/200 |
|
|
[考研] 349求調(diào)劑 +4 | 李木子啊哈哈 2026-03-25 | 4/200 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工(085600)296求調(diào)劑 +9 | 稻妻小編 2026-03-26 | 9/450 |
|
|
[考研] 299求調(diào)劑 +4 | 15188958825 2026-03-25 | 4/200 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 各位老師您好:本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
[考研] 070300,一志愿北航320求調(diào)劑 +3 | Jerry0216 2026-03-22 | 5/250 |
|