常规波束形成概述

  • 波束形成是将基阵各阵元输出进行一定的处理,从而形成基阵空间指向性的技术。一个波束形成器可以对空间某方位的信号有响应,而抑制其他方位的信号,所以又被称为空域滤波。一般分为常规波束形成(CBF)和自适应波束形成(ABF)。
  • 常规波束形成较为基础,一般又可以根据实现的途径不同,分为时延波束形成和频移波束形成。

时延波束形成

  • 在各阵元间添加时延使波束控制于不同方位的方法成为时延波束形成。设定一个目标方向,当信号来源方向在该目标方向时,响应最大;信号来源于其他方向时,响应较小。将这种响应大小随空间角度的变化作图,就是该波束形成器的波束图。

  • 如图,对于均匀直线阵,相邻阵元间存在时延。对于任意角度theta,相邻阵元时延差为:

  • 其中d为阵元间距,theta为信号来源方向与基阵法线所成夹角。以阵元0为参考阵元,其表达式为:

  • 引出各阵元的方向向量,各阵元相差一个时延,以五元阵为例,依次写出为一个行向量:

  • 各阵元的信号即为:

  • 为了使目标方向形成最大响应,得到最优加权向量,对各阵元补偿指定方向的相位,构成加权向量。

  • 输出信号为:

    所以,响应归一化后只需要:

    这就是各方向的响应。

  • 运行结果:

频域波束形成

频域波束形成是对各阵元信号的谱线进行频移再形成。以M元直线阵、N点信号为例:

  • 对各阵元信号进行目标方向的相位补偿

  • 对该补偿后的信号进行快速傅里叶变换,得到了M*N维矩阵的频谱

  • 对任意角度求取响应:任一角度上每个阵元信号所在谱线进行该角度的频移,如图,时域时延对应的频域乘以该部分。

  • 将各阵元某一谱线求和,对各谱线求和,得到该角度的响应。

运行结果如下:

两种波束形成方法比较

  • 研究两者形成过程可以看出

时域波束形成对整个信号进行时延,直接波束形成;
频域波束形成只选取了信号所在频率进行波束形成。
这里我并没有给信号加噪声,二者表现不明显,但是从过程中可以推测出,如果加上噪声,频域波束形成无疑会筛选掉非信号频带的噪声,波束形成效果会更好!

此处附上代码,信号处理新人,望大家指正:

%% 时延波束形成
clc;clear;close all;
c=1500;f0=800;
d=c/f0/5;fs=40000;
N=2^10;M=5;
t=(0:N-1)/fs;
% st=rand(1,M);
xt=cos(2*pi*f0*t);
theta=-90:0.1:90;theta0=30;
for i=1:Ma(i,:)=exp(1j*2*pi*(i-1)*d*sind(theta));w(1,i)=exp(1j*2*pi*(i-1)*d*sind(theta0));
end
BP=conj(w)*a/M;
figure(1)
subplot(121)
plot(theta,real(BP));title('时域波束形成波束图');xlabel('角度(°)');ylabel('幅值');
subplot(122)
polarplot(theta*pi/180,real(BP));title('时域波束形成极坐标图');
%% 频域波束形成
clc;clear;close all;
c=1500;f0=800;fs=40000;%水下声速、中心频率、采样频率
ome=2*pi*f0;%角频率
d=c/f0/2;%阵元间距
N=2^10;%采样点数
M=5;%阵元个数
t1=(0:2*N-1)/fs;
s0=cos(ome*t1);%参考信号
theta0=30;%波束形成目标角度
theta=-60:0.1:60;%刻画角度坐标
f=(0:N-1)*fs/N;%刻画频率坐标
for i=1:Mtao(i)=round((i-1)*d*sind(theta0)/c*fs);%各阵元在目标方向时延差取整,在整个信号上截取时延后的信号s(i,:)=s0(1,tao(i)+1:tao(i)+N);%生成各阵元时域信号sf(i,:)=fft(s(i,:),N);%转换到频域
end
figure(1)
plot(f,real(sf));legend('阵元1','阵元2','阵元3','阵元4','阵元5');fl=790;fh=810;%筛选出信号所在频率谱线进行波束形成
n1=round(fl/fs*N);n2=round(fh/fs*N);
for i=1:length(theta)for j=n1:n2for k=1:Mn11(i)=(k-1)*d*sind(theta(i))/c*fs/N;%对任意角度在该阵元上的时延转换为时域上平移的点sf1(k)=sf(k,j)*exp(-1j*2*pi*j*n11(i));%时域时延,频域相移ends1(j)=sum(sf1);%该谱线上各阵元求和ends2(i)=sum(abs(s1(j))^2);%该角度上各谱线求和
end
figure(2)
subplot(121)
plot(theta,s2/max(s2));title('频域波束形成波束图');xlabel('角度(°)');ylabel('幅值');
subplot(122)
polarplot(theta*pi/180,s2/max(s2));title('频域波束形成极坐标图');

