| 查看: 258 | 回复: 0 | |||
[交流]
【求助】请教ode45中加入噪声的问题
|
|
大家好！ 我想用ode45解一带噪声的微分方程组，就是在某一分量上架入噪声，也就是将函数的的第二行换成%%%后的语句。但是matalb能够开始计算却算不出结果，程序如下，请大家指点，谢谢！！ function s=output; tstart = 0; % 时间初始值 tstep = 1e-1; % 时间步长 wholetimes = 3e2; % 总的循环次数 options = odeset('RelTol',1e-4,'AbsTol',1e-5); tspan=tstart:tstep:wholetimes*tstep; [T,Y] = ode45(@Lorenz,tspan,[2 3 7],options); s=[T,Y]; function dY=Lorenz(t,Y); a=10; b=8/3; r=28; dY=zeros(3,1) dY=[a*(Y(2)-Y(1)) -Y(1)*Y(3)+r*Y(1)-Y(2) %%%%%%%%%%%%%%%-Y(1)*Y(3)+r*Y(1)-awgn(Y(2),2);%%对Y2加入信噪比为2的噪声 Y(1)*Y(2)-b*Y(3)]; |
回复此楼» 猜你喜欢
医学类期刊求推荐
已经有5人回复
-
生活琐事由它去
已经有4人回复
-
提交了我也来说说感想
已经有12人回复
-
青B发送上会通知了吗
已经有9人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有6人回复
-
论文撤稿了
已经有8人回复
-
化学专业申博
已经有4人回复
-
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！
已经有5人回复
-
26/27申博自荐
已经有9人回复
-
博士申请
已经有3人回复
