一种多突变结构任意截面复功率的计算方法技术

本发明专利技术属于微波无源器件CAD技术领域，具体涉及一种多突变结构任意截面复功率的计算方法，基于模式匹配法。本发明专利技术避开了电场与磁场的繁琐求解，减小了计算机资源消耗与计算时间；直接利用各级突变的散射矩阵获得任意截面的复功率，对模式匹配法应用的正确性提供了一种简单而有效的验证方式。可以快速准确地得到多级突变结构内部任意截面的复功率值，是一种快速高效、便于操作的计算方法。

A new method for calculating the complex power of arbitrary cross section

The invention belongs to the technical field of microwave passive device CAD, in particular relates to a method for calculating the complex power of any cross section of a multi mutation structure. The invention avoids the complicated solving electric field and magnetic field, reduce the consumption of computer resources and computing time; directly using the scattering matrix at all levels of mutation of complex power obtain arbitrary cross section, provides a simple and effective way to verify the correctness of the pattern matching method. It can quickly and accurately obtain the complex power value of any cross section of the multi-level catastrophe structure.

【技术实现步骤摘要】

本专利技术属于微波无源器件CAD
，具体涉及一种多突变结构任意截面复功率的计算方法(基于模式匹配法)。
技术介绍
微波无源器件在微波通信领域占有重要地位，被广泛地应用于现代通讯和雷达装置的微波毫米波系统中。典型的微波无源器件有高功率微波窗、波导滤波器、功率分配器、模式变换器、阻抗匹配器、定向耦合器、移相器、环行器、隔离器和衰减器等。这些无源器件一般都是由不同的规则几何体组成，包含一些渐变结构，阶梯结构等。微波无源器件的设计与传输特性的分析是一个重要的课题。包含多级不连续性结构的微波无源器件示意图如图1所示。整个微波无源器件由多个不同的均匀波导段构成，每个均匀波导段称为一个分区，共有n+1个分区，各分区从左到右依次编号1,2,…,n+1。不同均匀波导段的连接处构成了一个突变，共有n个突变，各突变从左到右依次编号1,2,…,n。模式匹配法是一种基于严格场理论的全波分析方法，具有仿真速度快、求解精度高、占用计算机内存资源小等优势，被广泛应用于微波无源器件的仿真设计与分析中。模式匹配法的主要思想是将各均匀波导段中的电磁场进行模式展开，在相邻波导段连接的不连续处，通过横向电场连续、横向磁场连续、或复功率守恒来实现两侧波导之间的模式耦合。这就涉及到如何在模式匹配法的框架下进行复功率的计算问题。目前，尚没有文献给出具体的基于模式匹配法中的复功率的计算方法。复功率Pc定义为：其中，和分别为归一化电场与磁场，一般情况下为复数；S0为任意截面；代表在横截面S0上的面积分；Pc的实部为有功功率，虚部为无功功率。利用HFSS、CST等三维电磁软件进行微波无源器件仿真和设计时，在求解出电场与磁场分布的基础上，可以利用软件的后处理功能计算坡印延矢量在横截面S0上的面积分，以得到单频点的复功率。
技术实现思路
本专利技术的目的是为了解决现有技术中对模式匹配法复功率的计算不明确的问题，提供了一种多突变结构任意截面复功率的计算方法。该计算方法基于模式匹配法，避开了电场与磁场的繁琐求解，直接利用各级突变的散射矩阵获得任意截面的复功率，对模式匹配法应用的正确性提供了一种简单而有效的验证方式。本专利技术技术方案为：一种多突变结构任意截面复功率的计算方法，包括以下步骤：S1、整个微波无源器件由多个不同的均匀波导段构成，每个均匀波导段称为一个分区，共有n+1个分区，各分区从左到右依次编号1,2,…,n+1，不同均匀波导段的连接处构成了一个突变，共有n个突变，各突变从左到右依次编号1,2,…,n。根据多突变结构所需求解复功率的截面S0的位置，将截面S0左侧的所有突变(共p个突变)划分为一个突变群，将截面S0右侧的所有突变(共n-p个突变)划分为一个突变群。共两个突变群，其划分如图2所示。对于截面S0刚好位于某突变上的情况，既可以把该突变划分到截面S0右侧的突变群中，也可以将其划分到截面S0左侧的突变群中。若截面S0左侧不包含突变结构，则左侧突变群可视为与截面S0左侧端口重合的一个截面；若截面S0右侧不包含突变结构，则右侧突变群可视为与截面S0右侧端口重合的一个截面。S2、运用模式匹配法给出每一个突变对应的广义散射矩阵，根据步骤S1划分的两个突变群，将截面S0左侧突变群包含的所有突变进行级联，得到级联矩阵SL；将截面S0右侧突变群包含的所有突变进行级联，得到级联矩阵SR。在模式匹配法中，各个均匀波导段内的横向电磁场采用归一化全波展开，即有：其中，ηi,r为区域i(1≤i≤n+1)中第r(1≤r≤Ri)个模式的波阻抗；Ri为区域i的模式展开数；与分别代表区域i中第r个模式的横向电场模矢量与横向磁场模矢量；分别代表区域i中沿z轴正向和负向的第r个模式的归一化功率幅度系数。在区域i中的任意截面Si上，各模式间满足归一化与正交性条件，即有：其中，a和b分别代表区域i中的模式a与模式b。故有，任意截面S0处的横向电场横向磁场其中，η0r为截面S0所在区域的模式r(1≤r≤N)的波阻抗；N为截面S0所在区域的模式展开数；与分别代表截面S0所在区域的模式r的横向电场模矢量与横向磁场模矢量；分别代表沿z轴正向和负向穿过截面S0的模式r的归一化功率幅度系数。且在截面S0上有：模式匹配法给出每一个突变对应的广义散射矩阵，突变1～n的广义散射矩阵分别为S1～Sn。第i(1≤i≤n)个突变的广义散射矩阵为阶的方阵，矩阵Si,11的维数为矩阵Si,12的维数为矩阵Si,21的维数为矩阵Si,22的维数为其中，为突变i左侧区域的模式展开数，为突变i右侧区域的模式展开数。由于突变i右侧的区域就是突变i+1左侧的区域，故该模式匹配过程为本领域常识，这里不再赘述。将截面S0左侧的p个突变的广义散射矩阵S1～Sp利用矩阵级联理论进行级联，可以得到截面S0左侧的级联矩阵同理，将截面S0右侧的n-p个突变的广义散射矩阵Sp+1～Sn利用矩阵级联理论进行级联，可以得到截面S0右侧的级联矩阵其中，矩阵的维数为M×M；矩阵的维数为M×N；矩阵的维数为N×M；矩阵的维数为N×N；矩阵的维数为N×N；矩阵的维数为N×K；矩阵的维数为K×N；矩阵的维数为K×K。其中，代表左侧突变群的左侧区域选取的模式数，代表两个等效突变群之间的过渡区域(即任意截面S0所在区域)选取的模式数，代表右侧等效突变群的右侧区域选取的模式数。该矩阵级联过程为本领域常识，这里不再赘述。对于截面S0左侧不包含突变结构的情形，截面S0左侧的级联矩阵可以视为OL、IL分别代表阶的全零矩阵和单位矩阵；对于截面S0右侧不包含突变结构的情形，截面S0右侧的级联矩阵可以视为OR、IR分别代表阶的全零矩阵和单位矩阵。S3、将步骤S1中划分的两个突变群分别等效为两个等效突变结，并将步骤S2中得到的级联矩阵SL作为截面S0左侧等效突变结的广义散射矩阵；将步骤S2中得到的级联矩阵SR作为截面S0右侧等效突变结的广义散射矩阵，构造出多突变结构的等效模型。等效模型仅包含两个等效突变结，分为三个区，从左到右依次编号一区、二区、三区。实际上，一区即尚未等效的原多级模型中的1区，二区即尚未等效的原多级模型中的p+1区，三区即为尚未等效的原多级模型中的n+1区。任意多突变模型的等效模型分析图如图3所示。其中，A1、A2、A3分别代表一区、二区、三区各模式的归一化功率幅度系数向量，元素个数与相应区域所选择的模式数量保持一致，分别为M、N、K；上标+、-分别代表沿z轴正向和负向的模式，上标L、R分别代表某传输区域的最左侧和最右侧；分别代表沿z轴正向和负向穿过截面S0的各模式的归一化功率幅度系数向量；D1为对角线元素组成的M阶对角矩阵，γ1m为一区中第m(1≤m≤M)个模式的传播常数，h1为一区长度；D2为对角线元素组成的N阶对角矩阵，γ2q为二区中第q(1≤q≤N)个模式的传播常数，h2为二区长度；D3为对角线元素组成的K阶对角矩阵γ3k为一区中第k(1≤k≤K)个模式的传播常数，h3为三区长度；左侧等效突变结与截面S0的距离为l1，相应的传输矩阵为DL，是由对角元素组成的N阶对角矩阵；右侧等效突变结与截面S0的距离为l2，相应的传输矩阵为DR，是由对角元素组成的N阶对角矩阵。对于截面S0左侧不包含突变结构的情形，满足h1＝0；对于截面S0右侧不包含突变结构本文档来自技高网...

