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

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

24小時(shí)熱門(mén)版塊排行榜

返回列表

下一頁(yè)

lylpiao

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 6742
紅花: 1
帖子: 399
在線: 234小時(shí)
蟲(chóng)號(hào): 1758002
注冊(cè): 2012-04-16
性別: MM
專(zhuān)業(yè): 微/納機(jī)械系統(tǒng)

[求助] 非線性二階微分方程組求解，matlab 已有1人參與

大家好，最近在計(jì)算過(guò)程中遇見(jiàn)一個(gè)矩陣求解問(wèn)題，方程組為非線性二階微分方程組，我不會(huì)求解，請(qǐng)大家指點(diǎn)下，問(wèn)題在圖片里面~~謝謝~

)8F_VRMW$Z)KN)A)~EZNYPH.jpg

回復(fù)此樓

把一切交給時(shí)間！

1樓 2014-03-25 15:58:00

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

回帖支持 ( 顯示支持度最高的前 50 名 )

george1990

金蟲(chóng) (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 1155.3
紅花: 2
帖子: 130
在線: 75.3小時(shí)
蟲(chóng)號(hào): 2142478
注冊(cè): 2012-11-22
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lylpiao: 金幣+10, ★有幫助, 謝謝 2014-03-30 11:47:01

這種問(wèn)題你可以參考機(jī)械動(dòng)力學(xué)里面的4階龍格庫(kù)塔解法它對(duì)于解這種方程很ok的

贊一下(1人)

回復(fù)此樓

天天多學(xué)一點(diǎn)

5樓2014-03-30 01:11:47

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

george1990

金蟲(chóng) (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 1155.3
紅花: 2
帖子: 130
在線: 75.3小時(shí)
蟲(chóng)號(hào): 2142478
注冊(cè): 2012-11-22
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)

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

★ ★ ★ ★ ★
lylpiao: 金幣+5, 哈，那我明天也去試一下，有結(jié)果再告訴你 2014-04-01 23:11:10

這個(gè)我當(dāng)時(shí)試了的消不消去都沒(méi)什么問(wèn)題兩種結(jié)果是一樣的至于多方程的情況我沒(méi)試過(guò)

贊一下(1人)

回復(fù)此樓

天天多學(xué)一點(diǎn)

11樓2014-04-01 17:37:47

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

普通回帖

WWJ1991

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 34 (小學(xué)生)
金幣: 2417.5
散金: 66
紅花: 10
帖子: 337
在線: 195.5小時(shí)
蟲(chóng)號(hào): 2423760
注冊(cè): 2013-04-18
專(zhuān)業(yè): 機(jī)構(gòu)學(xué)與機(jī)器人

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

感謝參與，應(yīng)助指數(shù) +1

參考：http://zhidao.baidu.com/link?url ... le8kHjgTQAnHcUx-EMa。
百度知道給的是一階的求法，你換成二階就行了�？梢詫⒚恳浑A微分的分量寫(xiě)在一起，組成矩陣形式。

贊一下

回復(fù)此樓

2樓2014-03-25 16:48:07

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

daiyouming

新蟲(chóng) (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 1.5
帖子: 3
在線: 3.2小時(shí)
蟲(chóng)號(hào): 2716773
注冊(cè): 2013-10-11

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

感謝參與，應(yīng)助指數(shù) +1

引用回帖:

2樓: Originally posted by WWJ1991 at 2014-03-25 16:48:07
參考：http://zhidao.baidu.com/link?url=n-6qzdYoFHcIxmVsFUFbKrknOkbLdZstCBx1IgbFeDhTH3rdU8gqbP1xQ5hPe-FAZ6Wle8kHjgTQAnHcUx-EMa。
百度知道給的是一階的求法，你換成二階就行了。可以將每一階微分的分量寫(xiě)在 ...

您好，請(qǐng)問(wèn)你的意思是變成24X24矩陣嗎,然后一個(gè)個(gè)一起求解，有這種參考案例嗎:

贊一下

回復(fù)此樓

3樓2014-03-25 18:39:05

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

WWJ1991

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 34 (小學(xué)生)
金幣: 2417.5
散金: 66
紅花: 10
帖子: 337
在線: 195.5小時(shí)
蟲(chóng)號(hào): 2423760
注冊(cè): 2013-04-18
專(zhuān)業(yè): 機(jī)構(gòu)學(xué)與機(jī)器人

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

★ ★ ★ ★ ★
lylpiao: 金幣+5, ★有幫助, 謝謝 2014-03-30 11:47:20

引用回帖:

3樓: Originally posted by daiyouming at 2014-03-25 18:39:05
您好，請(qǐng)問(wèn)你的意思是變成24X24矩陣嗎,然后一個(gè)個(gè)一起求解，有這種參考案例嗎:

:...

不是，還是12X12的。二階微分用d2y表示啊。

贊一下

回復(fù)此樓

4樓2014-03-25 20:44:17

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

lylpiao

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 6742
紅花: 1
帖子: 399
在線: 234小時(shí)
蟲(chóng)號(hào): 1758002
注冊(cè): 2012-04-16
性別: MM
專(zhuān)業(yè): 微/納機(jī)械系統(tǒng)

ben_ladeng: 求助不能申請(qǐng)程序強(qiáng)帖，應(yīng)助和精彩講座可以申請(qǐng)程序強(qiáng)帖。歡迎發(fā)表高質(zhì)量作品 2014-03-31 06:57:26

引用回帖:

5樓: Originally posted by george1990 at 2014-03-30 01:11:47
這種問(wèn)題你可以參考機(jī)械動(dòng)力學(xué)里面的4階龍格庫(kù)塔解法它對(duì)于解這種方程很ok的

請(qǐng)問(wèn)下，該如何求解啊，送你一朵鮮花哈

H{W}]M5NY@[C48~~W@42B}6.jpg

贊一下

回復(fù)此樓

把一切交給時(shí)間！

6樓2014-03-30 11:43:44

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

george1990

金蟲(chóng) (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 1155.3
紅花: 2
帖子: 130
在線: 75.3小時(shí)
蟲(chóng)號(hào): 2142478
注冊(cè): 2012-11-22
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)

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

function dy=lyl(t,y)
dy=zeros(4,1);
dy(1)=y(2);
dy(3)=y(4);
dy(2)=(30*(300*sin(t)-(0.5*y(2)+10*y(4)+1000*y(1)+600*y(4).*y(3)+300*y(2).*y(3))))/2999;
dy(4)=(100*cos(t)-(dy(2)+y(4)+10*y(2)+2*y(3)+3*y(2).*y(3)+20*y(1).*y(3)))/30;

close all;clc;
y0=[0 0 0 0]; %x初值和x導(dǎo)數(shù)的初值
tspan=[0 3];
options = odeset('RelTol',1e-8,'AbsTol',1e-2);
[t y]=ode23s(@lyl,tspan,y0,options);
%plot(t,y(:,2),'r');
我這里設(shè)置的初值是隨意的導(dǎo)致計(jì)算到2.7s就不能計(jì)算了你看看你那邊的初值是多少因?yàn)闆](méi)給這個(gè) 我是隨便設(shè)置的

