CN116663476A - A High-Order High-Temperature Superconducting Filter and Its Circuit Topology Rapid Design Method - Google Patents
A High-Order High-Temperature Superconducting Filter and Its Circuit Topology Rapid Design Method Download PDFInfo
- Publication number
- CN116663476A CN116663476A CN202310605221.0A CN202310605221A CN116663476A CN 116663476 A CN116663476 A CN 116663476A CN 202310605221 A CN202310605221 A CN 202310605221A CN 116663476 A CN116663476 A CN 116663476A
- Authority
- CN
- China
- Prior art keywords
- circuit
- coupling
- circuit topology
- filter
- order
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000013461 design Methods 0.000 title claims abstract description 38
- 238000010168 coupling process Methods 0.000 claims abstract description 90
- 238000005859 coupling reaction Methods 0.000 claims abstract description 90
- 230000008878 coupling Effects 0.000 claims abstract description 89
- 239000011159 matrix material Substances 0.000 claims abstract description 66
- 238000005457 optimization Methods 0.000 claims abstract description 46
- 238000004088 simulation Methods 0.000 claims abstract description 36
- 238000013528 artificial neural network Methods 0.000 claims abstract description 27
- 230000004044 response Effects 0.000 claims description 18
- 238000003062 neural network model Methods 0.000 claims description 7
- 239000000758 substrate Substances 0.000 claims description 3
- 230000001808 coupling effect Effects 0.000 claims 1
- 238000000151 deposition Methods 0.000 claims 1
- 239000010409 thin film Substances 0.000 claims 1
- 238000012549 training Methods 0.000 description 17
- 230000008569 process Effects 0.000 description 7
- 239000013598 vector Substances 0.000 description 7
- 239000000284 extract Substances 0.000 description 6
- 230000015572 biosynthetic process Effects 0.000 description 5
- 238000012892 rational function Methods 0.000 description 5
- 238000003786 synthesis reaction Methods 0.000 description 5
- 238000003780 insertion Methods 0.000 description 4
- 230000037431 insertion Effects 0.000 description 4
- 230000009466 transformation Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 229910002244 LaAlO3 Inorganic materials 0.000 description 2
- 238000006880 cross-coupling reaction Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 210000002569 neuron Anatomy 0.000 description 2
- 229910021521 yttrium barium copper oxide Inorganic materials 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 210000002364 input neuron Anatomy 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 210000004205 output neuron Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/39—Circuit design at the physical level
- G06F30/392—Floor-planning or layout, e.g. partitioning or placement
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/60—Superconducting electric elements or equipment; Power systems integrating superconducting elements or equipment
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Architecture (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Superconductor Devices And Manufacturing Methods Thereof (AREA)
Abstract
Description
技术领域technical field
本发明属于射频微波电路设计领域以及微波通信领域,特别涉及一种L波段、带宽在3%以内的高阶高温超导滤波器电路拓扑结构的高效、快速设计方法。The invention belongs to the field of radio frequency microwave circuit design and the field of microwave communication, in particular to an efficient and rapid design method for a circuit topology of a high-order high-temperature superconducting filter with an L-band and a bandwidth within 3%.
背景技术Background technique
滤波器是电子系统的关键组成部分,微波滤波器选择所需的有用微波信号并且抑制带外的干扰。随着微波系统的持续发展,滤波器的性能要求越来越高。但是在已经报道的工作中,传统的微带滤波器通常为2阶或者4阶,因为插入损耗随着滤波器阶数的增加而显著增加,导致高阶微带滤波器插入损耗太大而无法使用,这是因为,基于常规材料设计的微带谐振器Q值较低,当特别是滤波器的相对带宽小于5%时,插入损耗显著增加。而通过超导技术实现的微带谐振器Q值可以到达1万以上,滤波器可以实现高阶,例如8阶,10阶或者更高,同时具备很低的插入损耗以及更高的选择性。Filters are a key component of electronic systems. Microwave filters select desired useful microwave signals and suppress out-of-band interference. With the continuous development of microwave systems, the performance requirements of filters are getting higher and higher. However, in the reported work, the traditional microstrip filter is usually 2nd or 4th order, because the insertion loss increases significantly with the increase of the filter order, resulting in the high-order microstrip filter insertion loss is too large to be able to This is because the Q value of microstrip resonators designed based on conventional materials is low, and when the relative bandwidth of the filter is less than 5%, the insertion loss increases significantly. The Q value of the microstrip resonator realized by superconducting technology can reach more than 10,000, and the filter can achieve high order, such as 8th order, 10th order or higher, and has very low insertion loss and higher selectivity.
高阶超导滤波器的设计,需要通过电磁仿真进行长时间的电路拓扑结构参数优化,以达到理想的性能。然而,滤波器拓扑电路结构优化高度依赖于设计人员经验。特别是随着滤波器阶数的增加,电路拓扑结构越来越复杂,传统电磁仿真优化所花费的精力和时间越来越长,甚至多达数月,这在实际工程应用中难以接受。目前迫切需要一种高效、准确的高阶高温超导滤波器拓扑结构快速设计的方法。The design of high-order superconducting filters requires long-term optimization of circuit topology parameters through electromagnetic simulation to achieve ideal performance. However, the optimization of the filter topology is highly dependent on the designer's experience. Especially as the filter order increases, the circuit topology becomes more and more complex, and the traditional electromagnetic simulation optimization takes longer and longer energy and time, even up to several months, which is unacceptable in practical engineering applications. There is an urgent need for an efficient and accurate method for rapid design of high-order HTS filter topology.
减少滤波器拓扑结构优化时间的有效方法是寻找一种更加高效的模型来替代耗时的电磁仿真过程。本发明将sonnet电磁仿真软件与人工神经网络方法相结合,提出一种基于两次优化的,高效的高阶超导滤波器电路拓扑结构设计方法,解决高阶复杂超导滤波器电路设计效率低的问题。An effective way to reduce filter topology optimization time is to find a more efficient model to replace the time-consuming electromagnetic simulation process. The present invention combines the Sonnet electromagnetic simulation software with the artificial neural network method, proposes a high-efficiency high-order superconducting filter circuit topology design method based on two optimizations, and solves the problem of low design efficiency of high-order complex superconducting filter circuits .
发明内容Contents of the invention
本发明的目的是提出一种高阶高温超导滤波器及其电路拓扑结构快速设计方法,解决高阶高温超导滤波器(阶数为8阶或者10阶以上)电路拓扑结构设计过程中,依靠设计The purpose of the present invention is to propose a high-order high-temperature superconducting filter and its circuit topology rapid design method to solve the problem of high-order high-temperature superconducting filter (the order is 8 or 10 or more) circuit topology design process. rely on design
师的经验进行电磁仿真设计速度慢、效率低的问题。The problem of slow speed and low efficiency of electromagnetic simulation design based on the experience of the teacher.
高阶高温超导滤波器由50欧姆特性阻抗的输入馈线和输出馈线,以及多个半波长谐振器组成;谐振器的形式可以是单螺旋、双螺旋。多个谐振器的组合构成了高温超导滤波器的基础电路拓扑结构。The high-order high-temperature superconducting filter is composed of an input feeder and an output feeder with a characteristic impedance of 50 ohms, and multiple half-wavelength resonators; the form of the resonator can be single helix or double helix. The combination of multiple resonators forms the basic circuit topology of HTS filters.
一种高阶高温超导滤波器的电路拓扑结构快速设计方法基于两次优化过程,将sonnet仿真软件与人工神经网络相结合,用更加高效的模型替代一部分电磁仿真,极大的提升高阶高温超导滤波器的设计效率。A rapid circuit topology design method for high-order high-temperature superconducting filters is based on two optimization processes, combining sonnet simulation software with artificial neural networks, replacing part of the electromagnetic simulation with a more efficient model, and greatly improving the high-order high-temperature Design efficiency of superconducting filters.
为实现本发明的目的,采用以下技术方案:高阶高温超导滤波器的电路拓扑结构快速设计方法将sonnet电磁仿真与神经网络结合起来,进行建模,通过两次优化过程,对滤波器的电路拓扑结构进行调整,从而快速、高效地寻找到最优的电路拓扑结构参数。In order to realize the purpose of the present invention, adopt following technical scheme: the circuit topology fast design method of high-order high-temperature superconducting filter combines sonnet electromagnetic simulation and neural network, carries out modeling, by two optimization processes, to the filter The circuit topology is adjusted to quickly and efficiently find the optimal circuit topology parameters.
所述方法的流程如图1所示。首先根据滤波器综合设计理论得到理论计算的理想耦合矩阵,依据理想耦合矩阵构建初始的滤波器电路拓扑结构,并对拓扑结构定义几何变量。之后对该初始拓扑结构调用sonnet进行电磁仿真,依据sonnet电磁仿真结果提取实际电路的电路耦合矩阵。比较理想耦合矩阵与实际电路耦合矩阵之间的差距进行第一次优化,通过调整电路拓扑结构几何变量,使实际路耦合矩阵与理想耦合矩阵之间的差异最小化。之后构建神经网络模型,训练的样本来自于用sonnet进行电路拓扑结构电磁仿真得到的结果响应,并对电路拓扑结构参数的取值范围进行预定义。之后对训练充分的神经网络模型进行第二次优化,直到满足预设的条件。The process flow of the method is shown in FIG. 1 . Firstly, the theoretically calculated ideal coupling matrix is obtained according to the filter synthesis design theory, and the initial filter circuit topology is constructed according to the ideal coupling matrix, and geometric variables are defined for the topology. Then call Sonnet for electromagnetic simulation on the initial topology structure, and extract the circuit coupling matrix of the actual circuit according to the results of Sonnet electromagnetic simulation. The first optimization is performed by comparing the gap between the ideal coupling matrix and the actual circuit coupling matrix. By adjusting the geometric variables of the circuit topology, the difference between the actual coupling matrix and the ideal coupling matrix is minimized. Afterwards, the neural network model is constructed, and the training samples come from the result response obtained from the electromagnetic simulation of the circuit topology using Sonnet, and the value range of the circuit topology parameters is predefined. Afterwards, the fully trained neural network model is optimized for the second time until the preset conditions are met.
所述优化过程分为两个步骤:第一次优化(初步优化)以及第二次优化。The optimization process is divided into two steps: the first optimization (preliminary optimization) and the second optimization.
所述第一次优化中,利用矢量拟合方法和矩阵相似变换,从sonnet电磁仿真响应中提取出真实滤波器电路拓扑结构的电路耦合矩阵。将提取的真实耦合矩阵与理论计算得到的理想耦合矩阵进行比较,进行电路拓扑结构参数的第一次优化,电路耦合系数是关键优化参数,电路耦合系数与设计电路结构一一对应,In the first optimization, the circuit coupling matrix of the real filter circuit topology is extracted from the sonnet electromagnetic simulation response by using the vector fitting method and matrix similarity transformation. Compare the extracted real coupling matrix with the theoretically calculated ideal coupling matrix to optimize the circuit topology parameters for the first time. The circuit coupling coefficient is the key optimization parameter, and the circuit coupling coefficient corresponds to the designed circuit structure one by one.
电路耦合系数的调整可以通过相应谐振之间间距改变来实现的,对其中差距最大的矩阵元素优先调整,直到他们之间的差异最小,对提取的电路耦合矩阵完成第一次优化。The adjustment of the circuit coupling coefficient can be realized by changing the distance between the corresponding resonances, and the matrix elements with the largest gap are preferentially adjusted until the difference between them is the smallest, and the first optimization of the extracted circuit coupling matrix is completed.
所述第二次优化,利用sonnet电磁仿真结果产生的样本,建立并训练人工神经网络模型,以提供准确和快速的电磁响应预测,执行通带内S 11≤-20 dB的条件,对电路几何变量进行迭代,直到条件满足。The second optimization, using the samples generated by sonnet electromagnetic simulation results, establishes and trains the artificial neural network model to provide accurate and fast electromagnetic response prediction, and implements the condition of S 11≤-20 dB in the passband, for circuit geometry Variables are iterated until a condition is met.
所述神经网络如图2所示。它由一个三层神经网络,以及一个从耦合矩阵到滤波器频率响应的变换组成。选取滤波器电路拓扑结构的几何变量x的向量作为神经网络的输入,以耦合矩阵[M]为输出。该方法的依据是所有的的耦合矩阵元,包括非相邻耦合,都是由滤波器电路的拓扑结构决定的。The neural network is shown in Figure 2. It consists of a three-layer neural network, and a transformation from the coupling matrix to the filter frequency response. The vector of the geometric variable x of the filter circuit topology is selected as the input of the neural network, and the coupling matrix [ M ] is output. The basis of this approach is that all coupling matrix elements, including non-adjacent couplings, are determined by the topology of the filter circuit.
所述滤波器电路拓扑结构几何变量的定义如下,几何变量是一个l维向量:The definition of the geometry variable of the filter circuit topology is as follows, the geometry variable is an l- dimensional vector:
耦合元素m ij可以表示为l维向量空间中的拓扑结构几何变量的函数。利用人工神经网络学习几何变量x与耦合矩阵m ij之间的关系:Coupling element m ij can be expressed as a function of topological geometric variables in l- dimensional vector space. Use the artificial neural network to learn the relationship between the geometric variable x and the coupling matrix m ij :
在得到耦合矩阵[M]后可以通过理论公式计算出该电路拓扑结构的[S]参数,公式如下所示:After obtaining the coupling matrix [ M ], the [ S ] parameter of the circuit topology can be calculated through a theoretical formula, and the formula is as follows:
其中[M]是从滤波器拓扑结构中提取的实际电路耦合矩阵,[I]是(n+2)×(n+2)单位矩阵,s是复频率s=jw,[R]是一个全零矩阵,除了R 11=R n+2,n+2= 1。where [ M ] is the actual circuit coupling matrix extracted from the filter topology, [ I ] is the (n+2)×(n+2) identity matrix, s is the complex frequency s = jw , and [ R ] is a full Zero matrix, except that R 11 = R n+2, where n+2 = 1.
所述耦合矩阵的提取主要包含以下主要步骤:通sonnet电磁仿真得到[S]参数;对[S]参数进行相位去嵌入;将[S]参数转化为导纳参数[Y]参数;采用矢量拟合法提取[Y]参数有理函数的极点以及留数;求出N+2阶横向矩阵的矩阵元值;将横向矩阵通过相似变换变为实际滤波器的耦合结构,求出实际耦合矩阵;最后通过耦合矩阵求出[S]参数。The extraction of the coupling matrix mainly includes the following main steps: obtain the [ S ] parameter through sonnet electromagnetic simulation; carry out phase de-embedding to the [ S ] parameter; convert the [ S ] parameter into the admittance parameter [ Y ] parameter; Legally extract the poles and residues of the [ Y ] parameter rational function; find the matrix element values of the N+2 order transverse matrix; convert the transverse matrix into the coupling structure of the actual filter through similar transformation, and find the actual coupling matrix; finally pass Find the [ S ] parameter from the coupling matrix.
所述散射参数[S]转换为导纳参数[Y],采用如下公式:The scattering parameter [ S ] is converted into the admittance parameter [ Y ] using the following formula:
利用矢量拟合算法将上述[Y]参数曲线进行拟合,得到一个有理函数。有理函数的数学表示如下:Use the vector fitting algorithm to fit the above [ Y ] parameter curve to obtain a rational function. The mathematical representation of a rational function is as follows:
其中分子和分母的最高次数等于滤波器的阶数。滤波器的导纳参数y 11、y 12、y 21和y 22共享相同的根,应该在一次运行中进行拟合。以矩阵形式重写表达式上面的有理函数:where the highest degree of the numerator and denominator is equal to the order of the filter. The admittance parameters y 11 , y 12 , y 21 and y 22 of the filter share the same root and should be fitted in one run. Rewrite the rational function above the expression in matrix form:
由于一个典型的双端口微波滤波器可以等同于一个并行连接的横向网络。其导纳可以写成一个并联连接的横向阵列:Since a typical two-port microwave filter can be equivalent to a horizontal network connected in parallel. Its admittance can be written as a transverse array connected in parallel:
其中,M Sk和M Lk分别表示谐振器k与源/负载端口之间的耦合,M SL为源与负载之间的耦合,B k为谐振器k的自耦合。将实际滤波器的[S]参数转化而来的导纳Y (s),与横向网络导纳阵列Y N (s)的有理函数进行比较,可以得到N+2阶横向矩阵的每个矩阵元的值:Among them, M Sk and M Lk represent the coupling between resonator k and the source/load port, M SL is the coupling between source and load, and B k is the self-coupling of resonator k . Comparing the admittance Y (s) transformed from the [ S ] parameter of the actual filter with the rational function of the transverse network admittance array Y N (s) , each matrix element of the N+2 order transverse matrix can be obtained value of:
如上面的矩[M p ],它表示一个N+2耦合矩阵,谐振器之间不存在耦合,只与输入输出端口之间存在耦合。但是实际滤波器电路的耦合结构与它不同,相邻谐振器之间都会有耦合,为了将该横向耦合结构转化为实际的耦合结构,矩阵[M p ]进行相似变换旋转操作,该旋转操作的数学表示如下所示:As in the above moment [ M p ], it represents an N+2 coupling matrix, there is no coupling between resonators, and there is only coupling between the input and output ports. However, the coupling structure of the actual filter circuit is different from it, and there will be coupling between adjacent resonators. In order to transform the transverse coupling structure into the actual coupling structure, the matrix [ M p ] performs a similar transformation rotation operation, and the rotation operation The mathematical representation looks like this:
其中,R θ (i,j) 是旋转矩阵。通过选择合适的旋转轴[i,j]、旋转角度θ,可以消除横向矩阵[M p]中的非耦合矩阵元,将其转化为滤波器的实际耦合结构。where R θ (i, j) is the rotation matrix. By selecting the appropriate rotation axis [ i, j ] and rotation angle θ , the non-coupling matrix elements in the transverse matrix [ M p ] can be eliminated and transformed into the actual coupling structure of the filter.
所述选择滤波器拓扑结构几何变量组成的向量x作为模型的输入,调用sonnet仿真软件,得到电路拓扑结构的仿真响应曲线,用该仿真结果去训练该神经网络。在训练过程中,使用反向传播算法对神经网络中的连接权值和偏差参数进行更新,使训练误差最小如下:The vector x composed of geometric variables of the selected filter topology is used as the input of the model, and the sonnet simulation software is called to obtain the simulation response curve of the circuit topology, and the simulation result is used to train the neural network. During the training process, the connection weights and bias parameters in the neural network are updated using the backpropagation algorithm to minimize the training error as follows:
其中为神经网络输出的耦合矩阵元值,m ij 为sonnet仿真响应提取出来的训练耦合矩阵元值。E k是第k个样本的训练误差(平方误差),N是滤波器的阶数。E为所有训练样本的均方误差,n为训练样本的总数。在一个训练良好的神经网络中,输出值/>被期望等于训练值m ij 。in is the coupling matrix element value output by the neural network, and m ij is the training coupling matrix element value extracted from the sonnet simulation response. E k is the training error (squared error) of the kth sample, and N is the order of the filter. E is the mean square error of all training samples, and n is the total number of training samples. In a well-trained neural network, the output value /> is expected to be equal to the training value m ij .
高温超导滤波器一般采用500 nm厚的YBCO薄膜,沉积在0.5 mm厚的LaAlO3衬底上,相对介电常数为23.75。High-temperature superconducting filters generally use a 500 nm thick YBCO film deposited on a 0.5 mm thick LaAlO3 substrate with a relative permittivity of 23.75.
本发明的特点及有益效果:Features and beneficial effects of the present invention:
本发明提出了一种高阶高温超导滤波器及其电路拓扑结构快速设计方法,它采用两次优化的方法,结合sonnet电磁仿真软件以及人工神经网络的特点,极大提升了设计效率,并且适用于高阶复杂耦合的滤波器(包含非相耦合)快速。本发明适用于L波段以及更高的频率的高阶、以及复杂电路拓扑结构的超导滤波器设计。可以快速地得到滤波器电路拓扑结构的几何参数。经实例验证该方法得到的电路响应与电磁仿真软件的仿真结果一致。The present invention proposes a high-order high-temperature superconducting filter and a rapid design method for its circuit topology. It adopts a two-time optimization method, combined with the characteristics of sonnet electromagnetic simulation software and artificial neural network, which greatly improves the design efficiency, and Filters for high-order complex couplings (including non-phase couplings) are fast. The invention is applicable to the superconducting filter design of high-order and complex circuit topological structure in L band and higher frequency. The geometric parameters of the filter circuit topology can be quickly obtained. It is verified by an example that the circuit response obtained by this method is consistent with the simulation result of the electromagnetic simulation software.
附图说明Description of drawings
图1本发明的高阶高温超导滤波器电路拓扑结构快速设计方法流程图;Fig. 1 is a flow chart of the rapid design method of high-order high-temperature superconducting filter circuit topology of the present invention;
图2本发明的基于sonnet电磁仿真的神经网络训练模型;The neural network training model based on sonnet electromagnetic simulation of Fig. 2 of the present invention;
图3本发明的8阶高温超导滤波器理论综合得到的理想耦合矩阵;The ideal coupling matrix obtained by the theoretical synthesis of the 8th-order high-temperature superconducting filter of the present invention of Fig. 3;
图4本发明的8阶高温超导滤波器sonnet仿真接的提取出实际电路的电路耦合矩阵;The 8-order high-temperature superconducting filter sonnet simulation connection of Fig. 4 of the present invention extracts the circuit coupling matrix of the actual circuit;
图5本发明的8阶高温超导滤波器电路拓扑结构几何参数定义;Fig. 5 definition of geometric parameter of circuit topology structure of 8th-order high-temperature superconducting filter of the present invention;
图6本发明的8阶高温超导滤波器初始化电路仿真结果以及第一次优化结果;Fig. 6 is the 8th-order high-temperature superconducting filter initialization circuit simulation result and the first optimization result of the present invention;
图7本发明的8阶高温超导滤波器最终参数预测曲线与电磁仿真曲线比较;Fig. 7 compares the final parameter prediction curve of the 8th-order high-temperature superconducting filter of the present invention with the electromagnetic simulation curve;
图8本发明的高温超导滤波器示意图;Fig. 8 is a schematic diagram of a high temperature superconducting filter of the present invention;
图9本发明的8阶和高温超导滤波器实测结果;Fig. 9 8th-order and high-temperature superconducting filter actual measurement results of the present invention;
图10本发明的10阶高温超导滤波器电路拓扑结构几何参数定义;Fig. 10 Definition of geometrical parameters of the circuit topology of the 10th-order high-temperature superconducting filter of the present invention;
图11本发明的10阶高温超导滤波器理论综合得到的理想耦合矩阵;Fig. 11 is the ideal coupling matrix obtained by theoretical synthesis of the 10th-order high-temperature superconducting filter of the present invention;
图12本发明的10阶高温超导滤波器初始化电路仿真结果以及第一次优化结果;Fig. 12 is the 10th-order high-temperature superconducting filter initialization circuit simulation result and the first optimization result of the present invention;
图13本发明的10阶高温超导滤波器最终参数预测曲线与电磁仿真曲线比较;Fig. 13 compares the final parameter prediction curve of the 10th-order high-temperature superconducting filter of the present invention with the electromagnetic simulation curve;
图14本发明的10阶高温超导滤波器实测结果。Fig. 14 is the measured result of the 10-order high-temperature superconducting filter of the present invention.
具体实施方式Detailed ways
以下结合附图和具体实施例对本发明进行详细描述,但不作为对本发明的限定。The present invention will be described in detail below in conjunction with the accompanying drawings and specific embodiments, but not as a limitation of the present invention.
本发明提出了一种高阶高温超导滤波器及其电路拓扑结构快速设计方法。采用常规的电磁仿真,依靠设计师的经验设计电路的拓扑结构参数,得到满足要求的S参数响应曲线通常需要数周的时间。本发明基于两次优化过程,从sonnet电磁仿真与神经网络方法,可用于滤波器响应与几何变量的函数的快速预测。为了证明该方法的有效性,下面结合附图和具体实施例对本发明做进一步地描述。The invention proposes a high-order high-temperature superconducting filter and a rapid design method for its circuit topology. Using conventional electromagnetic simulation, relying on the designer's experience to design the topology parameters of the circuit, it usually takes several weeks to obtain the S-parameter response curve that meets the requirements. The present invention is based on two optimization processes, and can be used for fast prediction of filter response and function of geometric variables from sonnet electromagnetic simulation and neural network methods. In order to prove the effectiveness of the method, the present invention will be further described below in conjunction with the accompanying drawings and specific embodiments.
实施例1:8阶高温超导滤波器及其电路拓扑结构快速设计方法;Embodiment 1: 8-order high-temperature superconducting filter and its circuit topology rapid design method;
8阶高温超导滤波器由50欧姆特性阻抗的输入馈线和输出馈线,以及8个半波长谐振器组成。The 8th-order high-temperature superconducting filter consists of an input feeder and an output feeder with a characteristic impedance of 50 ohms, and 8 half-wavelength resonators.
该高阶滤波器中心设定频率和带宽分别为f 0= 508 MHz和BW = 8 MHz。The center frequency and bandwidth of this high-order filter are set to f 0 = 508 MHz and BW = 8 MHz, respectively.
步骤一,高阶超导滤波器电路拓扑结构初始值确定:根据指标要求通过滤波器综合理论得出理想耦合矩阵如图3所示。依据理想耦合矩阵,确定滤波器电路拓扑的初始拓扑结构,如图5所示,选取6个电路拓扑结构几何变量作为输入变量,即x= {q 1、L 1、g 1、g 2、g 3、g 4}。q 1是馈线的长度。L 1表示第一谐振器的参数。g 1、g 2、g 3、g 4是相邻谐振器的间距。Step 1, the initial value of the circuit topology of the high-order superconducting filter is determined: according to the requirements of the index, the ideal coupling matrix is obtained through the filter synthesis theory, as shown in Figure 3. According to the ideal coupling matrix, determine the initial topology of the filter circuit topology, as shown in Figure 5, select six geometric variables of the circuit topology as input variables, namely x = { q 1 , L 1 , g 1 , g 2 , g 3 , g 4 }. q 1 is the length of the feeder. L 1 represents the parameters of the first resonator. g 1 , g 2 , g 3 , g 4 are the distances between adjacent resonators.
步骤二,滤波器电路拓扑结构矩阵提取:根据前文的方法提取电路拓扑结构的实际电路耦合矩阵。Step 2, filter circuit topology matrix extraction: extract the actual circuit coupling matrix of the circuit topology according to the above method.
步骤三,电路拓扑结构几何参数初次优化:对步骤一中得到理想耦合矩阵与步骤二中提取的实际电路耦合矩阵进行比较,可以通过调整电路耦合系数优化电路耦合矩阵。Step 3, the initial optimization of the geometric parameters of the circuit topology: compare the ideal coupling matrix obtained in step 1 with the actual circuit coupling matrix extracted in step 2, and optimize the circuit coupling matrix by adjusting the circuit coupling coefficient.
电路耦合系数的调整通过相应谐振之间间距改变来实现的,对其中差距最大的矩阵元素优先调整,直到他们之间的差异最小,对提取的电路耦合矩阵完成第一次优化后如图4所示The adjustment of the circuit coupling coefficient is realized by changing the spacing between the corresponding resonances. The matrix elements with the largest gap are adjusted preferentially until the difference between them is the smallest. After the first optimization of the extracted circuit coupling matrix is completed, as shown in Figure 4 Show
步骤四,确定样本训练的范围,并通过sonnet仿真响应结果生成训练样本。本例中参数的初始值集为{3.32、5.04、0.98、1.58、1.46},参数的范围见表一。由sonnet软件生成81个样本,并在训练范围内以网格形式随机分散。Step 4, determine the range of sample training, and generate training samples through sonnet simulation response results. The initial value set of parameters in this example is {3.32, 5.04, 0.98, 1.58, 1.46}, and the range of parameters is shown in Table 1. 81 samples were generated by Sonnet software and randomly scattered in a grid form within the training range.
表1 8阶高温超导滤波器电路拓扑结构几何参数范围Table 1 Geometric parameter range of 8th-order HTS filter circuit topology
步骤五,提取所有样本耦合矩阵。Step five, extract all sample coupling matrices.
步骤六,根据步骤四和五得到的样本训练神经网络:采用隐层神经网络中包含11个神经元的三层神经网络学习几何变量与耦合矩阵之间的关系。在所有的样本中,有65个样本用于训练,16个样本用于测试。Step 6, training the neural network based on the samples obtained in steps 4 and 5: using a three-layer neural network including 11 neurons in the hidden layer neural network to learn the relationship between the geometric variables and the coupling matrix. Among all samples, 65 samples are used for training and 16 samples are used for testing.
步骤七:利用神经网络进行第二次优化:对充分训练的神经网络,使用指标条件max(S 11)≤-20 dB(504~512 MHz)进行第二次优化。几何参数被进一步调整。最优解为x={3.92、4.88、0.92、1.56、1.44、1.64}。Step 7: Use the neural network for the second optimization: For the fully trained neural network, use the index condition max (S 11) ≤ -20 dB (504~512 MHz) for the second optimization. Geometry parameters are further tuned. The optimal solution is x = {3.92, 4.88, 0.92, 1.56, 1.44, 1.64}.
步骤八:将最终得到的电路拓扑结构几何参数带入sonnet软件进行电磁仿真,与神经网络预测结果对比。Step 8: Bring the geometric parameters of the final circuit topology into Sonnet software for electromagnetic simulation, and compare with the prediction results of the neural network.
在本例中,如表2所示,训练神经网络模型的总时间为10 h,然而,一旦模型经过良好的训练,只需要不到0.1秒就可以预测神经网络模型的响应。然后将所求得的电路拓扑结构参数用sonnet进行仿真,如图7所示,仿真响应与模型响应吻合较好,满足设计规范。该高温超导滤波器采用500 nm厚的YBCO薄膜,沉积在0.5mm厚的LaAlO3衬底上,相对介电常数为23.75。该滤波器示意图如图8。该滤波器在65 K下由安捷伦E5072A网络分析仪测量,测量的响应如图9所示。In this example, as shown in Table 2, the total time to train the neural network model is 10 h, however, once the model is well trained, it takes less than 0.1 seconds to predict the response of the neural network model. Then the obtained circuit topology parameters are simulated by Sonnet, as shown in Figure 7, the simulation response is in good agreement with the model response and meets the design specifications. The high-temperature superconducting filter uses a 500 nm thick YBCO film deposited on a 0.5 mm thick LaAlO3 substrate with a relative permittivity of 23.75. A schematic diagram of the filter is shown in Figure 8. The filter was measured with an Agilent E5072A network analyzer at 65 K, and the measured response is shown in Figure 9.
表2 8节滤波器设计中本发明方法与常规电磁仿真所用时间比较Table 2 Comparison of the time used by the method of the present invention and conventional electromagnetic simulation in the filter design of section 8
实施例2: 10阶非相邻耦合高温超导滤波器拓扑结构快速优化模型;Embodiment 2: 10-order non-adjacent coupled high-temperature superconducting filter topology fast optimization model;
10阶高温超导滤波器由50欧姆特性阻抗的输入馈线和输出馈线,以及10个半波长谐振器组成。The 10th-order high-temperature superconducting filter consists of an input feeder and an output feeder with a characteristic impedance of 50 ohms, and 10 half-wavelength resonators.
为了提高通带的选择性,将滤波器阶数提升到10阶,并在谐振器2和谐振器9之间添加负交叉耦合,引入传输零点。如图10所示,设置参数v 1和h 1来控制负交叉耦合。该滤波器的中心频率为f 0= 1.025 GHz,带宽为 0.02 GHz。In order to improve the selectivity of the passband, the order of the filter is increased to 10, and a negative cross-coupling is added between resonator 2 and resonator 9 to introduce a transmission zero. As shown in Figure 10, set the parameters v1 and h1 to control the negative cross-coupling. The filter has a center frequency of f 0 = 1.025 GHz and a bandwidth of 0.02 GHz.
其快速优化的步骤与实例1相同。Its quick optimization steps are the same as Example 1.
首先依据滤波器综合理论得到理想的耦合矩阵如图11所示,建立滤波器的初始电路拓扑结构,定义9个电路拓扑结构几何变量x= {q 1、L 1、L 2、h 1、g 1、g 2、g 3、g 4、g 5}为输入向量,如图10所示。变量h 1用于调整谐振器2和9之间的耦合强度。通过sonnet仿真,得到频率响应,提取出实际电路的耦合矩阵,并与理想的耦合矩阵进行比较从而对滤波器结构参数进行初步优化,得到结果最大回波损耗至-13 dB,不能满足应用要求,第一次优化的结果如图12所示。Firstly, according to the filter synthesis theory, the ideal coupling matrix is obtained as shown in Figure 11, and the initial circuit topology of the filter is established, and nine circuit topology geometric variables x = { q 1 , L 1 , L 2 , h 1 , g 1 , g 2 , g 3 , g 4 , g 5 } are input vectors, as shown in Figure 10. The variable h1 is used to adjust the coupling strength between resonators 2 and 9. Through sonnet simulation, the frequency response is obtained, the coupling matrix of the actual circuit is extracted, and compared with the ideal coupling matrix, the filter structure parameters are initially optimized. The result is that the maximum return loss is -13 dB, which cannot meet the application requirements. The results of the first optimization are shown in Figure 12.
继续进行第二次优化,根据第一次优化的结果确定样本的范围,并通过sonnet得到所有样本的响应曲线,并提取耦合矩阵,进行神经网络的训练。训练数据的范围定义见表四。在此范围内,共随机分布了80个训练样本和20个测试样本。训练一个包含9个输入神经元、18个隐层神经元和144个输出神经元的三层神经网络来映射输入输出关系。在本例中,生成训练和测试样本,并训练神经网络模型需要21小时。加入指标要求,在通带1.015 -1.035 GHz内执行S11<-20 dB的条件。对几何变量的每次迭代,使用神经网络模型得到的结果只需0.06 s。该模型最终给出的10阶高温超导滤波器的拓扑结构几何参数优化值为x={1.36、1.16、1.36、1.14、1.3、0.92、1.76、2.2、4.2},通过电磁仿真验证,如图13所示,电磁仿真结果与本发明优化的结果一致。最终用与例1相同的高温超导材料制作该滤波器示意图如图8所示,实测性能如图14所示。如表所3所示,采用本发明的方法,10阶高温超导滤波器电路拓扑结构设计时间极大缩短,效率明显提升。Continue to carry out the second optimization, determine the range of samples according to the results of the first optimization, and obtain the response curves of all samples through Sonnet, and extract the coupling matrix to train the neural network. The range definition of the training data is shown in Table 4. Within this range, a total of 80 training samples and 20 testing samples are randomly distributed. Train a three-layer neural network with 9 input neurons, 18 hidden layer neurons and 144 output neurons to map the input-output relationship. In this example, it took 21 hours to generate training and test samples and train the neural network model. Add the index requirement, implement the condition of S11<-20 dB in the passband 1.015-1.035 GHz. For each iteration of the geometric variables, the results obtained using the neural network model only need 0.06 s. The optimal geometric parameters of the topological structure of the 10th-order high-temperature superconducting filter given by the model are x = {1.36, 1.16, 1.36, 1.14, 1.3, 0.92, 1.76, 2.2, 4.2}, which are verified by electromagnetic simulation, as shown in the figure As shown in 13, the electromagnetic simulation result is consistent with the optimized result of the present invention. Finally, the schematic diagram of the filter made of the same high-temperature superconducting material as in Example 1 is shown in Figure 8, and the measured performance is shown in Figure 14. As shown in Table 3, by adopting the method of the present invention, the design time of the topological structure of the 10-order high-temperature superconducting filter circuit is greatly shortened, and the efficiency is obviously improved.
表3 8节滤波器设计中本发明方法与常规电磁仿真所用时间比较Table 3 Comparison of the time used by the method of the present invention and the conventional electromagnetic simulation in the design of the 8-section filter
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。The above is only a preferred embodiment of the present invention, but the scope of protection of the present invention is not limited thereto, any person familiar with the technical field within the technical scope disclosed in the present invention, according to the technical solution of the present invention Any equivalent replacement or change of the inventive concepts thereof shall fall within the protection scope of the present invention.
Claims (9)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310605221.0A CN116663476A (en) | 2023-05-26 | 2023-05-26 | A High-Order High-Temperature Superconducting Filter and Its Circuit Topology Rapid Design Method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310605221.0A CN116663476A (en) | 2023-05-26 | 2023-05-26 | A High-Order High-Temperature Superconducting Filter and Its Circuit Topology Rapid Design Method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116663476A true CN116663476A (en) | 2023-08-29 |
Family
ID=87711158
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310605221.0A Pending CN116663476A (en) | 2023-05-26 | 2023-05-26 | A High-Order High-Temperature Superconducting Filter and Its Circuit Topology Rapid Design Method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116663476A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117634151A (en) * | 2023-11-07 | 2024-03-01 | 广州大学 | Two-layer optimization algorithm applied to microwave filter coupling matrix transformation |
CN118627572A (en) * | 2024-08-12 | 2024-09-10 | 国网山东省电力公司莱芜供电公司 | A neural network optimization method, system, device and medium for microwave filter |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3086470A1 (en) * | 2015-04-21 | 2016-10-26 | Airbus Defence and Space GmbH | Method for creating a coupling matrix for tuning filters and device for tuning filters |
CN115526107A (en) * | 2022-10-11 | 2022-12-27 | 东南大学 | Physical driving machine learning method for super-surface intelligent design |
CN116111984A (en) * | 2022-12-06 | 2023-05-12 | 中国电子科技集团公司信息科学研究院 | Filter design optimization method and device, filter, equipment and medium |
-
2023
- 2023-05-26 CN CN202310605221.0A patent/CN116663476A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3086470A1 (en) * | 2015-04-21 | 2016-10-26 | Airbus Defence and Space GmbH | Method for creating a coupling matrix for tuning filters and device for tuning filters |
CN115526107A (en) * | 2022-10-11 | 2022-12-27 | 东南大学 | Physical driving machine learning method for super-surface intelligent design |
CN116111984A (en) * | 2022-12-06 | 2023-05-12 | 中国电子科技集团公司信息科学研究院 | Filter design optimization method and device, filter, equipment and medium |
Non-Patent Citations (9)
Title |
---|
RICHARD J. CAMERON ET AL: "Microwave Filters For Communication Systems", vol. 2, 31 December 2018, JOHN WILEY & SONS,INC., pages: 250 - 264 * |
刘元甫: "基于耦合矩阵诊断的滤波器设计方法研究", 中国优秀硕士学位论文全文数据库 信息科技辑, no. 3, 15 March 2016 (2016-03-15), pages 135 - 1187 * |
吴章敏;包晓军;刘远曦;李琳;: "一种基于场路联合设计SIW滤波器的快速方法", 无线通信技术, no. 04, 15 December 2019 (2019-12-15), pages 27 - 30 * |
吴群: "微波工程技术", vol. 1, 30 September 2005, 哈尔滨工业大学出版社, pages: 164 * |
左涛等: "高温超导滤波器研制", 南开大学学报(自然科学版), vol. 39, no. 4, 31 August 2006 (2006-08-31), pages 70 * |
曹良足: "微波滤波器的设计与优化及案例分析", vol. 1, 31 December 2020, 西安电子科技大学出版社, pages: 92 * |
杜斌等: "工业智能化基础知识——建模、仿真、案例及代码", vol. 1, 30 September 2021, 冶金工业出版社, pages: 70 * |
杨虹等: "宽阻带高温超导带通滤波器的设计", 电子元件与材料, vol. 33, no. 1, 31 January 2014 (2014-01-31), pages 31 * |
田洪亮等: "基于灰狼算法的复合左右手超宽带高温超导滤波器设计", 2023年全国微波毫米波会议, 14 May 2023 (2023-05-14), pages 2 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117634151A (en) * | 2023-11-07 | 2024-03-01 | 广州大学 | Two-layer optimization algorithm applied to microwave filter coupling matrix transformation |
CN118627572A (en) * | 2024-08-12 | 2024-09-10 | 国网山东省电力公司莱芜供电公司 | A neural network optimization method, system, device and medium for microwave filter |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN116663476A (en) | A High-Order High-Temperature Superconducting Filter and Its Circuit Topology Rapid Design Method | |
CN113128119B (en) | Filter Reverse Design and Optimization Method Based on Deep Learning | |
CN109783905B (en) | Intelligent debugging method of microwave cavity filter based on particle swarm optimization algorithm | |
CN108170950A (en) | Multilayer Frequency-Selective Surfaces absorbing material modeling optimization method based on neural network | |
CN104811156B (en) | A kind of matching process for harmonic intensified broadband impedance matching circuit | |
Das et al. | Lossy coupling matrix synthesis approach for the realization of negative group delay response | |
Ohira et al. | Automated microstrip bandpass filter design using feedforward and inverse models of neural network | |
Rahola | Optimization of matching circuits for antennas | |
Lu et al. | High order narrow band superconducting filter design based on neural networks and extracted coupling matrix | |
CN115481591A (en) | A Neural Network-Based Extraction Method of Microwave Filter Coupling Matrix | |
CN110083958A (en) | A kind of negative group delay circuitry of three frequencies and design method based on T-connection and coaxial cable | |
Guney et al. | Input resistance calculation for circular microstrip antennas using adaptive neuro-fuzzy inference system | |
CN207098030U (en) | A kind of group delay modulator | |
CN116205082A (en) | A Design Method of Asymmetric Dual-Passband Filter Based on Frequency-variable Coupling Structure | |
CN113435147B (en) | Lumped parameter type negative group delay circuit and chip | |
CN212783739U (en) | T-type coupler-based band-pass negative group delay circuit | |
CN115859798A (en) | Design method and device of N-frequency unequal power divider based on BP neural network | |
Chen et al. | An Optimization Method for Wideband Bandpass Filters Design Using Autoencoder and Grey Wolf Optimizer | |
CN111934638A (en) | IC-shaped low-loss negative group delay circuit based on coupled microstrip line and its realization method | |
CN111490744A (en) | I L-shaped low-loss negative group delay circuit based on coupling microstrip line and design method | |
Shang et al. | Design and Fabrication of Ultra-Narrow Band High-Temperature Superconducting Filter with 1 MHz Bandwith | |
CN107508025A (en) | A kind of group delay modulator | |
Zhang et al. | A novel wideband power divider with multi-transmission zeros and poles | |
CN111611758B (en) | A Ti-type low-loss negative group delay circuit and its design method | |
CN210041761U (en) | Matching circuit for improving performance of broadband power amplifier |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |