存在激励幅度误差的方向图区间计算

2023-12-15 22:20

本文主要是介绍存在激励幅度误差的方向图区间计算,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

概述

阵列通常需要通过加权来实现特定的方向图特性,比如低副瓣(SLL,Low Sidelobe levels),高定向性(high directivity)和调零控制(null control),当阵列单元间存在来自机加工精度、装配精度、幅度漂移、器件不一致等因素所造成的幅度激励误差时,会导致方向图产生畸变,这种畸变主要包括法线方向增益、副瓣电平以及调零深度的变化。

幅相激励分析方法

有很多方法可以评估由于激励误差所导致的方向图畸变,最简单的方法是将激励误差建模为符合一定统计分布的随机变量,后通过蒙特卡洛仿真分析来实现激励误差所产生的影响;还有一种方法同样是将激励误差建模为随机变量,不过是通过理论推导来得到方向图期望作为畸变评估标准。

这两种方法的分析存在缺陷,在某些幅相失配条件下,方向图畸变可能会远超统计方法所给出的期望估计。因此学界将一种名为the interval arithmetic,IA技术引入,用来构建对方向图畸变更紧凑的估计边界,该方法能够处理不确定的被一定边界约束的误差,因此能够得到更可靠更有鲁棒性的阵列方向图估计。

目前该方法的相关研究主要是通过设计更合理的单元误差区域,来收紧方向图畸变的估计范围。从Cartesian IA到circular IA以及Taylor expansion-based方法等。这里参考文献内容给出一种混合基的IA方法,得到的方向图估计具有正定矩阵形式,更方便实际工程中的计算和应用。

理论分析

考虑一个均匀线阵ULA,包含N个等向同性天线单元,每个单元的幅度激励和相位激励分别为a_n和theta_n,n = 0,1,2,...,N-1。阵列辐射方向图(阵因子方向图)F(theta)可以表示为

其中的theta_n可由下式表述,

k为波束,d_n为第n个单元与参考点得距离,同时theta(没有n)为波束指向。由此公式可以得到功率方向图公式如下所示

假定第n个天线单元得激励幅度区间是以a_mid_n为圆心,以a_rad_n为半径得区域,那么第n个单元得激励幅度为

从而得到功率方向图为

上式第一项为功率方向图得中心值,第2项为由于幅度误差导致得1阶畸变,第3项为由于幅度误差导致得2阶畸变,1阶畸变有不等式

2阶畸变有不等式

分别利用1阶和2阶畸变项,可以构建P_inf下界项和P_sup上届项。功率方向图上限是中值项加上1阶和2阶畸变项,功率方向图下限是中值项减去1阶项。

仿真分析

根据前述推导,针对包含10个单元得ULA进行仿真分析,工作频率17GHz,最大扫描90°,单元间距半波长,幅度误差比为3%,仿真得到得波束指向0°,方向图得均值、上边界和下边界如图1所示。可以看到均值被包含在上下边界之中。

图1 波束指向零度时得均值、上下界估计(无幅度加权)

当存在20dB副瓣得切比雪夫幅度加权时,仿真得到得波束指向0°方向图得均值、上边界和下边界如图2所示。均值功率方向图同样被包含在上下边界之中。

图2 波束指向零度时得均值、上下界估计(20dB切比雪夫幅度加权)

小结

根据文献中给出得IA方法,仿真了存在幅度激励误差得10单元ULA阵列方向图畸变上下界,该方法可用于评估相控阵阵列工程设计。

这篇关于存在激励幅度误差的方向图区间计算的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/498117

相关文章

hdu 1754 I Hate It(线段树,单点更新,区间最值)

题意是求一个线段中的最大数。 线段树的模板题,试用了一下交大的模板。效率有点略低。 代码: #include <stdio.h>#include <string.h>#define TREE_SIZE (1 << (20))//const int TREE_SIZE = 200000 + 10;int max(int a, int b){return a > b ? a :

poj 1113 凸包+简单几何计算

题意: 给N个平面上的点,现在要在离点外L米处建城墙,使得城墙把所有点都包含进去且城墙的长度最短。 解析: 韬哥出的某次训练赛上A出的第一道计算几何,算是大水题吧。 用convexhull算法把凸包求出来,然后加加减减就A了。 计算见下图: 好久没玩画图了啊好开心。 代码: #include <iostream>#include <cstdio>#inclu

uva 1342 欧拉定理(计算几何模板)

题意: 给几个点,把这几个点用直线连起来,求这些直线把平面分成了几个。 解析: 欧拉定理: 顶点数 + 面数 - 边数= 2。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#inc

uva 11178 计算集合模板题

题意: 求三角形行三个角三等分点射线交出的内三角形坐标。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vector>#include <

XTU 1237 计算几何

题面: Magic Triangle Problem Description: Huangriq is a respectful acmer in ACM team of XTU because he brought the best place in regional contest in history of XTU. Huangriq works in a big compa

hdu4267区间统计

题意:给一些数,有两种操作,一种是在[a,b] 区间内,对(i - a)% k == 0 的加value,另一种操作是询问某个位置的值。 import java.io.BufferedInputStream;import java.io.BufferedReader;import java.io.IOException;import java.io.InputStream;import

hdu4417区间统计

给你一个数列{An},然后有m次查询,每次查询一段区间 [l,r] <= h 的值的个数。 import java.io.BufferedInputStream;import java.io.BufferedReader;import java.io.IOException;import java.io.InputStream;import java.io.InputStreamRead

hdu3333区间统计

题目大意:求一个区间内不重复数字的和,例如1 1 1 3,区间[1,4]的和为4。 import java.io.BufferedInputStream;import java.io.BufferedReader;import java.io.IOException;import java.io.InputStream;import java.io.InputStreamReader;

easyui同时验证账户格式和ajax是否存在

accountName: {validator: function (value, param) {if (!/^[a-zA-Z][a-zA-Z0-9_]{3,15}$/i.test(value)) {$.fn.validatebox.defaults.rules.accountName.message = '账户名称不合法(字母开头,允许4-16字节,允许字母数字下划线)';return fal

音视频入门基础:WAV专题(10)——FFmpeg源码中计算WAV音频文件每个packet的pts、dts的实现

一、引言 从文章《音视频入门基础:WAV专题(6)——通过FFprobe显示WAV音频文件每个数据包的信息》中我们可以知道,通过FFprobe命令可以打印WAV音频文件每个packet(也称为数据包或多媒体包)的信息,这些信息包含该packet的pts、dts: 打印出来的“pts”实际是AVPacket结构体中的成员变量pts,是以AVStream->time_base为单位的显