fig6_14.m

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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Figure 6.14
% Xiaomin Lu 6/26/00
% updated by K. Bell 11/19/00, 7/25/01, 11/12/01
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear all
close all

N = 10;
n = (-(N-1)/2:(N-1)/2)';
Vm = ones(N,1);
u = -1:1/1000:1;

k1 = 1;
for INR = 10.^([10 20]/10)
   for ui = [0.3 0.5]
      Vi = exp(j*n*pi*ui);
      Sn = INR*Vi*Vi' + eye(N);
      W = inv(Sn)*Vm/(Vm'*inv(Sn)*Vm);
      beam(k1,:) = W'*exp(j*n*pi*u);
      k1 = k1 + 1;
   end
end

beam = 20*log10(abs(beam));
subplot(2,2,1)
plot(u,beam(1,:))
xlabel('{\itu}')
ylabel('Beam pattern (dB)')
%title('\sigma_I^2 = 10dB, ui = 0.3')
axis([-1 1 -30 0])
line([0.3 0.3],[-30 -2.5])
grid

subplot(2,2,2)
plot(u,beam(2,:))
xlabel('{\itu}')
ylabel('Beam pattern (dB)')
%title('\sigma_I^2 = 10dB, ui = 0.5')
axis([-1 1 -30 0])
line([0.5 0.5],[-30 -2.5])
grid

subplot(2,2,3)
plot(u,beam(3,:))
xlabel('{\itu}')
ylabel('Beam pattern (dB)')
%title('\sigma_I^2 = 20dB, ui = 0.3')
axis([-1 1 -30 0])
grid
line([0.3 0.3],[-30 -2.5])

subplot(2,2,4)
plot(u,beam(4,:))
xlabel('{\itu}')
ylabel('Beam pattern (dB)')
%title('\sigma_I^2 = 20dB, ui = 0.5')
axis([-1 1 -30 0])
grid
line([0.5 0.5],[-30 -2.5])