-
随着现代通信系统[1]的发展,如码分多址[2](code division multiple access, CDMA)技术和正交频分多路复用[3](orthogonal frequency division multiplexing, OFDM)技术的应用,使得通信信号变得越来越复杂,导致微波元件中记忆效应[4]越来越突出。这种记忆效应主要表现在交调分量的幅度和相位的不对称性及频率依赖性。交调分量的这些特性引起了微波工程师的注意,在面对三阶调制分量(the third intermodulation, IM3)的频率依赖性和非对称性时,工程技术人员感到困惑,因为以前的理论认为IM3分量的下边带和上边带的振幅是对称的[5]。
已有研究对于这些微波元件的交调分量的幅度不对称性和频率依赖性进行了解释。文献[6-7]将这种频率依赖性和幅度不对称性定义为微波元件的记忆效应。文献[8]对IM3分量上下边带的不对称性影响进行了深入研究,并提出了相应的数学表达式。文献[9-10]利用时变增益调制函数(time varying gain modulation function, TVGMF)描述了长时基带记忆效应。然而,这些研究仅关注了放大器的记忆效应,而对于混频器的记忆效应则少有研究。文献[11]利用时变幂级数和沃尔特拉(Volterra)级数对CMOS混频器的IM3分量上下边带相位的不对称性进行了分析。文献[12]利用维纳(Wiener)模型表征了IM3分量上下边带相位的不对称性特征。
近年来,混频器的行为模型研究取得了进步。它们主要包括转换矩阵(conversion matrix, CM)模型[13]、散射参数混频器模型[14]、静态X参数混频器模型[15]、广义沃尔特拉级数(generalized volterra series, GVS)模型[16]以及多盒行为混频器模型(multi-box behavioral mixer, MBBM)[12]等。在这些模型中,GVS模型和MBBM模型可以表征混频器的记忆效应。动态X参数模型[4]虽然可以表征微波元件的记忆效应,但该模型目前仅用于表征放大器的记忆效应。文献[17]表明虽然测量和分析交调(intermadulation, IM)分量的振幅是双音测试中最常见的,但IM分量的相位测量和分析也引起了关注,并被证明在某些情况下其在记忆效应方面比IM分量的振幅更加敏感。GVS模型和MBBM模型虽然能够表征混频器的记忆效应,但不能直观地描述混频器交调分量的上下边带的相位具有相反的变化趋势这一特征。除此之外,这些模型也不能解释交调分量上下边带的相位具有相反变化趋势的原因。
本文的目的是通过直观的描述来解释IM3分量上下边带相位不对称性的原因,并对这些造成相位不对称性的因素进行区分与识别。本研究是在时变调制函数(time varying modulation function, TVMF)和MBBM的基础上进行的。结合上述两种理论方法对三阶交调分量上下边带相位不对称的原因做出了统一解释。该项工作的难点在于对混频器IM3分量相对相位的准确测量,这是由于在初始时刻难以确定的情况下,比较不同频率的信号相位没有意义。传统的测量方法[18-20]是在测试系统中加入一个参考通路,但这就增加了测量的复杂度。在本文中,由于R&S新一代矢量网络分析仪的内置源是相位相参的,所以在测试时,不再需要外接参考链路,只需利用未知直通校准即可完成IM3分量相对相位的测试,这是传统方法不具有的优势。最后,在仿真分析部分通过先进设计系统(advanced design system, ADS)验证了理论的正确性,在实验分析部分借助R&S公司新一代矢量网络分析仪强大的交调分量测量能力,对多个商用混频器的三阶交调分量的相位进行了测量,从而展示了这种相位不对称性现象。
-
文献[12]描述的多盒行为混频器模型结构如图1所示。在混频器的射频(radio frequency, RF)输入端,利用放大器的Wiener模型来描述混频器的非线性失真和记忆效应。在该模型中,放大器的Wiener模型首先按功能被分为两个部分:线性滤波器和加权多项式。这两部分被分别用来表征混频器射频端口的记忆效应和非线性特性。其次,将非线性部分再次分解为幅度失真(AM-AM)和相位失真(AM-PM)两个部分。对于混频器的本振(local oscillator, LO)端口,由于本振信号没有调制,忽略相位失真(AM-PM),只使用一个实数多项式来表征本振端口的幅度失真。通过这些修改,文献[12]中的多盒行为混频器模型的功能模块包括:用于表征记忆效应的线性滤波器模块,用于表征中频(intermediate frequency, IF)相位失真的相移多项式模块,以及用于表征射频和本振端口幅度失真的实数多项式模块。最后,将这些模块和一个理想混频器进行级联,从而得到了多盒混频器行为模型。
多盒行为混频器模型可以表述为:
$$ \begin{split} &\qquad\qquad\qquad\qquad {V_{{\rm{IF}}}}(t) =\\ &\left[ \begin{gathered} \sum\limits_{i = 1}^l {{r_i}\int\limits_{ - \infty }^{{\text{ + }}\infty } {\int\limits_{ - \infty }^{{\text{ + }}\infty } \cdots {\int\limits_{ - \infty }^{{\text{ + }}\infty } h } } (\tau )} \prod\limits_{i = 1}^n {{v_{{\rm{RF}}}}(t - {\tau _i})} {\rm{d}}{\tau _1} {\rm{d}}{\tau _2} \cdots {\rm{d}}{\tau _n} \times \\ \sum\limits_{j = 1}^J {{l_j}{{({v_{{\rm{LO}}}})}^j}} \\ \end{gathered} \right] \end{split}$$ (1) 式中,
$ {r_i} $ 和$ {l_i} $ 为实非线性系数;$ h(\tau ) $ 为具有记忆效应的线性滤波器。该模型无法直观描述混频器IM3分量上下边带相位不对称性机理。因此,需对IM3分量相位不对称性二阶记忆效应机理进行直观分析。 -
根据文献[12],放大器的Winner模型可以用来描述混频器的记忆效应和非线性。在文献[21]中,TVGMF用于识别和区分放大器的IM3相位不对称性的原因。首先对放大器的TVGMF进行了修正,得到TVMF。最后,基于多盒混频器模型的模块化构建思路,利用TVMF来替换多盒混频器模型中代表Winner模型功能的模块,从而用于分析混频器IM3分量上下边带相位不对称性的来源。这种相位不对称性表现在IM3分量的上边带和下边带随着双音间距的变化具有相反的变化趋势。一种新的基于TVMF多盒行为混频器特征描述结构如图2所示。
-
将等幅双音信号定义为:
$$ {V_{{\rm{RF}}}}(t) = {V_0}\cos ({\omega _1}t + {\varphi _1}) + {V_0}\cos ({\omega _2}t + {\varphi _2}) $$ (2) 式中,
$ {V_0} $ 、$ {\omega _1} $ 和$ {\varphi _1} $ 分别代表第一个输入信号的幅度、频率和初相;$ {V_0} $ 、$ {\omega _2} $ 和$ {\varphi _2} $ 分别代表第二个输入信号的幅度、频率和初相。利用三角恒等式,可以将双音输入信号重新表述为式(3),它由基带频率为$ {\omega _m} $ 的包络信号$ a(t) $ 以及频率为$ {\omega _c} $ 的载波信号组成,$ {\varphi _m} $ 和$ {\varphi _c} $ 分别表示包络信号和载波信号的初始相位:$$ {V_{{\rm{RF}}}}(t) = a(t)\cos ({\omega _c}t + {\varphi _c}) $$ (3) 其中:
$$ a(t) = 2{V_0}\cos ({\omega _m}t + {\varphi _m}) $$ (4) $$ {\omega _m} = \frac{{{\omega _2} - {\omega _1}}}{2} $$ (5) $$ {\omega _c} = \frac{{{\omega _2} + {\omega _1}}}{2} $$ (6) $$ {\varphi _m} = \frac{{{\varphi _2} - {\varphi _1}}}{2} $$ (7) $$ {\varphi _c} = \frac{{{\varphi _2} + {\varphi _1}}}{2} $$ (8) 二阶TVMF定义为:
$$ H\left( t \right){\text{ = }}{c_0} + {h_2}(t) \otimes V_{{\rm{RF}}}^2(t) $$ (9) 式中,
$ {h_2}(t) $ 为与二阶输入信号$ V_{{\rm{RF}}}^2(t) $ 相互作用的脉冲响应函数;符号$ \otimes $ 代表卷积;$ {c_0} $ 表示射频输入信号的损耗或增益。因此,对于射频端口,非线性分量定义为:$$ {S_{{\rm{out}}\_{\rm{ME}}}}(t) = [{c_0} + {h_2}(t) \otimes V_{{\rm{RF}}}^2(t)] {V_{{\rm{RF}}}}(t) $$ (10) 利用式(9)中的调制函数
$ H(t) $ 与式(3)中的双音信号相互作用可得:$$ \begin{array}{*{20}{c}} {{S_{{\rm{out}}\_{\rm{ME}}}}(t)} =\\ { \Big[ {H(0)V_0^2 + {c_0} + \displaystyle\sum\limits_{\scriptstyle - 1 \leqslant p \leqslant 2, - 1 \leqslant q \leqslant 2\atop \scriptstyle\;\;\;\;\;p \in Z,q \in Z}^{1 \leqslant \left| p \right| + \left| q \right| \leqslant 2;1 \leqslant p + q \leqslant 2} {{h_{p{\omega _1} + q{\omega _2}}}(t)} } \Big] {V_{{\rm{RF}}}}(t)} \end{array}$$ (11) 式中,
$ {h_{p{\omega _1} + q{\omega _2}}}(t) $ 为频率$ p{\omega _1} + q{\omega _2} $ 处的非线性调制函数。同时这些非线性调制函数也表示了双音信号的二阶交调分量,如下所示:$$ \begin{split} & \;\;\;\; {h_{2{\omega _m}}}(t) = h(t) \otimes V_0^2\cos [2({\omega _m}t + {\varphi _m})] = \\ & V_0^2\left| {H(2{\omega _m})} \right|\cos \left[ {2({\omega _m}t + {\varphi _m}) + \angle H(2{\omega _m})} \right] \end{split} $$ (12) $$ \begin{split} & \;\;\;\;{h_{2{\omega _c}}}(t) = h(t) \otimes V_0^2\cos [2({\omega _c}t + {\varphi _c})]= \\ & V_0^2\left| {H(2{\omega _c})} \right|\cos \left[ {2({\omega _c}t + {\varphi _c}) + \angle H(2{\omega _c})} \right] \end{split} $$ (13) $$ \begin{split} & \;\;\;\; {h_{2{\omega _2}}}(t) = h(t) \otimes \frac{{V_0^2}}{2}\cos [2({\omega _2}t + {\varphi _2})]= \\ & \frac{{V_0^2}}{2}\left| {H(2{\omega _2})} \right|\cos \left[ {2({\omega _2}t + {\varphi _2}) + \angle H(2{\omega _2})} \right] \end{split} $$ (14) $$ \begin{split} & \;\;\;\; h{(t)_{2{\omega _1}}} = h(t) \otimes \frac{{V_0^2}}{2}\cos [2({\omega _1}t + {\varphi _1})] = \\ & \frac{{V_0^2}}{2}\left| {H(2{\omega _1})} \right|\cos \left[ {2({\omega _1}t + {\varphi _1}) + \angle H(2{\omega _1})} \right] \end{split} $$ (15) $ h(t) $ 的傅里叶变换被定义为:$$ H(\omega ) = \left| {H(\omega )} \right|\exp [j(\angle H(\omega ))] $$ (16) 本振信号被定义为:
$$ {V_{{\rm{LO}}}}(t) = {V_{{\rm{LO}}}}\cos ({\omega _{{\rm{LO}}}}t + {\varphi _{{\rm{LO}}}}) $$ (17) 式中,
$ {V_{{\rm{LO}}}} $ 、$ {\omega _{{\rm{LO}}}} $ 和$ {\varphi _{{\rm{LO}}}} $ 分别代表本振信号的幅度、频率和初相。本振信号的展开多项式为:$$ {S_{{\rm{LO}}}}(t) = \sum\limits_{m = 1}^l {{l_m}{{\left[ {{V_{{\rm{LO}}}}(t)} \right]}^m}} $$ (18) 根据图2,理想混频器的输出为:
$$ \begin{split} & \;\;\;\;\;\;\;\;\;\;\;\; {S_{{\rm{out}}}}(t) = {S_{{\rm{out}}\_{\rm{ME}}}}(t) {S_{{\rm{LO}}}}(t)=\;\; \\ & \left\{ \begin{gathered} \left[ {{c_0}{\text{ + }}H(0)V_0^2} \right] {V_{{\rm{RF}}}}(t){V_{{\rm{LO}}}}(t){\text{ + }} \\ \sum\limits_{\begin{subarray}{l} - 1 \leqslant p \leqslant 2, - 1 \leqslant q \leqslant 2 \\ \;\;\;\;\;p \in Z,q \in Z \end{subarray}} ^{1 \leqslant \left| p \right| + \left| q \right| \leqslant 2;1 \leqslant p + q \leqslant 2} {S{{(t)}_{{\rm{out}}\_{\rm{p}}{\omega _1} + q{\omega _2}\_{\rm{ME}}}}} \\ \end{gathered} \right\} \end{split} $$ (19) -
根据式(19),在基带频率处,基带信号与时变调制函数相互作用的分量如式(20)所示,其中最后两项表现为对IM3分量的贡献:
$$ \begin{split} &\qquad\qquad\;\;\;\;\;\;\;\;\;\;\;\;\; {S_{{\rm{out}}\_2{\omega _m}\_{\rm{ME}}}}(t) = \\ & \sum\limits_{\begin{subarray}{l} m = 2k - 1 \\ \;\;k \in {Z^ + } \end{subarray}} ^l {\left[ \begin{gathered} \frac{{C_m^{\tfrac{{m + 1}}{2}}V_{{\rm{LO}}}^m{l_m}V_0^3\left| {H(2{\omega _m})} \right|}}{{{2^m}}} \times \\ \left\{ \begin{gathered} \cos [({\omega _2} - {\omega _{{\rm{LO}}}})t + ({\varphi _2} - {\varphi _{{\rm{LO}}}}) + \\ \;\;\;\; \angle H(2{\omega _m})]+ \\ \cos [({\omega _1} - {\omega _{{\rm{LO}}}})t + ({\varphi _1} - {\varphi _{{\rm{LO}}}}) - \\ \;\;\;\;\;\;\; \angle H(2{\omega _m})]+ \\ \cos [(2{\omega _2} - {\omega _1} - {\omega _{{\rm{LO}}}})t + (2{\varphi _2} - {\varphi _1} - \\ \;\;\;\;\;\;\; {\varphi _{{\rm{LO}}}}){\text{ + }}\angle H(2{\omega _m})] + \\ \cos [(2{\omega _1} - {\omega _2} - {\omega _{{\rm{LO}}}})t + (2{\varphi _1} - {\varphi _2} - \\ \;\;\;\;\;\;\; {\varphi _{{\rm{LO}}}}) - \angle H(2{\omega _m})] \\ \end{gathered} \right\} \\ \end{gathered} \right]} \end{split} $$ (20) 式(20)中,
$ H(2{\omega _m}) $ 为$ {h_{2{\omega _m}}}(t) $ 的傅里叶变换系数。需要注意的是,在最后两项中,基带调制函数的相位$ \angle H(2{\omega _m}) $ 被叠加到IM3分量的上边带的贡献项的相位中,相反,IM3分量的下边带的贡献项的相位中却减去了基带调制函数的相位$ \angle H(2{\omega _m}) $ 。因此,从该式的最后两项中可以得知这会导致IM3分量上边带相位与下边带相位具有相反的变化趋势。 -
式(21)和式(22)中,
$ H(2{\omega _1}) $ 和$ H(2{\omega _2}) $ 分别为$ {h_{2{\omega _1}}}(t) $ 和$ {h_{2{\omega _2}}}(t) $ 的傅里叶变化系数。式(21)和式(22)的最后一项显示了由二次谐波调制对IM3分量的上边带和下边带的贡献。在这两项中,应该注意的是,$ \angle H(2{\omega _1}) $ 被叠加到IM3分量的下边带的相位上,同时$ \angle H(2{\omega _2}) $ 也被叠加到IM3分量的上边带的相位上。根据式(21)和式(22),当双音间距足够小时,
$ 2{\omega _1} $ 与$ 2{\omega _2} $ 接近,即$ \angle H(2{\omega _1}) $ 近似等于$ \angle H(2{\omega _2}) $ 。因此,可以得知由二次谐波调制产生的IM3分量上下边带的相位具有相同的变化趋势。这也说明在窄带情况下,IM3分量上下边带的相位具有相反的变化趋势是通过式(20)所示的基带调制实现的。基带调制机制不同于窄带情况下的二次谐波调制。在窄带情况下,基带调制机制很容易与二次谐波调制区别。$$ \begin{split} &\qquad\qquad \;\;\;\;\;\;\;\;\;\;\;\;\; {S_{{\rm{out}}\_2{\omega _1}\_{\rm{ME}}}}(t)= \\ & \sum\limits_{\begin{subarray}{l} m = 2k - 1 \\ k \in {Z^ + } \end{subarray} }^l {\left[ \begin{gathered} \frac{{C_m^{\tfrac{{m + 1}}{2}}V_{{\rm{LO}}}^m{l_m}V_0^3}}{{{2^{m{\text{ + }}1}}}}\left| {H(2{\omega _1})} \right| \times \\ \left\{ \begin{gathered} \cos [({\omega _1} - {\omega _{{\rm{LO}}}})t + ({\varphi _1} - {\varphi _{{\rm{LO}}}})+ \\ \;\;\;\;\; \angle H(2{\omega _1})]+ \\ \cos [(3{\omega _1} - {\omega _{{\rm{LO}}}})t + (3{\varphi _1} - {\varphi _{{\rm{LO}}}}) + \\ \;\;\;\;\;\;\; \angle H(2{\omega _1})] + \\ \cos [(2{\omega _1}{\text{ + }}{\omega _2} - {\omega _{{\rm{LO}}}})t + (2{\varphi _1}{\text{ + }}{\varphi _2} - \\ \;\;\;\;\;\;\; {\varphi _{{\rm{LO}}}}) + \angle H(2{\omega _1})]+ \\ \cos [(2{\omega _1} - {\omega _2} - {\omega _{{\rm{LO}}}})t + (2{\varphi _1} - {\varphi _2} - \\ \;\;\;\;\;\;\; {\varphi _{{\rm{LO}}}}) + \angle H(2{\omega _1})] \\ \end{gathered} \right\} \\ \end{gathered} \right]} \end{split} $$ (21) $$ \begin{split} &\qquad\qquad \;\;\;\;\;\;\;\;\;\;\;\;\; {S_{{\rm{out}}\_2{\omega _2}\_{\rm{ME}}}}(t) = \\ & \sum\limits_{\begin{subarray}{l} m = 2k - 1 \\ k \in {Z^ + } \end{subarray}} ^l {\left[ \begin{gathered} \frac{{C_m^{\tfrac{{m + 1}}{2}}V_{{\rm{LO}}}^m{l_m}V_0^3}}{{{2^{m{\text{ + }}1}}}}\left| {H(2{\omega _2})} \right| \times \\ \left\{ \begin{gathered} \cos [({\omega _2} - {\omega _{{\rm{LO}}}})t + ({\varphi _2} - {\varphi _{{\rm{LO}}}}) + \\ \;\;\;\; \angle H(2{\omega _2})] + \\ \cos [(3{\omega _2} - {\omega _{{\rm{LO}}}})t + (3{\varphi _2} - {\varphi _{{\rm{LO}}}}) + \\ \;\;\;\;\;\;\; \angle H(2{\omega _2})] + \\ \cos [(2{\omega _2}{\text{ + }}{\omega _1} - {\omega _{{\rm{LO}}}})t + (2{\varphi _2}{\text{ + }}{\varphi _1} - \\ \;\;\;\;\;\;\; {\varphi _{{\rm{LO}}}}) + \angle H(2{\omega _2})] + \\ \cos [(2{\omega _2} - {\omega _1} - {\omega _{{\rm{LO}}}})t + (2{\varphi _2} - {\varphi _1} - \\ \;\;\;\;\;\;\; {\varphi _{{\rm{LO}}}}) + \angle H(2{\omega _2})] \\ \end{gathered} \right\} \\ \end{gathered} \right]} \end{split} $$ (22) 当双音间距较宽时,
$ 2{\omega _1} $ 和$ 2{\omega _2} $ 差异变大,相位$ \angle H(2{\omega _1}) $ 不再近似等于$ \angle H(2{\omega _2}) $ 。此时,IM3分量上下边带的相位的反向趋势的来源可能包括二次谐波调制,如当$ \angle H(2{\omega _1}) $ 和$ \angle H(2{\omega _1}) $ 随双音间距的变化具有相反的变化趋势,或这两项具有相反的符号时,IM3分量上下边带的相位也具有相反的变化趋势。
Study on the Mechanisms of Phase Asymmetry of Intermodulation Products in Mixers
-
摘要: 基于二阶记忆效应机制研究了混频器交调分量的相位不对称性机理。该研究是基于时变调制函数(TVMF)和多盒行为混频器模型(MBBM)提出的,这种相位不对称性表现在三阶交调分量的上下边带的相位相对于双音间距具有反向变化趋势。 研究表明,基带调制和二次谐波调制是导致三阶交调分量上下边带的相位不对称的原因。但是,基带调制产生的三阶交调分量的上下边带相位具有相反变化趋势的机制与二次谐波调制造成的三阶交调分量相位具有相反变化趋势的机制不同,可以以此加以区分。 最后,利用ADS对混频器的三阶交调分量上下边带的相位进行仿真,并结合实测结果验证了上述分析。Abstract: The phase asymmetry of the mixers’ intermodulation (IM) products is studied in this paper based on the second-order memory mechanisms. This study is developed based on the time-varying modulation function (TVMF) and multi-box behavioral mixer (MBBM) model. This phase asymmetry manifests in the phase where the lower and upper third intermodulation (IM3) products have inverse changing trends vs. tone spacing. The research shows that the inverse trends of the phase of the upper and lower IM3 products are caused by baseband modulation and second harmonics modulations. However, the upper and lower IM3 products’ phase inverse trend mechanism created by baseband modulation is different from that created by second harmonics modulations, and thus may be readily distinguished. Finally, the phases of the upper and lower IM3 products of the mixer are simulated based on advanced design system (ADS), and the above analysis is verified with the measured results.
-
[1] ZHAO J, LIU J, YANG L, et al. Future 5G-oriented system for urban rail transit: Opportunities and challenges[J]. China Communications, 2021, 18(2): 1-12. doi: 10.23919/JCC.2021.02.001 [2] YAN Y, MA M. Novel frequency-domain oversampling receiver for CP MC-CDMA systems[J]. IEEE Communications Letters, 2015, 19(4): 661-664. doi: 10.1109/LCOMM.2015.2391103 [3] ARRANO H F, AZURDIA-MEZA C A. OFDM: Today and in the future of next generation wireless communications[C]//IEEE Central America and Panama Student Conference. Guatemala: IEEE, 2016: 1-6. [4] VERSPECHT J, HORN J, BETTS L, et al. Extension of X-parameters to include long-term dynamic memory effects[C]//International Microwave Symposium Digest. Boston: IEEE, 2009: 741-744. [5] CARVALHO N B, PEDRO J C. A comprehensive explanation of distortion sideband asymmetries[J]. IEEE Transactions on Microwave Theory & Techniques, 2002, 50(9): 2090-2101. doi: 10.1109/TMTT.2002.802321 [6] BOSCH W, GATTI G. Measurement and simulation of memory effects in predistortion linearizers[J]. IEEE Transactions on Microwave Theory & Techniques, 1989, 37(12): 1885-1890. [7] APARIN V, PERSICO C. Effect of out-of-band terminations on intermodulation distortion in common-emitter circuits[C]//1999 IEEE MTT-S International Microwave Symposium Digest. Anaheim: IEEE, 1999: 977-980. [8] CARVALHO N B, PEDRO J C. Two-tone IMD asymmetry in microwave power amplifiers[C]//2000 IEEE MTT-S International Microwave Symposium Digest. Boston: IEEE, 2000: 445-448. [9] ASBECK, KOBAYASHI P M, IWAMOTO H, et al. Augmented behavioral characterization for modeling the nonlinear response of power amplifiers[C]//2002 IEEE MTT-S International Microwave Symposium Digest. Seattle: IEEE, 2002: 135-138. [10] DRAXLER P, LANGMORE I, HUNT P, et al. Time domain characterization of power amplifiers with memory effects[C]//2003 IEEE Microwave Symposium Digest, MTT-S International. Philadelphia: IEEE, 2003: 803-806. [11] SAKIAN P, MAHMOUDI R, ROERMUND A V. A circuit-level analysis of third order intermodulation mechanisms in CMOS mixers using time-invariant power and Volterra series[C]//2011 41st European Microwave Conference. Manchester: IEEE, 2011: 595-598. [12] OZGUN A, NESIMOGLU T, DEMIR S. A multi-box behavioural mixer model and its validation using measurements[C]//2019 IEEE International Microwave Symposium. Los Angeles: IEEE, 2019: 607-610. [13] STIGLITZ M R, CALLAHAN J. Microwave mixers[M]. 2nd ed. Norwood, MA: Artech House Press, 1993: 21-25. [14] WILLIAMS D F, NDAGIJIMANA F, REMLEY K A, et al. Scattering-parameter models and representations for microwave mixers[J]. IEEE Transactions on Microwave Theory & Techniques, 2005, 53(1): 314-321. [15] ROOT E D, VERSPECHT J, HORN J, et al. X-Parameters: Characterization, modeling, and design of nonlinear RF and microwave components[M]. Cambridge: Cambridge University Press, 2016. [16] VANDERSTEEN G, VORST A V. Nonlinear model and measurement method for microwave mixers[J]. IEEE Microwave & Guided Wave Letters, 1993, 3(5): 133-135. [17] ALIZADEH M, RNNOW D. A two-tone test for characterizing nonlinear dynamic effects of radio frequency amplifiers in different amplitude regions[J]. Measurement, 2016, 89: 273-279. doi: 10.1016/j.measurement.2016.04.027 [18] SHEN Y. Computation-based phase measurement of RF power-amplifier intermodulation products[J]. IEEE Transactions on Instrumentation & Measurement, 2011, 60(8): 2934-2941. [19] WALKER A, STEER M, GARD K G. A vector intermodulation analyzer applied to behavioral modeling of nonlinear amplifiers with memory[J]. IEEE Transactions on Microwave Theory & Techniques, 2006, 54(5): 1991-1999. [20] LEE S Y, LEE Y S, JEONG Y H. A novel phase measurement technique for IM3 components in RF power amplifiers[J]. IEEE Transactions on Microwave Theory & Techniques, 2006, 54(1): 451-457. [21] REMLEY K A, WILLIAMS D F, SCHREURS M P, et al. Simplifying and interpreting two-tone measurements[J]. IEEE Transactions on Microwave Theory & Techniques, 2004, 52(11): 2576-2584.