how to add pink noise to force

lasha dalaqishvili さんによって質問されました 2018 年 1 月 9 日

最新アクティビティ Hristo Zhivomirov さんによって回答されました 2018 年 10 月 7 日

13 ビュー (過去 30 日間)

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 2018 年 1 月 9 日

What spectral dharacteristics define ‘pink noise’ in your application?

lasha dalaqishvili 2018 年 1 月 9 日

I don't have one in this function, I coudn't add it. pink noise (flicker noise) there is a program in mathlab. This is pink noise generator https://www.mathworks.com/matlabcentral/fileexchange/42919-pink--red--blue-and-violet-noise-generation-with-matlab-implementation?focused=6636118&tab=function

lasha dalaqishvili 2018 年 1 月 11 日

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

サインイン to comment.

1 件の回答

リンク

回答者: Hristo Zhivomirov 2018 年 10 月 7 日

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

0 件のコメント

サインイン to comment.

サインインしてこの質問に回答します。

現在この質問をフォロー中です

lasha dalaqishvili (view profile)