% question three point three
function v=whitenoise(p,N)
p=0.038;
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')
%產生在第三組數據下的u（n）序列
a1=-1.9114;
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.9114，a2=0.95')
ylabel('幅度')
xlabel('n')
%產生在第三組數據下估計的r（n）序列
r=xcov(u,1024,'coeff');
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];%估計ak
ex1=(-r(2)*r(1)+r(2)*r(3))/den
ex2=-(r(1)*r(3)-r(2)^2)/den
title('AR模型下估計的r（n）  n取1024 a1=-1.9114，a2=0.95')
ylabel('幅度')
xlabel('n')