zzstephen

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 10
帖子: 5
在線: 1.3小時(shí)
蟲號(hào): 1478687
注冊(cè): 2011-11-06
專業(yè): 常微分方程與動(dòng)力系統(tǒng)

[求助] 【緊急求助】帶有微分平方項(xiàng)的ODE求解

dx/(1+(dx)^2) + x = 0.1sin(2t)

想用 Matlab 解，自己弄了很久了

dsolve 解不出來，又不知道怎么轉(zhuǎn)化成ode45 能解的格式，請(qǐng)大家?guī)蛶兔Γ馕鼋饣蛘邤?shù)值解的方法都好

回復(fù)此樓

» 猜你喜歡

0703化學(xué)求調(diào)劑，各位老師看看我�。�！已經(jīng)有3人回復(fù)
339求調(diào)劑已經(jīng)有5人回復(fù)
287求調(diào)劑已經(jīng)有8人回復(fù)
一志愿北京化工大學(xué) 070300 學(xué)碩 336分求調(diào)劑已經(jīng)有8人回復(fù)
085602 化工專碩 338分求調(diào)劑已經(jīng)有5人回復(fù)
0856，材料與化工321分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿北京化工大學(xué)材料與化工 264分各科過A區(qū)國家線已經(jīng)有10人回復(fù)
266求調(diào)劑已經(jīng)有8人回復(fù)
265求調(diào)劑已經(jīng)有3人回復(fù)
286求調(diào)劑已經(jīng)有6人回復(fù)

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

幫忙求解一道微分方程的題目已經(jīng)有7人回復(fù)
LABVIEW四階龍格庫塔法求解一階微分方程組已經(jīng)有12人回復(fù)
用MATLAB ode45求解2階微分方程已經(jīng)有4人回復(fù)
【求助】matlab 求解微分方程中的未知參數(shù) 已經(jīng)有20人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問題已經(jīng)有21人回復(fù)
【求助】用pdepe求解偏微分方程已經(jīng)有10人回復(fù)
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已經(jīng)有22人回復(fù)
【求助完畢】用ode逆向求解微分方程問題已經(jīng)有4人回復(fù)
【求助】積分微分方程matlab求解已經(jīng)有6人回復(fù)
【求助】Python 中求解微分方程 odeint 的用法已經(jīng)有11人回復(fù)
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復(fù)
【求助】求教這種微分方程的詳細(xì)求解過程已經(jīng)有14人回復(fù)

1樓 2011-11-06 12:41:23

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

guojunhang

金蟲 (正式寫手)

應(yīng)助: 7 (幼兒園)
金幣: 3820.6
散金: 32
紅花: 3
帖子: 516
在線: 174.2小時(shí)
蟲號(hào): 1438713
注冊(cè): 2011-10-12
性別: MM
專業(yè): 金屬材料的力學(xué)行為

★
soliton923(金幣+1): 謝謝參與討論~~~ 2011-11-06 18:01:45

現(xiàn)在感覺這個(gè)問題屬于“提法不當(dāng)”，但是又看不出來。
我認(rèn)為你這個(gè)方程的形式要稍微變換一下，再寫出來。
變成一般的微分方程，大家好討論
比如說
f(t,x,dx/dt,....) = 0
的形式。

贊一下(1人)

回復(fù)此樓

2樓2011-11-06 14:40:34

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

yimeng123456

鐵蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 87
帖子: 60
在線: 22.6小時(shí)
蟲號(hào): 1478997
注冊(cè): 2011-11-06
專業(yè): 常微分方程與動(dòng)力系統(tǒng)

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

★
soliton923(金幣+1): 謝謝參與討論~~~ 2011-11-06 18:01:55

換成微分方程組唄令dx=y ,再對(duì)y求導(dǎo)

贊一下(1人)

回復(fù)此樓

3樓2011-11-06 17:37:47

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

meveny

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 317.2
帖子: 40
在線: 23.5小時(shí)
蟲號(hào): 1159404
注冊(cè): 2010-11-30
性別: MM
專業(yè): 常微分方程與動(dòng)力系統(tǒng)

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

我試了下，貌似沒法直接用MATLAB的ode45求解，因?yàn)榍蟛怀鰀x.是否可以先用符號(hào)函數(shù)求解工具（maple?）解二次方程求出dx,再用matlab的ode45做？

贊一下

回復(fù)此樓

4樓2011-11-06 19:48:41

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

