%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