| 3 | 1/1 | 返回列表 |
| 查看: 1680 | 回復(fù): 2 | ||
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)各位大神看看原因,謝謝了。 |
至尊木蟲 (著名寫手)
驃騎將軍

銀蟲 (初入文壇)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
[考研] 352分 化工與材料 +3 | 海納百川Ly 2026-03-27 | 3/150 |
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
[考博] 26申博 +3 | 加油沖! 2026-03-26 | 3/150 |
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[碩博家園] 北京林業(yè)大學(xué)碩導(dǎo)招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
[考研] 一志愿陜師大生物學(xué)071000,298分,求調(diào)劑 +5 | SYA! 2026-03-23 | 5/250 |
|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績338 +4 | YYYYYNNNNN 2026-03-26 | 4/200 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +4 | 李軼男003 2026-03-20 | 4/200 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
[考研] 284求調(diào)劑 +15 | Zhao anqi 2026-03-22 | 15/750 |
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生(材料、化學(xué)、電化學(xué),環(huán)境) +4 | 我愛學(xué)電池 2026-03-23 | 4/200 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +5 | WW.' 2026-03-23 | 7/350 |
|
|
[考博] 26申博自薦 +3 | whh869393 2026-03-24 | 3/150 |
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|