波束形成中的主瓣宽度

阵列信号处理相关基础知识及主瓣宽度

      • 导向矢量
      • 阵列方向图
      • 确知波束形成
      • 普通波束形成
      • 主瓣宽度
        • 确知波束形成主瓣宽度
        • 普通波束形成主瓣宽度

在讨论主瓣宽度之前,首先得了解导向矢量、波束形成、阵列方向图的概念,这些是阵列信号处理中最基础的知识。

导向矢量

在这里插入图片描述

图1. 均匀线性阵接收信号模型

对于均匀线阵,俯仰角 θ \theta θ的定义域通常为 θ ∈ ( − 9 0 ∘ , 9 0 ∘ ) \theta \in (-90^{\circ},90^{\circ}) θ(90,90)。设阵列参考点为 ο \omicron ο,即左起第一个阵元。由几何关系我们可以知道,第 m m m个阵元相对于参考点的波程差为 ( m − 1 ) d s i n θ (m-1)d\rm{sin}\theta (m1)dsinθ,因此我们可以得到第 m m m个阵元相对于参考点的时延 τ m \tau_m τm
τ m = ( m − 1 ) d s i n θ c \tau_m=\frac{(m-1)d\rm{sin}\theta}{c} τm=c(m1)dsinθ
利用上式,均匀线阵的导向矢量可以表示为:
a ( θ ) = [ 1 , e j 2 π d s i n θ λ , . . . , e j 2 π ( M − 1 ) d s i n θ λ ] T \bm{a}(\theta)=[1,e^{j\frac{2\pi d\rm{sin}\theta}{\lambda}},...,e^{j\frac{2\pi (M-1)d\rm{sin}\theta}{\lambda}}]^{T} a(θ)=[1,ejλ2πdsinθ,...,ejλ2π(M1)dsinθ]T
在均匀线性阵中,要求相邻阵元间距 d ≤ λ / 2 d\leq \lambda/2 dλ/2,否则会造成相位混叠,进而影响单脉冲测角。
由导向矢量,可以得到来波方向为 θ \theta θ的信号 s ( t ) s(t) s(t)的阵列输出为:
y = a ( θ ) s ( t ) y=\bm{a}(\theta)s(t) y=a(θ)s(t)
一个包含 M M M个阵元的阵列,在 N N N个快拍下,接收到的信号可以排列为一个 N × M N\times M N×M的矩阵,有没有一种方法,将这 M M M个阵元接收到的信号进行加权求和得到一个还是 N × 1 N\times 1 N×1的向量呢?这就是波束形成,不同的权值计算方法对应不同的波束形成。在了解波束形成之前,首先要了解阵列方向图的概念。

阵列方向图

首先要知道,方向图是一条曲线,横坐标是角度,纵坐标是强度/增益。方向图是指给定权矢量对不同角度信号的阵列响应,它表明了阵列输出与来波方向之间的关系。阵列方向图能够反映波束形成的效果,它与来波方向 θ \theta θ有关,定义为:
F ( θ ) = W H a ( θ ) F(\theta)=W^Ha(\theta) F(θ)=WHa(θ)
W = [ w 1 , w 2 , ⋯   , w M ] T W=[w_1,w_2,\cdots,w_M]^T W=[w1,w2,,wM]T
通常情况下,取模的平方,再对其做归一化处理,得到幅度方向图增益为
G ( θ ) = ∣ F ( θ ) ∣ m a x ∣ F ( θ ) ∣ G(\theta)=\frac{|F(\theta)|}{max|F(\theta)|} G(θ)=maxF(θ)F(θ)
功率方向图增益为
G ( θ ) = ∣ F ( θ ) ∣ 2 m a x ∣ F ( θ ) ∣ 2 G(\theta)=\frac{|F(\theta)|^2}{max|F(\theta)|^2} G(θ)=maxF(θ)2F(θ)2
通常以分贝为单位,幅度方向图为:
G ( θ ) = 10 l g ∣ F ( θ ) ∣ m a x ∣ F ( θ ) ∣ G(\theta)=10lg\frac{|F(\theta)|}{max|F(\theta)|} G(θ)=10lgmaxF(θ)F(θ)
功率方向图为:
G ( θ ) = 20 l g ∣ F ( θ ) ∣ m a x ∣ F ( θ ) ∣ G(\theta)=20lg\frac{|F(\theta)|}{max|F(\theta)|} G(θ)=20lgmaxF(θ)F(θ)
一般情况下,我们约定阵列发现方向为 0 ∘ 0^{\circ} 0,顺时针方向为正角度方向,逆时针方向为负角度方向。即 θ ∈ [ − 9 0 ∘ , 9 0 ∘ ] \theta\in[-90^{\circ},90^{\circ}] θ[90,90]。一般情况下,方向图分为两类,一是静态方向图,即为阵列输出的直接相加,其方向图的最大增益出现再阵列的法线方向,即 θ = 0 ∘ \theta=0^{\circ} θ=0;二是带有指向的方向图,其信号的指向通过加权相位的控制来实现。