【技术保护点】
一种多突变结构任意截面复功率的计算方法，包括以下步骤：S1、整个微波无源器件由多个不同的均匀波导段构成，每个均匀波导段称为一个分区，共有n+1个分区，各分区从左到右依次编号1,2,…,n+1，不同均匀波导段的连接处构成了一个突变，共有n个突变，各突变从左到右依次编号1,2,…,n；根据多突变结构所需求解复功率的截面S0的位置，将截面S0左侧的所有共p个突变划分为一个突变群，将截面S0右侧的所有共n‑p个突变划分为一个突变群，共计两个突变群；对于截面S0刚好位于某突变上的情况，既可以把该突变划分到截面S0右侧的突变群中，也可以将其划分到截面S0左侧的突变群中；若截面S0左侧不包含突变结构，则左侧突变群可视为与截面S0左侧端口重合的一个截面；若截面S0右侧不包含突变结构，则右侧突变群可视为与截面S0右侧端口重合的一个截面；S2、运用模式匹配法给出每一个突变对应的广义散射矩阵，根据步骤S1划分的两个突变群，将截面S0左侧突变群包含的所有突变进行级联，得到级联矩阵SL；将截面S0右侧突变群包含的所有突变进行级联，得到级联矩阵SR。在模式匹配法中，各个均匀波导段内的横向电磁场采用归一化全波展开，即有：其中，ηi,r为区域i(1≤i≤n+1)中第r(1≤r≤Ri)个模式的波阻抗；Ri为区域i的模式展开数；与分别代表区域i中第r个模式的横向电场模矢量与横向磁场模矢量；分别代表区域i中沿z轴正向和负向的第r个模式的归一化功率幅度系数；在区域i中的任意截面Si上，各模式间满足归一化与正交性条件，即有：其中，a和b分别代表区域i中的模式a与模式b；故有，任意截面S0处的横向电场横向磁场其中，η0r为截面S0所在区域的模式r(1≤r≤N)的波阻抗，N为截面S0所在区域的模式展开数，与分别代表截面S0所在区域的模式r的横向电场模矢量与横向磁场模矢量，分别代表沿z轴正向和负向穿过截面S0的模式r的归一化功率幅度系数，且在截面S0上有：模式匹配法给出每一个突变对应的广义散射矩阵，突变1～n的广义散射矩阵分别为S1～Sn；第i(1≤i≤n)个突变的广义散射矩阵为阶的方阵，矩阵Si,11的维数为矩阵Si,12的维数为矩阵Si,21的维数为矩阵Si,22的维数为其中，为突变i左侧区域的模式展开数，为突变i右侧区域的模式展开数；由于突变i右侧的区域就是突变i+1左侧的区域，故将截面S0左侧的p个突变的广义散射矩阵S1～Sp利用矩阵级联理论进行级联，可以得到截面S0左侧的级联矩阵同理，将截面S0右侧的n‑p个突变的广义散射矩阵Sp+1～Sn利用矩阵级联理论进行级联，可以得到截面S0右侧的级联矩阵其中，矩阵的维数为M×M，矩阵的维数为M×N，矩阵的维数为N×M，矩阵的维数为N×N，矩阵的维数为N×N，矩阵的维数为N×K，矩阵的维数为K×N，矩阵的维数为K×K；其中，代表左侧突变群的左侧区域选取的模式数，代表两个等效突变群之间的过渡区域，即任意截面S0所在区域选取的模式数，代表右侧等效突变群的右侧区域选取的模式数；对于截面S0左侧不包含突变结构的情形，截面S0左侧的级联矩阵可以视为OL、IL分别代表阶的全零矩阵和单位矩阵；对于截面S0右侧不包含突变结构的情形，截面S0右侧的级联矩阵可以视为OR、IR分别代表阶的全零矩阵和单位矩阵；S3、将步骤S1中划分的两个突变群分别等效为两个等效突变结，并将步骤S2中得到的级联矩阵SL作为截面S0左侧等效突变结的广义散射矩阵；将步骤S2中得到的级联矩阵SR作为截面S0右侧等效突变结的广义散射矩阵，构造出多突变结构的等效模型；等效模型仅包含两个等效突变结，分为三个区，从左到右依次编号一区、二区和三区；其中，A1、A2、A3分别代表一区、二区、三区各模式的归一化功率幅度系数向量，元素个数与相应区域所选择的模式数量保持一致，分别为M、N、K；上标+、‑分别代表沿z轴正向和负向的模式，上标L、R分别代表某传输区域的最左侧和最右侧；分别代表沿z轴正向和负向穿过截面S0的各模式的归一化功率幅度系数向量；D1为对角线元素组成的M阶对角矩阵，γ1m为一区中第m(1≤m≤M)个模式的传播常数，h1为一区长度；D2为对角线元素组成的N阶对角矩阵，γ2q为二区中第q(1≤q≤N)个模式的传播常数，h2为二区长度；D3为对角线元素组成的K阶对角矩阵γ3k为一区中第k(1≤k≤K)个模式的传播常数，h3为三区长度；左侧等效突变结与截面S0的距离为l1，相应的传输矩阵为DL，是由对角元素组成的N阶对角矩阵；右侧等效突变结与截面S0的距离为l2，相应的传输矩阵为DR，是由对角元素组成的N阶对角矩阵；对于截面S0左侧不包含突变结构的情形，满足h1＝0；对于截面S0右侧不包含突变结构的情形，满足h3＝0；S4、根据模式匹配法和电磁传输理论，得到截面S0上的归一化功率...

