mohushanmian.m

对波束形成中的非模糊扇面和模糊扇面作了仿真研究

close all
clear all
f0=10000;c=1523;lamda=c/f0;
d=(lamda/2)*10/5;
N=7;
stheta=-1:0.0001:1;
theta=asin(stheta);
fai=2*pi*d*stheta/lamda;
R=sin(N*fai/2)./(N*sin(fai/2));
figure
plot(theta*180/pi,R)
axis([-90 90 -0.2 1.1]);
grid

figure
fai8=90*pi/180;%扫描角度范围
ai=sin(fai8)/8;%根据主波束覆盖的范围在sin（thea）轴上计算最小波束控制角
fai0=[-8*ai -7*ai -6*ai -5*ai -4*ai -3*ai -2*ai -ai 0 ai  2*ai  3*ai  4*ai  5*ai  6*ai  7*ai  8*ai];%sin（thea）轴上各主波束的分布
theta0=asin(fai0);
beida=2*pi*d*sin(theta0)/lamda ;
R1=sin(N*(fai-beida(9))/2)./(N*sin((fai-beida(9))/2));
plot(theta*180/pi,R1,'g')
hold on
R1=sin(N*(fai-beida(8))/2)./(N*sin((fai-beida(8))/2));
plot(theta*180/pi,R1,'y')
R1=sin(N*(fai-beida(10))/2)./(N*sin((fai-beida(10))/2));
plot(theta*180/pi,R1,'b')
R1=sin(N*(fai-beida(7))/2)./(N*sin((fai-beida(7))/2));
plot(theta*180/pi,R1,'r')
R1=sin(N*(fai-beida(11))/2)./(N*sin((fai-beida(11))/2));
plot(theta*180/pi,R1,'k')

R1=sin(N*(fai-beida(12))/2)./(N*sin((fai-beida(12))/2));
plot(theta*180/pi,R1,'c')

R1=sin(N*(fai-beida(6))/2)./(N*sin((fai-beida(6))/2));
plot(theta*180/pi,R1,'m')
plot(-17.45,0:0.001:1)
plot(17.45,0:0.001:1)
hold off

figure
R=sin(N*(fai-beida(9))/2)./(N*sin((fai-beida(9))/2));
polar(-(theta-pi/2),R,'g')
hold on
R=sin(N*(fai-beida(8))/2)./(N*sin((fai-beida(8))/2));
polar(-(theta-pi/2),R,'y')
R=sin(N*(fai-beida(10))/2)./(N*sin((fai-beida(10))/2));
polar(-(theta-pi/2),R,'b')
R=sin(N*(fai-beida(7))/2)./(N*sin((fai-beida(7))/2));
polar(-(theta-pi/2),R,'r')
R=sin(N*(fai-beida(11))/2)./(N*sin((fai-beida(11))/2));
polar(-(theta-pi/2),R,'k')
R=sin(N*(fai-beida(6))/2)./(N*sin((fai-beida(6))/2));
polar(-(theta-pi/2),R,'c')
R=sin(N*(fai-beida(12))/2)./(N*sin((fai-beida(12))/2));
polar(-(theta-pi/2),R,'r')

hold off
theta0*180/pi