 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 1298 | 回復(fù): 0 | ||
月下行云新蟲(chóng) (小有名氣)
|
[求助]
matlab解DAE遇到的問(wèn)題
|
|
各位同事： 我在解DAE時(shí)，用示例1是可以調(diào)試出結(jié)果的，但是實(shí)例2卻總是 Error using daeic12 (line 77) This DAE appears to be of index greater than 1. Error in ode15s (line 311) [y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,... Error in sunli_Nsolve (line 12) [t,x]=ode15s(@odefun,[0 10],x0,options);% 各位知道為什么嗎？ 示例1程序 function dx = rigid(t,x,k,l,m) dx= zeros(4,1); dx(1)=-0.2*x(1)+k*x(2)*x(3)+l*0.3*x(1)*x(2); dx(2)=2*x(1)*x(2)-5*x(2)*x(3)-2*m*x(2)*x(2); dx(3)=x(1)+x(2)+x(3)-1; dx(4)=-x(4)-1+x(3);%微分方程組 end M=[1 0 0 0;0 1 0 0;0 0 0 0;0 0 0 0;];%質(zhì)量矩陣 options=odeset('mass',M);%對(duì)以DAE問(wèn)題，mass屬性必須設(shè)置 x0=[0.8;0.1;0.1;0.1];%初值 k=1;l=1;m=1; [t,x]=ode15s(@rigid,[0 20],x0,options,k,l,m);%這里好像不能使用ode45 figure('numbertitle','off','name','DAE demo—by Matlabsky') plot(t,x) legend('x1(t)','x2(t)','x3(t)') 實(shí)例2程序： function dx = odefun(t,x) dx=zeros(6,1); l=1.33;k0=19212;k2=0.000560; c0=180;c1=10600;c2=59830;d1=500;d2=3000;tao=17; rb=96;Dw=820;ut=76;Dsw=20;hec=2000; dx(1)=1/c1*((hec-d1)*Dw+(l*x(2)-hec)*Dsw+(d1-l*x(2))*x(6)+k0*rb);%pm dx(2)=1/c2*((hec-d2)*Dw+(l*x(2)-hec)*Dsw+(d2-l*x(2))*x(6)+k0*rb);%hm dx(3)=k2*(x(6)*x(5)-x(5)*hec); dx(4)=x(1)-(x(6)-Dsw)^2/(74554*(x(1)-x(4))^1.3); dx(5)=(l*x(2)*(x(6)-Dsw)+Dsw*hec)/x(6); dx(6)=ut*4877.6*x(4)*(-854+13.47*x(4)+0.4*x(5)-0.003*x(4)*x(5));%微分方程組 end M=[1 0 0 0 0 0;0 1 0 0 0 0;0 0 0 0 0 0;0 0 0 0 0 0;0 0 0 0 0 0;0 0 0 0 0 0];%質(zhì)量矩陣 options=odeset('mass',M);%對(duì)以DAE問(wèn)題，mass屬性必須設(shè)置 x0=[25;2800;950;20;2500;850] [t,x]=ode15s(@odefun,[0 10],x0,options);%這里好像不能使用ode45 figure('numbertitle','off','name','DAE demo—by Matlabsky') plot(t,x) legend('x1(t)','x2(t)','x3(t)') 非常疑惑，基本一樣，就加了幾個(gè)參數(shù) |
回復(fù)此樓» 猜你喜歡
求助
已經(jīng)有4人回復(fù)
-
0805材料320求調(diào)劑
已經(jīng)有3人回復(fù)
-
297求調(diào)劑
已經(jīng)有11人回復(fù)
-
268求調(diào)劑
已經(jīng)有9人回復(fù)
-
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑
已經(jīng)有10人回復(fù)
-
278求調(diào)劑
已經(jīng)有7人回復(fù)
-
26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑
已經(jīng)有5人回復(fù)
-
280求調(diào)劑
已經(jīng)有9人回復(fù)
-
279求調(diào)劑
已經(jīng)有5人回復(fù)
-
286求調(diào)劑
已經(jīng)有5人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
求各位大神解決一個(gè)有關(guān)催化的問(wèn)題
已經(jīng)有5人回復(fù)
- 求大神 幫忙設(shè)計(jì)一下下面這個(gè)東西怎樣合成
已經(jīng)有8人回復(fù)
- matlab中遇到的畫(huà)圖小問(wèn)題，程序沒(méi)問(wèn)題。標(biāo)題要長(zhǎng)@@@@@@@@@@@@@@@@@@@@@@
已經(jīng)有8人回復(fù)
- matlab編程中遇到的問(wèn)題
已經(jīng)有4人回復(fù)
- Matlab: 用bvp解二階微分方程（邊界值）遇到困難，請(qǐng)高手指點(diǎn)
已經(jīng)有5人回復(fù)
- 采用matlab求解微分方程遇到問(wèn)題，請(qǐng)幫忙分析一下錯(cuò)誤-重謝
已經(jīng)有3人回復(fù)
- 請(qǐng)教matlab反應(yīng)動(dòng)力學(xué)參數(shù)估計(jì)遇到的問(wèn)題，謝謝
已經(jīng)有15人回復(fù)
- 最近用matlab采集數(shù)據(jù)遇到一個(gè)問(wèn)題，求助！
已經(jīng)有6人回復(fù)
- Matlab中RBF建模，仿真是遇到的問(wèn)題
已經(jīng)有6人回復(fù)
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問(wèn)題
已經(jīng)有21人回復(fù)
- 【求助】修改MATLAB支持向量機(jī)核函數(shù)程序遇到的疑問(wèn)？求助
已經(jīng)有5人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲(chóng)開(kāi)始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 月下行云 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 278求調(diào)劑 +7 | 煙火先于春 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12 | yangfz 2026-03-17 | 12/600 |
|
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 一志愿中國(guó)石油大學(xué)（華東） 本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-17 | 3/150 |
|
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
|
[考研] 330求調(diào)劑 +4 | 小材化本科 2026-03-18 | 4/200 |
|
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計(jì)算機(jī)專業(yè)調(diào)劑） +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +5 | 陽(yáng)陽(yáng)哇塞 2026-03-14 | 10/500 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學(xué)666 2026-03-16 | 7/350 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
|
[論文投稿] 有沒(méi)有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語(yǔ) 2026-03-15 | 5/250 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|