 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 877 | 回復(fù): 0 | ||
lijiawei2013金蟲(chóng) (小有名氣)
|
[求助]
龍格庫(kù)塔計(jì)算結(jié)果問(wèn)題
|
|
我在matlab用龍格庫(kù)塔oed45計(jì)算下面微分方程組時(shí)，為什么隨著求解區(qū)間的增大，它的計(jì)算結(jié)果為NaN？ function dx=fun1(t,x) k12=4.85*10^9;afa=pi/9;r1=0.48745;r2=0.1078,T1=2000;k1x=2*10^9;k1y=2*10^9 k2x=2.5*10^9;k2y=2.5*10^9;I1=152.16;I2=0.69;m1=1096.04;m2=101.7 dx=zeros(12,1) dx(1)=x(2) dx(2)=(k12*((x(2)-x(1))*cos(afa)+(x(3)-x(7))*sin(afa)+r1*x(5)-r2*x(11))*cos(afa)-k1x*x(1))/m1 dx(3)=x(4) dx(4)=(-k12*((x(2)-x(1))*cos(afa)+(x(3)-x(7))*sin(afa)+r1*x(5)-r2*x(11))*sin(afa)-k1y*x(3))/m1 dx(5)=x(6) dx(6)=(T1+r1*k12*((x(2)-x(1))*cos(afa)+(x(3)-x(7))*sin(afa)+r1*x(5)-r2*x(11)))/I1 dx(7)=x(8) dx(8)=(-k12*((x(2)-x(1))*cos(afa)+(x(3)-x(7))*sin(afa)+r1*x(5)-r2*x(11))*cos(afa)-k2x*x(7))/m2 dx(9)=x(10) dx(10)=(k12*((x(2)-x(1))*cos(afa)+(x(3)-x(7))*sin(afa)+r1*x(5)-r2*x(11))*sin(afa)-k2y*x(10))/m2 dx(11)=x(12) dx(12)=(-r2*k12*((x(2)-x(1))*cos(afa)+(x(3)-x(7))*sin(afa)+r1*x(5)-r2*x(11)))/I2 end |
回復(fù)此樓» 猜你喜歡
考研調(diào)劑
已經(jīng)有3人回復(fù)
-
289求調(diào)劑
已經(jīng)有5人回復(fù)
-
材料292調(diào)劑
已經(jīng)有12人回復(fù)
-
化學(xué)調(diào)劑
已經(jīng)有4人回復(fù)
-
292求調(diào)劑
已經(jīng)有14人回復(fù)
-
322求調(diào)劑
已經(jīng)有4人回復(fù)
-
26申博
已經(jīng)有3人回復(fù)
-
307求調(diào)劑
已經(jīng)有9人回復(fù)
-
283求調(diào)劑（080500）
已經(jīng)有4人回復(fù)
-
287求調(diào)劑
已經(jīng)有10人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
關(guān)于用1stpro擬合數(shù)值的問(wèn)題。
已經(jīng)有4人回復(fù)
- 動(dòng)力學(xué)方程組參數(shù)求解問(wèn)題
已經(jīng)有9人回復(fù)
- 求助��！對(duì)于調(diào)用龍格庫(kù)塔法如何對(duì)結(jié)果判斷反饋
已經(jīng)有6人回復(fù)
- Matlab中ode45出現(xiàn)大誤差，關(guān)于用ODE45解偏微分方程的疑問(wèn)，有待解決，大家共同討論！
已經(jīng)有4人回復(fù)
- 編程有點(diǎn)問(wèn)題，求助
已經(jīng)有5人回復(fù)
- 請(qǐng)問(wèn)什么是分段龍格庫(kù)塔的方法？需要用matlab軟件采用這種方法求二階微分方程的解
已經(jīng)有3人回復(fù)
- 求助：matlab解17個(gè)未知數(shù)，17個(gè)方程。
已經(jīng)有4人回復(fù)
- 微分方程組求參數(shù)問(wèn)題，求高人指教，感謝
已經(jīng)有11人回復(fù)
- 龍格-庫(kù)塔（runger-kuta）求解微分方程
已經(jīng)有3人回復(fù)
- Matlab用龍格庫(kù)塔法求二階非線(xiàn)性方程
已經(jīng)有13人回復(fù)
- 龍格庫(kù)塔法 matlab
已經(jīng)有6人回復(fù)
- 4階龍格庫(kù)塔方法疑問(wèn)
已經(jīng)有6人回復(fù)
- 如何求解描述振動(dòng)的二階微分方程
已經(jīng)有21人回復(fù)
- matlab-常微分方程參數(shù)估計(jì)
已經(jīng)有12人回復(fù)
- LABVIEW四階龍格庫(kù)塔法求解一階微分方程組
已經(jīng)有12人回復(fù)
- 【求助】Matlab中利用四階龍格-庫(kù)塔法求解微分方程！�。�！
已經(jīng)有9人回復(fù)
- 【求助】請(qǐng)教二階微分方程的特解（含初始條件）
已經(jīng)有8人回復(fù)
- 【求助】四階龍格-庫(kù)塔方法！�。�！
已經(jīng)有17人回復(fù)
- 【求助】這個(gè)怪異公式是如何推導(dǎo)出來(lái)的？
已經(jīng)有9人回復(fù)
- 【求助】使用Matlab預(yù)估動(dòng)力學(xué)方程問(wèn)題
已經(jīng)有13人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲(chóng)開(kāi)始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 lijiawei2013 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|---|---|---|---|---|
|
[考研] 287求調(diào)劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
|
[考研] 085602 化工專(zhuān)碩 338分 求調(diào)劑 +6 | 路癡小琪 2026-03-27 | 6/300 |
|
|
|
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +4 | 0857求調(diào)劑. 2026-03-24 | 4/200 |
|
|
|
[考研] 一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 丹青奶蓋 2026-03-26 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +17 | 碩星赴 2026-03-23 | 17/850 |
|
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(qǐng)(333分) +9 | holy J 2026-03-21 | 9/450 |
|
|
|
[考研] 0856求調(diào)劑 +8 | zhn03 2026-03-25 | 9/450 |
|
|
|
[考研] 291 求調(diào)劑 +7 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 8/400 |
|
|
|
[考研] 299求調(diào)劑 +4 | 15188958825 2026-03-25 | 4/200 |
|
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
|
[考研] 318求調(diào)劑 +5 | plum李子 2026-03-21 | 8/400 |
|
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
|
[考研] 一志愿南大，0703化學(xué)，分?jǐn)?shù)336，求調(diào)劑 +3 | 收到VS 2026-03-21 | 3/150 |
|