matlab 求副瓣电平,一种阵列天线一维方向图的副瓣电平计算方法与流程
本发明涉及阵列天线分析技术领域,尤其是一种阵列天线一维方向图的副瓣电平计算方法。
背景技术:
在计算阵列天线一维方向图的副瓣电平时,传统算法采用两种途径。第一种途径是先找到方向图最大值,作为主瓣峰值,然后向左右两边搜索,得到左右两边的第一零点。搜索左边零点以左区间的最大值,作为左边最大副瓣电平;搜索右边零点以右区间最大值,作为右边最大副瓣电平。取左右最大副瓣电平值的较大者,作为整个一维方向图的副瓣电平。第二种途径是搜索一维方向图的各个局部峰值,从方向图的最左边向右搜索,记录所有峰值的幅值和位置。其中最大峰值必然对应主瓣;第二峰值则对应最高旁瓣,其幅值为方向图的副瓣电平。
由于阵列天线的幅相加权问题,某些时候阵列天线的一维方向图上会出现“肩膀”,如图3所示,即方向图上出现拐点。这是由于方向图主瓣和第一副瓣甚至第二副瓣过于接近而融合在一起,主瓣和副瓣之间、相邻副瓣之间没有零陷出现。反观图2,则是标准的一维方向图,主瓣和副瓣、副瓣之间均正常分离而出现零陷。
在计算一维方向图副瓣电平的时候,如果主瓣存在拐点,则应该把拐点视为一个副瓣,并将其电平值纳入统计之中,因为拐点是主瓣和副瓣融合导致。传统算法在计算副瓣电平的时候,只是搜索方向图峰值来确定最高副瓣,因而不能识别拐点。所以目前需要一种新的副瓣电平计算方法,既能处理常规的方向图,也能处理带有拐点的方向图,并进一步降低运算量。
技术实现要素:
本发明所要解决的技术问题在于,提供一种阵列天线一维方向图的副瓣电平计算方法,能够很好处理方向图主瓣和第一副瓣融合在一起的情况,使副瓣电平计算得更加可靠和准确。
为解决上述技术问题,本发明提供一种阵列天线一维方向图的副瓣电平计算方法,包括如下步骤:
(1)针对阵列天线一维方向图E,搜索方向图最大值,即主瓣峰值,记为E(p),此处p表示主瓣峰值的坐标;为方向图E建立对应的差分数组D,E与D长度都为N,相同位置的方向性值和差分值具有对应关系,将E(p)对应的差分值D(p)置为0;
(2)搜索主瓣峰值左边第一个零点或者拐点,其坐标记为lC,然后搜索lC点左边区间(1~lC)中的最大值,作为主瓣左边区域的副瓣电平值,记为SL;
(3)搜索主瓣峰值右边第一个零点或者拐点,其坐标记为rC,然后搜索rC点右边区间(rC~N)中的最大值,作为主瓣右边区域的副瓣电平值,记为SR;
(4)将SL和SR进行比较,选择其中较大的一个,作为整个一维方向图的副瓣电平,作为结果进行输出。
优选的,步骤(2)中,所述的左边第一零点或者拐点搜索,其步骤如下:
(21)计算主瓣峰值左边相邻点的差分,计算方法如下:
D(l1)=E(l1)-E(p)
此处l1表示主瓣峰值左边第一个离散坐标;
(22)计算主瓣峰值左边第二个点的差分:
D(l2)=E(l2)-E(l1)
此处l2表示l1左边的相邻坐标;
(23)将D(l2)与D(l1)进行比较;如果D(l2)小于等于D(l1),则将l2和l1均减1,然后进入第(22)步,继续左向搜索;如果D(l2)大于D(l1),则停止左向搜索,将此时的l2记为lC。
优选的,步骤(3)中,所述的右边第一零点或者拐点搜索,其步骤如下:
(31)计算主瓣峰值右边相邻点的差分,计算方法如下:
D(r1)=E(r1)-E(p)
此处r1表示主瓣峰值右边第一个离散坐标;
(32)计算峰值右边第二个点的差分:
D(r2)=E(r2)-E(r1)
此处r2表示r1右边的相邻坐标。
(33)将D(r2)与D(r1)进行比较3;如果D(r2)小于等于D(r1),则将r2和r1均加1,然后进入第(32)步,继续右向搜索;如果D(r2)大于D(r1),则停止右向搜索,将此时的r2记为rC。
本发明的有益效果为:本发明在计算阵列天线一维方向图的幅瓣电平时,所采用的双向差分算法能基于最小的运算量,精确地找到零点或拐点的位置,进一步计算出方向图的幅瓣电平;算法实现的复杂度较低,计算量不大,需要的内存消耗很小;与传统算法相比,本发明能够很好处理方向图主瓣和第一副瓣融合在一起的情况,即识别主瓣中的拐点,使副瓣电平计算得更加可靠和准确。
附图说明
图1为本发明的方法流程示意图。
图2为采用传统方法和双向差分算法对标准的一维方向图进行副瓣电平计算的结果示意图。
图3为采用传统方法和双向差分算法对带有拐点的一维方向图进行副瓣电平计算的结果。
具体实施方式
如图1所示,本发明提出的阵列天线一维方向图副瓣电平计算方法,即双向差分算法。首先找到方向图的主瓣峰值,然后分别向左和向右搜索差分变大的位置,将此两处作为左右两边的零点或拐点,二者之间即为方向图主瓣。搜索主瓣左右两边区间的最大值,作为整个一维方向图的副瓣电平。该算法的特点是根据方向图相邻点之间的差分来判断方向性的变化,不仅可以找到峰值和零点,还可以发现拐点。
具体步骤如下:
步骤1):搜索一维方向图主瓣峰值
针对阵列天线一维方向图E,搜索方向图最大值,即主瓣峰值,记为E(p),此处p表示主瓣峰值的坐标。为方向图E建立对应的差分数组D,E与D长度都为N,相同位置的方向性值和差分值具有对应关系。将E(p)对应的差分值D(p)置为0。
步骤2):左向搜索零点或者拐点
(1)计算主瓣峰值左边相邻点的差分,计算方法如下:
D(l1)=E(l1)-E(p)
此处l1表示主瓣峰值左边第一个离散坐标。由于方向图主瓣的值向左不断下降,所以D(l1)为负数。
(2)计算主瓣峰值左边第二个方向图值的差分:
D(l2)=E(l2)-E(l1)
此处l2表示l1左边的相邻坐标。
(3)将D(l2)与D(l1)进行比较。如果D(l2)小于等于D(l1),则将l2和l1均减1,然后进入第(2)步,继续左向搜索;如果D(l2)大于D(l1),则说明l2这个位置是一个零点或者拐点,此时可以停止左向搜索,进入下一步。
如果l2处不是拐点或者零点,则仍然处于方向图主瓣的加速下降区,D(l2)和D(l1)均为负数,且D(l2)的绝对值大于等于D(l1)的绝对值。此时必有D(l2)小于等于D(l1)。
在遇到拐点的时候,D(l2)和D(l1)均为负数,只不过D(l2)的绝对值会小于D(l1)的绝对值。在遇到零点的时候,D(l1)为负数,D(l2)有可能为正数,即方向图由加速下降突然转变为上升。这两种情况下,均必有D(l2)大于D(l1)。
步骤3):计算主瓣左边副瓣电平
将此时的l2记为lC,表示主瓣峰值左边的拐点或零点。搜索lC点左边方向图的最大值,作为主瓣左边区域的副瓣电平值,记为SL。搜索范围为:E(l),0≤l≤lC。
步骤4):右向搜索零点或者拐点
(1)计算主瓣峰值右边相邻点的差分,计算方法如下:
D(r1)=E(r1)-E(p)
此处r1表示主瓣峰值右边第一个离散坐标。由于方向图主瓣的值向右不断下降,所以D(r1)为负数。
(2)计算峰值右边第二个方向图值的差分:
D(r2)=E(r2)-E(r1)
此处r2表示r1右边的相邻坐标。
(3)将D(r2)与D(r1)进行比较。如果D(r2)小于等于D(r1),则将r2和r1均加1,然后进入第(2)步,继续右向搜索;如果D(r2)大于D(r1),则说明r2这个位置是一个零点或者拐点,此时可以停止右向搜索,进入下一步。
如果r2处不是拐点或者零点,则仍然处于方向图主瓣的加速下降区,D(r2)和D(r1)均为负数,且D(r2)的绝对值大于等于D(r1)的绝对值。此时必有D(r2)小于等于D(r1)。
在遇到拐点的时候,D(r2)和D(r1)均为负数,只不过D(r2)的绝对值会小于D(r1)的绝对值。在遇到零点的时候,D(r1)为负数,D(r2)有可能为正数,即方向图由加速下降突然转变为上升。这两种情况下,均必有D(r2)大于D(r1)。
步骤5):计算主瓣右边副瓣电平
将此时的r2记为rC,表示主瓣峰值右边的拐点或零点。搜索rC点右边方向图的最大值,作为主瓣右边区域的副瓣电平值,记为SR。搜索范围为:E(r),rC≤r≤N。
步骤6):结果输出
由以上分析可知,0~lC为左边副瓣区,lC~rC为主瓣区,rC~N为右边副瓣区。SL为左边副瓣区最高副瓣电平,SR为右边副瓣区最高副瓣电平。将SL和SR进行比较,选择其中较大的一个,作为整个一维方向图的副瓣电平,记为SALL,作为结果输出。
在将SALL作为结果输出时,也可以输出其他一些重要参数,包括差分数组D,主瓣零点宽度rC-lC等,对分析阵列天线方向性具有参考价值。
实施例1:本例进行阵列天线一维方向图的副瓣电平计算。天线方向图如图2所示,是一个标准的一维方向图,主瓣指向为θ0=90°,角度θ从-90°到90°。不论是本发明提出的双向差分算法,还是传统的副瓣电平计算方法,都能很准确地计算出副瓣电平为-30dB。运算量统计表明,双向差分算法的运算量小于传统算法。
实施例2:本例进行阵列天线一维方向图的副瓣电平计算。天线方向图如图3所示,主瓣指向为θ0=90°,角度θ从-90°到90°。方向图中存在主瓣和第一、第二副瓣融合在一起的情况,所以主瓣出现了拐点,表现为两级“肩膀”。传统算法由于不能识别拐点,所以算出副瓣电平为-30dB;双向差分算法能识别拐点,故能准确算出副瓣电平为-14dB。运算量统计表明,双向差分算法的运算量小于传统算法。
尽管本发明就优选实施方式进行了示意和描述,但本领域的技术人员应当理解,只要不超出本发明的权利要求所限定的范围,可以对本发明进行各种变化和修改。
matlab 求副瓣电平,一种阵列天线一维方向图的副瓣电平计算方法与流程相关推荐
- 振子天线三维方向图 matlab仿真,1阵列天线方向图的MATLAB实现
阵列天线方向图的MATLAB实现 第25卷第1期 2010年6月西藏大学学报(自然科学版)JOURNALOFTIBETUNIVERSnYV01.25No.1Jun.2010 阵列天线方向图的MATLA ...
- 微分方程求解matlab冲激信号,用Matlab求冲激响应的几种方法
·实验技术与方法· 用 Matlab 求冲激响应的几种方法 贺富堂,应柏青,张 锋,孙 敏 ( 西安交通大学 电工电子教学实验中心,陕西 西安 710049) 摘 要: 用 Matlab 软件求冲激响 ...
- 相控阵天线(一):直线阵列天线特性和阵列因子(方向图乘积定理、波束扫描、含python代码)
目录 方向图乘积定理 阵列因子方向图 波束扫描 阵列方向图和单元方向图 方向图乘积定理的python代码示例 方向图乘积定理 任意形式单元天线构成的直线阵如下图所示: 阵中第n个单元的远区辐射场可表示 ...
- 天线阵列计算远场方向图
天线阵列计算远场方向图 通过对一维或二维相位梯度阵列的自定义排布,以此计算该系统的远场电磁波辐射特性. 简介 阵列天线顾名思义,是指由若干个天线排列组成的辐射器.和单元天线不同,它可以形成更多.更复杂 ...
- 什么是天线的方向图?
天线方向图又叫辐射方向图 ( radiation pattern ) .远场方向图 ( far-field pattern ) .从方向图上面不能得到天线增益,由方向图得到的是方向系数.天线增益=方向 ...
- 混合波束成形| 通过天线空间方向图理解波束成形的物理意义
波束成形的物理意义 如图, 是在各种教材中经常看到的天线方向图. 上图表示的是当前天线经过波束成形后在空间中 指向 30度方向 (一般考虑的方向是0-180度).这里解释一下天线方向: 不失一般性的, ...
- 泰勒阵列天线综合与matlab,阵列天线综合之切比雪夫低副瓣阵列设计Matlab
在<自适应天线与相控阵>这门课中,我了解到了关于理想低副瓣阵列设计的一些方法,其中切比雪夫等副瓣阵列设计方法是一种基础的方法,故将其设计流程写成maltab程序供以后学习使用.在此分享一下 ...
- 用matlab计算稳态误差,利用Matlab求稳态误差的两种方法.
利用Matlab求稳态误差的两种方法 摘要:稳态误差是系统控制精度或抗扰动能力的一种度量,它是稳态性能的一个重要指标.本文介绍利用Matlab的控制系统工具箱和Simulink工具箱求取系统误差稳态的 ...
- 【matlab 圆周率计算】matlab 求圆周率的两种算法实现比较
%author:laidefa %data:2014-09-19 %丘德诺夫斯基公式求圆周率 function mpi=qdnfsj(m) i=m; s=13591409; for n=1:i A=( ...
- matlab模式识别提取特征向量,一种基于小波特征向量提取的手机检测方法与流程...
本发明涉及到手机检测领域,尤其涉及到一种基于小波特征向量提取的手机检测方法. 背景技术: 随着保密要求的不断提高,很多场合严禁携带手机.录音笔.录像机等电子产品,亟需一种设备可以检测出该类电子产品.目 ...
最新文章
- 边缘计算 — 困境与破局
- maven 常用插件3
- 设置函数环境——setfenv
- Vue中集成高德地图API实现定位与自定义样式信息窗体
- secondarynamenode异常
- win 2016 ssh_如何使用BPF将SSH会话转换为结构化事件
- ERROR 1130: Host '192.168.1.3' is not allowed to c
- 通过fiddler和逍遥模拟器模拟抓包android手机
- 在你的代码中使用Boost智能指针
- 操作键盘事件源码解析(常用的鼠标事件、 键盘事件对象之keyCode属性)
- windows和android结合,Android和Windows 10可以很好地合作的10种方式
- 机器学习实战——3.4 示例:使用决策树预测隐形眼镜类型
- 生活美学 | 8种咖啡冲煮器具分别有什么特点
- 有趣吧–中国最大的在线制作漫画社区 开站了
- 访问学者博士后面签后的几种情况?
- DTI-ATS入门(1):DTI综述
- DELPHI获取网卡MAC地址
- [pwn]格式化字符串:0ctf 2015 login writeup
- u盘上1T linux系统部署到4T硬盘,bios使用Legace启动(实践可用)
- A Novel Proof-of-Reputation Consensus for Storage Allocation in Edge Blockchain Systems 精读笔记(三)
热门文章
- Spring源码解析(一)
- 谷粒商城:秒杀系统设计与编写
- 公文流转 php,河北金力集团公文流转系统----节选
- 程序员考证书,有用吗?
- python rgb565_读取RGB565格式的图像
- 小趣味:js编写斗地主规则(完整代码)
- Java读取计算 PPT,Word,excel的页数
- 宇视服务器硬件如何安装,宇视科技无需后端平台与服务器支撑 即可形成小型人脸识别方案...
- gflags 调试内存_gflags 检查内存越界 | 学步园
- google gflags 库完全使用