版塊導(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í)熱門版塊排行榜

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

返回列表

monitor2885

至尊木蟲 (職業(yè)作家)

應(yīng)助: 6 (幼兒園)
貴賓: 0.339
金幣: 39466
帖子: 4994
在線: 2174.2小時(shí)
蟲號(hào): 701577

[交流] 【求助】ode45方程中的參數(shù)

function testfunction
global u1 u2
u1=[1;2;3;4;5];
u2=[6;7;8;9;10];
for i=1:5
[t,y]=ode45(@fun,[0 20],[1 2]);
z1(i)=y(:,1);
z2(i)=y(:,2);
end

function dy=fun(t,y)
global u1 u2
dy=zeros(2,1);
dy(1)=u1(i)*y(2);
dy(2)=u2(i)*(1-y(1)^2)*y(2)-y(1);

u1和u2是微分方程的參數(shù)，是5*1的數(shù)組。我要利用ode45函數(shù)運(yùn)行5次，即for=1:5，第一次運(yùn)行，參數(shù)u1和u2分別是1和6，即u1(1)和u2(1)；第二次運(yùn)行，參數(shù)u1和u2分別是2和7，即u1(2)和u2(2)……以此類推。每次運(yùn)行的結(jié)果y1和y2都儲(chǔ)存到z1和z2里面。請(qǐng)問(wèn)，上面的程序需要哪些改動(dòng)？謝謝！

回復(fù)此樓

» 猜你喜歡

請(qǐng)問(wèn)共同通訊和共同一作的認(rèn)可度問(wèn)題已經(jīng)有7人回復(fù)
285求調(diào)劑已經(jīng)有7人回復(fù)
一志愿山東大學(xué)，085600，344 已經(jīng)有4人回復(fù)
一志愿北京科技材料科學(xué)與工程288分，求調(diào)劑已經(jīng)有12人回復(fù)
280求調(diào)劑已經(jīng)有7人回復(fù)
315求調(diào)劑已經(jīng)有5人回復(fù)
26考研調(diào)劑已經(jīng)有5人回復(fù)
301求調(diào)劑已經(jīng)有7人回復(fù)
求調(diào)劑22408 288分已經(jīng)有4人回復(fù)
求調(diào)劑已經(jīng)有4人回復(fù)

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

如何求解描述振動(dòng)的二階微分方程已經(jīng)有21人回復(fù)
matlab擬合方程參數(shù)時(shí)初值的選擇已經(jīng)有15人回復(fù)
matlab解微分方程已經(jīng)有10人回復(fù)
matlab 求指點(diǎn) 動(dòng)力學(xué)方程擬合過(guò)程中導(dǎo)數(shù)的獲取已經(jīng)有10人回復(fù)
matlab-常微分方程參數(shù)估計(jì) 已經(jīng)有12人回復(fù)
如何用MATLAB 實(shí)現(xiàn)化學(xué)反應(yīng)方程式（寫程序代碼）? 已經(jīng)有20人回復(fù)
再次問(wèn)如何顯示function里面的變量（附論文和相應(yīng)的程序）已經(jīng)有8人回復(fù)
幫幫忙看看我的難題已經(jīng)有15人回復(fù)
【求助】催化反應(yīng)動(dòng)力學(xué) 已經(jīng)有5人回復(fù)
【求助完畢】用ode逆向求解微分方程問(wèn)題已經(jīng)有4人回復(fù)
【求助】matlab怎么求解偏微分方程組啊，先謝謝了已經(jīng)有13人回復(fù)
【求助】使用Matlab預(yù)估動(dòng)力學(xué)方程問(wèn)題已經(jīng)有13人回復(fù)

» 搶金幣啦！回帖就可以得到:

查看全部散金貼

1樓 2011-02-07 15:00:54

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

xiegangmai

版主 (職業(yè)作家)