确知波束形成

下面以静态方向图为例,它也被称之为确知波束形成,即各个阵元信号进行均匀加权, W = [ 1 , 1 , ⋯   , 1 ] T W=[1,1,\cdots,1]^T W=[1,1,,1]T
对于确知波束形成,可以得到方向图功率增益为
G ( θ ) = 1 N 2 ∣ ∑ n = 1 N e j 2 π λ d ( n − 1 ) s i n θ ∣ 2 = ∣ s i n ( π d N s i n θ / λ ) N s i n ( π d s i n θ / λ ) ∣ 2 G(\theta)=\frac{1}{N^2}|\sum_{n=1}^{N}e^{j\frac{2\pi}{\lambda}d(n-1)sin\theta}|^2=|\frac{sin(\pi dNsin \theta /\lambda)}{Nsin(\pi d sin\theta /\lambda)} |^2 G(θ)=N21n=1Nejλ2πd(n1)sinθ2=Nsin(πdsinθ/λ)sin(πdNsinθ/λ)2
阵元数 M = 16 M=16 M=16,阵元间距为 1 / 2 λ 1/2\lambda 1/2λ的归一化功率方向图:

在这里插入图片描述

图1. 均匀加权ULA阵的功率增益方向图

可以看到,当权值都为1的时候,波束形成后,在0°方向形成峰值,其他方向进行了抑制。

普通波束形成

上面的确知波束形成,权重 W = [ 1 , 1 , ⋯   , 1 ] T W=[1,1,\cdots,1]^T W=[1,1,,1]T,使得阵列在0°方向形成了最大增益,然而在实际情况中,感兴趣的目标信号不一定在法线方向。因此有必要对权值进行调整,使得对阵元接收数据进行加权求和后,主波束总指向期望信号的方向,因为方向图会进行归一化,所以这个权值没必要是大于1的幅度,我们直接取与一个复指数就 e j φ e^{j\varphi} ejφ行了。可以证明,当权重 W W W等于期望方向 θ 0 \theta_0 θ0的导向矢量时,权系数能够行号对各个阵元接收到的信号进行相位补偿,使得相加前的各阵元信号相位相同。
W = a ( θ 0 ) W=a(\theta_0) W=a(θ0)
此时在期望方向获得的增益为:
F ( θ 0 ) = a H ( θ 0 ) a ( θ 0 ) = M F(\theta_0)=a^{H}(\theta_0)a(\theta_0)=M F(θ0)=aH(θ0)a(θ0)=M
增益值刚好等于阵元数 M M M,所以阵元数越多,获得的增益也就越强,后面我们对波束宽度的分析还可以知道,当阵元越多的时候,波束的宽度会越小,即阵列的方向性越强,空间分辨能力越强。如果把一个线阵看作是一个线性孔径的等间距采样,它等效为:阵元数越多的时候,天线的孔径越大。
使用此系数的波束形成器被称之为普通波束形成器(CBF)。
下面是一个指向10°,阵元数 M = 16 M=16 M=16,阵元间距为 1 / 2 λ 1/2\lambda 1/2λ的归一化功率方向图:
在这里插入图片描述

图2.普通波束形成ULA阵的方向图

主瓣宽度

我们除了想知道主瓣指向的角度外,还希望得到空间的分辨能力。我们定义主波束宽度为方向图功率相对最高增益下降3dB处的宽度,即为半功率增益

确知波束形成主瓣宽度

