blp.m

这主要是线性预测算法的一些程序

%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%BLP算法%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%

Xb=Xt([2:M],:);


Wblp=(((((Xb*Xb')/N)^-1)*Xb*Xt(1,:)'/N)').';

t=[M-1:-1:0]';
theta1=[-90:0.1:90];%基阵接受信号方向
atheta1=exp(-j*2*pi*f*t*d*sin(theta1/180*pi)/c);
for i=1:size(theta1,2)
    Pblp(i)=1/((abs(atheta1(:,i)'*[1;-Wblp]))^2);
end

plot(theta1,10*log10(real(Pblp)/max(real(Pblp))));