1164185188

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 636.4
散金: 10
帖子: 32
在線: 19.8小時
蟲號: 3571957
注冊: 2014-12-02
性別: GG
專業(yè): 微/納機械系統(tǒng)

[求助] 請問如何用MATLAB求解四階微分方程并做出未知函數(shù)圖像已有1人參與

$y''''=1/(1-y)^2+1/(1-y)^4$

回復(fù)此樓

» 猜你喜歡

326求調(diào)劑已經(jīng)有8人回復(fù)
333求調(diào)劑已經(jīng)有7人回復(fù)
一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑已經(jīng)有3人回復(fù)
考研化學(xué)學(xué)碩調(diào)劑，一志愿985 已經(jīng)有7人回復(fù)
化學(xué)工程321分求調(diào)劑已經(jīng)有21人回復(fù)
0703化學(xué)調(diào)劑已經(jīng)有15人回復(fù)
0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷已經(jīng)有14人回復(fù)
326求調(diào)劑已經(jīng)有4人回復(fù)
311求調(diào)劑已經(jīng)有13人回復(fù)
311求調(diào)劑已經(jīng)有7人回復(fù)

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

MATLA與CAD 已經(jīng)有0人回復(fù)

1樓 2016-05-15 20:29:27

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

回帖置頂 ( 共有1個 )

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
1164185188: 金幣+40, ★有幫助 2016-05-16 07:49:23

引用回帖:

4樓: Originally posted by 1164185188 at 2016-05-15 20:54:05
我是這么編程的，出現(xiàn)這種狀況，請問各位是什么問題？
dsolve('D4y=1/(1-y)^2+1/(1-y)^4','y(0)=0','y(1)=0','Dy(0)=0','Dy(1)=0','x')

ans =

月只藍

Dsolve是求解析解的函數(shù)，不適用。原方程具有非線性，大部分情況下，是沒有解析解的。只能求數(shù)值解。
原方程四階的，通過變量替換，可降階成含有4個一階方程的常微分方程組，然后用MATLAB ode函數(shù)解之。

二階降階的例子，參見：http://www.gaoyang168.com/bbs/viewthread.php?tid=8979085&target=self&page=1
四階的做類似處理即可。

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

1164185188

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

5樓2016-05-15 21:52:49

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

普通回帖

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
1164185188: 金幣+10, ★有幫助 2016-05-15 20:36:52

再補上4個初值條件，可以求得數(shù)值解。
MATLAB的ode 系列函數(shù)。

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

1164185188

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

2樓2016-05-15 20:31:00

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

1164185188

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 636.4
散金: 10
帖子: 32
在線: 19.8小時
蟲號: 3571957
注冊: 2014-12-02
性別: GG
專業(yè): 微/納機械系統(tǒng)

送紅花一朵

引用回帖:

2樓: Originally posted by 月只藍 at 2016-05-15 20:31:00
再補上4個初值條件，可以求得數(shù)值解。
MATLAB的ode 系列函數(shù)。

你好，如果加上這四個初值條件，請問該如何編程呢？
f(0)=0, f(1)=0,f'(0)=0,f’（1）=0

贊一下

回復(fù)此樓

3樓2016-05-15 20:39:54

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

1164185188

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 636.4
散金: 10
帖子: 32
在線: 19.8小時
蟲號: 3571957
注冊: 2014-12-02
性別: GG
專業(yè): 微/納機械系統(tǒng)

我是這么編程的，出現(xiàn)這種狀況，請問各位是什么問題？
dsolve('D4y=1/(1-y)^2+1/(1-y)^4','y(0)=0','y(1)=0','Dy(0)=0','Dy(1)=0','x')

ans =

[ empty sym ] @月只藍

贊一下

回復(fù)此樓

4樓2016-05-15 20:54:05

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

1164185188

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 636.4
散金: 10
帖子: 32
在線: 19.8小時
蟲號: 3571957
注冊: 2014-12-02
性別: GG
專業(yè): 微/納機械系統(tǒng)

引用回帖:

5樓: Originally posted by 月只藍 at 2016-05-15 21:52:49
Dsolve是求解析解的函數(shù)，不適用。原方程具有非線性，大部分情況下，是沒有解析解的。只能求數(shù)值解。
原方程四階的，通過變量替換，可降階成含有4個一階方程的常微分方程組，然后用MATLAB ode函數(shù)解之。

二階降 ...

非常感謝！

回復(fù)此樓

6樓2016-05-16 07:50:09

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

1164185188

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 636.4
散金: 10
帖子: 32
在線: 19.8小時
蟲號: 3571957
注冊: 2014-12-02
性別: GG
專業(yè): 微/納機械系統(tǒng)

引用回帖:

再次感謝，已經(jīng)做出來了！

贊一下

回復(fù)此樓

7樓2016-05-16 09:03:51

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

1164185188

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 636.4
散金: 10
帖子: 32
在線: 19.8小時
蟲號: 3571957
注冊: 2014-12-02
性別: GG
專業(yè): 微/納機械系統(tǒng)

送紅花一朵

引用回帖:

你好，麻煩再問一下，我的初值條件是這四個，'y(0)=0','y(1)=0','Dy(0)=0','Dy(1)=0'，請問用ode45編程時，這個邊界條件如何編程，是這樣嗎？怎么老感覺不對呢？
clc,clear,close all
tspan=0:0.1:1;
y0=[0,0,0,0];
[x, y]=ode45(@myodes,tspan,y0);

贊一下

回復(fù)此樓

8樓2016-05-16 14:18:42

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

月只藍

主管區(qū)長 (職業(yè)作家)

專家經(jīng)驗: +1059
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線: 3291.4小時
蟲號: 1122189
注冊: 2010-10-14
專業(yè): 宇宙學(xué)
管轄: 計算模擬區(qū)

引用回帖:

8樓: Originally posted by 1164185188 at 2016-05-16 14:18:42
你好，麻煩再問一下，我的初值條件是這四個，'y(0)=0','y(1)=0','Dy(0)=0','Dy(1)=0'，請問用ode45編程時，這個邊界條件如何編程，是這樣嗎？怎么老感覺不對呢？
clc,clear,close all
tspan=0:0.1:1;
y0=;
=od ...

'y(0)=0','y(1)=0','Dy(0)=0','Dy(1)=0' 這種不是初值問題了，而是邊值問題。y(0)=0，y(1)=0，這算一個邊值條件，針對y的；Dy(0)=0，Dy(1)=0，這是第二個邊值條件，針對y'的。原方程高達4階，而邊值條件僅僅給到了一階，需要補充上二階和三階的邊值條件，否則無法定解。
5樓說的降階的方法，是適用的，但代碼是針對初值問題來設(shè)計的，對于邊值問題，ode函數(shù)不適用，而需要用 bvp4c 函數(shù)。
可在MATLAB主界面中doc bvp4c ，查看相關(guān)代碼例子。

贊一下

回復(fù)此樓

MATLAB、MS小問題、普通問題請發(fā)帖求助！時間精力有限，恕不接受無償私信求助。

9樓2016-05-16 14:44:26

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

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

[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3	邱gl 2026-03-20	7/350	2026-03-21 19:05 by 15709483992
[考研] 材料工程專碩 348分求調(diào)劑 +3	冬辭. 2026-03-17	5/250	2026-03-21 18:47 by 學(xué)員8dgXkO
[考研] 材料 271求調(diào)劑 +5	展信悅_ 2026-03-21	5/250	2026-03-21 17:29 by 學(xué)員8dgXkO
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +3	沿岸有貝殼6 2026-03-21	3/150	2026-03-21 10:35 by 暮云清寒
[考研] 求調(diào)劑 +6	Mqqqqqq 2026-03-19	6/300	2026-03-21 08:04 by JourneyLucky
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6	77chaselx 2026-03-20	6/300	2026-03-21 07:24 by JourneyLucky
[考研] 一志愿山大07化學(xué) 332分四六級已過本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-16	3/150	2026-03-21 03:59 by JourneyLucky
[考研] 301求調(diào)劑 +10	yy要上岸呀 2026-03-17	10/500	2026-03-21 03:14 by JourneyLucky
[考研] 354求調(diào)劑 +5	Tyoumou 2026-03-18	8/400	2026-03-21 00:35 by JourneyLucky
[考研] 317求調(diào)劑 +5	申子申申 2026-03-19	9/450	2026-03-20 22:26 by JourneyLucky
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 289求調(diào)劑 +6	懷瑾握瑜l 2026-03-20	6/300	2026-03-20 20:30 by 學(xué)員8dgXkO
[考研] 353求調(diào)劑 +3	拉鉤不許變 2026-03-20	3/150	2026-03-20 19:56 by JourneyLucky
[考研] 工科材料085601 279求調(diào)劑 +7	困于星晨 2026-03-17	9/450	2026-03-20 17:38 by 無懈可擊111
[考研] 招收調(diào)劑碩士 +4	lidianxing 2026-03-19	12/600	2026-03-20 12:25 by lidianxing
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82