版塊導航: 正在加載中...

登錄注冊

應《網絡安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網跟帖服務。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學院2026年研究生招生接收調劑公告

返回列表

CQU_wjq

鐵桿木蟲 (著名寫手)

應助: 20 (小學生)
金幣: 5137
紅花: 1
帖子: 1634
在線: 221.5小時
蟲號: 2333035
注冊: 2013-03-09
性別: GG
專業(yè): 金屬材料表面科學與工程

[求助] Matlab解微分方程組已有1人參與

想用Matlab解一組微分方程組，具體如下：
①先編輯了一個.m文件：
function dy=eq1(t,y)
dy=zeros(2,1);
dy(1)=-4/((1-y(2))^2)*((0.8+0.1*t^2)/(0.8+0.2*t^2)+Q/(sqrt(3))/(1-y(2))/(1-y(2)));
dy(2)=0.4*(y(1)-1)*(-0.86*y(2)^2+0.345*y(2)+0.515)/(0.16+0.04*t^2)/(y(2)-(1-y(1)/(2.571-y(2)-y(2)*log(1-y(2)))));
②邊界條件：
t=1時，y（1）=0，y（2）=0.5；
③關于方程中的Q：
Q為一不確定常數(shù)，滿足條件：t=0，y（1）=0；

如何采用龍格-庫塔法解上述方程呢，先求出Q？求指導。

回復此樓

» 猜你喜歡

» 本主題相關價值貼推薦，對您同樣有幫助:

matlab 解非線性偏微分方程組已經有6人回復
非線性二階微分方程組求解，matlab 已經有12人回復
利用matlab求一階變系數(shù)微分方程組的數(shù)值解。已經有7人回復
matlab 非線性微分方程求解已經有3人回復
matlab如何求解一個非線性微分方程組已經有8人回復
求助，用matlab解微分方程組，希望高手能給予指點已經有17人回復
如何用Matlab求時變微分方程的解？已經有5人回復
matlab解偏微分方程求助已經有12人回復
請教用matlab求解一個非線性偏微分方程組的數(shù)值解已經有8人回復
matlab中有關微分方程求解的幾個小問題已經有4人回復
matlab有限差分法求解復雜的偏微分方程組已經有19人回復
求助MATLAB解方程組已經有6人回復
matlab解微分方程已經有10人回復
Matlab求解偏微分方程組已經有22人回復
matlab數(shù)值求解邊界條件微分方程組已經有7人回復
高金求助matlab解微分方程組已經有12人回復
matlab解微分方程組已經有15人回復
matlab解方程組已經有8人回復
【求助】matlab 求解微分方程中的未知參數(shù) 已經有20人回復
【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題已經有21人回復
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經有13人回復
【求助】常微分方程組的解析解已經有3人回復
【求助】用MATLAB 求微分方程已經有18人回復

逆風奔跑。

1樓 2014-12-25 20:31:39

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

CelestialCYJ

木蟲 (小有名氣)

應助: 33 (小學生)
金幣: 2482.6
紅花: 12
帖子: 144
在線: 43.3小時
蟲號: 3277506
注冊: 2014-06-16
專業(yè): 結構工程

【答案】應助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
CQU_wjq: 金幣+10, ★★★★★最佳答案 2015-01-19 16:19:57

tspan=[1 0],y0=[0 0.5],eq1要改成eq1(t,y,flag,Q),然后用ode45函數(shù)，Q通過不斷的迭代，最終找到滿足y1（t=0）=0的Q，就可以把這組微分方程解出來了，至于Q怎么尋優(yōu)，你可以用各種方法，比如MATLAB自帶的fminsearch函數(shù)，等等

贊一下

回復此樓

2樓2014-12-27 22:53:01

已閱回復此樓關注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關版塊跳轉我要訂閱樓主 CQU_wjq 的主題更新