% question three point two
function v=whitenoise(p,N)
p=0.0322;
N=2048;
a=sqrt(p);
v=rand(1,N);
v=v/std(v);
v=v-mean(v);
v=a*v;
subplot(3,1,1);
plot(v);
title('白噪聲v（n）  n取2048')
ylabel('幅度')
xlabel('n')
%產(chǎn)生在第二組數(shù)據(jù)下的u（n）序列
a1=-1.5955;
a2=0.95;
u1=0;
u2=0;
u=[u1,u2];
for n=3:1:2048
    u(n)=-a1*u(n-1)-a2*u(n-2)+v(n);
end
subplot(3,1,2);
plot(u);
title('AR模型下輸出u（n）  n取2048 a1=-1.5955，a2=0.95')
ylabel('幅度')
xlabel('n')
%產(chǎn)生在第二組數(shù)據(jù)下估計(jì)的r（n）序列
r=xcorr(u,1024,'coeff');%auto-corrlation function
k=0:1023;
r=r(1025:2048);
subplot(3,1,3)
plot(k,r);
r(1)
r(2)
r(3)
den=[r(1)^2-r(2)^2];%估計(jì)ak
ex1=(-r(2)*r(1)+r(2)*r(3))/den
ex2=-(r(1)*r(3)-r(2)^2)/den
title('AR模型下估計(jì)的r（n）  n取1024 a1=-1.5955，a2=0.95')
ylabel('幅度')
xlabel('n')