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

返回列表

monitor2885

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

隊(duì)長

應(yīng)助: 6 (幼兒園)
貴賓: 0.339
金幣: 39405.5
散金: 147
紅花: 22
帖子: 4990
在線: 2172小時(shí)
蟲號(hào): 701577
注冊: 2009-02-15
性別: GG
專業(yè): 生態(tài)系統(tǒng)生態(tài)學(xué)

[求助] 隨機(jī)數(shù)如何傳入ode方程

u1=-5:0.1:5; %u1的范圍
p1=normpdf(u1,0.3,1); %均值0.3，方差1
u2=-6:0.1:6; %u2的范圍
p2=normpdf(u2,0.5,2); %均值0.5，方差2
plot(u1,p1);hold on;
plot(u2,p2);hold off;

function testplot
global u1 u2
t=0:1:50;
[t,y]=ode45(@ivpodefun,t,[0 2]);
plot(t,y(:,1));hold on
plot(t,y(:,2));hold off

function dydt=ivpodefun(t,y)
global u1 u2
dydt=zeros(2,1);
dydt(1)=u1*y(2);
dydt(2)=u2*(1-y(1)^2)*y(2)-y(1);

我想在參數(shù) u1，u2 的正態(tài)分布函數(shù)中，用randn各隨機(jī)取10個(gè)數(shù)，然后帶入下面的ode方程求解。做一個(gè)不確定性分析。請(qǐng)問，如何編程，謝謝！

回復(fù)此樓

» 猜你喜歡

288求調(diào)劑已經(jīng)有5人回復(fù)
281求調(diào)劑（0805）已經(jīng)有19人回復(fù)
材料考研調(diào)劑已經(jīng)有3人回復(fù)
346求調(diào)劑[0856] 已經(jīng)有6人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有9人回復(fù)
328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷已經(jīng)有12人回復(fù)
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑已經(jīng)有5人回復(fù)
085601材料工程專碩求調(diào)劑已經(jīng)有7人回復(fù)
321求調(diào)劑已經(jīng)有6人回復(fù)
一志愿中海洋材料工程專碩330分求調(diào)劑已經(jīng)有7人回復(fù)

Retirement

1樓 2011-07-21 06:54:00

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

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

請(qǐng)教matlab中如何繪制該圖？已經(jīng)有7人回復(fù)
fortran 隨機(jī)數(shù) 每次運(yùn)行都是一樣的！已經(jīng)有14人回復(fù)
用MATLAB ode45求解2階微分方程已經(jīng)有4人回復(fù)
聚甲基丙烯酸甲酯如何粉碎用什么機(jī)器已經(jīng)有3人回復(fù)
【求助】急求各位大俠，關(guān)于用ode45解微分方程已經(jīng)有8人回復(fù)
【求助】如何用匯編語言產(chǎn)生隨機(jī)數(shù) 已經(jīng)有5人回復(fù)
【求助完畢】用ode逆向求解微分方程問題已經(jīng)有4人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

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