| 查看: 265 | 回复: 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)]; |
回复此楼» 猜你喜欢
评委有多少概率知道其他专家手中有哪些人的本子？
已经有4人回复
-
E0414, 我的本子有没有希望？
已经有17人回复
-
青A35岁以下通知答辩了吗
已经有4人回复
-
小城的小雨
已经有3人回复
-
看《给阿ma的情书》有感
已经有5人回复
-
国自然申请五篇代表作大比拼，感觉这个是最重要的
已经有4人回复
-
雷雨
已经有3人回复