贊一下

回復(fù)此樓

天天多學(xué)一點(diǎn)

7樓2014-03-31 18:46:58

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

george1990

金蟲(chóng) (小有名氣)

應(yīng)助: 11 (小學(xué)生)
金幣: 1155.3
紅花: 2
帖子: 130
在線: 75.3小時(shí)
蟲(chóng)號(hào): 2142478
注冊(cè): 2012-11-22
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lylpiao: 金幣+5, ★有幫助, 謝謝~ 2014-04-01 10:37:09
lylpiao: 金幣+5, ★有幫助, 謝謝~ 2014-04-01 10:40:34

額第一段程序發(fā)錯(cuò)了應(yīng)該是這樣的：
function dy=lyl(t,y)
dy=zeros(4,1);
dy(1)=y(2);
dy(3)=y(4);
dy(2)=(300*sin(t)-(dy(4)+0.5*y(2)+10*y(4)+1000*y(1)+600*y(4).*y(3)+300*y(2).*y(3)))/100;
dy(4)=(100*cos(t)-(dy(2)+y(4)+10*y(2)+2*y(3)+3*y(2).*y(3)+20*y(1).*y(3)))/30;

贊一下

回復(fù)此樓

天天多學(xué)一點(diǎn)

8樓2014-03-31 18:53:11

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

lylpiao

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 6742
紅花: 1
帖子: 399
在線: 234小時(shí)
蟲(chóng)號(hào): 1758002
注冊(cè): 2012-04-16
性別: MM
專(zhuān)業(yè): 微/納機(jī)械系統(tǒng)

引用回帖:

