0.8corrcoef.m

本程序实现随机数的产生和在相同调制方式下功率谱比较

%－－随机产生一信号序列其相邻信号自相关系数为0.2
n=randn(1,200);                          %－－n为高斯白噪声                        
a=zeros(1,200);
a(1,1)=n(1,1);
for i=2:200                              %－－信号序列用高斯白噪声序列表示
    a(1,i)=n(1,i)+0.2*a(1,i-1);
end
c=zeros(1,200);                        
for i=1:199                              %－－信号序列进行移位操作
    c(1,i)=a(1,i+1);
end
c(1,200)=0;
d=corrcoef(a,c);                         %－－验证a,c的互相关系数即相邻信号的互相关系数
f=mean(a);                               %－－验证信号序列的均值


%－－随机产生一信号序列其相邻信号自相关系数为0.8
n1=randn(1,200);                         %－－n1为高斯白噪声
a1=zeros(1,200);
a1(1,1)=n1(1,1);
for i=2:200                              %－－信号序列用高斯白噪声序列表示
    a1(1,i)=n1(1,i)+0.8*a1(1,i-1);
end
c1=zeros(1,200);
for i=1:199                              %－－信号序列进行移位操作
    c1(1,i)=a1(1,i+1);
end
c1(1,200)=0;
d1=corrcoef(a1,c1);                      %－－验证a1,c1的互相关系数即相邻信号的互相关系数
f1=mean(a1);                             %－－验证信号序列的均值

%－－打印所产生的两列信号
b=1:200;
figure(1)
subplot(211)
plot(b,a,'r:.')
grid on
title('series whose corrcoef is 0.2')
subplot(212)
plot(b,a1,'r:.')
grid on
title('series whose corrcoef is 0.8');