一种基于PGP的超宽带天线的设计方法与流程

本发明属于微波电磁领域，具体涉及一种基于渐进高斯过程的超宽带天线的设计方法。

背景技术：

近年来，超宽带天线引起了人们对无线通信和传感应用的关注。由于具有宽带宽，简单结构和全向辐射方向图的特性，超宽带天线已被应用到无线通信。高斯过程是一个随机过程，近几年被广泛应用在微波器件的建模中，并且取得了较好的结果。但是对于复杂的天线往往需要大批量的训练样本去描述在工作频率范围内天线结构中不同几何和结构参数的变化，训练样本往往是通过电磁场数值计算或者测量得到的，当需要的训练样本数量过多时，电磁场的计算工作量会增大，通常会有训练输入样本获取困难，样本数量不够，导致模型的泛化能力不强的问题。

在现有技术中，常利用神经网络(ann)、支持向量机(svm)、高斯过程(gp)作为替代模型代替电磁仿真软件对微波器件进行优化设计。通过从电磁仿真软件获取训练样本对模型进行建模，利用智能优化算法对训练好的模型进行优化。

对于神经网络而言，如果遇到复杂的微波器件，往往需要大批量的训练样本去描述在工作频率范围内微波电路结构中不同几何和结构参数的变化，使得神经网络模型的建立变得非常有难度。而svm和gp也会面临因为训练样本不够，导致模型不精确，泛化能力不够的问题。

技术实现要素：

发明目的：为解决现有替代模型对复杂天线的训练样本获取困难的难点，本发明的目的在于提出一种渐进高斯过程(progressivegaussianprocess，pgp)的超宽带天线的设计方法，应用于超宽带天线中，通过增加训练样本的数量，来提高模型泛化能力。

本发明所采用的技术方案是：一种基于pgp的超宽带天线的设计方法，包括以下步骤：

s000：采用训练数据库对gp模型进行迭代训练，得到pgp模型；

s100：采用pso算法对pgp模型寻优，得到满足指标的超宽带天线尺寸参数。

进一步的，所述pgp模型的建模步骤包括：

s010：建立训练数据库；

s020：采用训练数据库中的训练样本对gp模型进行训练，直至收敛，得到已训练的gp模型；

s030：采用测试样本对已训练的gp模型和hfss模型进行预测和求解，得到对应的两种结果，采用式(1)计算得到这两种结果的相关系数r；判断相关系数r是否大于阈值，若大于，则当前gp模型即为pgp模型，并执行s100使用粒子群算法对pgp模型寻优，得到满足指标的天线尺寸；否则，进行s040；

式中，n代表频率点的个数，yi代表hfss模拟值，代表模型的预测值，代表模型的预测平均值，代表hfss平均值；

s040：随机产生若干组粒子，通过已训练的gp模型得出预测值，并计算得到其预测方差值，判断该预测方差值是否大于方差阈值，若大于，则将该预测值对应的样本带入hfss模型中进行仿真，得到的仿真结果放入s010的数据库中，并记录更新数据库的次数，执行s020。

进一步的，所述训练数据库中的训练样本的获取步骤包括：

利用正交试验设计法，选取部分正交表获取训练输入样本；

将训练输入样本送入hfss模型进行仿真，得到的仿真结果作为训练输出样本；

所述训练输入样本和训练输出样本构成训练数据库中的训练样本。

进一步的，所述s100具体包括以下步骤：

s110：随机生成初始种群，作为pgp模型的输入，通过pgp模型得到相应的输出；

s120：继承s110的输出，并根据超宽带天线的设计指标设定对应的适应度函数；

s130：通过粒子群算法更新初始种群，得到满足适应度函数的最优解；

s140：将s130中得到的最优解带入高频结构仿真进行仿真，验证是否满足设计指标。

进一步的，所述适应度函数可定义为：

fit＝min(∑|yi-vi|)

其中，vi为超宽带天线fi处要求的电压驻波比vswr的幅值，yi为fi处pgp模型的输出。

有益效果：本发明利用了高斯过程的预测方差这一特性筛选出对新样本获取有效的数据点，通过不断更新数据库来达到对模型精度的要求，提高泛化能力。本发明可以快速获得pgp模型的训练数据，克服了gp模型对复杂天线数据获取困难的问题。pgp模型训练完成以后，利用其对超宽带天线进行优化，得到满足指标的超宽带天线尺寸，pgp测试输出都较好的贴合hfss精确数值，可见本发明方法构造的渐进高斯过程确实有较高的稳定性。

附图说明：

图1是pgp的建模流程图；

图2是e型槽口的超宽带天线结构图；

图3是e型槽口的超宽带天线hfss仿真模型；

图4是优化后的e型槽口的超宽带天线的仿真结果。

具体实施方式

本发明的一种基于pgp的超宽带天线的设计方法，包括以下步骤：

s000：渐进高斯过程的建模，现对该过程详细表述如下：

s010：训练样本的获取：

利用正交试验设计法，选取部分正交表获取训练输入样本，将训练输入样本送入高频结构仿真(highfrequencystructuresimulator，hfss)中进行仿真，将高频结构仿真(highfrequencystructuresimulator，hfss)的仿真结果作为训练输出样本，将训练输入样本和训练输出样本当作初始数据库s。

s020：渐进高斯过程的构造：

s021：随机产生一组粒子通过传统gp模型和hfss进行预测和求解，对两种结果求相关系数r，如果相关系数r大于阈值，则进行s100，利用粒子群算法对gp模型进行寻优，不断迭代获取全局最优解；否则进行s022；

相关系数r的计算公式为：

