24小时热门版块排行榜    

查看: 275  |  回复: 0

monitor2885

至尊木虫 (职业作家)

队长

[求助] 随机数如何传入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);

我想在参数 u1,u2 的正态分布函数中,用randn各随机取10个数,然后带入下面的ode方程求解。做一个不确定性分析。请问,如何编程,谢谢!
回复此楼

» 猜你喜欢

Retirement
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 monitor2885 的主题更新
信息提示
请填处理意见