| 查看: 227 | 回复: 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人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有10人回复
-
版面费该交吗
已经有12人回复
-
面上可以超过30页吧？
已经有5人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有13人回复
-
“人文社科而论，许多学术研究还没有达到民国时期的水平”
已经有5人回复
-
什么是人一生最重要的？
已经有4人回复