对于确知波束形成,可以得到方向图功率增益为
G ( θ ) = 1 N 2 ∣ ∑ n = 1 N e j 2 π λ d ( n − 1 ) s i n θ ∣ 2 = ∣ s i n ( π d N s i n θ / λ ) N s i n ( π d s i n θ / λ ) ∣ 2 G(\theta)=\frac{1}{N^2}|\sum_{n=1}^{N}e^{j\frac{2\pi}{\lambda}d(n-1)sin\theta}|^2=|\frac{sin(\pi dNsin \theta /\lambda)}{Nsin(\pi d sin\theta /\lambda)} |^2 G(θ)=N21n=1Nejλ2πd(n1)sinθ2=Nsin(πdsinθ/λ)sin(πdNsinθ/λ)2
根据主波束宽度的定义,

G ( θ ) = ∣ s i n ( π d N s i n θ 0.5 / λ ) N s i n ( π d s i n θ 0.5 / λ ) ∣ 2 = 0.5 G(\theta)=|\frac{sin(\pi dNsin \theta_{0.5} /\lambda)}{Nsin(\pi d sin\theta_{0.5} /\lambda)} |^2=0.5 G(θ)=Nsin(πdsinθ0.5/λ)sin(πdNsinθ0.5/λ)2=0.5
通常主波束很窄, ∣ θ 0.5 ∣ |\theta_{0.5}| θ0.5非常小,因此 s i n ( π d s i n θ 0.5 / λ ) ≈ π d s i n θ 0.5 / λ sin(\pi d sin\theta_{0.5}/\lambda)\approx \pi d sin\theta_{0.5}/\lambda sin(πdsinθ0.5/λ)πdsinθ0.5/λ,因此可以将上式化为:
s i n ( π d N s i n θ 0.5 / λ ) N π d s i n θ 0.5 / λ = 1 2 \frac{sin(\pi dNsin \theta_{0.5} /\lambda)}{N\pi d sin\theta_{0.5} /\lambda} =\frac{1}{\sqrt{2}} Nπdsinθ0.5/λsin(πdNsinθ0.5/λ)=2 1
上式近似为sinc函数,由此可以得到
N π d s i n θ 0.5 λ = 1.39 \frac{N\pi dsin\theta_{0.5}}{\lambda}=1.39 λNπdsinθ0.5=1.39

s i n θ 0.5 = 1.39 λ N π d sin\theta_{0.5}=\frac{1.39\lambda}{N\pi d} sinθ0.5=Nπd1.39λ
又有 s i n θ 0.5 ≈ θ 0.5 sin\theta_{0.5}\approx\theta_{0.5} sinθ0.5θ0.5,得到
θ 0.5 = 1.39 λ N π d ( r a d ) = 1.39 × 180 λ π 2 N d ( ∘ ) ≈ 25.35 λ N d ( ∘ ) \theta_{0.5}=\frac{1.39\lambda}{N\pi d}(rad)=\frac{1.39\times180\lambda }{\pi^2 Nd}(^\circ)\approx\frac{25.35\lambda }{ Nd}(^\circ) θ0.5=Nπd1.39λ(rad)=π2Nd1.39×180λ()Nd25.35λ()
要注意的是,主瓣宽度是两倍的主波束最大角度,如下图所示:
在这里插入图片描述

图3.主瓣宽度示意图

可以求得主波束宽度为:
θ m b = 2 θ 0.5 = 25.35 × 2 λ N d ( ∘ ) = 50.7 λ N d ( ∘ ) \theta{mb}=2\theta_{0.5}=\frac{25.35\times2\lambda }{ Nd}(^\circ)=\frac{50.7\lambda }{ Nd}(^\circ) θmb=2θ0.5=Nd25.35×2λ()=Nd50.7λ()
对于常见的阵列结构,阵元间距 d = 1 / 2 λ d=1/2\lambda d=1/2λ的情况下,可以得到主波束宽度为:
θ m b = 2 θ 0.5 = 25.35 × 2 λ N d ( ∘ ) = 50.7 λ N 1 / 2 λ ( ∘ ) = 101.4 N ( ∘ ) \theta{mb}=2\theta_{0.5}=\frac{25.35\times2\lambda }{ Nd}(^\circ)=\frac{50.7\lambda }{ N1/2\lambda}(^\circ)=\frac{101.4 }{ N}(^\circ) θmb=2θ0.5=Nd25.35×2λ()=N1/2λ50.7λ()=N101.4()

普通波束形成主瓣宽度

