# how to add pink noise to force

17 views (last 30 days)
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);
lasha dalaqishvili on 11 Jan 2018
function y = pinknoise(N)
% function: y = pinknoise(N) % N - number of samples to be returned in a row vector % y - a row vector of pink (flicker) noise samples
% The function generates a sequence of pink (flicker) noise samples. % In terms of power at a constant bandwidth, pink noise falls off at 3 dB/oct, i.e. 10 dB/dec.
% difine the length of the vector and % ensure that the M is even, this will simplify the processing if rem(N, 2) M = N+1; else M = N; end

Hristo Zhivomirov on 7 Oct 2018
Hi Lasha!
Please, see my code Pink, Red, Blue and Violet Noise Generation with Matlab on Matlab File Exchange. You can use example2.m in order to add some pink noise to your initial signal.
All best,
Hristo Zhivomirow