 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (808)
- >考研 (136)
- >導(dǎo)師招生 (103)
- >蟲(chóng)友互識(shí) (52)
- >招聘信息布告欄 (12)
- >考博 (12)
- >碩博家園 (11)
- >文獻(xiàn)求助 (8)
- >基金申請(qǐng) (5)
- >教師之家 (3)
- >論文投稿 (3)
- >休閑灌水 (3)
- >職場(chǎng)人生 (2)
- >博后之家 (2)
- >論文道賀祈福 (2)
- >找工作 (2)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 732 | 回復(fù): 1 | ||
zhuang1017鐵蟲(chóng) (小有名氣)
|
[求助]
matlab 編程求助 已有1人參與
|
|
各位前輩，matlab編程求助，關(guān)于碰撞振動(dòng)的。 g(1)=0.02; g(2)=0.02; u=0.02; odefunhdb = @(t,y)[y(2); (-2*g(1)*y(2)-y(1)+2*u*g(2)*(y(4)-y(2))+u*(y(3)-y(1))+cos(t));... y(4); (-2*g(2)*(y(4)-y(2))-(y(3)-y(1)))];%把動(dòng)力學(xué)方程中的二階微分換成一階 h = 0.05; t = (0:h:100)'; y = zeros(length(t),4); y(1,: )=[0,0,0,0];%參數(shù)設(shè)置與初值 for n = 1:length(t)-1 k1 = odefunhdb(t(n),y(n,: )'); k2 = odefunhdb(t(n)+h/2,y(n,: )'+h/2*k1); k3 = odefunhdb(t(n)+h/2,y(n,: )'+h/2*k2); k4 = odefunhdb(t(n)+h,y(n,: )'+h*k3); y(n+1,: ) = y(n,: )+h/6*(k1+2*k2+2*k3+k4)'; end plot(t,y(:,1)) 這是沒(méi)有發(fā)生碰撞時(shí)的動(dòng)力學(xué)方程編程，現(xiàn)需要加入一個(gè)發(fā)生碰撞的條件，當(dāng)在循環(huán)過(guò)程中滿(mǎn)足y(3)-y(1)=0.8+1e-4時(shí)，有 y(2)+=0.9765*y(2)-+0.235*y(4)- y(4)+=1.1765*y(2)--0.1765*y(4)- 其中+表示變化后的，-表示變化之前，然后用此刻的y(1),y(2)+,y(3),y(4)+作為初值繼續(xù)在龍格庫(kù)塔計(jì)算中循環(huán)，最后需要得到所有的y(1) |
回復(fù)此樓» 猜你喜歡
計(jì)算與實(shí)驗(yàn)科學(xué)工程國(guó)際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有4人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有234人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院
已經(jīng)有28人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復(fù)
-
歡迎調(diào)劑到贛南師范大學(xué)智能制造與未來(lái)能源學(xué)院
已經(jīng)有26人回復(fù)
-
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
UJN物理學(xué)專(zhuān)業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有0人回復(fù)
chendequan
鐵蟲(chóng) (小有名氣)
- 應(yīng)助: 152 (高中生)
- 金幣: 918.6
- 紅花: 17
- 帖子: 262
- 在線: 208.6小時(shí)
- 蟲(chóng)號(hào): 994474
- 注冊(cè): 2010-04-11
- 專(zhuān)業(yè): 結(jié)構(gòu)工程
2樓2017-06-17 15:09:36
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zhuang1017 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 273求調(diào)劑 +19 | 李芷新1 2026-03-31 | 19/950 |
|
|---|---|---|---|---|
|
[考研] 環(huán)境工程297分求調(diào)劑一志愿杭高院 +10 | GENJIOW 2026-03-31 | 11/550 |
|
|
|
[考研] 085600 一志愿9 總分351 求調(diào)劑學(xué)校 +7 | czhcz 2026-03-31 | 9/450 |
|
|
|
[考研] 生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績(jī)338 +8 | YYYYYNNNNN 2026-03-26 | 10/500 |
|
|
|
[碩博家園] 考研調(diào)劑 +5 | 駱駝男人 2026-04-01 | 5/250 |
|
|
|
[考研] 336求調(diào)劑 +3 | kiyy 2026-04-01 | 3/150 |
|
|
|
[考研] 085601英二數(shù)二求調(diào)劑 總分325 +4 | 余航航 2026-03-31 | 4/200 |
|
|
|
[考研] 322求調(diào)劑：一志愿湖南大學(xué) 材料與化工（085600），已過(guò)六級(jí)。 +10 | XX小鄧 2026-03-29 | 10/500 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +13 | 辰沐5211314 2026-03-26 | 13/650 |
|
|
|
[考研] 266分，求材料冶金能源化工等調(diào)劑 +8 | 哇呼哼呼哼 2026-03-27 | 10/500 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +16 | kikiki7 2026-03-30 | 16/800 |
|
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +12 | Li李魚(yú) 2026-03-26 | 12/600 |
|
|
|
[考研] 抱歉 +3 | 田洪有 2026-03-30 | 3/150 |
|
|
|
[考研] 342求調(diào)劑 +4 | 加油a李zs 2026-03-26 | 4/200 |
|
|
|
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4 | 幸運(yùn)的三木 2026-03-29 | 5/250 |
|
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
|
[考研] 求調(diào)劑 +7 | 爭(zhēng)取九點(diǎn)睡 2026-03-28 | 8/400 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
|
[考研] 330一志愿中國(guó)海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
|
[考研] 275求調(diào)劑 +10 | jjjjjjjjjjl 2026-03-27 | 10/500 |
|