仿真EPI: 3
應(yīng)助: 157 (高中生)
貴賓: 5.735
金幣: 40609
帖子: 4046
在線: 1330.2小時(shí)
蟲號(hào): 653607

monitor2885(金幣+5): 運(yùn)行對(duì)了，怎么把每次運(yùn)行的結(jié)果y1和y2都儲(chǔ)存到z1和z2里面？ 2011-02-08 03:32:41

引用回帖:

Originally posted by monitor2885 at 2011-02-07 15:00:54:
function testfunction
global u1 u2
u1=[1;2;3;4;5];
u2=[6;7;8;9;10];
for i=1:5
[t,y]=ode45(@fun,[0 20],[1 2]);
z1(i)=y(:,1);
z2(i)=y(:,2);
end

function dy=fun(t,y)
global u1 u2
dy=zeros ...

用參數(shù)來(lái)實(shí)現(xiàn)吧。

CODE:

function testfunction

u1 = [ 1; 2; 3; 4; 5 ];

u2 = [ 6; 7; 8; 9; 10 ];

for i = 1 : 5

    [ t, y ] = ode45(@( t, y )fun( t, y, u1( i ), u2( i ) ), [ 0 20 ], [ 1 2 ] );

%     z1( i ) = y( :, 1 );

%     z2( i ) = y( :, 2 );

    plot( t, y );

    hold on

end

function dy = fun( t, y, u1, u2 )

dy = zeros( 2, 1 );

dy( 1 ) = u1 * y( 2 );

dy( 2 ) = u2 * ( 1 - y( 1 ) ^ 2 ) * y( 2 ) - y( 1 );

贊一下(1人)

回復(fù)此樓

2樓2011-02-07 18:39:46

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

monitor2885

至尊木蟲 (職業(yè)作家)

應(yīng)助: 6 (幼兒園)
貴賓: 0.339
金幣: 39466
帖子: 4994
在線: 2174.2小時(shí)
蟲號(hào): 701577

引用回帖:

Originally posted by xiegangmai at 2011-02-07 18:39:46:

用參數(shù)來(lái)實(shí)現(xiàn)吧。

[code]
function testfunction

u1 = [ 1; 2; 3; 4; 5 ];
u2 = [ 6; 7; 8; 9; 10 ];
for i = 1 : 5
[ t, y ] = ode45(@( t, y )fun( t, y, u1( i ), u2( i ) ), [ 0 20 ], [ 1 ...

運(yùn)行對(duì)了，怎么把每次運(yùn)行的結(jié)果y1和y2都儲(chǔ)存到z1和z2里面？

贊一下

回復(fù)此樓

3樓2011-02-09 02:38:00

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

xiegangmai

版主 (職業(yè)作家)

仿真EPI: 3
應(yīng)助: 157 (高中生)
貴賓: 5.735
金幣: 40609
帖子: 4046
在線: 1330.2小時(shí)
蟲號(hào): 653607

monitor2885(金幣+5): cell是個(gè)做什么的函數(shù)呢 2011-02-10 02:30:53

引用回帖:

Originally posted by monitor2885 at 2011-02-09 02:38:00:

運(yùn)行對(duì)了，怎么把每次運(yùn)行的結(jié)果y1和y2都儲(chǔ)存到z1和z2里面？

因?yàn)槊看斡?jì)算出的結(jié)果中數(shù)組維度大小不一樣，可以用cell或struct實(shí)現(xiàn)。
例如用cell數(shù)據(jù)：

CODE:

function testfunction

u1 = [ 1; 2; 3; 4; 5 ];

u2 = [ 6; 7; 8; 9; 10 ];

z1 = cell( 5, 1 );

z2 = cell( 5, 1 );

for i = 1 : 5

    [ t, y ] = ode45(@( t, y )fun( t, y, u1( i ), u2( i ) ), [ 0 20 ], [ 1 2 ] );

    z1{ i, 1 } = y( :, 1 );

    z2{ i, 1 } = y( :, 2 );

    plot( t, y );

    hold on