普通波束形成的情况下,波束指向期望方向,这时,主瓣宽度相对于确知波束形成而言相对较大,我们也可以将确知波束形成当作是普通波束形成的一种特例,即波束指向0°时为确知波束形成,当指向其他的期望方向时为普通波束形成。当天线阵列在进行扫描时,波束就会偏离法线方向,假设波束指向偏离法线方向的角度为 θ Δ \theta_{\Delta} θΔ,此时阵列在空间上有效的接收信号的孔径为 N d c o s θ Δ Ndcos\theta_{\Delta} NdcosθΔ,经推导可以得到主波束的宽度为:
θ m b s = θ m b c o s θ Δ = 50.7 λ N d c o s θ Δ ( ∘ ) \theta_{mbs}=\frac{\theta{mb}}{cos\theta_{\Delta}}=\frac{50.7\lambda }{ Ndcos\theta_{\Delta}}(^\circ) θmbs=cosθΔθmb=NdcosθΔ50.7λ()
对于常见的阵列结构,阵元间距 d = 1 / 2 λ d=1/2\lambda d=1/2λ的情况下,可以得到主波束宽度为:
θ m b s = 2 θ 0.5 = 50.7 λ N 1 / 2 λ c o s θ Δ ( ∘ ) = 101.4 N c o s θ Δ ( ∘ ) \theta{mbs}=2\theta_{0.5}=\frac{50.7\lambda }{ N1/2\lambda cos\theta_{\Delta}}(^\circ)=\frac{101.4 }{ Ncos\theta_{\Delta}}(^\circ) θmbs=2θ0.5=N1/2λcosθΔ50.7λ()=NcosθΔ101.4()
下面是一些常见的实例:
波束指向0°的时候,阵元数为16,主波束宽度为6.34°
波束指向5°的时候,阵元数为16,主波束宽度为6.36°
波束指向0°的时候,阵元数为8,主波束宽度为12.68°
波束指向5°的时候,阵元数为8,主波束宽度为12.72°

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:/a/152148.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

编译智能合约以及前端交互工具库(Web3项目一实战之三)

我们已然在上一篇 Web3项目灵魂所在之智能合约编写(Web3项目一实战之二) ,为项目写好了智能合约代码。 但身为开发人员的我们,深知高级编程语言所编写出来的代码,都是需要经过编译,而后外部方能正常调用。很显然,使用solidity这门新的高级编程语言编写出来的智能合约,也…

【机器学习】线性回归算法:原理、公式推导、损失函数、似然函数、梯度下降

1. 概念简述 线性回归是通过一个或多个自变量与因变量之间进行建模的回归分析,其特点为一个或多个称为回归系数的模型参数的线性组合。如下图所示,样本点为历史数据,回归曲线要能最贴切的模拟样本点的趋势,将误差降到最小。 2. 线…

基于旗鱼算法优化概率神经网络PNN的分类预测 - 附代码

基于旗鱼算法优化概率神经网络PNN的分类预测 - 附代码 文章目录 基于旗鱼算法优化概率神经网络PNN的分类预测 - 附代码1.PNN网络概述2.变压器故障诊街系统相关背景2.1 模型建立 3.基于旗鱼优化的PNN网络5.测试结果6.参考文献7.Matlab代码 摘要:针对PNN神经网络的光滑…

云表|低代码软件开发“外挂”,新时代的黑科技

随着技术的日新月异,现代企业对于软件开发的需求愈加迫切,传统的软件开发方式已然无法满足快速迭代和创新的需求。在这种背景下,低代码开发平台如破茧而出,应运而生。这种平台通过提供可视化的开发工具和预构建的组件,…

NVIDIA安装

电脑显卡类型 两种方法: 选择对应的版本 产品系列下载Notebooks,这样产品才会出现Laptop的GPU(Laptop是代表笔记本)。 下载完之后双击安装,更改下载路径后,选择默认的下载即可。 卸载 如果之后要卸载…

Spring6(三):面向切面AOP

文章目录 4. 面向切面:AOP4.1 场景模拟4.1.1 声明接口4.1.2 创建实现类4.1.3 创建带日志功能的实现类4.1.4 提出问题 4.2 代理模式4.2.1 概念4.2.2 静态代理4.2.3 动态代理4.2.4 测试 4.3 AOP概念4.3.1 相关术语①横切关注点②通知(增强)③切…