【技术特征摘要】
1.一种多突变结构任意截面复功率的计算方法，包括以下步骤：S1、整个微波无源器件由多个不同的均匀波导段构成，每个均匀波导段称为一个分区，共有n+1个分区，各分区从左到右依次编号1,2,…,n+1，不同均匀波导段的连接处构成了一个突变，共有n个突变，各突变从左到右依次编号1,2,…,n；根据多突变结构所需求解复功率的截面S0的位置，将截面S0左侧的所有共p个突变划分为一个突变群，将截面S0右侧的所有共n-p个突变划分为一个突变群，共计两个突变群；对于截面S0刚好位于某突变上的情况，既可以把该突变划分到截面S0右侧的突变群中，也可以将其划分到截面S0左侧的突变群中；若截面S0左侧不包含突变结构，则左侧突变群可视为与截面S0左侧端口重合的一个截面；若截面S0右侧不包含突变结构，则右侧突变群可视为与截面S0右侧端口重合的一个截面；S2、运用模式匹配法给出每一个突变对应的广义散射矩阵，根据步骤S1划分的两个突变群，将截面S0左侧突变群包含的所有突变进行级联，得到级联矩阵SL；将截面S0右侧突变群包含的所有突变进行级联，得到级联矩阵SR。在模式匹配法中，各个均匀波导段内的横向电磁场采用归一化全波展开，即有：其中，ηi,r为区域i(1≤i≤n+1)中第r(1≤r≤Ri)个模式的波阻抗；Ri为区域i的模式展开数；与分别代表区域i中第r个模式的横向电场模矢量与横向磁场模矢量；分别代表区域i中沿z轴正向和负向的第r个模式的归一化功率幅度系数；在区域i中的任意截面Si上，各模式间满足归一化与正交性条件，即有：其中，a和b分别代表区域i中的模式a与模式b；故有，任意截面S0处的横向电场横向磁场其中，η0r为截面S0所在区域的模式r(1≤r≤N)的波阻抗，N为截面S0所在区域的模式展开数，与分别代表截面S0所在区域的模式r的横向电场模矢量与横向磁场模矢量，分别代表沿z轴正向和负向穿过截面S0的模式r的归一化功率幅度系数，且在截面S0上有：模式匹配法给出每一个突变对应的广义散射矩阵，突变1～n的广义散射矩阵分别为S1～Sn；第i(1≤i≤n)个突变的广义散射矩阵为阶的方阵，矩阵Si,11的维数为矩阵Si,12的维数为矩阵Si,21的维数为矩阵Si,22的维数为其中，为突变i左侧区域的模式展开数，为突变i右侧区域的模式展开数；由于突变i右侧的区域就是突变i+1左侧的区域，故将截面S0左侧的p个突变的广义散射矩阵S1～Sp利用矩阵级联理论进行级联，可以得到截面S0左侧的级联矩阵同理，将截面S0右侧的n-p个突变的广义散射矩阵Sp+1～Sn利用矩阵级联理论进行级联，可以得到截面S0右侧的级联矩阵其中，矩阵的维数为M×M，矩阵的维数为M×N，矩阵的维数为N×M，矩阵的维数为N×N，矩阵的维数为N×N，矩阵的维数为N×K，矩阵的维数为K×N，矩阵的维数为K×K；其中，代表左侧突变群的左侧区域选取的模式数，代表两个等效突变群之间的过渡区域，即任意截面S0所在区域选取的模式数，代表右侧等效突变群的右侧区域选取的模式数；对于截面S0左侧不包含突变结构的情形，截面S0左侧的级联矩阵可以视为OL、IL分别代表阶的全零矩阵和单位矩阵；对于截面S0右侧不包含突变结构的情形，截面S0右侧的级联矩阵可以视为OR、IR分别代表阶的全零矩阵和单位矩阵；S3、将步骤S1中划分的两个突变群分别等效为两个等效突变结，并将步骤S2中得到的级联矩阵SL作为截面S0左侧等效突变结的广义散射矩阵；将步骤S2中得到的级联矩阵SR作为截面S0右侧等效突变结的广义散射矩阵，构造出多突变结构的等效模型；等效模型仅包含两个等效突变结，分为三个区，从左到右依次编号一区、二区和三区；其中，A1、A2、A3分别代表一区、二区、三区各模式的归一化功率幅度系数向量，元素个数与相应区域所选择的模式数量保持一致，分别为M、N、K；上标+、-分别代表沿z轴正向和负向的模式，上标L、R分别代表某传输区域的最左侧和最右侧；分别代表沿z轴正向和负向穿过截面S0的各模式的归一化功率幅度系数向量；D1为对角线元素组成的M阶对角矩阵，γ1m为一区中第m(1≤m≤M)个模式的传播常数，h1为一区长度；D2为对角线元素组成的N阶对角矩阵，γ2q为二区中第q(1≤q≤N)个模式的传播常数，h2为二区长度；D3为对角线元素组成的K阶对角矩阵γ3k为一区中第k(1≤k≤K)个模式的传播常数，h3为三区长度；左侧等效突变结与截面S0的距离为l1，相应的传输矩阵为DL，是由对角元素组成的N阶对角矩阵；右侧等效突变结与截面S0的距离为l2，相应的传输矩阵为DR，是由对角元素组成的N阶对角矩阵；对于截面S0左侧不包含突变结构的情形，满足h1＝0；对于截面S0右侧不包含突变结构的情形，满足h3＝0；S4、根据模式匹配法和电磁传输理论，得到截面S0上的归一化功率幅度系数与端口输入条件之间的关系；根据等效模型分析图，结合模式匹配法与电磁传输理论，有：A1R+=D1·A1L+A3L-=D3·A3R----(6)]]>A2R+=D2·A2L+A2L-=D2·A2R----(7)]]>A0+=DL·A2L+A0-=DR·A2R----(8)]]>左侧等效突变结的广义散射矩阵SL满足：A1R-A2L+=S11LS12LS21LS22LA1R+A2L----(9)]]>右侧等效突变结的广义散射矩阵SR满足：A2R-A3L+=S11RS12RS21RS22RA2R+A3L----(10)]]>由9式和10式：A2L+=S21L·A1R++S22L·A2L----(11)]]>A2R-=S11R·A2R++S12R·A3L----(12)]]>将7式分别代入11式和12式，可得：D2-1·A2R+=S21L·A1R++S22L·D2·A2R----(13)]]>D2-1·A2L-=S11R·D2·A2L++S12R·A3L----(14)]]>将12式代入13式，可得：D2-1·A2R+=S21L·A1R++S22L·D2·(S11R·A2R++S12R·A3L-)---(15)]]>将11式代入14式，可得：D2-1·A2L-=S11R·D2·(S21L·A1R++S22L·A2L-)+S12R·A3L----(16)]]>整理15式与16式：A2R+=(I-D2·S22L·D2·S11R)-1·D2·(S21L·A1R++S22L·D2·S12R·A3L-)---(17)]]>A2L-=(I-D2·S11R·D2·S22L)-1·D2·(S11R·...

【专利技术属性】
技术研发人员：朱小芳，郝毅亮，胡权，胡玉禄，李斌，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川;51