微塵、夢想

木蟲 (知名作家)

程序強帖: 6
應助: 2 (幼兒園)
貴賓: 0.353
金幣: 4757.9
散金: 3089
紅花: 31
沙發(fā): 247
帖子: 8788
在線: 1125小時
蟲號: 1203290
注冊: 2011-02-14
專業(yè): 制造系統(tǒng)與自動化

[求助] 隱式歐拉法求解一階常微分方程

算法：

MATLAB程序：

CODE:

Format long;

N=(b-a)/h;

y=zeros(N+1,1);

y(1)=y0;

x=a:h:b;

var=findsym(f);

for i=2:N+1

   fx=Funval(f,var(1),x(i));

   gx=y(i-1)+h*fx-varvec(2);

   y(i)=NewtonRoot(gx,-10,10,eps);%用牛頓法得出下步的迭代值

end

format short;

其中，F(xiàn)unval是求函數(shù)f在點(var(1),x(i))處的值

PS：顯示法求解常微分方程的算法都用C語言實現(xiàn)，但隱式法不知道怎么實現(xiàn)，主要是不知道如何求后一個點的導數(shù)值，希望高人給指點一下，呵呵……

回復此樓

» 猜你喜歡

被我言中：新模板不強調(diào)格式了，假專家開始管格式了已經(jīng)有3人回復
085600材料與化工求調(diào)劑已經(jīng)有5人回復
材料與化工求調(diào)劑已經(jīng)有5人回復
277調(diào)劑已經(jīng)有4人回復
290求調(diào)劑已經(jīng)有4人回復
材料工程專碩調(diào)劑已經(jīng)有3人回復
26考研求調(diào)劑已經(jīng)有6人回復
求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學。調(diào)劑到任何專業(yè)都可以已經(jīng)有5人回復
311求調(diào)劑已經(jīng)有7人回復
303求調(diào)劑已經(jīng)有3人回復

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

大腸桿菌超聲破碎參數(shù)，求解！��！已經(jīng)有15人回復
畢設，數(shù)學相關專業(yè)，求指導！已經(jīng)有18人回復
奇怪的質粒電泳圖，求解。。已經(jīng)有3人回復
求解MALDI-TOF MS圖譜已經(jīng)有18人回復
【求助】ANSYS或abaqus 等CAE軟件的迭代求解的初值問題已經(jīng)有7人回復
【求助】求解隱函數(shù) 已經(jīng)有6人回復
【求助】求助一個一階常微分方程的初值問題已經(jīng)有16人回復
【求助】一個二階常微分方程的解法（重獎金幣�。� 已經(jīng)有16人回復
【求助】一個常微分方程的通解已經(jīng)有6人回復

任風云變幻，我笑對人生！

1樓 2012-02-23 12:20:33

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

sudo

木蟲 (正式寫手)

程序強帖: 16
應助: 6 (幼兒園)
金幣: 1297.6
散金: 1486
紅花: 20
帖子: 588
在線: 641小時
蟲號: 1211394
注冊: 2011-02-24
性別: GG
專業(yè): 文藝美學

【答案】應助回帖

感謝參與，應助指數(shù) +1
余澤成(金幣+1): 謝謝參與應助！ 2012-02-25 04:34:01
微塵、夢想(金幣+30): ★★★★★最佳答案謝謝，sudo教會了我很多東西呀，呵呵…… 2012-02-27 13:57:13

隱式歐拉法迭代式右邊f(xié)(x[n+1], y[n+1])里面的y[n+1]，其實也只是用估計的方法得出來的

一般可以用歐拉法（或者牛頓法和其他別的什么方法）得到它的估計值，所以這個過程實際上是：

CODE:

yt[n+1] = y[n] + h * f(x[n], y[n]);

y[n+1] = y[n] + h * f(x[n+1], yt[n+1]);

梯形法什么的類似，迭代式中f里面的y[n+1]，無論怎樣都是需要某種方法先估計出來的

贊一下

回復此樓

2樓2012-02-23 13:25:46

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

sudo

木蟲 (正式寫手)

程序強帖: 16
應助: 6 (幼兒園)
金幣: 1297.6
散金: 1486
紅花: 20
帖子: 588
在線: 641小時
蟲號: 1211394
注冊: 2011-02-24
性別: GG
專業(yè): 文藝美學

