1.   引言

无源互调^[1-3](Passive InterModulation, PIM)是指在大功率条件下，当输入两个或者两个以上载波时，由于微波部件的非线性导致载波信号相互调制，产生新的载波频率落入接收通带内造成干扰的现象。正是由于微波部件的非线性，通过微波部件的载波频率出现了多种不同的组合频率，即不同阶数的互调。由于无源互调的功率随阶数的增加而减小，因此较低阶互调如3阶和5阶更容易引起干扰。然而，当通过微波部件的载波信号功率较大时，一些高阶互调的功率也会很大，足以引起严重干扰。在通信系统中，为了应对日趋紧张的频谱资源，提升系统通信容量，多载波或多系统共存现象将进一步增加。这就使得同一传输信道内会有多个不同频率的信号，这些基频信号由于各种非线性特性的影响而产生新的频率分量。如果这些互调信号恰好落在接收频段内，就会对通信系统产生干扰，进而降低传输系统的接收灵敏度，严重时导致通信系统崩溃。尤其，随着卫星通信技术的快速发展，空间通信系统对大功率发射和高灵敏度接收也提出了更高要求。因此，无论对于地面移动通信还是空间卫星通信，无源互调都是制约通信系统性能提高的亟待解决的关键问题之一。

目前，国际上普遍认为无源互调是由通信系统中无源微波部件的非线性^[4-10]特性引起的。有关无源互调的研究也多从其产生机理、预测分析、抑制方法、测量技术等方面展开。在无源互调预测分析方面，常用的方法有：等效电路法、幂级数法、MMicroscope法、Volterra级数法等。其中，代表性的研究有：西班牙的Vicente等人^[5]提出使用等效电路法分析预测矩形波导法兰连接结的无源互调。在此基础上，西安交通大学的Chen等人^[6]使用等效电路模型方法对同轴连接器的无源互调进行了预测分析。西安电子科技大学的张世全等人^[8]基于幂级数方法进行低阶无源互调功率预测高阶无源互调功率的研究，空军工程大学的王海宁等人^[9]提出使用IM Microscope 方法预测无源互调功率。以上这些预测研究普遍认为，在输入条件不变的情况下，对于固定阶数的无源互调功率为一定值。然而，实际研究中发现，无源互调功率并非一个定值，而是随时间变化的。北京理工大学的田露^[10]基于Volterra级数法建立的无源互调预测模型中考虑到了时间因素，并进行了无源互调对消抑制技术研究，但Volterra级数法的准确性受到模型参数数量的影响，当参数数量增多时，计算过程也会变得极其复杂。同时，这些预测方法的共同点是，都需要通过已有数据建立相应的主观模型，并根据这些人为设定的主观模型进行预测。而混沌方法不用预先建立主观模型，只是根据数据本身的特性规律进行预测，这样就解决了由于人为主观建模导致模型的不准确，以至造成预测的误差，达到完全用客观数据预测时的可信度。

混沌^[11-13]是研究系统非线性的一种有效方法，它是指在确定性系统中出现的一种貌似无规则的、类似随机的现象，是非线性系统的一种新的存在形式，并且已在多个领域得到广泛应用。其中，混沌电路的研究实现了对复杂非线性现象和混沌理论从数学模型到硬件电路研究的转变，为混沌通信的研究和应用奠定了基础。利用混沌信号的类随机性、连续宽带功率谱特性，研究人员将混沌应用于通信系统中的多路接入技术、数据加密、编码方法、调制技术等^[14,15]多个层次，获得比传统通信更好的通信性能。由于产生无源互调的根源是微波部件的非线性，因此可以将该微波部件看作一个非线性系统来处理。基于此，本文拟从非线性系统的整体角度出发，研究无源互调功率与时间的关系。以通信系统中最常用的同轴连接器作为待研究的典型微波部件，以通信系统中最常见的3阶无源互调为分析对象，通过实验方法得到其3阶无源互调功率随时间的变化关系，获得的实验数据可以看作无源互调功率的时间序列。基于获得的无源互调功率时间序列，分别从定性和定量角度对其进行混沌特性识别。在验证了无源互调功率时间序列具有混沌特性的基础上，基于最大Lyapunov指数^[16-18]方法对无源互调功率时间序列进行预测分析研究。

