稀疏阵的方向图仿真，给出了不同条件下的稀疏阵列的方向图

稀疏阵的方向图仿真，给出了不同条件下的稀疏阵列的方向图

%8阵元均匀线阵方向图，来波方向为0度
clc;
clear all;
close all;
imag=sqrt（-1）;
element_num=9;%满阵阵元数为9
element_num1=39;%稀布阵对应的满阵阵元为39
d_lamda=1/2;%阵元间距d与波长lamda的关系
theta=linspace（-pi/2pi/22000）;
theta0=0;%来波方向
x=random（‘unif‘00.51element_num）;
%阵元数为9的均匀直线阵
w=exp（imag*2*pi*d_lamda*sin（theta0）*[0:element_num-1]‘）;
%阵元数为9的稀布阵，对应阵元39的直线阵
d_lamda1=[01.52.54.588.5101419];
w1=exp（imag*2*pi*sin（theta0）*d_lamda1‘）;
%阵元数为39的均匀直线阵
w2=exp（imag*2*pi*d_lamda*sin（theta0）*[0:element_num1-1]‘）;
%阵元距离随机分布

w4=exp（imag*2*pi*d_lamda*sin（theta0）*（[0:element_num-1]+x）‘）;
for  j=1:length（theta）
   a=exp（imag*2*pi*d_lamda*sin（theta（j））*[0:element_num-1]‘）;
   a1=exp（imag*2*pi*sin（theta（j））*d_lamda1‘）;
   a2=exp（imag*2*pi*d_lamda*sin（theta（j））*[0:element_num1-1]‘）;
   a4=exp（imag*2*pi*d_lamda*sin（theta（j））*（[0:element_num-1]+x）‘）;
   p（j）=w‘*a;
   p1（j）=w1‘*a1;
   p2（j）=w2‘*a2;
   p4（j）=w4‘*a4;
end
figure（1）;
plot（theta*180/pidb（abs（p）/max（abs（p））））grid on
xlabel（‘theta‘）
ylabel（‘amplitude/db‘）
title（‘9阵元均匀线阵方向图与9元稀布阵方向图‘）
hold on
plot（theta*180/pidb（abs（p1）/max（abs（p1）））‘r‘）;
hold on
plot（theta*180/pidb（abs（p4）/max（abs（p4）））‘g‘）;

figure（2）;
plot（theta*180/pidb（abs（p2）/max（abs（p2））））grid on;
xlabel（‘theta‘）
ylabel（‘amplitude/db‘）
title（‘39阵元均匀线阵方向图与其对应的9元稀布阵方向图‘）
hold on
plot（theta*180/pidb（abs（p1）/max（abs（p1）））‘g‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1499  2012-11-08 17:31  稀布阵方向图.m