Методами узагальненог матриц розсювання i часткових областей розв'язано задачу розсЮ-вання хвилi - Н10 на структурi, утвореног лан-цюжком зв'язаних хрестоподiбних резонаторiв, як частково заповнюють позамежний хвиле-вiд по висотi. Отримана багатомодова модель, яка верифшована експериментом. Використання новог моделi в методi ттелектуального синтезу i оптимiзацiг конструкций багатоланкових ХДР-фiльтрiв дозволяв враховувати довЫьну кшь-тсть хвиль в позамежному хвилеводi, що знижув похибку проектування
Ключовi слова: НВЧ-фшьтр, хвилеводно^е-лектричт резонатори, хрестоподiбний хвилеви),
узагальнена матриця розсювання
□-□
Методами обобщенной матрицы рассеяния и частичных областей решена задача рассеяния волны - Н10 на многозвенной структуре, образованной цепочкой связанных крестообразных резонаторов, частично заполняющих запредельный волновод по высоте. Полученная многомо-довая модель верифицирована экспериментом. Использование новой модели в методе интеллектуального синтеза и оптимизации конструкций многозвенных ВДР-фильтров, позволяет учитывать произвольное число волн в запредельном волноводе, что снижает погрешность проектирования
Ключевые слова: СВЧ-фильтр, волново-дно-диэлектрические резонаторы, крестообразный волновод, обобщенная матрица рассеяния -□ □-
UDC 007:159.955:519.768:621.372.852:621.372.413
I DOI: 10.15587/1729-4061.2016.70340 |
DEVELOPMENT OF A MATHEMATICAL MODEL OF MICROWAVE FILTER BASED ON THE PARTIALLY FILLED CROSS-SHAPED WAVEGUIDE-DIELECTRIC RESONATORS
D. Mamedov
Research Assistant* E-mail: [email protected] A. Yushchenko
Candidate of Physical and Mathematical Sciences, Professor* E-mail: [email protected] *Department "Systems of Information" National Technical University «Kharkiv Polytechnic Institute» Bagaliya str., 21, Kharkiv, Ukraine, 61002
1. Introduction
Ultra-high frequency range is traditionally used for designing information channels in space technology, computer networks and mobile providers' networks, also in radars and defense complexes. The need for the super wide-band and super fast transmission of information led to transferring ultramodern telecommunication technologies into the UHF range [1, 2]. The characteristics of radio-relay stations substantially depend on the electrical parameters of multitier band pass filters, installed in their transceivers. Among the known microwave filters, it is expedient to consider the designs on the basis of the partially filled waveguide-dielectric resonators (WDR) with the use of leucosapphire or quartz, due to such high quality indicators as: its own high quality, rare spectrum parasitic oscillations, high level of delivered power [3]. The advantages of these filters also include technological feasibility in the UHF-range.
2. Analysis of scientific literature and the problem statement
The presence of both strict multimode electrodynamic models and efficient methods of optimization is neces-
sary for solving the problem of optimization of the filters designs. Among the latter, preferable are contemporary methods of the artificial intelligence, such as genetic algorithms [4], artificial neural networks [5] and, in particular, expert systems. The key idea of the proposed method of intellectual synthesis and optimization of multitier WDR-filters is the logical analysis of frequency response, which is calculated on the basis of the mathematical model [6]. At the same time, this method, when applied to partially-filled cross-shaped resonators, relies on approximate calculation, which in program implementation takes into account an arbitrary number of waves at the cross point of a regular waveguide and the cut-off waveguide, and limited - in the cross-shaped [7]. Therefore, to reduce an error of the calculations, there appears a necessity to develop more efficient mathematical models, on the basis of the methods of the generalized scattering matrix (GSM) [8, 9], with the use of the method of transverse resonance [10] and numerical-analytical PRM in combination with the method of decomposition [11], allowing for the creation of numerical algorithms, which potentially consider an arbitrary number of modes on all scatterers. Given approaches are advisable to take as the basis of the electrodynamic model of WDR filters with cross-shaped cross-section.
©
3. The purpose and objectives of the study
The purpose of this study is an error reduction in the known solution of the problem of the H10 wave scattering on the chain of partially-filled WDR with cross-shaped cross-section, by developing of a new multimode model on the basis of the GSM method.
To achieve the set goal, the following tasks are to be solved:
- to solve the problem of the scattering by the GSM method, using physically substantiated assumptions;
- to implement numerically the obtained solution and to study the convergence of calculations;
- to conduct experimental verification of the calculated data.
4. Mathematical model of microwave filter with cross-shaped dielectric resonators, partially filled on height
Let us examine a planar junction (Fig. 1) of a cross-shaped layered (port 1) and hollow rectangular (port 2) waveguides. Horizontal dotted line y = 0 corresponds to ideal electrical wall in the plane of symmetry, z=0 - ideal magnetic wall. The planes of symmetry correspond to the symmetry of the basic H10 mode of uniform rectangular waveguide, considered as excitatory. Considering the symmetry, let us confine ourselves to the examination of modes in one fourth of a waveguide. To find the S-matrix of a planar junction, we will use the method of partial regions, assuming the infinite conductivity of waveguide walls and imperfection of a dielectric. Let us designate
j(0_
{a«}* and b ^{b«}*
|kz|2 < maxx2 -k0,
(1)
¿0) -<2) a
port 1 port 2
b0)
O z
Fig. 1. Planar junction of rectangular and cross-shaped layered waveguides
In Fig. 1 max x is the maximum transverse wave number of modes of the basis of a rectangular waveguide.
The field in the considered waveguide junction is represented in the form
r E > Ni Za« n_1 r E ^ HM n N1 +Zbn1) n_1 'E ^ f5[n1-) , n
< H > N2 Zai2' n_1 ' E H(2,-) n n2 +Zbn2) n_1 ' E ^ I hH
, z < 0 (port 1),
, z > 0(port 2),
(2)
where
b(;>_Ya WsM+Ya pw
n ¿-I p np ¿-I p np '
p_1 p_1
i_ 1,2, n_ 1,2,...,Nj,
(3)
snpj), i,j _ 1,2,n _ 1,2,...,N,, p _ 1,2,...,N
are the elements of S-matrix.
Here additional indices "+" and "-" determine the propagation (or damping) of wave along the axis z or in the opposite direction, accordingly. The indices N determine the number of modes of projection basis in the port number i.
Conditions of tangential components of electrical and magnetic fields at the cross point must be fulfilled:
[Ef-xs2 (x,y)E[2)_0, z _0, (x,y)eS1, [Hf-H [2) _ 0, z _ 0, (x,y)eS2,
(4)
the vectors of the amplitudes of the incident and scattered waves in the port i _ 1,2. The ratio between the number of the considered modes for the best convergence is determined by the following rule: at the master frequency, the basis of the frequency-dependent cross-shaped layered waveguide must be joined by all propagating modes and those damped modes, for which
where S1 and S2 are the considered quarters of cross sections of cross-shaped and rectangular waveguides, respectively.
xs2 (x,y)_
[1, (x,y )eS2J |0, (x,y) i S2
(5)
where max x is the maximum transverse wave number of modes of the basis of a rectangular waveguide.
n
is the characteristic function of the set S2.
The first equation in (4) is multiplied vectorially on the right by the system of functions h™, k_ 1,2,...,N1. We project the result on the axis z and integrate by the section of the comprehensive cross-shaped waveguide. Let us vectorially on the left multiply the second equation in (4) by the system of functions ek2), k _ 1,2,...,N2. We project the result on the axis o z and integrate by the section of the rectangular waveguide. As a result we obtain
J[E[1) xh(k1)]■ zds- J[E(t2) xh(k1)]■ zds _ 0, k _ 1,2,...,N1,
S, S2
J^k^H«-h(2))|zds_0, k_ 1,2,...,N2,
(6)
where tangential components of the fields are considered only in the joint's plane z _ 0.
Using the representations (2), we obtain
X Î^Ca^+Sa^jlfÊ ni-x h «].zds-
n=1 V p=1 N, / N1
P=1
-XIX a^f + IXaP2^
n=1V P=1
p=i
xJ [Ën2'+) xh«]■ zds = Xan1' J[Ëni,+) xh«]■ zds +
S2 n=1 S1
+Xan2)J[Ë n2,-)x h «l zds, k = 1,2,...,N1,
n=1 s2
XX i]XaP1^)+XaP2)Snp2)]j[ek2)x H ^^ds -
n=1 V P=1 P=1 ) s2
-XI X a(1)S(21) + X a(2)S(:
Z-| ^ P nP ^ p n
n=1V P=1
P=1
x J [42) x ii M], zds = ^an1) J [e k2) x fi nu)] ■ zds +
S2 n=1 s2
+Xan2) J [ek2) x H (n2'-)] ■ Zds, k = 1,2,...,N2.
M
g(2) ) ' h,2P+1,2q I
= 4(-1)P+qn I2P+1X
Nh,2P+1,2q I a m=0
a« + Aht
ki^L
|kJ ra^
(7)
xMimP (3,31 k<1l,2nq,o,h1
1 M.
a<IL e - Aht^
k I
MimP (3,3. )i(s)i kyi.)m,2na,o,hi
2p +1
1,2 k(2)
a(2)+a(2)
e,P h,P 1
ra^o^
2 k(2)
R(2)- B(2) K kx,P
Be,P Bh,P |kz| W
2p +1
r(3)_ r(3) k2 kX(P Be,P - Rh,P|kJ ra^o
By using the property of the orthogonality of the modes of the waveguides
q3 2b
A (2)k(2) - A (2) A__
Ae>Pkx>P Ah P |kz| ra^
i(c)( e2^,-«3,h2 )+
I(C)V k«!2^,-b/2,h( ]
I(s)i k<2L,2na,y2,h2 )-
ëi1,f)=eni), ffni,,±)=±hni)
ëifn2;±)=^^^hn2),
n n Vw(2) n
we obtain
(8) -
(9)
R(2) k(2)+ R(2) _
Re,Pkx,P + RhP|kz| ra^
R(3)k(3) + R(3) K._
Re-kx,P + RhP|kz| ra^o
i(s)| ky2m ,2nq,-y3,h2 )-
i(s)i ky3m - bAb
(12)
^a^ + ^aP2^^
V P=. N.
P=1 N.
Dk -
Mkn =
-^>/wn^[^aPl)Sn2Pl) + ^aP2)si2P2)
n=1 V P=1 P=1
= -ani)Dk + X>/W(iTan2)Mkn, k = 1,2,...,N
n=1
■ N, [ N,
VWk^X: ^a^+^aP2^'
n=1 V P=1 P=1
+:!:aP1)SkkP1) + ^a(P2)sk2P2) =
(1,2)
M„
P=1
P=1
^ -4J[en2)xh«]-
M,,„ = 4
zds,
k = 1,2,...,N1, n = 1,2,...,N2.
M(h (1),eei+i,2q )=
4(-l)P+q n I 2q X
I b X
' e,2P+1,2q
A?l+ Ail
k2 kxm
ra^o^
M<CP(a,ai )I3C)[ kyi,)m,2nq,o,hi
(10)
2P +1M
A(i) k(i) - A(i) k ' Aemkxm Ahm |kz| ra^
M<mP (a, a, )l3s)[ kyim,2nq,o,hi
q3 2b
1 2 k(2)
A(2) + A (2) k kx,P
Ae,P + Ah,P I, I
P P |kz| ra^
I^ ky2L,^,y2,h2 )-
(2) _R(2) kz kl2P
R(2) - Rl
Re,P Rh,P
^>/Wkii^ani)Mnk + fê, k = 1,2,...,N2.
where designation for the integral relation is introduced
4c)( ky2m ,2bq,-y3,h2 )-
R(3) _R(3) k2 kX3P
Re,P - Rh,P|kJ ra^o
ra^o^
I3c)| kyt^,-W2,h:
2p +1
2
A (2)k(2) - A (2) __
Ae>Pkx>P AhP|kz| ra^
i3s)V ky2m ,2nq,y2,h2 )-
(11)
R(2) k(2)+ R(2) A Re>Pkx>P + RhP|k^ ra^
Let us write out these integrals for the cases of both polarizations of the waves of a rectangular waveguide.
R(3) k(3) + R(3) K
Re,Pkx,P + Rh,P|kz| ra^o
^ ky2L,^,-y3,h2 )-
is)i ky3m ,2nq,-b2,h3
(13)
4
where
I(3c) (a, P, d, h) _ Je-1 Ima|h cos (ay) cos (p (y + d)) dy _
h, a _ p _ 0,
h . . sin (a(2h + d)) sin (ad) —cos (ad) +-------—'-
4a
2
a _ p ^ 0, sin (ah -p(h + d)) 2(a-p) + sin (ah + P(h + d)) psin(pd) + 2(a + p) + a2-p2 , a^p.
4a
(14)
S(2,2) _ (i + W21/2MtD-1MWf )-1 x
x(l - W21/2MtD-1MW21/2 ). (19)
From the third equation in the (17) we find S(1,1):
S(1,1) _ D-1MW21/2S(2i) -1. (20)
Substituting this expression into the fourth equation of the system (17), we obtain
S(2,1) _ 2(I + WfMTD-1MWf )-1 W21/2Mt. (21)
Let us define
X _ MW21/2 is the matrix of size N1 x N2, (22)
then
h - Ima h
re
I(s) (a,p,d,h)_ J-sin (ay) sin (p(y + d))dy _
0, p _ 0,
- p cos (p(h + d)) + ^ [sin (p(h + d)) - sin (pd)], a_0, p^0,
s(1,1) _ D-1XS(21)-1 = = 2D-1X (I + XTD-1X)-1 XT -1,
(23)
S(1,2)_D-1X(I + S(2,2) I_
h . \ sin (a(2h + d)) sin (ad) —cos (ad)----'--'-—-—'-
2
a _ p ^ 0,
4a
4a
(15)
1
sin (ah -p(h + d)) sin (ah + p(h + d)) a sin (pd) 2(a-p) 2(a + p) + a2-p2
ap
_ D-1X (I + (I + XTD-1X) (2I - (I + XTD-1X _ 2D-1X (I + XTD-1X)-1,
S(21)_ 2(I + XTD-1X)-1 XT,
S(2,2) _ (I+XTD-1X) (I - XTD-1X).
(24)
(25)
(26)
Let us write down the system (10) in the matrix form
(DS(11) - MWfS(21))a<1) + + (DS(12) -MW21/2S(2,2))e(2) _ -De(1) + MWfa<2), (wfMTS(11) + S(21))a(1) +
(16)
W12MTS<12) + S<2'2))a(2) _ W21/2MTa(1) + a(2).
By considering separately cases a(1) _ 0 and a(2) _ 0 , we obtain
DS(1,2) - MWfS<2'2)_ MW12, W212M tS(1,2)+S(22)_ I, DS(1,1) - MWfs'21 _ -D, W,1/2MtS(1,1) + S(21)_ w12MT.
From the first equation in (17) we find S(1,2):
S(12)_ D-1MW,1/2 (I + S(2,2)).
(17)
(18)
Substituting this expression into the second equation of the system (17), we obtain
By defining
A = I + XT D-1X, (27)
we obtain for S(2,2)
AS(2,2)_ I - XTD-1X (28)
and for the remaining blocks of the scattering matrix
AS(2,1) _ 2XT, (29)
S(1,2) _ 2D-1XA-1 _ D-1S(21)T, (30)
S(11)_ S(1,2)XT -1, (31)
For obtaining the formula (30), the symmetry of the matrix A was considered
AT _ A. (32)
Thus, for finding S-matrix it suffices to invert the matrix A of the size N2 x N2. Let us note that the elements of S-matrix do not fulfill the condition of reciprocity (S(2,1) ^ S(12)T) because the basis of a cross-shaped waveguide was not normalized according to the power (D ^ 1).
Since in reality the values of the protrusions of a cross-shaped waveguide are sufficiently low, then as the approximation of the solution of the task, let us examine a symmet-
_e
rical planar joint of the rectangular waveguide c x b with a rectangular waveguide axb of less width (a < c) and identical height. Fig. 2 displays a half of such a joint. Horizontal dot-and-dash line x = 0 corresponds to an ideal magnetic plane of symmetry.
AS(2,2) = I - XTX, AS(21) = 2XT, S(12) = 2XA-1 = S(21)T. S(1,1) = S(1,2)XT -1, where
X = W1-12 MWf, A = I + XTX,
M = {Mkn is the matrix of moments, here
a/2 b/2
Mkn = 4 J dx J[e<2)x hk1)]-zdy,
0 0
k = 1,2,...,N1, n = 1,2,...,N2.
M(h(1) e(2) lvl\uh,2p+1,2q>eh,2r+1,2s
a/2 b/2
= 4 J dx J[e<22 r+1,2s x hh1)2p+1,2q ]■ zdy =
0 0
,p+q+r+s n2b
= S„
(-1)p
Nh,2p+1,2qNh,2r+1,2s
M«(a,c) + 4£M«(a,c)
M (h(1) e(2) ) = 0
(41)
M(h(e12p+1,2q,ei
.(2)
,2p+1,2q>v-h,2r+1,2s
a/2 b/2
= s
4 J dx J[ eh,2 r+1,2s x ^¿p+Uq } zdy = 0 0
2 (_1)p+q+r+s n2
Ne,2p+1,2qNh,2r+1,2s
^ M^c)-^ M^c)
(42)
t(hC2F
Mlh«1 e(2)
a/2 b/2
4 J dx J[ x hei+1,2q } zdy =
Fig. 2. Symmetrical planar joint of two rectangular waveguides of identical height
For calculating the scattering matrix of this joint, the approach, described above, is used, with the only difference being that the ratio (9) and orthonormality condition is fulfilled for the bases of both ports. Scattering matrices in this case are calculated by the formulas
(_1)p+q+r+s n 2b = S -x
qs n N
e,2p+1,2q e,2r+1,2s
(43)
(33)
(34)
(35)
(36)
(37)
(38)
Because both joined rectangular waveguides have an identical height, the modes with different indices n along y do not interact with each other. Therefore, the efficiency of the calculation of the S-matrix of a planar junction can be increased if we examine separately the tasks of scattering for each value of the index n _ 0,2,..., followed by the assignment of found coefficients of partial S-matrices to the corresponding elements of complete S-matrix.
In contrast to the scattering matrix S, expressing amplitudes of the scattered modes through the amplitudes of the falling modes (3), the matrix of transfer T expresses the amplitudes of modes (falling and scattered) in the port 1 through the amplitudes of modes in the port 2.
fb(1) = T(11)a (2) + T(1,2)a(2), I a(1) = T(21)a(2) + T(2,2)b(2).
(44)
The blocks of the matrix of transfer can be calculated through the blocks of the scattering matrix according to the formulas
(39)
Let us give the values of relation integrals between the modes of different polarization
tM = S<1,2) _ S(1,1) [S(2,1)]-1 S(2,2) T(12 ) = S(1,1)[S(2,1)]-1,
T(2,1) = _
T(2,2) =
[S(21)]-1 S<22),
[S«21)]-1.
(45)
(46)
(47)
(48)
(40)
5. Results of calculation and experiment
Numerical algorithm was implemented programmatical-ly based on the developed mathematical model, and a study of convergence with the purpose of establishing a necessary
number of waves, considered in the cut-off waveguide, was also carried out; the diagrams of the convergence of mathematical model are represented in Fig. 3.
Fig. 3 shows that in order to ensure less than 0,2 % error in the calculation, it suffices to consider 5-7 modes in the cut-off cross-shaped waveguide for the selected parameters of the problem.
Mathematical calculation of the designs of one-tier filters and their comparison to the experimental prototypes was performed for the purpose of the verification of the developed mathematical model; the results of comparison are in Table 1; diagrams of frequency responses of one-tier filters are represented in Fig. 4,5.
Table 1
Study of experimental prototypes
Parameters of one-tier filters First Second
Calculated frequency, GHz 20,77 20,73
Experimental frequency, GHz 20,76 20,71
Width of the regular waveguide, mm 11,0 11,0
Width of the cut-off cross-shaped waveguide ("protrusion" - 0.3 mm), mm 4,5 4,5
Height of the waveguides, mm 5,5 5,5
Length of the cut-off waveguide, mm 15,0 15,0
Width of resonator, mm 4,5 4,5
Height of resonator, mm 3,2 2,9
Length of resonator, mm 2,4 3,0
Dielectric permittivity of the resonator 11,5 11,5
§ U c <-i 0.5 mm
•
•
i - 4 mm
11,0 mn 4,5 mm t
i 15.0 mm
90
19
20
23
21 22 Frequency, GHz
Fig. 4. Calculated frequency response of the first prototype
24
Conducted experimental study showed that the divergence between the calculated and measured frequencies did not exceed an error of measurements, which is a confirmation of the adequacy of the calculations, fulfilled on the base of the developed model.
The obtained mathematical model of microwave filter with cross-shaped WDR makes it possible to consider an arbitrary number of waves in the cut-off waveguide, and so to increase the accuracy of calculations. This model is implemented in the original method of intellectual synthesis and optimization of designs of multitier microwave filters [12]. The results of the study are planned to apply in the development of a new model of microwave filter with resonators of different classes.
0,8 0,7 0,6 0,5 0,4 0,3 0,2 0,1 0
1 0.4 mill
3.6 nui 1.2 mm
• 2.2 mm
10
15 N
20
25
Fig. 3. Diagrams of convergence of the microwave filter for the frequencies 73,5 and 83,5 GHz, where 5f is the frequency change in per cents, N is the number of calculated waves in the cut-off waveguide
21 22 Frequency, GHz Fig. 5. Calculated frequency response of the second prototype
6. Conclusions
1. Mathematical model of microwave filter with the cross-shaped cross section of a waveguide, partially filled with dielectric by height, was developed. This model is the solution to the problem of H10 wave scattering along the chain of partially filled WDR. The solution to the problem was carried out by applying the GSM and partial regions methods. It was shown that the applying of a matrix apparatus ultimately led to the system of equations, taking into consideration both the propagating and damped waves. Therefore, with the programmatic implementation of the model, there appears a possibility to consider an arbitrary number of waves in the cut-off waveguide, which increases the accuracy of calculations when designing filters of this type.
2. Numerical algorithm was programmatically implemented and the convergence of calculations was studied, as a result of which we established that in order to ensure the error in the calculation less than 0.2 %, it was sufficient to consider not more than 5-7 modes in the cut-off cross-shaped waveguide.
3. The experimental study of the calculated designs of one-tier filters was performed. This research showed that the divergence between the calculated and measured frequencies did not exceed one percent that is comparable to an error of measurement and is the confirmation of the adequacy of the calculations, fulfilled on the basis of the developed model.
30
References
1. Vishnevsky, V. Millimeter range as an industrial reality. The IEEE 802.15.3c standard and WirelessHD specification [Text] / V. Vishnevsky, S. Frolov, I. Shahnovich // Journal Electronica: Nauka, Technologiya, Biznes. - 2010. - Vol. 3. - P. 70-78.
2. Shao-Qiu Xiao, Ming-Tuo Zhou Millimeter wave technology in wireless PAN, LAN, and MAN [Text] / S.-Q. Xiao, M.-T. Zhou. -CRC Press, 2008. - 448 p.
3. Yushchenko, A. G. High Unloaded Qs MM Wave WDR Filters Designing [Text] / A. G. Yushchenko // International Journal of Infrared and Millimeter Waves. - 2001. - Vol. 22, Issue 12. - P. 1831-1836.
4. Thabet, R. Design of Metallic Cylindrical Waveguide Bandpass Filters Using Genetic Algorithm Optimization [Text] / R. Thabet, M. L. Riabi // Progress in Electromagnetics Research Symposium Abstracts, 2009. - P. 785-786.
5. Yahia, M. Ridged Waveguide Filter Optimization Using the Neural Networks and a Modified Simplex Method [Text] / M. Yahia, J. W. Tao, H. Benzina and M. N. Abdelkrim // International Journal of Innovation, Management and Technology. - 2010. - Vol. 1, Issue 3. - P. 259-263.
6. Yushchenko, A. G. Intellectual CAD for Three-Tier Wide Band WDR Filters [Text] / A. G. Yushchenko, D. B. Mamedov, D. M. Zaytsev // Wireless Engineering and Technology. - 2012. - Vol. 3. - P. 30-35. doi: 10.4236/wet.2012.31005
7. Ilchenko, M. E. Waveguide-Dielectric Filters Based on cross Shaped Waveguides [Text] / M. E. Ilchenko, A. G. Yushchenko, S. F. Shibalkin, V. V. Popov // International Conference on Mil-limeter and Submillimeter Waves and Applications. - 1994. -Vol. 2250. - P. 571-572.
8. Arndt, F. Theory and design of low-insertion loss fin-line filters [Text] / F. Arndt, J. Bornemann, D. Grauerholz, R. Vahldieck // IEEE Transactions on Microwave Theory and Techniques. - 1982. - Vol. 30, Issue 2. - P. 155-163. doi: 10.1109/tmtt.1982.1131041
9. Sieverding, T. Mode-matching CAD of rectangular or circular multiaperture narrow-wall couplers [Text] / T. Sieverding, U. Papziner, F. Arndt // IEEE Transactions on Microwave Theory and Techniques. - 2002. - Vol. 45, Issue 7. - P. 1034-1040. doi: 10.1109/22.598438
10. Kirilenko, A. Decomposition Approach to Multilayer Circuits Electromagnetic Modeling [Text] / A. Kirilenko, D. Kulik, Yu. Park-homenko, L. Rud, V. Tkachenko // Proc. of MMET. 2000 Conf. Kharkov, 2000. - P. 21-26.
11. Steshenko, S. A. Metod chastichnyh oblastej s uchetom osobennostej vo vnutrennih zadachah s proizvol'nymi kusochno-koor-dinatnymi granicami. Chast' 2. Plosko-poperechnye soedinenija i "in-line" ob'ekty [Text] / S. A. Steshenko, S. A. Prikolotin, A. A. Kirilenko, D. Ju. Kulik, L. A. Rud', S. L. Senkevich // Radiofizika i elektronika. - 2013. - Vol. 4 (18), Issue 3. - P. 13-21.
12. Yushchenko, A. G. Evolutionary design of seven-tier LM-mode filters optimized with original knowledge-based CAD system [Text] / A. G. Yushchenko, D. B. Mamedov // Visnyk NTU «HPY» Tehnika ta elektrofizyka vysokyh naprug. - 2014. - Vol. 21. -P. 159-171.