flp.m

来自「这主要是线性预测算法的一些程序」· M 代码 · 共 19 行

%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%FLP算法%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%
for i=1:M-1
    Xf(i,:)=Xt(M-i,:);
end

Wflp=(((((Xf*Xf')/N)^-1)*Xf*Xt(M,:)'/N)').';

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

end

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