dingd

鐵桿木蟲 (職業(yè)作家)

應(yīng)助: 1641 (講師)
金幣: 15037.3
散金: 101
紅花: 234
帖子: 3410
在線: 1223.7小時(shí)
蟲號(hào): 291104
注冊(cè): 2006-10-28

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

zzstephen(金幣+1): 多謝提供解題思路！后來用matlab 的 ode15i 可以解，對(duì)比之后跟你的答案相同。 2011-11-07 00:37:59

有點(diǎn)類似隱式微分方程了，用1stOpt應(yīng)該可以求解，假設(shè)x初值為0.1，t范圍=【0,2*pi】：

CODE:

Variable t=[0,2*pi], x=0.1;

Plot x;

ODEFunction x' = (0.1*sin(2*t)-x)*(1+(x')^2);

t x(t)
0 0.1
0.20943951 0.0850711745007771
0.41887902 0.0801305233227716
0.62831853 0.0813026887426947
0.83775804 0.0846100642268956
1.0471976 0.0864196821522688
1.2566371 0.0839373785627877
1.4660766 0.0755683106928769
1.6755161 0.0610888133870998
1.8849556 0.0416396783046324
2.0943951 0.019519065516734
2.3038346 -0.00226168631502244
2.5132741 -0.0206005920207564
2.7227136 -0.0329220995968774
2.9321531 -0.0376176957124688
3.1415927 -0.0342741951426507
3.3510322 -0.0237389947743921
3.5604717 -0.00803680859667298
3.7699112 0.00991370306354871
3.9793507 0.026808362307764
4.1887902 0.039582546783106
4.3982297 0.0459562014148314
4.6076692 0.0447762626323625
4.8171087 0.0361645432606872
5.0265482 0.0215109045350385
5.2359878 0.00329051787674805
5.4454273 -0.0153456621450357
5.6548668 -0.0311675068120555
5.8643063 -0.0414765431671204
6.0737458 -0.0445541890163028
6.2831853 -0.0398986264011418

贊一下

回復(fù)此樓

5樓2011-11-06 21:43:51

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

zzstephen

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 10
帖子: 5
在線: 1.3小時(shí)
蟲號(hào): 1478687
注冊(cè): 2011-11-06
專業(yè): 常微分方程與動(dòng)力系統(tǒng)

引用回帖:

2樓: Originally posted by guojunhang at 2011-11-06 14:40:34:
現(xiàn)在感覺這個(gè)問題屬于“提法不當(dāng)”，但是又看不出來。
我認(rèn)為你這個(gè)方程的形式要稍微變換一下，再寫出來。
變成一般的微分方程，大家好討論
比如說
f(t,x,dx/dt,....) = 0
的形式。

恩，這個(gè)問題我也困擾了很久，后來發(fā)現(xiàn)這個(gè)是隱式的微分方程要用ode15i 解

贊一下

回復(fù)此樓

6樓2011-11-07 00:39:23

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

zzstephen

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 10
帖子: 5
在線: 1.3小時(shí)
蟲號(hào): 1478687
注冊(cè): 2011-11-06
專業(yè): 常微分方程與動(dòng)力系統(tǒng)

引用回帖:

試過了，可能水平不夠沒弄出來，呵呵。謝謝回復(fù)！

贊一下

回復(fù)此樓

7樓2011-11-07 00:40:53

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

zzstephen

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 10
帖子: 5
在線: 1.3小時(shí)
蟲號(hào): 1478687
注冊(cè): 2011-11-06
專業(yè): 常微分方程與動(dòng)力系統(tǒng)

引用回帖:

4樓: Originally posted by meveny at 2011-11-06 19:48:41:
我試了下，貌似沒法直接用MATLAB的ode45求解，因?yàn)榍蟛怀鰀x.是否可以先用符號(hào)函數(shù)求解工具（maple?）解二次方程求出dx,再用matlab的ode45做？

多謝牛人！根據(jù)你的思路后來發(fā)現(xiàn) matlab ode15i 可以解，答案跟你貼出來的一樣

贊一下

回復(fù)此樓

8樓2011-11-07 00:42:08

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

zzstephen

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 10
帖子: 5
在線: 1.3小時(shí)
蟲號(hào): 1478687
注冊(cè): 2011-11-06
專業(yè): 常微分方程與動(dòng)力系統(tǒng)

引用回帖:

7樓: Originally posted by zzstephen at 2011-11-07 00:40:53:
試過了，可能水平不夠沒弄出來，呵呵。謝謝回復(fù)！

matlab 代碼，遇到同樣問題的筒子可以瞄一眼

f = @(t,x,dx) dx + (x-0.1*sin(2*t))*(1+dx^2);
t0 = 0;
x0 = 0.1;
[x0mod,dx0mod] = decic(f,t0,x0,[],0,[]);
[t,x] = ode15i(f,[0 2*pi],x0mod,dx0mod);
plot(t,x,'o-')

贊一下(1人)

回復(fù)此樓

9樓2011-11-07 00:45:57

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

zhangjob

捐助貴賓 (著名寫手)

應(yīng)助: 52 (初中生)
金幣: 1047.4
散金: 2216
紅花: 67
帖子: 1935
在線: 128.7小時(shí)
蟲號(hào): 1001348
注冊(cè): 2010-04-20
性別: GG
專業(yè): 動(dòng)力學(xué)與控制

我這有ode45包的模版

贊一下

回復(fù)此樓

以前不懂，看貼總是不回，一直沒賺到金幣，也沒成為泰斗精英；現(xiàn)在我明白了，回貼賺金幣，不回白不回，回了也白回，于是，我就把這句話復(fù)制下來，遇貼就回，撈經(jīng)

10樓2012-05-25 13:49:17

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

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

[考研] 一志愿華東理工大學(xué)081700，初試分?jǐn)?shù)271 +6	kotoko_ik 2026-03-23	7/350	2026-03-27 12:29 by 惠州彭于晏
[考研] 322求調(diào)劑 +3	我真的很想學(xué)習(xí) 2026-03-23	3/150	2026-03-27 12:05 by 不吃魚的貓
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +6	曼111 2026-03-24	7/350	2026-03-27 11:58 by 不吃魚的貓
[考研] 一志愿鄭大085600，310分求調(diào)劑 +5	李瀟可 2026-03-26	5/250	2026-03-27 11:14 by 不吃魚的貓
[考研] 0703化學(xué)338求調(diào)劑！ +6	Zuhui0306 2026-03-26	7/350	2026-03-27 10:35 by shangxh
[考研] 安徽大學(xué)專碩生物與醫(yī)藥專業(yè)(086000)324分，英語已過四六級(jí)，六級(jí)521，求調(diào)劑 +3	美味可樂雞翅 2026-03-26	3/150	2026-03-27 08:46 by hxsm
[考研] 286求調(diào)劑 +13	Faune 2026-03-21	13/650	2026-03-26 19:52 by peike
[考研] 資源與環(huán)境調(diào)劑申請(qǐng)(333分) +9	holy J 2026-03-21	9/450	2026-03-26 15:47 by 161765490
[考研] 347求調(diào)劑 +4	L when 2026-03-25	4/200	2026-03-25 13:37 by cocolv
[考研] 0854電子信息求調(diào)劑 324 +4	Promise-jyl 2026-03-23	4/200	2026-03-25 11:36 by Sugarlight
[考研] 311求調(diào)劑 +3	冬十三 2026-03-24	3/150	2026-03-24 21:31 by peike
[考研] 300求調(diào)劑，材料科學(xué)英一數(shù)二 +5	leaflight 2026-03-24	5/250	2026-03-24 16:25 by laoshidan
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7	jhybd 2026-03-23	12/600	2026-03-24 09:03 by jhybd
[考研] 333求調(diào)劑 +3	ALULU4408 2026-03-23	3/150	2026-03-23 19:04 by macy2011
[考研] 生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑 +6	zxts12 2026-03-21	9/450	2026-03-23 18:37 by macy2011
[論文投稿] 急發(fā)核心期刊論文 +3	賢達(dá)問津 2026-03-23	5/250	2026-03-23 17:13 by 妹子不好惹
[考研] 306求調(diào)劑 +5	來好運(yùn)來來來 2026-03-22	5/250	2026-03-22 16:17 by BruceLiu320
[考研] 求調(diào)劑 +4	要好好無聊 2026-03-21	4/200	2026-03-21 18:57 by 學(xué)員8dgXkO
[考研] 297求調(diào)劑 +3	喜歡還是不甘心 2026-03-20	3/150	2026-03-21 18:33 by 學(xué)員8dgXkO
[考研] 330求調(diào)劑0854 +3	assdll 2026-03-21	3/150	2026-03-21 13:01 by 搏擊518