2.   PIM功率时间序列

由于通过微波部件的载波信号都是随时间变化的，因此由这些载波信号经过微波部件的非线性调制而产生的不同阶数的PIM信号也是随时间变化的。研究发现，当输入条件固定时，对确定阶数的PIM信号，PIM功率并不是一个恒定值，而是随时间变化的。因此，PIM功率具有时间特性。如果将随时间变化的PIM功率看作时间序列，则该时间序列不仅包含该微波部件过去的PIM信息，还包含了PIM演化的大量信息。

为了获得PIM功率时间序列，本文搭建了一个S频段的PIM测试系统(如图1所示)。该测试系统的发射通带为2.16～2.21 GHz，灵敏度达到–160 dBm/Hz，可对3～15阶的PIM功率进行测量。待测件端口的最大功率可达200 W，可以满足不同输入功率条件下的PIM特性研究。鉴于同轴连接器是通信系统中最常用的典型微波部件，因此，本文以DIN(7/16)型同轴连接器作为产生PIM的源，以通信系统中最为关注的3阶PIM为研究对象，研究PIM功率与时间的依赖关系。实验过程中的输入信号频率分别为2.16 GHz和2.21 GHz，输入功率均为20 W，实验测试数据为3阶PIM功率。测试过程中，记录采集前5000 s的3阶PIM功率，测试结果如图2所示，该结果即为PIM功率的时间序列。从图2可以看出，PIM功率的变化具有时变性，同时也具有一定的随机性和不确定性。如果能够证明该时间序列具有混沌特性，则可以运用混沌方法对PIM功率进行预测分析。

图 1  PIM测试系统及待测样件

下载: 全尺寸图片 幻灯片

图 2  3阶PIM功率时间序列

下载: 全尺寸图片 幻灯片

3.   PIM功率时间序列混沌特性识别

使用混沌方法对PIM功率时间序列进行分析研究之前，需要先判断该时间序列是否具有混沌特性。判断一个系统的动态行为是否混沌，就是判断该系统是否具有混沌吸引子。定性分析方法主要是根据时间序列在时域或频域内表现出的特殊性质对序列的主要特征进行粗略分析。常用的有相图法、功率谱法、庞加莱截面法和代替数据法等。定量分析的方法主要是对描述混沌系统的重要特性指标包括关联维数、最大Lyapunov指数和Kolmogorov熵等特性指标定量分析，从而进行混沌识别。本文分别选用相图法和最大Lyapunov指数法，从定性和定量角度对PIM功率时间序列进行混沌识别。

3.1   确定嵌入维数和延迟时间

3.1.1   相空间重构理论

通常对于大多数实际工程问题来说，并不知道其具体的数学模型，能获得的信息也比较有限，一般只能获得1维的时间序列信号。为了从一个变量时间序列获得其他变量的信息，有研究人员提出时间延迟坐标的思想，重构出观测到的动力学系统的相空间。由于相空间重构^[19-21]的结果可以保留原系统的动力学性质和几何性质，因此，相空间重构是非线性时间序列分析的重要步骤。其具体原理如下。

假设时间序列为{x(t_i), i=1, 2, ···,N}，其采样时间间隔为$\Delta t$。通过时间延迟可构成m维的向量：X(t_i)=(x(t_i), x(t+τ), ···, x(t_i+(m–1)τ)), i=1, 2, ···, M，其中m表示嵌入维数，τ表示时间延迟，X(t_i)为m维相空间中的相点，M为相点个数，且有M=N–(m–1)τ 得到的集合{X(t_i), i=1, 2, ···, M}描述了系统在相空间中的演化轨迹。只要m, τ选择合适，重构的相空间与原系统就具有相同的拓扑性质。因此，如何确定时间序列恰当的m和τ，对时间序列的研究至关重要。

3.1.2   C-C方法确定嵌入维数和延迟时间

关于延迟时间τ和最佳嵌入维数m的确定方法，国内外研究人员已经进行了大量的研究。本文选用C-C方法对PIM功率时间序列进行分析。C-C方法是由Kim等人^[21]提出的，其优点在于可通过关联积分对时间延迟τ和延迟时间窗口τ_W进行同时估算。具体步骤为