恩，如果不要求通用的方法，只需要針對特定函數(shù)的話，也可以先化簡迭代式，或者可能可以直接導出y[n+1]的公式什么的~

贊一下(1人)

回復此樓

3樓2012-02-23 13:29:11

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

微塵、夢想

木蟲 (知名作家)

程序強帖: 6
應助: 2 (幼兒園)
貴賓: 0.353
金幣: 4757.9
散金: 3089
紅花: 31
沙發(fā): 247
帖子: 8788
在線: 1125小時
蟲號: 1203290
注冊: 2011-02-14
專業(yè): 制造系統(tǒng)與自動化

引用回帖:

3樓: Originally posted by sudo at 2012-02-23 13:29:11:
恩，如果不要求通用的方法，只需要針對特定函數(shù)的話，也可以先化簡迭代式，或者可能可以直接導出y[n+1]的公式什么的~

那個MATLAB程序就是用牛頓法得到的下一個值，不過，不知道牛頓法的算法是什么呀，我只知道用牛頓法解方程的公式：
x[n+1]=x[n]-f(x[n])/f`(x[n])
這是一種切線法，從一端向一個方向逼近方程的根，可是跟這個聯(lián)系不到一塊呀

贊一下

回復此樓

任風云變幻，我笑對人生！

4樓2012-02-23 13:38:29

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

sudo

木蟲 (正式寫手)

程序強帖: 16
應助: 6 (幼兒園)
金幣: 1297.6
散金: 1486
紅花: 20
帖子: 588
在線: 641小時
蟲號: 1211394
注冊: 2011-02-24
性別: GG
專業(yè): 文藝美學

余澤成(金幣+2): 謝謝參與應助！ 2012-02-25 04:34:19

引用回帖:

4樓: Originally posted by 微塵、夢想 at 2012-02-23 13:38:29:
那個MATLAB程序就是用牛頓法得到的下一個值，不過，不知道牛頓法的算法是什么呀，我只知道用牛頓法解方程的公式：
x[n+1]=x[n]-f(x[n])/f`(x[n])
這是一種切線法，從一端向一個方向逼近方程的根，可是跟這個聯(lián) ...

牛頓法是用來做什么的呢？一般是用來求非線性方程的根的：