end

function dy = fun( t, y, u1, u2 )

dy = zeros( 2, 1 );

dy( 1 ) = u1 * y( 2 );

dy( 2 ) = u2 * ( 1 - y( 1 ) ^ 2 ) * y( 2 ) - y( 1 );

贊一下

回復(fù)此樓

4樓2011-02-09 09:21:06

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

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

返回列表

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

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

[考研] 285求調(diào)劑 +14	AZMK 2026-04-02	14/700	2026-04-02 15:54 by 上九天攬?jiān)拢ê?/a>
[考研] 初試301，代碼085701環(huán)境工程，本碩一致，四六級(jí)已過(guò)，有二區(qū)一作，共發(fā)表5篇論文 +6	axibli 2026-04-01	6/300	2026-04-02 13:42 by Ecowxq666！
[考研] 化學(xué)308分調(diào)劑 +14	你好明天你好 2026-03-30	15/750	2026-04-02 13:39 by 嘻嘻笑22
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-04-02	4/200	2026-04-02 13:03 by yulian1987
[考研] 一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑 +9	雙馬尾痞老板2 2026-04-01	9/450	2026-04-02 12:14 by oooqiao
[考研] 271求調(diào)劑 +15	勒布朗@ 2026-03-31	20/1000	2026-04-02 11:24 by Sammy2
[考研] 321求調(diào)劑 +9	y-yh 2026-04-01	10/500	2026-04-02 10:14 by 不吃魚的貓
[考研] 各位老師好，我的一志愿為北京科技大學(xué)085601材料專碩 +13	Koxui 2026-03-28	13/650	2026-04-02 09:35 by 噠噠噠呱呱呱
[考研] 一志愿346上海大學(xué)生物學(xué) +3	上海大學(xué)346調(diào)劑 2026-04-01	3/150	2026-04-02 08:36 by w蟲蟲123
[考研] 求調(diào)劑，一志愿南京師范大學(xué)計(jì)算機(jī)專碩，初試373，六級(jí)通過(guò)， +3	計(jì)算機(jī)追夢(mèng)人 2026-04-01	3/150	2026-04-02 07:57 by fxue1114
[考研] 070300化學(xué)279求調(diào)劑 +15	哈哈哈^_^ 2026-03-31	17/850	2026-04-01 21:37 by 給你你注意休息
[考研] 材料專業(yè)調(diào)劑 +5	啦啦啦哭 2026-03-31	6/300	2026-04-01 16:48 by JourneyLucky
[考研] 309求調(diào)劑 +19	誰(shuí)不是少年 2026-03-29	19/950	2026-04-01 15:47 by jp9609
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +7	12.27 2026-03-27	7/350	2026-03-31 21:59 by lbsjt
[考研] 070300化學(xué)354求調(diào)劑 +15	101次希望 2026-03-28	15/750	2026-03-31 17:58 by jp9609
[考研] 286求調(diào)劑 +5	丟掉懶惰 2026-03-27	8/400	2026-03-31 11:27 by Delta2012
[考研] 22408 359分調(diào)劑 +4	Qshers 2026-03-27	8/400	2026-03-31 08:53 by Qshers
[考研] 085701環(huán)境工程求調(diào)劑 +11	多久上課 2026-03-27	12/600	2026-03-30 21:21 by 研究僧導(dǎo)導(dǎo)
[有機(jī)交流] 考研調(diào)劑 +8	watb 2026-03-26	8/400	2026-03-30 18:40 by 544594351
[考研] 342求調(diào)劑 +4	加油a李zs 2026-03-26	4/200	2026-03-30 16:39 by 晶體之美

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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

monitor2885

[交流] 【求助】ode45方程中的參數(shù)

» 猜你喜歡

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

» 搶金幣啦！回帖就可以得到: 查看全部散金貼

xiegangmai

monitor2885

xiegangmai

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

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

» 搶金幣啦！回帖就可以得到:

查看全部散金貼