gaussnoise.m

来自「基于AR模型的现代谱估计」· M 代码 · 共 10 行

function y=gaussnoise(m,q,N);

%GAUSSNOISE(m,q,N) generate the N length gause distribution noise,
%with the mean m and the std q.

randn('state', sum(100*clock));  %initialize the state of randn
y=randn(N,1);
y=y/std(y);  %normalize to unit std value.
y=y-mean(y);  %shift to 0 mean value
y=m+q*y;   %generate (m,q) gauss noise with length N