homework5_1.m

ar模型 bt算法 ls rls lms music等数字信号处理的源代码

clear
clc
N=25;
n=1:N;
w0=input('please inputw0(between 0 and 1.4*pi)');
M=input('please input M(tne order of AR model)');
u=exp(j*pi*n-j*pi)+exp(j*w0*n-j*0.79*pi)+sqrt(0.1)*randn(1,N);
bfH=zeros(1,N-M);
for n=1:N-M
    bfH(n)=u(n+M);
end
AfH=zeros(M,N-M);
for m=1:M
    for n=1:N-M
        AfH(m,n)=u(n-m+M);
    end
end
bbH=zeros(1,N-M);
for n=1:N-M
    bbH(n)=u(n);
end
AbH=zeros(M,N-M);
for m=1:M
    for n=1:N-M
        AbH(m,n)=u(m+n);
    end
end
bH=[bfH conj(bbH)];
AH=[AfH conj(AbH)];
[U,S,V]=svd(AH');

n1=rank(S);
w=zeros(M,1)
for n=1:n1
    w=w+V(:,n)*(((V(:,n))'*AH*bH'))/(S(n,n))^2;
end
theta=linspace(0,2*pi,1001);
a=-conj(w);
P0=ones(1,1001);
for n=1:M
    P0=P0+a(n)*exp(-j*n*theta);
end
PAR=1./(abs(P0)).^2;
plot(theta,10*log10(PAR))
xlabel('w')
ylabel('PAR(dB)')
hold on