 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術科研互動平臺 » 計算模擬區(qū) » 計算模擬 » 軟件使用與操作 » 求解二元微分方程組
| 1 | 1/1 | 返回列表 |
| 查看: 354 | 回復: 0 | ||
可樂不提神新蟲 (小有名氣)
|
[求助]
求解二元微分方程組
|
|
請看我寫的主函數： clc; clear; tp=0.1; n=5; ttp=tp/n; tspan=[0:ttp:tp]; y0=[0,4.5,0,4.5]'; [t y]=ode45(@myfunn,tspan,y0); plot(4.5*t/0.075,y(:,4)/4.5); 下面是myfunn子函數的代碼： function dy=myfunn(t,y) dy=zeros(4,1); Hf=0.075; M2=80; k=353.33*353.33*M2; c=0.2*2*(k*M2)^0.5; Mf=M2/50; A=Mf/(2700*Hf); M1=250*Mf; ee=y(1)/Hf; Y=0； if y(1)>=Y p=3.5e6+0.6e6*ee/(1-ee)^2; else p=0 ; end dy(1)=y(2); dy(2)=c/M1*(y(4)-y(2))+k/M1*(y(3)-y(1))-A*p/M1+9.8; dy(3)=y(4); dy(4)=-c/M2*(y(4)-y(2))-k/M2*(y(3)-y(1))+9.8; end if y(1)>=Y Y=y(1); else Y=Y; end 我的問題來了，如下： 在這個微分方程中，p值是隨y(1)的變化而變化的。假如在第t步計算中，y(1)>Y則p=3.5e6+0.6e6*ee/(1-ee)^2;當y(1)<Y時，p=0。然后計算(t+1)步的y(1)值。把(t+1)步的y(1)和Y對比。當y(1)>Y,則用y(1)替換Y。我寫的代碼中每步的Y初始值都是0開始，但是每一次y(1)與Y對比之后，Y的值有可能發(fā)生變化。而且我的代碼中y(1)在每一次調用子函數的過程中都是不變的。 請問怎么解決？ |
回復此樓» 猜你喜歡
中德博士后2026項目5月30日截止- 國內博士后入站、德國兩年全職
已經有0人回復
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內博后入站德國兩年全職
已經有0人回復
-
物理學I論文潤色/翻譯怎么收費?
已經有240人回復
-
【新加坡】納米電子器件項目組有“聯合培養(yǎng)博士生”名額
已經有2人回復
-
重慶交大材料學院26年碩士招生復試實施細則已出。歡迎加入光子學微結構與器件課題組。
已經有0人回復
-
北工大宋曉艷教授團隊招收2026年9月入學博士【學博和專博名額均有】！
已經有1人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收10人
已經有0人回復
-
天津城建大學理學院化學、環(huán)境工程、材料物理等方向接收調劑，10人
已經有0人回復
-
計算與實驗科學工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學換能器分會S8
已經有0人回復
-
0702一志愿吉大B區(qū)求調劑有論文
已經有2人回復
相關版塊跳轉
我要訂閱樓主 可樂不提神 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 279求調劑 +6 | j的立方 2026-03-29 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 生物技術與工程 +7 | 1294608413 2026-03-25 | 8/400 |
|
|
|
[考研] 一志愿哈爾濱工業(yè)大學材料與化工方向336分 +11 | 辰沐5211314 2026-03-26 | 11/550 |
|
|
[考研] 290求調劑 +3 | dfffsar 2026-03-29 | 3/150 |
|
|
|
[考研] 調劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|
|
[考研] 289求調劑 +5 | BrightLL 2026-03-29 | 5/250 |
|
|
|
[考研] 070305高分子化學與物理 304分求調劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
|
[考研] 279求調劑 +4 | 蝶舞輕繞 2026-03-29 | 4/200 |
|
|
|
[考研] 295材料工程專碩求調劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|
|
[考研] 一志愿華北電力大學能動專碩，293，求調劑 +3 | 15537177284 2026-03-23 | 5/250 |
|
|
|
[考研] 265求調劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
|
[考研] 一志愿上海理工能源動力（085800）310分求調劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 安徽大學專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級，六級521，求調劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
|
[考研] 085601 材料工程 313分 求調劑 +5 | Ong3 2026-03-27 | 5/250 |
|
|
|
[考研] 324求調劑 +5 | hanamiko 2026-03-26 | 5/250 |
|
|
|
[考研] 343求調劑 +4 | 贈我一本書 2026-03-23 | 4/200 |
|
|
|
[考研] 290分調劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
|
[考研] 【2026考研調劑】制藥工程 284分 求相關專業(yè)調劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
|
[考研] 一志愿山東大學藥學學碩求調劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|