依据重建的相空间中的相点，

获得嵌入时间序列的关联积分为

其中，m为嵌入维数，N是时间序列的数据个数，r为计算中所取的搜索半径，θ为Heaviside函数

${d_{ij}} = {\left\| {{x_i} - {x_j}} \right\|_{\left( \infty \right)}}$表示${\infty - }$范数。该关联积分表示相空间中任意两点之间距离小于r的概率。此处，点与点之间的距离用矢量之差的${\infty - }$范数表示。

定义序列x={x_i}的检验统计量为

将时间序列平均分成t个互不相交的子序列，即有

根据式(4)，有

${S_1}\left( {m,r,t} \right) \sim t$反映了时间序列的自相关特性。选择最大和最小的两个半径r，并定义差量

$\Delta {S_1}\left( {m,t} \right)$表示${S_1}\left( {m,r,t} \right) \sim t$对所有半径r的最大偏差。进而得到计算公式^[16]

因此，最优时间延迟τ，即为$\overline {{S_1}} \left( t \right) \sim t$的第1个零点(见图3(a))或$\Delta \overline {{S_1}} \left( t \right) \sim t$的第1个局部极小值点(见图3(b))，其表示重构吸引子轨道在相空间完全展开。综合考虑$\overline {{S_1}} \left( t \right)$与$\Delta \overline {{S_1}} \left( t \right)$，定义

图 3  C-C方法计算的最优时间延迟和最优延迟时间窗口

下载: 全尺寸图片 幻灯片

${S_{1\_{\rm{cor}}}}(t) \sim t$的全局最小点即为最优延迟时间窗口τ_W(见图3(c))，其表示了平均轨道周期的最优估计。

本文使用C-C方法对获得的同轴连接器的3阶PIM功率时间序列进行分析，其计算过程如图3所示。从图3(a)和图3(b)可以看出，第1个零点和第1个极小值点的位置都在时间点为τ=6，因此可知最佳延迟时间τ为6。从图3(c)可以看到，全局极小值在时间为20处取得，因此最优延迟时间窗口τ_W为20。根据公式${\tau _W} = \left( {m - 1} \right)\tau$，可得嵌入维数$m \approx 4$。

3.2   混沌特性识别

3.2.1   相图法

相图法将通过绘制相空间重构后点的轨迹分布规律作为混沌特性判断的依据。如果所分析的时间序列具有混沌特性，则系统的相空间轨迹通常表现为在有限空间内不断伸长和折叠形成的回复性永不相交的非周期运动，该运动不同于毫无规律的随机运动，但也不是周期函数的重复性运动，即相图出现混沌吸引子。由于该方法具有简单、直观的特性，已广泛用于系统混沌特性的识别中。图4为基于本文研究的同轴连接器的3阶PIM功率时间序列而建立的相空间重构的3个分量构成的3维相图。从图中可以看出，相点的轨迹并不是杂乱无章的，存在混沌吸引子，从而表明了PIM功率时间序列具有混沌特性。

图 4  PIM功率时间序列3维相空间图

下载: 全尺寸图片 幻灯片

3.2.2   最大Lyapunov指数法

Lyapunov指数表征了混沌系统相空间中两条相邻轨迹之间收敛或者发散的平均指数率，是宏观层次上对混沌吸引子的一个表述，也是定量判断一个系统是否具有混沌特性的重要指标。如果所研究的时间序列的最大Lyapunov指数大于0，则证明该时间序列具有混沌特性。并且，Lyapunov指数越大，说明混沌特性越明显，即混沌程度越高。目前，用于计算最大Lyapunov指数的常用方法有定义法、wolf法、正交法和小数据量法。本文选用小数据量法进行求解，其计算原理如下。

假定已知混沌时间序列重构后的相空间由式(1)所示，通过欧氏距离最小值求得给定点上各个点的最相邻点

