| 3 | 1/1 | 返回列表 |
| 查看: 662 | 回復(fù): 2 | ||
huangshuang1木蟲 (著名寫手)
|
[求助]
Matlab ode45求解非線性方程 已有1人參與
|
|
function f=mmm(t,y) f=zeros(4,1); f(1)=y(2); f(2)=(50*pi-y(2))/((50*pi-y(2))^2+595.8566)*1.858*10^6-1.52*10^5*(0.0625*y(1)-0.149*y(3))-1.43*10^7*(0.0625*y(1)-0.149*y(3))^2-2.86*10^6*(0.0625*y(1)-0.149*y(3))^3-32.36*(0.0625*y(2)-0.149*y(4)); f(3)=y(4); f(4)=-86.994-2.437*10^4*(0.149*y(3)-0.0625*y(1))-2.29*10^6*(0.149*y(3)-0.0625*y(1))^2-4.58*10^5*(0.149*y(3)-0.0625*y(1))^3-5.185*(0.149*y(4)-0.0625*y(2)); y0=[0,50*pi,0,50*pi]; [t,y]=ode45(@mmm,[0:0.002:20],y0); plot(t,y(:,2),t,y(:,4),':') 用ODE45一直處于busy狀態(tài),然后我換成ODE23S,可是結(jié)果振蕩得很厲害,不收斂。 請問各位,這個該怎么解?? |
木蟲 (著名寫手)
木蟲 (小有名氣)
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +11 | 新時代材料 2026-03-27 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 本科新能源科學(xué)與工程,一志愿華理能動285求調(diào)劑 +3 | AZMK 2026-03-27 | 5/250 |
|
|
[考研] 329求調(diào)劑 +6 | 星野? 2026-03-26 | 6/300 |
|
|
[考研] 0856求調(diào)劑 +11 | zhn03 2026-03-25 | 12/600 |
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|
|
[考研] 調(diào)劑 +4 | 13853210211 2026-03-24 | 4/200 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|