this is the function for the oscillator with the additional force of white noise (awgn). and can I somehow add pink noise to the force? instead of awgn
function dy=x(t,y) global w0 alfa f gama teta b %k1 = wgn(1000,1,100); dy=zeros(2,1);
dy(1)=y(2); %dy(2)=-teta*y(2)+w0^2*(1-alfa*y(1)^2)*y(1)+f*cos(gama*t); %dy(2)=-teta*y(2)+w0*y(1)-alfa*y(1)^3+f*cos(gama*t)+wgn(1000,1,0);
dy(2)=-teta*y(2)+w0^2*(1-alfa*y(1)^2)*y(1)+f*cos(gama*t)+b*awgn(dy(2),0.01);
3 件のコメント
Star Strider (view profile)
このコメントへの直接リンク
https://jp.mathworks.com/matlabcentral/answers/376201-how-to-add-pink-noise-to-force#comment_522728
lasha dalaqishvili (view profile)
このコメントへの直接リンク
https://jp.mathworks.com/matlabcentral/answers/376201-how-to-add-pink-noise-to-force#comment_523076
lasha dalaqishvili (view profile)
このコメントへの直接リンク
https://jp.mathworks.com/matlabcentral/answers/376201-how-to-add-pink-noise-to-force#comment_523838
サインイン to comment.