www.pudn.com > matlabzixiangguan.rar > zuoye2.m

% 
%设法产生零均值、相邻信号（相隔距离为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)))