對于f(x)=0
迭代式x[n+1] = x[n] - f(x[n])/f`(x[n])
迭代的結果就是f(x)=0的數(shù)值解

=============================================

對于y[n+1] = y[n] + h * f(x[n+1], y[n+1])

CODE:

   fx=Funval(f,var(1),x(i)); %里面的x(i)就相當于x[n+1]，var(1)是一個符號變量（其實是等下要求解的y[n+1]），把它們代入到f中之后，得到一個符號表達式

   gx=y(i-1)+h*fx-varvec(2); %其實這里就是方程g(x)=0左邊的g(x)，就是y[n] + h*fx - y[n+1] = 0這個方程，里面y[n]是上一步知道的值，fx是含有變量y[n+1]的式子，外面的y[n+1]也視為變量

   y(i)=NewtonRoot(gx,-10,10,eps); %這就是用牛頓法來解非線性方程了

也就是說，這段MATLAB程序里面涉及了符號的運算

贊一下

回復此樓

5樓2012-02-23 15:34:57

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

sudo

木蟲 (正式寫手)

程序強帖: 16
應助: 6 (幼兒園)
金幣: 1297.6
散金: 1486
紅花: 20
帖子: 588
在線: 641小時
蟲號: 1211394
注冊: 2011-02-24
性別: GG
專業(yè): 文藝美學

余澤成(金幣+1): 鼓勵交流！ 2012-02-25 04:34:33

如果它的方法要通用的話，在NewtonRoot還必須對gx進行一次符號求導（用diff函數(shù)），所以這個方法其實比較蛋疼

還不如直接用前向歐拉給估計了...

贊一下

回復此樓

6樓2012-02-23 15:41:38

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

微塵、夢想

木蟲 (知名作家)

程序強帖: 6
應助: 2 (幼兒園)
貴賓: 0.353
金幣: 4757.9
散金: 3089
紅花: 31
沙發(fā): 247
帖子: 8788
在線: 1125小時
蟲號: 1203290
注冊: 2011-02-14
專業(yè): 制造系統(tǒng)與自動化

引用回帖:

6樓: Originally posted by sudo at 2012-02-23 15:41:38:
如果它的方法要通用的話，在NewtonRoot還必須對gx進行一次符號求導（用diff函數(shù)），所以這個方法其實比較蛋疼

還不如直接用前向歐拉給估計了...

嗯，我也就是想把它搞明白它到底是如何工作的，解常微分方程的顯示解法，我都搞明白了，從最簡單的歐拉法到四階龍格-庫塔法，不過我在網(wǎng)上找到一篇關于算法分析的文章，結論是隱式算法的穩(wěn)定性普遍比顯式算法高，且只有隱式歐拉方法才是絕對穩(wěn)定的。
也沒仔細看，反正用四階龍格-庫塔法的精度已經(jīng)非常高了，呵呵……
這是文章的地址：
http://wenku.baidu.com/view/8bcb52492e3f5727a4e96207.html

贊一下

回復此樓

任風云變幻，我笑對人生！

7樓2012-02-23 16:50:00

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

sudo

木蟲 (正式寫手)

程序強帖: 16
應助: 6 (幼兒園)
金幣: 1297.6
散金: 1486
紅花: 20
帖子: 588
在線: 641小時
蟲號: 1211394
注冊: 2011-02-24
性別: GG
專業(yè): 文藝美學

余澤成(金幣+1): 鼓勵討論！ 2012-02-25 04:34:54

引用回帖:

7樓: Originally posted by 微塵、夢想 at 2012-02-23 16:50:00:
嗯，我也就是想把它搞明白它到底是如何工作的，解常微分方程的顯示解法，我都搞明白了，從最簡單的歐拉法到四階龍格-庫塔法，不過我在網(wǎng)上找到一篇關于算法分析的文章，結論是隱式算法的穩(wěn)定性普遍比顯式算法高， ...

嗯，隱式歐拉的確是有它的意義的。。。不過你好像誤會了我的意思

我說的是估計y[n+1]那里不用牛頓法而用顯式歐拉法來做會方便點

贊一下

回復此樓

8樓2012-02-23 17:19:55

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

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學。調(diào)劑到任何專業(yè)都可以 +4	小丁想進步 2026-03-11	5/250	2026-03-17 16:05 by 外星文明
[考研] 304求調(diào)劑 +7	小熊joy 2026-03-14	7/350	2026-03-17 08:53 by 霧散后相遇lc
[考研] 化學工程321分求調(diào)劑 +10	大米飯！ 2026-03-15	13/650	2026-03-17 06:57 by laoshidan
[考研] [導師推薦]西南科技大學國防/材料導師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	4/200	2026-03-16 20:13 by Wangjingyue
[考研] 機械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[考研] 085600調(diào)劑 +5	漾漾123sun 2026-03-12	6/300	2026-03-16 15:58 by 漾漾123sun
[考研] 0703 物理化學調(diào)劑 +3	我可以上岸的對?/a> 2026-03-13	5/250	2026-03-16 10:50 by 我可以上岸的對?/a>
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 288求調(diào)劑 +4	奇點0314 2026-03-14	4/200	2026-03-14 23:04 by JourneyLucky
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6	邱gl 2026-03-12	7/350	2026-03-13 23:24 by 邱gl
[考研] 341求調(diào)劑 +3	番茄頭--- 2026-03-10	3/150	2026-03-13 23:07 by JourneyLucky
[考研] 285 求調(diào)劑資源與環(huán)境一志愿北京化工大學 +3	未名考生 2026-03-10	3/150	2026-03-13 23:04 by JourneyLucky
[考研] 308求調(diào)劑 +5	是Lupa啊 2026-03-11	5/250	2026-03-13 22:13 by JourneyLucky
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4	寶貴熱 2026-03-13	4/200	2026-03-13 22:01 by 星空星月
[考研] 【0856】化學工程（085602）313 分，本科學科評估A類院�；瘜W工程與工藝，誠求調(diào)劑 +7	小劉快快上岸 2026-03-11	7/350	2026-03-13 16:06 by ruiyingmiao
[考研] 考研調(diào)劑 +4	芬達46 2026-03-12	4/200	2026-03-13 16:04 by ruiyingmiao
[考研] 277求調(diào)劑 +4	anchor17 2026-03-12	4/200	2026-03-13 11:15 by 白夜悠長
[考研] 279求調(diào)劑 +3	莫xiao 2026-03-10	4/200	2026-03-11 08:06 by 斬魂滴兔子！