版塊導(dǎo)航: 正在加載中...

登錄注冊

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

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

prince1000

銀蟲 (小有名氣)

應(yīng)助: 2 (幼兒園)
金幣: 35.3
紅花: 1
帖子: 137
在線: 141.8小時
蟲號: 836061
注冊: 2009-08-28
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 再次懇求高手幫忙ode45解微分方程組的問題

function dy=diff(r,y)
dy=zeros(3,1)
dy(1)=-((371.04*cos(123.68*r)/(20+sin(123.68*r))+1/r+140/(20+sin(123.68*r))^2/r)/(1+140/(20+sin(123.68*r))^2))*y(1)-((137.2848*cos(123.68*r)/(20+sin(123.68*r))/r-1/r^2)/(1+140/(20+sin(123.68*r))^2))*y(2)+(140/(20+sin(123.68*r))^2/r^2/(1+140/(20+sin(123.68*r))^2))*y(2)
dy(2)=y(1)

用ode45（'diff',[r0,r],[y1,y2]）命令解以上一階常微分方程組，這個解區(qū)間[r0,r]應(yīng)該為[0 100],可是零是要出錯的，因為方程中有1/r,那這個r0怎么設(shè)定呢？是不是取一個很小的數(shù)，如1.D-50，可是我試了用[1.D-50 100]和[1.D-40 100]，算出來的y（1）差好幾個數(shù)量級，也就是說y（1）強烈的依賴于r0的取值，這是什么原因呢？解區(qū)間的開始的值會導(dǎo)致發(fā)散的情況，這個解區(qū)間的開始值應(yīng)怎么取呢？

回復(fù)此樓

» 猜你喜歡

328求調(diào)劑已經(jīng)有7人回復(fù)
340求調(diào)劑已經(jīng)有5人回復(fù)
一志愿南昌大學(xué)324求調(diào)劑已經(jīng)有4人回復(fù)
求調(diào)劑推薦材料 304 已經(jīng)有15人回復(fù)
291求調(diào)劑已經(jīng)有19人回復(fù)
331環(huán)境科學(xué)與工程求調(diào)劑已經(jīng)有3人回復(fù)
289求調(diào)劑已經(jīng)有8人回復(fù)
315分求調(diào)劑已經(jīng)有7人回復(fù)
295求調(diào)劑已經(jīng)有6人回復(fù)
291求調(diào)劑已經(jīng)有7人回復(fù)

1樓 2013-06-25 21:02:37

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

智能機器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

求高手解答Matlab求解微分方程的初始值和邊界值問題已經(jīng)有5人回復(fù)
一階常微分方程組，matlab進(jìn)行數(shù)值求解已經(jīng)有3人回復(fù)
如何用matlab求解矩陣系數(shù)的二階微分方程已經(jīng)有13人回復(fù)
matlab解微分方程已經(jīng)有10人回復(fù)
lsqnonlin函數(shù)擬合微分方程組參數(shù)擬合問題已經(jīng)有10人回復(fù)
matlab解微分方程組已經(jīng)有15人回復(fù)
一個微分方程組，求解已經(jīng)有18人回復(fù)
求助利用matlab的ODE45求解微分方程已經(jīng)有16人回復(fù)
【緊急求助】帶有微分平方項的ODE求解已經(jīng)有11人回復(fù)
用MATLAB ode45求解2階微分方程已經(jīng)有4人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題已經(jīng)有21人回復(fù)
【求助】急求各位大俠，關(guān)于用ode45解微分方程已經(jīng)有8人回復(fù)
【求助】請教非齊次常微分方程組的解析解法已經(jīng)有4人回復(fù)
【求助】ode45數(shù)值求解常微分方程總是錯，找不出原因已經(jīng)有3人回復(fù)
【求助完畢】用ode逆向求解微分方程問題已經(jīng)有4人回復(fù)
【求助】積分微分方程matlab求解已經(jīng)有6人回復(fù)
【求助】常微分方程組求解中系數(shù)與某變量值關(guān)聯(lián)的問題已經(jīng)有12人回復(fù)
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復(fù)
【求助】解常微分方程組已經(jīng)有7人回復(fù)
【求助】一個時變微分方程組的matlab實現(xiàn) 已經(jīng)有8人回復(fù)

點擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 prince1000 的主題更新