| 查看: 304 | 回复: 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)]; |
回复此楼» 猜你喜欢
各位基础老师们，看完中央《从严治党》的决定后，对学术界更加有信心了吧
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有13人回复
