原创力文档- 1、本文档共31页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
1.均匀线阵方向图
%8 阵元均匀线阵方向图，来波方向为 0 度
clc;
clear all;
close all;
imag=sqrt(-1);
element_num=8;% 阵元数为 8
d_lamda=1/2;% 阵元间距 d 与波长 lamda 的关系
theta=linspace(-pi/2,pi/2,200);
theta0=0;% 来波方向
w=exp(imag*2*pi*d_lamda*sin(theta0)*[0:element_num-1]');
for j=1:length(theta)
a=exp(imag*2*pi*d_lamda*sin(theta(j))*[0:element_num-1]');
p(j)=w'*a;
end
figure;
plot(theta,abs(p)),grid on
xlabel('theta/radian')
ylabel('amplitude')
title('8 阵元均匀线阵方向图 ')
当来波方向为 45 度时，仿真图如下：
8 阵元均匀线阵方向图如下，来波方向为 0 度， 20log （dB ）
随着阵元数的增加，波束宽度变窄，分辨力提高：仿真图如下：
2. 波束宽度与波达方向及阵元数的关系
clc
clear all
close all
ima=sqrt(-1);
element_num1=16; %阵元数
element_num2=128;
element_num3=1024;
lamda=0.03; %波长为 0.03 米
d=1/2*lamda; %阵元间距与波长的关系
theta=0:0.5:90;
for j=1:length(theta);
fai(j)=theta(j)*pi/180-asin(sin(theta(j)*pi/180)-lamda/(element_num1*d));
psi(j)=theta(j)*pi/180-asin(sin(theta(j)*pi/180)-lamda/(element_num2*d));
beta(j)=theta(j)*pi/180-asin(sin(theta(j)*pi/180)-lamda/(element_num3*d));
end
figure;
plot(theta,fai,'r',theta,psi,'b',theta,beta,'g'),grid on
xlabel('theta');
ylabel('Width in radians')
title(' 波束宽度与波达方向及阵元数的关系 ')
仿真图如下：
3. 当阵元间距 d / 2 时，会出现栅瓣，导致空间模糊。仿真图如下：
4. 类似于时域滤波，天线方向图是最优权的傅立叶变换
仿真程序和仿真图如下：
clc
clear all
close all
ima=sqrt(-1);
element_num=32; %阵元数
source_num=1; %信源数
d_lamda=1/2; %阵元间距与波长的关系
theta=linspace(-pi/2,pi/2,200);
theta0=0; % 来波方向
w=exp(ima*2*pi*d_lamda*sin(theta0)*[0:element_num-1]');
for j=1:length(theta);
a=exp(ima*2*pi*d_lamda*sin(theta(j))*[0:element_num-1]');
p(j)=w'*a;
end
figure;
subplot(1,2,1)
plot(theta,abs(p)),grid on
xlabel('theta/radian')
ylabel('amplitude')
title(' 按定义的方向图 ')
pfft=fftshift(fft(w,128));
subplot(1,2,2)
plot
文档评论（0）