白噪声及有色噪声序列的产生

%白噪声及有色噪声序列的产生

设ξ(k) 为均值为0，方差为1的高斯白噪声序列，e(k)为有色噪声序

列：

1

1

1

12

123

()

()()()()

()

10.50.2

()

1 1.50.70.1

C z

e k G z k k

D z

z z

k

z z z

ξξ

ξ

-

-

-

--

---

==

++

=

-++

高斯白噪声序列ξ(k)在Matlab中由rand()函数产生，程序如下：clear all; close all;

L=500; %仿真长度

d=[1 -1.5 0.7 0.1]; c=[1 0.5 0.2]; % 分子分母多项式系数

nd=length(d)-1 ;nc=length(c)-1; %阶次

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

ek=zeros(nd,1);

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

for k=1:L

e(k)=-d(2:nd+1)*ek+c*[xi(k);xik]; %产生有色噪声

%数据更新

for i=nd:-1:2

ek(i)=ek(i-1);

end

ek(1)=e(k);

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('有色噪声序列');