均匀直线阵的常规波束形成matlab相关推荐

  1. 【通信】盖氏圆盘方法(GDE)计算均匀直线阵(ULA)中信号源个数附matlab代码

    1 简介 信源数估计是空间谱估计中的关键技术,研究符合实际应用环境的稳健的信源数估计方法具有十分重要的现实意义.基于空间谱估计中用于估计信源数的传统盖氏圆盘法,该程序用于计算均匀直线阵(ULA)中信号 ...

  2. 六元均匀直线阵的各元间距为_实验二 均匀直线阵

    实验二 均匀直线阵 一.实验目的: 通过 MATLAB 编程,了解均匀直线阵的辐射特性,熟悉影响天线阵辐射的各种因素及其 产生的影响. 二.实验环境: MATLAB 软件 三.实验原理: 单个天线的方 ...

  3. 六元均匀直线阵的各元间距为_天线原理考试卷B.doc

    天线原理考试卷B 4.(8分)某矩形微带天线w=11.43cm,L=7.62cm,εr=2.25,工作于(0,1)模,求对应该模式的频率? 二.(11分)一个七元引向天线,反射器与有源振子的间距为0. ...

  4. 六元均匀直线阵的各元间距为_梦得杯—镀覆技术论文大赛丨微机电系统器件电镀镍厚度均匀性的模拟与改进【刘瑞 许文杰 袁妍妍】...

    微机电系统器件电镀镍厚度均匀性的模拟与改进 刘瑞,许文杰,袁妍妍(江苏科技大学材料科学与工程国家级实验教学示范中心) 作者简介:刘瑞,博士,副教授,主要研究方向为微纳米器件与材料. 文章全文 近年来, ...

  5. 六元均匀直线阵的各元间距为_给棉花地选购滴灌带时记住这几点,不再为棉花滴水时发愁...

    现在马上进入棉花大量采收阶段,好多农户已经把棉花地里的90管拆了拉回来了,而且不少滴灌带生产厂家也主动和农户联系2021年的滴灌带兑换事项,农户们也在相互打听哪家的质量好.价格合适. 选择滴灌带是棉花 ...

  6. 六元均匀直线阵的各元间距为_小间距led显示屏的封装方式有哪些?本文带你了解!...

    小间距LED显示屏产品(一般市场定义为点间距不大于2.5mm)因稳定性.可靠性.耐久性及易维护性开启了各类指挥调度中心的新时代,并不断满足其他中高端应用场景精细化.个性化的产品需求.而对于这些LED小 ...

  7. 六元均匀直线阵的各元间距为_地采暖安装有哪些要求?盘管间距误差多少?盘管固定有何要求?...

    文章导语:地采暖根据热源的不同,其实又分为了水地暖和电地暖.对于集中供热,在供热过程中基本上都是使用的水地暖.家居杂坛想在这篇文章中给大家详细介绍一下水地暖的具体做法的要求.个人发现很多地方的做法都很 ...

  8. 六元均匀直线阵的各元间距为_微波技术与天线复习题

    1 1 .一般来说,微波是指频率为 _________ 到 _________ 的电磁波. 一.填空题 1 . ( 2006 . 2007 . 2009 )波导中的电磁场模式为 ___________ ...

  9. 阵列信号处理之常规波束形成基础+matlab仿真(二)

    DOA估计的三种常见方法:CBF.MVDR(CAPON).MUSIC 统计阵列处理理论:入射波包含了期望信号.干扰信号和噪声,它们中的部分或者全部都可以看成空时随机过程的样本函数.利用信号.噪声和干扰 ...

最新文章

  1. 《C语言及程序设计》实践参考——乘法口诀表
  2. CSS clip:rect矩形剪裁功能
  3. 2018-2019-1 20189210 《LInux内核原理与分析》第四周作业
  4. n分解后的最大乘积java_【编程题】乘积最大的拆分
  5. 科技范足 多家上市川企获省科技进步奖
  6. copy-and-swap idiom详解和实现安全自我赋值
  7. Spring(Data-Rest)
  8. java Structs 介绍
  9. Html5固定链接生成二维码的技巧
  10. 怎样设置路由器禁用其他设备
  11. netlimiter 4 功能使用TCP UDP上行下行流量监控
  12. 修改NPM的默认安装路径
  13. 快速列出字母+数字的全部组合
  14. MacBook Pro 触控栏不能正常使用怎么解决
  15. 分子动力学模拟gro格式转换为 car
  16. 各大互联网大厂年终奖一览表,又是别人家的公司!
  17. matlab编写拉格朗日插值代码函数
  18. 一篇文章说清楚shared_ptr,unique_ptr的区别和关系
  19. CTF新近学习经验和疑难
  20. 成都大数据科技产业园5月开建 一期工程年底建成

热门文章

  1. MySQL之InnoDB的两个配置参数INNODB_FAST_SHUTDOWN和INNODB_FORCE_RECOVERY,用不好可能导致线上万数据丢失
  2. Android关于Color你所知道的和不知道的一切
  3. 直观地简单理解Tarjan算法(寻找有向图中的强连通图)
  4. 特征选择—过滤法(Filter)、嵌入法(Embedded)和包装法(Wrapper)
  5. Dockerfile 详解
  6. inspect python_inspect
  7. SpringBoot基础系列-使用日志
  8. spring-boot-maven-plugin插件坐标爆红
  9. java打印鸭子,python – 鸭子打字和(java)界面的概念
  10. AD学习笔记(二)——PCB布局布线