8樓: Originally posted by george1990 at 2014-03-31 18:53:11
額第一段程序發(fā)錯(cuò)了應(yīng)該是這樣的：
function dy=lyl(t,y)
dy=zeros(4,1);
dy(1)=y(2);
dy(3)=y(4);
dy(2)=(300*sin(t)-(dy(4)+0.5*y(2)+10*y(4)+1000*y(1)+600*y(4).*y(3)+300*y(2).*y(3)))/100;
dy(4)=(10 ...

你好，我有個(gè)疑惑，按照你的思路，一個(gè)方程里面同時(shí)含有二階量，dy(2)和dy(4)，這可以嗎
dy(2)=(300*sin(t)-(dy(4)+0.5*y(2)+10*y(4)+1000*y(1)+600*y(4).*y(3)+300*y(2).*y(3)))/100;
dy(4)=(100*cos(t)-(dy(2)+y(4)+10*y(2)+2*y(3)+3*y(2).*y(3)+20*y(1).*y(3)))/30;

贊一下

回復(fù)此樓

把一切交給時(shí)間！

9樓2014-04-01 10:32:53

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

lylpiao

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 6742
紅花: 1
帖子: 399
在線: 234小時(shí)
蟲(chóng)號(hào): 1758002
注冊(cè): 2012-04-16
性別: MM
專(zhuān)業(yè): 微/納機(jī)械系統(tǒng)

引用回帖:

9樓: Originally posted by lylpiao at 2014-04-01 10:32:53
你好，我有個(gè)疑惑，按照你的思路，一個(gè)方程里面同時(shí)含有二階量，dy(2)和dy(4)，這可以嗎
dy(2)=(300*sin(t)-(dy(4)+0.5*y(2)+10*y(4)+1000*y(1)+600*y(4).*y(3)+300*y(2).*y(3)))/100;
dy(4)=(100*cos(t)-(dy(2) ...

查閱相關(guān)資料：若兩個(gè)高階微分方程同時(shí)含有最高價(jià)導(dǎo)數(shù)項(xiàng)，需要先進(jìn)行相應(yīng)處理，先消去某個(gè)高階導(dǎo)數(shù)。但是方程組少的話計(jì)算還行，多的話消去就太麻煩了。所以還是不知道求解。不知道我理解的對(duì)不對(duì)呢

贊一下

回復(fù)此樓

把一切交給時(shí)間！

10樓2014-04-01 11:13:33

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

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

返回列表

下一頁(yè)

普通表情龍兔虎貓

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

[考研] 278求調(diào)劑 +6	煙火先于春 2026-03-17	6/300	2026-03-21 01:57 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3	研究生何瑤明 2026-03-18	3/150	2026-03-21 01:30 by JourneyLucky
[考研] 324分 085600材料化工求調(diào)劑 +4	llllkkkhh 2026-03-18	4/200	2026-03-21 01:24 by JourneyLucky
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +5	守候夕陽(yáng)CF 2026-03-18	5/250	2026-03-21 01:06 by JourneyLucky
[考研] 354求調(diào)劑 +5	Tyoumou 2026-03-18	8/400	2026-03-21 00:35 by JourneyLucky
[考研] 287求調(diào)劑 +7	晨昏線與星海 2026-03-19	8/400	2026-03-20 22:19 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4	晨昏線與星海 2026-03-19	4/200	2026-03-20 22:15 by JourneyLucky
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3	niko- 2026-03-19	6/300	2026-03-20 21:58 by luoyongfeng
[考研] 一志愿西南交通專(zhuān)碩材料355 本科雙非求調(diào)劑 +5	西南交通專(zhuān)材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 一志愿武理材料工程348求調(diào)劑 +3	￣^￣゜汗 2026-03-19	4/200	2026-03-20 21:01 by zhukairuo
[考研] 289求調(diào)劑 +6	懷瑾握瑜l 2026-03-20	6/300	2026-03-20 20:30 by 學(xué)員8dgXkO
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3	牛奶很忙 2026-03-20	4/200	2026-03-20 19:34 by JourneyLucky
[考研] 286分人工智能專(zhuān)業(yè)請(qǐng)求調(diào)劑愿意跨考！ +3	lemonzzn 2026-03-17	4/200	2026-03-20 11:04 by lemonzzn
[考研] 一志愿西安交通大學(xué)材料工程專(zhuān)業(yè) 282分求調(diào)劑 +5	楓橋ZL 2026-03-18	7/350	2026-03-19 14:52 by 功夫瘋狂
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 材料專(zhuān)碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 275求調(diào)劑 +4	太陽(yáng)花天天開(kāi)心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