其中，$\hat t{\rm{ = }}1{\rm{,}}2{\rm{,}} ··· {\rm{,}}M$，且$t \ne \hat t$, $X{\rm{(}}t{\rm{)}}$为相空间中的任一点，${T_m}$为时间序列的平均周期，${d_t}{\rm{(0)}}$表示一组最近邻域值在初始时的欧氏距离。

应用邻近点的发散的平均值来预估最大Lyapunov指数为

其中，$\Delta t$表示样本周期，i为离散时间，${d_t}{\rm{(}}i{\rm{)}}$表示在经过i段时间步长后两邻近对之间演化的距离值。结合Sato等人^[18]的估计式有

对式(13)进行化简

由式(14)可得直线的斜率值就是λ₁的数值。采用最小平方法，求出相应的回归直线的斜率值，即

其中，q为${d_t}{\rm{(}}i{\rm{)}}$的个数。

通过对本文获得的同轴连接器的3阶PIM功率时间序列，使用小数据量法进行分析(计算结果如图5所示)，可得其最大Lyapunov指数为0.0075。根据混沌序列的判定标准，由于分析的PIM功率的时间序列的最大Lyapunov指数大于0，因此，该时间序列具有混沌特性。

图 5  小数据量法计算的最大Lyapunov指数

下载: 全尺寸图片 幻灯片

4.   基于最大Lyapunov指数的PIM功率时间序列预测

经过前面的分析可知本文研究的同轴连接器的3阶PIM功率时间序列具有混沌特性。下面将结合获得的最大Lyapunov指数，对该时间序列进行进一步的PIM预测^[22,23]研究。

假定时间序列的最大Lyapunov指数为λ₁, Y_M为预测的中心点，Y_M的最近临近点为Y_K，两者之间的欧氏距离为${d_M}{\rm{(0)}}$，则有

Y_M和Y_K经一步演化后分别成为Y_M+1和Y_K+1，根据最大Lyapunov指数的物理含义，有

Y_M, Y_K, Y_K+1是相空间中的点，Y_M+1中只有1个要预测的PIM功率是未知的，设为Y_M+1,m，令

则有

图6为基于最大Lyapunov指数得到的PIM功率预测结果。图中对原时间序列的5000 s之后的100 s的PIM功率进行预测。从图中可以看出，理论预测值与实验值吻合很好，其最大误差为2.61%。

图 6  前100 s的预测结果

下载: 全尺寸图片 幻灯片

图7为对原时间序列的5000 s之后的200 s的PIM功率的预测结果。从图中可以看出，理论预测值与实验值在前137个点吻合很好，误差较小且稳定，而在138 s之后误差变大，最大误差为6.71%。这是因为混沌时间序列分析存在最大可预测尺度，其定义为

从式(20)可知，${\lambda _1}$越大，最大可预测尺度越小，即混沌特性越明显，可准确预测的点的个数就越少。通过式(20)计算得到本文分析PIM功率时间序列的最大可预测尺度为T_max=1/0.0075=133.33。也就是说，该PIM功率时间序列的预测只在前134个点具有较好的预测结果，该结论与图7的分析结果基本一致。

图 7  前200 s的预测结果

下载: 全尺寸图片 幻灯片

5.   结束语

无源互调是影响通信系统性能的重要因素，而微波部件自身的非线性特性是导致通信系统产生无源互调的根源。鉴于混沌是研究系统非线性的一种有效方法，并且在研究中发现，无源互调也具有混沌特征。因此，本文以通信系统中的典型微波部件——同轴连接器为对象，以通信系统中最常见的3阶无源互调为分析对象，验证微波部件无源互调的混沌特性，并基于混沌理论进行无源互调功率预测。首先，通过实验系统测试得到同轴连接器的3阶无源互调，获得无源互调功率时间序列。其次，通过相图法和最大Lyapunov指数法，分别从定性和定量角度验证了无源互调功率时间序列具有混沌特性。最后，基于最大Lyapunov指数法，对不同时间的无源互调功率进行预测。通过实验对比可知，在最大可预测尺度范围内，预测结果与实验结果吻合较好，最大误差为2.61%。在最大可预测尺度范围外，则会产生较大误差。本文所讨论的无源互调功率预测方法，对研究无源互调功率与时间的关系，以及为开展无源互调的抑制技术研究提供了一种新的思路。