fig5_26.m

来自「《最优阵列处理》一书第五章的MATLAB例程」· M 代码 · 共 23 行

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 5.26
% Wavenumber spectra for complex AR(2) process
% z_1 = 0.9exp(-j0.5pi), z_2 = 0.9exp(-j0.7pi)
% K. Bell 10/13/00
% updated by K. Bell 7/24/01, 10/23/01
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all
close all

psi = [-1:0.001:1]*pi;
nf = size(psi,2);

z1 = 0.9*exp(-j*pi*0.5);
z2 = 0.9*exp(-j*pi*0.7);

P = ones(1,nf)./(((abs(ones(1,nf)-z1*exp(-j*psi))).^2).*((abs(ones(1,nf)-z2*exp(-j*psi))).^2));

plot(psi/pi,10*log10(P))
xlabel('\psi/\pi','Fontsize',14)
ylabel(['{\itP}_{\itxx}(\psi)/\sigma^2_{\itu} (dB)'],'Fontsize',14)

grid on