zuoye2.m

用matlab产生自相关系数固定的序列,然后再观察其功率谱密度.

%
%设法产生零均值、相邻信号（相隔距离为1）的自相关系数分别为0.2及0.8的随机序列，并完成如下工作：
%a) 说明产生方法，对两序列进行比较，并验证序列参数的正确性（理论/仿真）
%b) 仿真给出两序列的自相关函数、功率谱密度，并画出图形，说明两者之间的关系。
%
clc;
clear;
close all;

p=0.8;                                   %设定自相关系数
q=1-p*p;                                 %以此来保证s(k)信号方差为1
s=zeros(1,1000);
a=normrnd(0,1,[1,1000]);                 %初始信号为均值为0方差为1
s(1)=a(1);  
b=normrnd(0,q,[1,1000]);                 %躁声均值为0 方差为q
for k=2:1:1000 
    s(k)=p*s(k-1)+b(k-1);
end                                      %得到所需要的序列
figure
stem(s)

for m=1:1000
 j=0;
   for k=m:1000
    j=j+s(k)*s(k-m+1);
end
  r(m)=j/(10000-m+1);
end                                      %求出自相关函数


for i=1:201
if  i<=100
 r1(i)=r(100-i+1);
else
   r1(i)=r(i-100);
end                                     %取部分范围(-100,100)内的值
end
figure
  plot(-100:100,r1(1:201))
  
f=fft(r(1:1000),1000);                  %求出功率谱密度
l=2*pi/1000:2*pi/1000:2*pi;
figure
plot(l,abs(f(1:1000)))