第三章(一)_优化阵列信号处理_

clc;clf;close all; %% 均匀线列阵阵元间距与栅瓣的关系 M = 10; d1_lambda = 1/2;%d/λ，阵元间距 d2_lambda = 1; d3_lambda = 3/2; d4_lambda = 2; theta0_1 = 0;%期望波束观察方向，单位：度 theta0_2 = 30; theta0_3 = 90; scanAngle = -90:0.1:90;%波束扫描范围 B1 = sin(pi*M*d1_lambda*(sin(scanAngle/180*pi)-sin(theta0_1/180*pi)))... ./(M*sin(pi*d1_lambda*(sin(scanAngle/180*pi)-sin(theta0_1/180*pi)))); B1_dB = 20*log10(abs(B1)); B2 = sin(pi*M*d2_lambda*(sin(scanAngle/180*pi)-sin(theta0_1/180*pi)))... ./(M*sin(pi*d2_lambda*(sin(scanAngle/180*pi)-sin(theta0_1/180*pi)))); B2_dB = 20*log10(abs(B2)); B3 = sin(pi*M*d3_lambda*(sin(scanAngle/180*pi)-sin(theta0_1/180*pi)))... ./(M*sin(pi*d3_lambda*(sin(scanAngle/180*pi)-sin(theta0_1/180*pi)))); B3_dB = 20*log10(abs(B3)); B4 = sin(pi*M*d4_lambda*(sin(scanAngle/180*pi)-sin(theta0_1/180*pi)))... ./(M*sin(pi*d4_lambda*(sin(scanAngle/180*pi)-sin(theta0_1/180*pi)))); B4_dB = 20*log10(abs(B4)); B5 = sin(pi*M*d1_lambda*(sin(scanAngle/180*pi)-sin(theta0_2/180*pi)))... ./(M*sin(pi*d1_lambda*(sin(scanAngle/180*pi)-sin(theta0_2/180*pi)))); B5_dB = 20*log10(abs(B5)); B6 = sin(pi*M*d2_lambda*(sin(scanAngle/180*pi)-sin(theta0_2/180*pi)))... ./(M*sin(pi*d2_lambda*(sin(scanAngle/180*pi)-sin(theta0_2/180*pi)))); B6_dB = 20*log10(abs(B6)); B7 = sin(pi*M*d3_lambda*(sin(scanAngle/180*pi)-sin(theta0_2/180*pi)))... ./(M*sin(pi*d3_lambda*(sin(scanAngle/180*pi)-sin(theta0_2/180*pi)))); B7_dB = 20*log10(abs(B7)); B8 = sin(pi*M*d4_lambda*(sin(scanAngle/180*pi)-sin(theta0_2/180*pi)))... ./(M*sin(pi*d4_lambda*(sin(scanAngle/180*pi)-sin(theta0_2/180*pi)))); B8_dB = 20*log10(abs(B8)); B9 = sin(pi*M*d1_lambda*(sin(scanAngle/180*pi)-sin(theta0_3/180*pi)))... ./(M*sin(pi*d1_lambda*(sin(scanAngle/180*pi)-sin(theta0_3/180*pi)))); B9_dB = 20*log10(abs(B9)); B10 = sin(pi*M*d2_lambda*(sin(scanAngle/180*pi)-sin(theta0_3/180*pi)))... ./(M*sin(pi*d2_lambda*(sin(scanAngle/180*pi)-sin(theta0_3/180*pi)))); B10_dB = 20*log10(abs(B10)); B11 = sin(pi*M*d3_lambda*(sin(scanAngle/180*pi)-sin(theta0_3/180*pi)))... ./(M*sin(pi*d3_lambda*(sin(scanAngle/180*pi)-sin(theta0_3/180*pi)))); B11_dB = 20*log10(abs(B11)); B12 = sin(pi*M*d4_lambda*(sin(scanAngle/180*pi)-sin(theta0_3/180*pi)))... ./(M*sin(pi*d4_lambda*(sin(scanAngle/180*pi)-sin(theta0_3/180*pi)))); B12_dB = 20*log10(abs(B12)); subplot(3,1,1); plot(sin(scanAngle/180*pi),B1_dB,sin(scanAngle/180*pi),B2_dB,... sin(scanAngle/180*pi),B3_dB,sin(scanAngle/180*pi),B4_dB); title('\theta_0=0^°'); ax = gca; ax.YLim = [-60 0]; legend('d=λ/2','d=λ','d=3λ/2','d=2λ'); subplot(3,1,2); plot(sin(scanAngle/180*pi),B5_dB,sin(scanAngle/180*pi),B6_dB,... sin(scanAngle/180*pi),B7_dB,sin(scanAngle/180*pi),B8_dB); title('\theta_0=30^°'); ax = gca; ax.YLim = [-60 0]; legend('d=λ/2','d=λ','d=3λ/2','d=2λ'); subplot(3,1,3); plot(sin(scanAngle/180*pi),B9_dB,sin(scanAngle/180*pi),B10_dB,... sin(scanAngle/180*pi),B11_dB,sin(scanAngle/180*pi),B12_dB); title('\theta_0=90^°'); ax = gca; ax.YLim = [-60 0]; legend('d=λ/2','d=λ','d=3λ/2','d=2λ');