| 查看: 276 | 回复: 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)]; |
回复此楼» 猜你喜欢
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有12人回复
-
面上提前没消息，有中的吗
已经有14人回复
-
这种情况还有戏吗
已经有12人回复
-
微信指数没变化，科研之友没阅读
已经有19人回复
-
准备明年的基金了
已经有7人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有24人回复
-
2026年国自然面上资助率
已经有21人回复
-
2027年申博
已经有3人回复
-
一个有机合成实验室都需要哪些设备？
已经有7人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
