波束形成Matlab程序

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 度，20lo g（dB） 随 着 阵 元 数 的 增 加 ， 波 束 宽 度 变 窄 ， 分 辨 力 提 高 ： 仿 真 图 如 下 ： 2. 波 束 宽 度 与波 达方向及阵元数的关系 clc clear all close all ima=sqrt(-1); element_nu m1=16; %阵元数 element_nu m2=128; element_nu m3=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_nu m1*d)); psi(j)=theta(j)*pi/180-asin(sin(theta(j)*pi/180)-lamda/(element_nu m2*d)); beta(j)=theta(j)*pi/180-asin(sin(theta(j)*pi/180)-lamda/(element_nu m3*d)); end figu re; plot(theta,fai,'r',theta,psi,'b',theta,beta,'g'),grid on x label('theta'); y label('Width in radians') title('波束宽度与波达方向及阵元数的关系') 仿真图如下： 3. 当 阵 元 间 距 时 ， 会 出 现 栅 瓣 ， 导 致 空 间 模 糊 。仿真图如下： / 2d 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...