chrome 浏览器个别字体模糊不清

特别是在虚拟机里,有些字体看不清,但是有些就可以,设置办法: chrome://settings/fonts 这里明显可以看到有些字体就是模糊的状态: 把这种模糊的字体换掉即可解决一部分问题。 另外,经过观察,…

Neuro-Oncology | IF:15.9 CUTTag和RNA-seq联合解析胶质母细胞瘤的耐药性

发表单位:德克萨斯大学圣安东尼奥分校 发表日期:2023年1月18日 期 刊:Neuro-Oncology(IF: 15.9) 研究技术:CUT&Tag-seq、RNA-seq、RT-qPCR(爱基百客均可提供) 2023年1月1…

如何在10亿级别用户中检查用户名是否存在?

题目 不知道大家有没有留意过,在使用一些app注册的时候,提示你用户名已经被占用了,需要更换一个,这是如何实现的呢?你可能想这不是很简单吗,去数据库里查一下有没有不就行了吗,那么假如用户数量…

【人工智能实验】A*算法求解8数码问题 golang

人工智能经典问题八数码求解 实际上是将求解转为寻找最优节点的问题,算法流程如下: 求非0元素的逆序数的和,判断是否有解将开始状态放到节点集,并设置访问标识位为true从节点集中取出h(x)g(x)最小的节点判断取出的节点的状态是不…

Redis - 订阅发布替换 Etcd 解决方案

为了减轻项目的中间件臃肿,由于我们项目本身就应用了 Redis,正好 Redis 的也具备订阅发布监听的特性,正好应对 Etcd 的功能,所以本次给大家讲解如何使用 Redis 消息订阅发布来替代 Etcd 的解决方案。接下来,我们先看 R…

linux之shell

一、是什么 Shell是一个由c语言编写的应用程序,它是用户使用 Linux 的桥梁。Shell 既是一种命令语言,又是一种程序设计语言 它连接了用户和Linux内核,让用户能够更加高效、安全、低成本地使用 Linux 内核 其本身并不是内核的一部分&#x…

Java实现自定义windows右键菜单

要添加Java应用程序到Windows桌面的右键菜单,可以按照以下步骤操作: 创建一个新的.reg文件,并在文本编辑器中打开它。 添加以下代码到.reg文件中,将名称和路径替换为您的Java应用程序的名称和路径。 Windows Registry Editor V…

虚拟化热添加技术在数据备份上的应用

虚拟化中的热添加技术主要是指:无需停止或中断虚拟机的情况下,在线添加物理资源(如硬盘、内存、CPU、网卡等)的技术。热添加技术也是相比物理机一个非常巨大的优势,其使得资源分配变得更加灵活。 虚拟化中的热添加技术…

SOP作业指导书系统如何帮助厂家实现数字化转型

SOP(Standard Operating Procedure,标准操作程序)电子作业操作手册的应用对于厂家实现数字化转型起着至关重要的作用。本文将探讨SOP电子作业操作手册如何帮助厂家实现数字化转型的重要性和优势。 首先,SOP作业指导书可以提高生产…

idea菜单栏任务栏放缩比例修改

在编辑自定义VM选项中增加 -Dide.ui.scale0.8 参数 Help -> Edit Custom VM Options

这家提供数据闭环完整链路的企业,已拿下多家头部主机厂定点

“BEV感知数据闭环”已经成为新一代自动驾驶系统的核心架构。 进入2023年,小鹏、理想、阿维塔、智己、华为问界等汽车品牌正在全力推动从高速NOA到城区NOA的升级。在这一过程当中,如何利用高效的算力支撑、完善的算法模型、大量有效的数据形成闭环&…

Ubuntu部署OpenStack踩坑指南:还要看系统版本?

正文共:1515 字 12 图,预估阅读时间:2 分钟 到目前为止,我对OpenStack还不太了解,只知道OpenStack本身是一个云管理平台(什么是OpenStack?)。那作为云管理平台,我能想到最…

解决网络编程中的EOF违反协议问题:requests库与SSL错误案例分析

1. 问题背景 近期,一个用户在使用requests库进行网络编程时遭遇到了一个不寻常的问题,涉及SSL错误,并提示错误消息为SSLError(SSLEOFError(8, uEOF occurred in violation of protocol (_ssl.c:661)),))。该用户表示已经采取了多种方法来解决…