波束形成中的主瓣宽度

news2025/2/25 23:18:36

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

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

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

导向矢量

在这里插入图片描述

图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°

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

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

【LeetCode刷题-滑动窗口】--1423.可获得的最大点数

1423.可获得的最大点数 思路: 数组cardPoints的长度为n,由于只能从开头和末尾拿k张卡牌,所以最后剩下的必然是连续的n-k张卡牌,可以通过求出剩余卡牌点数之和的最小值,来求出拿走卡牌点数之和的最大值 算法&#xff…

centos7升级python2到python3.6.8使用yum安装问题

背景 公司爬虫需要使用python3.6.8版本,因此升级centos的python2到python3.6.8。但是当使用yum安装包时 ,出现如下异常: [rootlocalhost bin]# yum install npm -y Loaded plugins: fastestmirror, product-id, search-disabled-repos, sub…

S25FL256S介绍及FPGA实现思路

本文介绍 S25FL256S 这款 FLASH 芯片,并进行 FPGA 读写控制的实现(编程思路及注意事项)。 文章目录 S25FL-S 介绍管脚功能说明SPI 时钟模式SDRDDR 工作模式FLASH存储阵列(地址空间映射)常用寄存器及相关指令Status Reg…

Google codelab WebGPU入门教程源码<4> - 使用Uniform类型对象给着色器传数据(源码)

对应的教程文章: https://codelabs.developers.google.com/your-first-webgpu-app?hlzh-cn#5 对应的源码执行效果: 对应的教程源码: 此处源码和教程本身提供的部分代码可能存在一点差异。 class Color4 {r: number;g: number;b: number;a: number;constructor(pr 1.0, …

如何选择一款车规级电源开关/驱动器芯片?NCV8405ASTT3G自保护低压侧驱动器

关于车规级芯片(Automotive Grade Chip),车规级芯片是专门用于汽车行业的芯片,具有高可靠性、高稳定性和低功耗等特点,以满足汽车电子系统的严格要求。这些芯片通常用于车载电子控制单元(ECU)和…

zabbix基本介绍 安装部署 页面访问

这里写目录标题 一、zabbix 监控1、zabbix 监控架构2、zabbix 监控报警渠道3、Zabbix 优点4、Zabbix 缺点5、Zabbix 监控系统监控对象6、Zabbix监控方式7、zabbix 架构1、Server2、数据库存储3、Web界面4、Proxy 代理服务器5、Agent监控代理6、数据流 8、Zabbix常用术语的含义1…

贝茄莱(BR)实时以太网通讯测试

测试目的 测试系统实时以太网通讯的稳定性,以及系统掉电、数据量大小、数据传输速度对通讯的影响,并测试系统对通讯状态参数的监控。 测试标准 以PLC系统间Profinet通讯作为测试对象,在数据传输周期为20ms,双向数据传输变量都不…

【机器学习8】采样

1 均匀分布随机数 均匀分布是指整个样本空间中的每一个样本点对应的概率(密度) 都是相等的。 根据样本空间是否连续, 又分为离散均匀分布和连续均匀分布。编程实现均匀分布随机数生成器一般可采用线性同余法(Linear Congruential…

android初集成flutter,遇到的问题

环境 studio版本:2022.1.1 flutter版本:2.8.0 电脑:mac flutter项目总是报错,编译不过 以 Resources Root 加载 记得设置dart:主工程和flutter项目都需要设置,否则不出现手机链接 下面这个样子就是好了&…

计算机视觉基础(9)——相机标定与对极几何

前言 本节我们将学习相机标定和对极几何两部分的内容。 在相机标定部分,我们将学习直接线性变换(Direct Linear Transform, DL),张正友标定法(Zhang’s Method)和 Perspective-n-Point (PnP) 这三种方法。 在对极几何部…

openbabel包的安装

今天在安装openbabel的时候,使用pip / 通过源码的方式进行安装openbabel,但是都并没有什么用,直接使用: conda install openbabel 就openbabel3.1.1成功了

最新完美版积分商城系统-奇偶商城系统源码+独立代理后台+附搭建教程

源码简介: 最新完美版积分商城系统,网购商城系统源码,是更新的奇偶商城系统源码,它拥有独立代理后台,而且内附搭建教程。 1.演示环境:Linux Centos7以上版本 宝塔 2.Nginx 1.18.0 PHP7.0 Mysql5.6 3…

Talk | UCSB博士生赵宣栋:生成式AI时代的水印技术

本期为TechBeat人工智能社区第546期线上Talk! 北京时间11月15日(周三)20:00,UC Santa Barbara博士生—赵宣栋的Talk已准时在TechBeat人工智能社区开播! 他与大家分享的主题是: “生成式AI时代的水印技术”,系统地介绍了他的团队在…

栈和队列的实现及相关面试题

栈和队列 栈概念与结构栈的功能栈的实现头文件Stack.h栈的结构体 Stack 源文件Stack.c初始化 void StackInit(Stack* ps)压栈 void StackPush(Stack* ps, STDataType data)出栈 void StackPop(Stack* ps)返回栈顶的值 STDataType StackTop(Stack* ps)返回栈中元素的个数 int St…

一个 不用充钱 也能让你变强的 VSCode 插件

今天给大家推荐一款不用充钱也能让你变强的 vscode 插件 通义灵码(TONGYI Lingma),可以称之为 copilot 的替代甜品 💪 什么是 通义灵码(TONGYI Lingma) 通义灵码(TONGYI Lingma)&am…

《向量数据库指南》——TruLens 用于语言模型应用跟踪和评估

TruLens 用于语言模型应用跟踪和评估 TruLens 是一个用于评估语言模型应用(如 RAG)的性能的开源库。通过 TruLens,我们还可以利用语言模型本身来评估输出、检索质量等。 构建语言模型应用时,多数人最关心的问题是 AI 幻觉(hallucination)。RAG 通过为语言模型提供检索上下文…

校招面试技巧:如何回答好一个问题

今天分享的主题是“如何回答好一个问题”。通过与同学们进行交流,发现一个比较普遍的现象,大家在回答问题时,时常会让面试官觉得自己在背八股,从而影响面试官的印象。那么我们应该怎么避免这种情况呢? 1.基础知识点 …

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

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

解决Python Requests库中处理重定向时的多重Cookie问题

问题背景 在更新至f73bda06e9版本后,用户发现某些请求会引发CookieConflictError。具体来说,任何设置饼干且重定向到设置相同饼干的页面的请求都会引发CookieConflictError。 重现步骤 1、更新Requests至上述版本。 2、从中国以外的任何地方向baidu…

Java 设计模式——中介者模式

目录 1.概述2.结构3.案例实现3.1.抽象中介类3.2.抽象同事类3.3.具体同事类3.4.具体中介类3.5.测试 4.优缺点5.使用场景 1.概述 (1)一般来说,同事类之间的关系是比较复杂的,多个同事类之间互相关联时,他们之间的关系会…