產生有色噪聲e(k) = x(k) + 0.5*x(k-1)。其中，x(k)為方差為1的白噪聲

clear all; close all;

clc

L=500; %模擬長度

c = [1 -0.5];

nc = length(c) - 1;

xik=zeros(nc,1); %白噪聲初值

xi=randn(L,1); %產生均值為0，方差為1的高斯白噪聲序列

for k=1:L

e(k)=c*[xi(k);xik]; %產生有色噪聲

%資料更新

for i=nc:-1:2

xik(i)=xik(i-1);

end

xik(1)=xi(k);

end

subplot(2,1,1);

plot(xi);

xlabel("k");ylabel("噪聲幅值");title("白噪聲序列");

subplot(2,1,2);

plot(e);

xlabel("k");ylabel("噪聲幅值");title("有色噪聲序列");

%測試功率譜

[y1,f1] = Spectrum_Calc(xi",512);

p1 = 1/L * y1.*conj(y1);

figure(2)

subplot(211)

plot(f1,p1)

[y2,f2] = Spectrum_Calc(e,512);

p2 = 1/L * y2.*conj(y2);

subplot(212)

plot(f2,p2)