%
%設法產生零均值、相鄰信號（相隔距離為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)))