 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2465 | 回復(fù): 3 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
glgjw銀蟲 (小有名氣)
|
[求助]
MATLAB求解二階常微分方程初值問題 已有1人參與
|
||
|
主程序如下： >> global m1 m2 k1 c1 c2 v0 >> m1=4080;m2=12.5e3;k1=7.35e6;c1=17.56e3;c2=0.225e3;v0=3.2/60 >> [t,y]=ode45(@doublemass,[0 5],[-0.0154;-0.0049;0;0]) 其中的doublemass的函數(shù)M文件定義如下： function dyy=doublemass(t,y) global m1 m2 k1 c1 c2 v0 dyy=[y(3);(-4.48e6*v0*t/(6-y(1))-(k1+4.48e6/(6-y(1)))*y(1)-(c1+c2)*y(2)+c2*y(4)+4.48e6*y(2)/(6-y(1)))/m1;y(4);(4.48e6*v0*t/(6-y(1))-m2*10+c2*y(3)+4.48e6*y(1)/(6-y(1))-c2*y(4)-4.48e6*y(2)/(6-y(1)))/m2] 但計算出的結(jié)果卻不對，前面兩個初值不是前面設(shè)置好的-0.0154；-0.0049，而是 -0.0000 -0.0000，如下圖： y = 1.0e+009 * -0.0000 -0.0000 0 0 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0001 -0.0000 -0.0000 -0.0000 0.0001 -0.0000 -0.0000 -0.0000 0.0001 -0.0000 -0.0000 -0.0000 0.0001 -0.0000 -0.0000 -0.0000 0.0002 -0.0000 -0.0001 -0.0000 0.0002 -0.0000 -0.0001 -0.0000 0.0002 -0.0000 -0.0001 -0.0000 0.0003 -0.0000 -0.0001 -0.0000 0.0003 -0.0000 -0.0001 -0.0000 0.0004 -0.0000 -0.0001 -0.0000 0.0005 -0.0000 -0.0002 -0.0000 0.0005 -0.0000 -0.0002 -0.0000 0.0006 -0.0000 -0.0002 -0.0000 0.0006 -0.0000 -0.0002 -0.0000 0.0007 -0.0000 -0.0002 -0.0000 0.0007 -0.0000 -0.0002 -0.0000 0.0007 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0010 -0.0001 -0.0003 -0.0000 0.0013 -0.0001 -0.0003 -0.0000 0.0017 -0.0001 -0.0004 -0.0000 0.0022 -0.0001 -0.0004 -0.0000 0.0030 -0.0001 -0.0004 -0.0000 0.0041 -0.0001 -0.0004 -0.0000 0.0054 -0.0001 -0.0004 -0.0000 0.0069 -0.0002 -0.0004 -0.0000 0.0094 -0.0002 -0.0004 -0.0001 0.0123 -0.0002 -0.0004 -0.0001 0.0157 -0.0002 -0.0004 -0.0001 0.0197 -0.0003 -0.0004 -0.0001 0.0261 -0.0003 -0.0004 -0.0001 0.0335 -0.0003 -0.0004 -0.0002 0.0420 -0.0003 -0.0004 -0.0002 0.0516 -0.0004 -0.0004 -0.0002 0.0678 -0.0004 -0.0004 -0.0003 0.0864 -0.0005 -0.0004 -0.0004 0.1076 -0.0005 -0.0004 -0.0004 0.1314 -0.0006 -0.0005 -0.0005 0.1585 -0.0007 -0.0005 -0.0006 0.1885 -0.0007 -0.0005 -0.0007 0.2214 -0.0008 -0.0005 -0.0008 0.2575 -0.0008 -0.0005 -0.0009 0.2967 -0.0009 -0.0005 -0.0010 0.3392 -0.0010 -0.0006 -0.0011 0.3851 -0.0010 -0.0006 -0.0013 0.4345 -0.0011 -0.0006 -0.0014 0.4875 -0.0012 -0.0006 -0.0016 0.5441 -0.0013 -0.0006 -0.0017 0.6045 -0.0013 -0.0006 -0.0019 0.6688 -0.0014 -0.0007 -0.0021 0.7372 -0.0015 -0.0007 -0.0023 0.8096 -0.0016 -0.0007 -0.0025 0.8863 -0.0017 -0.0007 -0.0027 0.9673 -0.0018 -0.0007 -0.0029 1.0528 -0.0019 -0.0007 -0.0032 1.1428 -0.0020 -0.0008 -0.0034 1.2376 -0.0021 -0.0008 -0.0037 1.3371 -0.0021 -0.0008 -0.0040 1.4416 -0.0023 -0.0008 -0.0043 1.5512 -0.0024 -0.0008 -0.0046 1.6659 -0.0025 -0.0009 -0.0049 1.7859 -0.0026 -0.0009 -0.0051 1.8795 -0.0027 -0.0009 -0.0054 1.9762 -0.0027 -0.0009 -0.0056 2.0761 -0.0028 -0.0009 -0.0059 2.1792 -0.0029 -0.0009 不知是何原因？ |
回復(fù)此樓» 猜你喜歡
085600材料與化工調(diào)劑
已經(jīng)有14人回復(fù)
-
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑
已經(jīng)有10人回復(fù)
-
材料科學(xué)與工程求調(diào)劑
已經(jīng)有9人回復(fù)
-
一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩
已經(jīng)有3人回復(fù)
-
材料科學(xué)與工程調(diào)劑
已經(jīng)有10人回復(fù)
-
274求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿食品科學(xué)與工程083200求調(diào)劑
已經(jīng)有3人回復(fù)
-
311求調(diào)劑一志愿合肥工業(yè)大學(xué)
已經(jīng)有3人回復(fù)
-
266分，求材料相關(guān)專業(yè)調(diào)劑
已經(jīng)有8人回復(fù)
-
材料求調(diào)劑 一志愿哈工大總分298分，前三科223分
已經(jīng)有10人回復(fù)
|
本帖內(nèi)容被屏蔽 |
3樓2016-02-25 10:00:29
wurongjun
專家顧問 (職業(yè)作家)
-
專家經(jīng)驗: +831 - 數(shù)學(xué)EPI: 9
- 應(yīng)助: 791 (博后)
- 貴賓: 0.308
- 金幣: 24609
- 散金: 310
- 紅花: 75
- 帖子: 3004
- 在線: 881.4小時
- 蟲號: 1368482
- 注冊: 2011-08-14
- 性別: GG
- 專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算
- 管轄: 數(shù)學(xué)
2樓2016-02-24 17:47:24
glgjw
銀蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 252.4
- 帖子: 232
- 在線: 93.2小時
- 蟲號: 2115194
- 注冊: 2012-11-08
- 性別: GG
- 專業(yè): 信息處理方法與技術(shù)
4樓2016-02-25 11:11:56
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 生物考研337分求調(diào)劑 +3 | cgxin 2026-03-30 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 求調(diào)劑，一志愿 南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +9 | @taotao 2026-03-29 | 9/450 |
|
|
[考研] 285求調(diào)劑 +6 | AZMK 2026-03-29 | 9/450 |
|
|
|
[考研] 一志愿鄭大材料工程290求調(diào)劑 +11 | Youth_ 2026-03-30 | 11/550 |
|
|
|
[考研] 26考研-291分-廈門大學(xué)（085601）-柔性電子學(xué)院材料工程專業(yè)求調(diào)劑 +5 | min3 2026-03-24 | 6/300 |
|
|
|
[考研] 310求調(diào)劑 +10 | 爭取九點睡 2026-03-30 | 10/500 |
|
|
|
[考研] 329求調(diào)劑 +8 | 星野? 2026-03-26 | 8/400 |
|
|
|
[考研] 070300一志愿211，312分求調(diào)劑院校 +4 | 小黃鴨寶 2026-03-30 | 4/200 |
|
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +10 | 景隨風(fēng) 2026-03-29 | 15/750 |
|
|
|
[考研] 求化學(xué)調(diào)劑 +11 | wulanna 2026-03-28 | 11/550 |
|
|
|
[考研] 294分080500材料科學(xué)與工程求調(diào)劑 +8 | 柳溪邊 2026-03-26 | 8/400 |
|
|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
|
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
|
[考研] 305求調(diào)劑 +8 | RuiFairyrui 2026-03-28 | 8/400 |
|
|
|
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調(diào)劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 300分，材料，求調(diào)劑，英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
|
[考研]
 材料考研調(diào)劑生
+3
|
黃粱一夢千年 2026-03-24 | 3/150 |
|
|
|
[考研] 292求調(diào)劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|