n代表频率点的个数，yi代表hfss模拟值，代表模型预测值，代表模型预测平均值，代表hfss平均值；在迭代过程中，均需要对相关系数r进行计算，故此处的模型可以为传统的gp模型也可以为pgp模型；

s022：随机产生若干组粒子，通过gp模型得出预测值，对整个样本中预测方差最大的样本筛选出来代入hfss中进行仿真，将得到的结果加入到数据库s中，记录更新数据库的次数为n，然后数据库s继续训练gp模型，用同一组测试样本继续测试该gp模型，计算其与hfss的相关系数r，直到r大于阈值即可。此过程形成的模型成为渐进高斯过程。

相关概念如下：

高斯过程描述了一种函数分布，它是无限数量的随机变量组成任何子集都符合高斯分布的集合，其性质可由均值函数u(x)＝e[y(x)]和协方差函数c(x,x')＝e[(y(x)-u(x))(y(x')-u(x'))]决定，故高斯过程可定义为：

f(x)～gp(u(x),c(x,x'))(2)

式中：x,x'∈x为任意d维矢量。

假设包含n个观测值的有限数据集d＝{(xⁱ,tⁱ),i＝1,...,n}作为高斯模型的训练样本，假设观测目标值t被加性噪声ε污染，ε服从正态分布，则模型可表示为：

tⁱ＝f(xⁱ)+εⁱ,i＝1,...,nr^d→r(3)

式中：xⁱ∈x表示训练输入矢量组成的d×n维训练输入矩阵；tⁱ表示相应的n个训练输出标量组成的训练输出矢量；ε为服从正态分布的随机变量，即：

n个训练输出t和n^*个测试输出t^*组成的联合高斯先验分布是：

式中：c(x,x^*)为训练输入样本与测试输入样本之间的n×n^*阶协方差矩阵；c(x^*,x^*)是测试输入样本本身的协方差矩阵。

在测试点x^*和训练集d已给定的前提下，贝叶斯预测概率的目的是计算出概率p(t^*|d,x^*)。根据贝叶斯后验概率公式得：

式中：t^*的期望和方差分别为：

包含了和x*中测试输入向量有关的测试输出的最有可能的值，协方差矩阵给出了相应的预测方差。

高斯过程的协方差函数必须满足mercer条件：对任意一点集都能够保证产生一个非负正定协方差矩阵。本发明选择ardmatern52协方差函数作为高斯过程的协方差函数。

其中：

为信号方差。高斯过程的均值函数和协方差函数的性质由一组超参数决定，这也是高斯过程唯一需要确定的参数。

其中，c1和c2是加速常数或称作学习因子；rand()是用于产生(0，1)之间的随机数；和分别为粒子i在第k次迭代中第d维的速度和位置；是单个粒子的个体极值位置，是全体粒子的全局极值的位置。

s100：对天线的优化设计：

s110：随机生成初始种群，作为渐进高斯过程模型的输入，通过渐进高斯过程模型得到相应的输出；

s120：优化过程中，根据设计指标设定对应的适应度函数；假设设计指标在fi(i＝1,2,…)处要求的电压驻波比vswr的幅值分别为vi(v＝1,2,…)，那么适应度函数为：

fit＝min(∑|yi-vi|)(13)

其中，yi为fi(＝1,2,…)处pgp模型的输出；

s130：通过粒子群算法更新初始种群，找到满足适应度函数的最优解；

s140：将s130中得到的最优解带入hfss中验证是否满足设计指标。

为使本发明的目的、技术方案和优点更加清楚，下面结合具体实施方式，进一步阐明本发明。

实施例：

为了验证本发明用于优化设计天线的可行性，本实施例采用pgp模型对图2所示的带e型槽的超宽带(ultar-wideband,uwb)天线进行设计，图中，w0和l为介质基板的长宽，优化设计指标为优化3.1ghz-10.6ghz频率范围内天线的电压驻波比vswr小于2，其中5ghz-6ghz频率范围内vswr大于2，e型槽口的长(ls2)和宽(ls3)以及其离介质基板顶端距离(lsc)和e的宽度(w4)决定了天线的带宽和带阻的频率段等性能，其取值范围与其余固定不变的尺寸参数如表1所示。

表一基于共轭波导的带e型槽口的超宽带天线的结构参数

本实施例的优化设计步骤如下：

如图3所示，建立hfss模型，介质层的材质使用的是pcb中最常用的玻璃纤维环氧树脂(fr4)，其相对介电常数εr＝2.65，耗角正切tanδ＝0.8，通过matlab调用得到各尺寸所对应的仿真结果，得到的电压驻波比vswr。

gp模型的输入样本为x[lsc,ls2,ls2,w4]和工作频率f，其中，工作频率f的范围为2～12ghz，步长为0.1ghz。训练样本采用部分组合正交实验设计获得，一共是9*101＝909个样本，将其记为数据库s，每次更新一个样本，更新数据库次数记为n，最后的pgp模型一共有1616个样本，更新数据库7次，粒子群的参数，c1＝c2＝2，相关系数r由0.8692增加到0.9697。

将训练好的pgp模型利用pso算法寻优，得到的尺寸参数代入hfss进行验证。

通过粒子群算法对网络进行训练，最终所得的渐进高斯过程模型的相关系数为0.9697，最终得到的尺寸参数分别为[7.4596,11.6115,3.2404,0.4127]，并用测试样本对网络进行测试，结果如图4所示，其中纵坐标是vswr，实线是hfss的仿真图，虚线为pgp的预测图，两条曲线有着较好的拟合性，满足设计指标要求。