| 3 | 1/1 | 返回列表 |
| 查看: 1682 | 回復(fù): 2 | |||
waixxi金蟲(chóng) (初入文壇)
|
[求助]
請(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)各位大神看看原因,謝謝了。 |
至尊木蟲(chóng) (著名寫(xiě)手)
驃騎將軍

銀蟲(chóng) (初入文壇)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 086000生物與醫(yī)藥調(diào)劑 +3 | Feisty。 2026-03-28 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 食品工程專碩一志愿中海洋309求調(diào)劑 +4 | 小張zxy張 2026-03-26 | 8/400 |
|
|
[考研] 壓國(guó)家一區(qū)線,求導(dǎo)師收留,有恩必謝! +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
[考研] 0856,材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么! 2026-03-26 | 7/350 |
|
|
[考博] 26申博 +3 | 加油沖啊! 2026-03-26 | 3/150 |
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
[考研] 292求調(diào)劑 +4 | 求求了收下我吧?/a> 2026-03-26 | 4/200 |
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
[考研] 297求調(diào)劑 +6 | 田洪有 2026-03-26 | 6/300 |
|
|
[考研] 081700 調(diào)劑 267分 +11 | 迷人的哈哈 2026-03-23 | 11/550 |
|
|
[考研] 環(huán)境專碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 材料/農(nóng)業(yè)專業(yè),07/08開(kāi)頭均可,過(guò)線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問(wèn)津 2026-03-23 | 5/250 |
|