IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES A PUBLICATION OF THE IEEE MICROWAVE THEORY AND TECHNIQUES SOCIETY JR MTT-S NOVEMBER 2018 VOLUME 66 NUMBER 11 IETMAB (ISSN 0018-9480) THIS ISSUE INCLUDES THE JOURNAL WITHIN A JOURNAL ON MICROWAVE SYSTEMS AND APPLICATIONS REGULAR PAPERS OF THE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES EM Theory and Analysis Techniques Application of Belevitch Theorem for Pole-Zero Analysis of Microwave Filters With Transmission Lines and Lumped Elements .............................................................................................. E. L. Tan and D. Y. Heh 4669 Closed-Form Solution of Rough Conductor Surface Impedance ............................................... D. N. Grujic 4677 Devices and Modeling A GaN HEMT Global Large-Signal Model Including Charge Trapping for Multibias Operation ....................... . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . G. P. Gibiino, A. Santarelli, and F Filicori 4684 15-Gb/s 50-cm Wireless Link Using a High-Power Compact III-V 84-GHz Transmitter ................................ . . ... . . .. ... .. . . .. . .. .. .. ... .... .. ... .. . .. . . ... ... . J. Wang, A. Al-Khalidi, L. Wang, R. Morariu, A. Ofiare, and E. Wasige 4698 Nonreciprocal Components Based on Switched Transmission Lines ......................................................... . ... .. .. ....... ..... .. ..... .. A. Nagulu, T. Dine, Z. Xiao, M. Tymchenko, D. L. Saunas, A. Alu, and H. Krishnaswamy 4706 Millimeter-Wave Double-Ridge Waveguide and Components ............. S. Manafi, M. Al-Tarifi, and D. S. Filipovic 4726 Design of Frequency Selective Surface-Based Hybrid Nanocomposite Absorber for Stealth Applications ............. . ..... .. .. .. ..... .... ... . . .. ... ... ... ... .... ... . V. K. Chakradhary, H. B. Baskey, R. Roshan, A. Pathik, and M. J. Akhtar 4737 Conversion Rules Between X-Parameters and Linearized Two-Port Network Parameters for Large-Signal Operating Conditions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . R. Essaadali, A. Jarndal, A. B. Kouki, and F M. Ghannouchi 4 7 45 Passive Circuits Theoretical Analysis of RF Pulse Termination in Nonlinear Transmission Lines .......................................... . . .. .. .. . .. .... ... .. . . ... .. ... .... .. ....... .... .. ... . . ... .......... M. Samizadeh Nikoo, S. M.-A. Hashemi, and F Farzaneh 4757 Integrated Full-Hemisphere Space-to-Frequency Mapping Antenna With CRLH Stripline Feed Network .............. . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . M. D. Enders, J. H. Choi, and J. K. Lee 4765 Multiport ln-Phase/Antiphase Power Dividing Network With Bandpass Response Based on Dielectric Resonator .... ... .. ...... ... .. ........ ...... .. ..... .. .. .. ... .. ... ..... ........... .. ...... ... .......... .... .. .. ..... ...... ... W. Yu and J.-X. Chen 4773 Generalized Synthesized Technique for the Design of Thickness Customizable High-Order Bandpass Frequency-Selective Surface ............................................................. K. Payne, K. Xu, and J. H. Choi 4783 (Contents Continued on Page 4667) +.IEEE (Contents Continued from Front Cover) Asymmetrical Impedance Inverter for Quasi-Optical Bandpass Filters With Transmission Lines of Fixed Length .... ...................................................................................................... P. K. Loo and G. Goussetis 4794 A New Balanced Bandpass Filter With Improved Performance on Right-Angled Isosceles Triangular Patch Resonator ..................................................... Q. Liu, J. Wang, L. Zhu, G. Zhang, F. Huang, and W. Wu 4803 Systematic Evaluation of Spikes Due to Interference Between Cascaded Filters ............................................ ................... A. Morini, G. Venanzoni, P. M. Iglesias, C. Ernst, N. Sidiropoulos, A. Di Donato, and M. Farina 4814 Dual-ModeCharacteristicsofHalf-ModeSIWRectangularCavityandApplicationstoDual-BandFiltersWithWidely Separated Passbands ...................................................................... K. Zhou, C.-X. Zhou, and W. Wu 4820 Hybrid and Monolithic RF Integrated Circuits Analysis and Design of N-Path RF Bandstop Filters Using Walsh-Function-Based Sequence Mixing .................. .................................................................................................... A. Agrawal and A. Natarajan 4830 Compact Series Power Combining Using Subquarter-Wavelength Baluns in Silicon Germanium at 120 GHz ......... ............................................................................................. S. Daneshgar and J. F. Buckwalter 4844 0.3–14 and 16–28 GHz Wide-Bandwidth Cryogenic MMIC Low-Noise Amplifiers ........................................ .............. E. Cha, N. Wadefalk, P.-Å. Nilsson, J. Schleeh, G. Moschetti, A. Pourkabirian, S. Tuzi, and J. Grahn 4860 A 1.8–3.8-GHz Power Amplifier With 40% Efficiency at 8-dB Power Back-Off ........................................... ............................................................................... P. Saad, R. Hou, R. Hellberg, and B. Berglund 4870 Instrumentation and Measurement Techniques Nonuniformly Distributed Electronic Impedance Synthesizer ................... Y. Zhao, S. Hemour, T. Liu, and K. Wu 4883 Jitter Sensitivity Analysis of the Superconducting Josephson Arbitrary Waveform Synthesizer ........................... .......................................... C. A. Donnelly, J. A. Brevik, P. D. Dresselhaus, P. F. Hopkins, and S. P. Benz 4898 JOURNAL WITHIN A JOURNAL ON MICROWAVE SYSTEMS AND APPLICATIONS JOURNALWITHINAJOURNALPAPERS Wireless Communication Systems Effect of Out-of-Band Blockers on the Required Linearity, Phase Noise, and Harmonic Rejection of SDR Receivers Without Input SAW Filter .............................................................. A. Rasekh and M. Sharif Bakhtiar 4913 An Ultralow-Power RF Wireless Receiver With RF Blocker Energy Recycling for IoT Applications ................... ....................................... O. Elsayed, M. Abouzied, V. Vaidya, K. Ravichandran,and E. Sánchez-Sinencio 4927 AFullyIntegrated300-MHzChannelBandwidth256QAMTransceiverWithSelf-InterferenceSuppressioninClosely Spaced Channels at 6.5-GHz Band ............................... Y. Zhang, N. Jiang, F. Huang, X. Tang, and X. You 4943 A Real-Time Architecture for Agile and FPGA-Based Concurrent Triple-Band All-Digital RF Transmission ......... ................ D. C. Dinis, R. Ma, S. Shinjo, K. Yamanaka, K. H. Teo, P. V. Orlik, A. S. R. Oliveira, and J. Vieira 4955 A 0.4-to-4-GHz All-Digital RF Transmitter Package With a Band-Selecting Interposer Combining Three Wideband CMOS Transmitters ............................................................................................................... ................ N.-C. Kuo, B. Yang, A. Wang, L. Kong, C. Wu, V. P. Srini, E. Alon, B. Nikolic´, and A. M. Niknejad 4967 Extraction of the Third-Order 3×3 MIMO Volterra Kernel Outputs Using Multitone Signals ........................... ......................................................................... Z. A. Khan, E. Zenteno, P. Händel, and M. Isaksson 4985 InstantaneousSampleIndexedMagnitude-SelectiveAffine Function-BasedBehavioralModelfor DigitalPredistortion of RF Power Amplifiers .......................................................................... Y. Li, W. Cao, and A. Zhu 5000 CompositeNeuralNetworkDigitalPredistortionModelforJointMitigationofCrosstalk, I/Q Imbalance,Nonlinearity in MIMO Transmitters ...................................................... P. Jaraut, M. Rawat, and F. M. Ghannouchi 5011 Wireless Power Transfer and RFID Systems Increasing the Range of Wireless Power Transmission to Stretchable Electronics .......................................... ................................................................. E. Siman-Tov, V. F.-G. Tseng, S. S. Bedair, and N. Lazarus 5021 Bootstrapped Rectifier–Antenna Co-Integration for Increased Sensitivity in Wirelessly-Powered Sensors ............... .......................................................................................... J. Kang, P. Chiang, and A. Natarajan 5031 Microwave Imaging and Radar Applications A Linear Synthetic Focusing Method for Microwave Imaging of 2-D Objects .............................................. ................................................................................ T. Gholipur, M. Nakhkash, and M. Zoofaghari 5042 W-BandMIMOFMCWRadarSystemWithSimultaneousTransmissionofOrthogonalWaveformsforHigh-Resolution Imaging .... S.-Y.Jeon, M.-H. Ka, S. Shin, M. Kim, S. Kim, S. Kim, J. Kim, A. Dewantari, J. Kim, and H.Chung 5051 (Contents Continued on Page 4668) (Contents Continued from Page 4667) An Efficient Algorithm for MIMO Cylindrical Millimeter-Wave Holographic 3-D Imaging .............................. ................................................................................. J. Gao, B. Deng, Y. Qin, H. Wang, and X. Li 5065 A Portable K-Band 3-D MIMO Radar With Nonuniformly Spaced Array for Short-Range Localization ............... .................................................................................................................. Z. Peng and C. Li 5075 Integration of SPDT Antenna Switch With CMOS Power Amplifier and LNA for FMICW Radar Front End ......... ...................................................................................... B. Kim, J. Jang, C.-Y. Kim, and S. Hong 5087 A C-Band FMCW SAR Transmitter With 2-GHz Bandwidth Using Injection-Locking and Synthetic Bandwidth Techniques .......................................................... S. Balon, K. Mouthaan, C.-H. Heng, and Z. N. Chen 5095 A Fundamental-and-Harmonic Dual-Frequency Doppler Radar System for Vital Signs Detection Enabling Radar Movement Self-Cancellation .................................................................. F. Zhu, K. Wang, and K. Wu 5106 Microwave Sensors and Biomedical Applications Novel Microwave Tomography System Using a Phased-Array Antenna ...................................................... ........................... Y. Abo Rahama, O. Al Aryani, U. Ahmed Din, M. Al Awar, A. Zakaria, and N. Qaddoumi 5119 An RF-Powered Crystal-Less Double-Mixing Receiver for Miniaturized Biomedical Implants ........................... ................................................................................... M. Cai, Z. Wang, Y. Luo, and S. Mirabbasi 5129 IEEE MICROWAVETHEORYAND TECHNIQUES SOCIETY TheMicrowaveTheoryandTechniquesSocietyisanorganization,withintheframeworkoftheIEEE,ofmemberswithprincipalprofessionalinterestsinthefieldofmicrowave theoryandtechniques.AllmembersoftheIEEEareeligibleformembershipintheSocietyuponpaymentoftheannualSocietymembershipfeeof$24.00andobtainelectronic access, plus $50.00 per year for print media. For information on joining, write to the IEEE at the address below. Member copies of Transactions/Journals are for personal useonly. ADMINISTRATIVECOMMITTEE T.BRAZIL, President D.SCHREURS, PresidentElect A.ZHU, Secretary A.ABUNJAILEH, Treasurer Electedmembers A.ABUNJAILEH M.BOZZI T.BRAZIL N.CARVALHO G.CHATTOPADHYAY W.CHE K.GHORBANI R.GUPTA R.HENDERSON P.KHANNA D.KISSINGER S.KOUL G.LYONS M.MADIHIAN J.NAVARRO D.PASQUET G.PONCHAK S.RAMAN J.E.RAYAS-SANCHEZ S.REISING A.SANADA D.SCHREURS M.STEER Editors-in-ChiefofthepublicationsoftheIEEEMicrowaveTheoryandTechniquesSociety IEEETransactionsonMicrowave IEEEMicrowave IEEEMicrowaveandWireless IEEETransactionsonTerahertz IEEEJournalofElectromagnetics,RF IEEEJournalonMultiscaleand TheoryandTechniques Magazine ComponentsLetters ScienceandTechnology andMicrowavesinMedicineandBiology MultiphysicsComputationalTechniques J.C.PEDRO L.PERREGRINI R.CAVERLEY N.SCOTTBARKER J.STAKE J.-C.CHIAO Q.H.LIU HonoraryLifeMembers DistinguishedLecturers PastPresidents J.BARR R.SPARKS 2016–2018 2017–2019 2018–2020 D.WILLIAMS(2017) T.ITOH P.STAECKER C.CAMPBELL P.ROBLIN W.Y.ALI-AHMAD S.BASTIOLI M.GARDILL K.WU(2016) T.NAGATSUMA N.SHINOHARA N.B.CARVALHO J.EVERARD T.LEE(2015) MTT-SChapterChairs Albuquerque:E.FARR Delhi/India:S.K.KOUL LongIsland/NewYork: Phoenix:C.SCOTT/S.ROCKWELL Switzerland:N.PARRAMORA Argentina:F.J.DIVRUNO Denver:T.SAMSON/ S.PADMANABHAN PikesPeak:K.HU Syracuse:M.C.TAYLOR Atlanta:W.WILLIAMS W.N.KEFAUVER LosAngeles,Coastal:H.-Y.PAN Poland:W.J.KRZYSZTOFIK Taegu:Y.-H.JEONG Austria:W.BOESCH EasternNorthCarolina:T.NICHOLS LosAngeles,Metro/SanFernando: Portugal:R.F.S.CALDEIRINHA Tainan:C.-L.YANG Bahia:M.TAVARESDEMELO Egypt:E.A.ELH.HASHEESH J.C.WEILER Princeton/CentralJersey:A.KATZ Taipei:Y.-J.E.CHEN Baltimore:R.C.PARYANI Finland:K.HANEDA/V.VIIKARI Macau:W.-W.CHOI Queensland:M.SHAHPARI Thailand:T.ANGKAEW Bangalore/India:V.V.SRINIVASAN FloridaWestCoast:J.WANG Madras/India:V.ABHAIKUMAR RiodeJaneiro:J.R.BERGMANN Tiblisi,Rep.ofGeorgia: Beijing:W.CHEN Foothills:M.CHERUBIN Malaysia:F.C.SEMAN Rochester:J.D.MAJKOWSKI K.TAVZARASHVILI Belarus:S.MALYSHEV France:A.GHIOTTO/P.DESCAMPS Malaysia,Penang:P.W.WONG Romania:T.PETRESCU Tokyo:M.NAKATSUGAWA Benelux:G.VANDENBOSCH Germany:P.KNOTT MexicanCouncil:R.LMIRANDA Russia,Moscow:V.A.KALOSHIN Toronto:G.V.ELEFTHERIADES Boston:A.ZAI Greece:S.KOULOURIDIS Milwaukee:S.S.HOLLAND Russia,Novosibirsk: Tucson:H.XIN/M.LI Bombay/India:A.D.JAGATIA/ Guadalajara:Z.BRITO Montreal:K.WU A.B.MARKHASIN Tunisia:N.BOULEJFEN Q.H.BAKIR Gujarat/India:M.B.MAHAJAN Morocco:M.ESSAAIDI Russia,Saratov/Penza: Turkey:O¨.ERGU¨L Brasilia/Amazon:M.V.A.NUNES Guangzhou:Q.-X.CHU Nagoya:T.SEKINE M.D.PROKHOROV TwinCities:C.FULLER Buenaventura:C.SEABURY Harbin:Q.WU Nanjing:W.HONG Russia,Tomsk:D.ZYKOV UK/RI:A.REZAZADEH Buffalo:M.R.GILLETTE Hawaii:A.SINGH NewHampshire:D.SHERWOOD SanDiego:T.E.BABAIAN Ukraine,East: K.V.ILYENKO Bulgaria:M.HRISTOV Hiroshima:K.OKUBO NewJerseyCoast:A.AGARWAL SantaClaraValley/SanFrancisco: Ukraine,Kiev:Y.PROKOPENKO Canada,Atlantic:C.J.WHITT HongKong:H.W.LAI/K.M.SHUM NewSouthWales:R.M.HASHMI O.E.LANEY Ukraine,Vinnytsya:O.O.KOVALYUK CedarRapids/CentralIowa:M.K.ROY Houston:S.A.LONG NewZealand:A.WILLIAMSON Seattle:D.HEO/M.P.ANANTRAM Ukraine,West:M.I.ANDRIYCHUK Central&SouthItaly:L.TARRICONE Hungary:L.NAGY NorthItaly:G.OLIVERI Seoul:J.-S.RIEH UnitedArabEmirates:N.K.MALLAT CentralNorthCarolina:F.SUCO Huntsville:B.J.WOLFSON NorthJersey:A.PODDAR SerbiaandMontenegro: UttarPradesh/India:A.R.HARISH CentralTexas:A.ALU Hyderabad/India:Y.K.VERMA NorthernAustralia:J.MAZIERSKA Z.MARINKOVIC Vancouver:D.G.MICHELSON Centro-NorteBrasil:A.P.LANARIBO Indonesia:M.ALAYDRUS NorthernCanada:A.K.IYER/ Shanghai:J.F.MAO Venezuela:J.B.PENA Chengdu:B.-Z.WANG Islamabad:H.CHEEMA M.DANESHMAN Singapore:X.CHEN Victoria:E.VINNAL Chicago:K.J.KACZMARSKI Israel:S.AUSTER NorthernNevada:B.S.RAWAT SouthAfrica:D.DEVILLIERS VirginiaMountain:G.WILLIAMS Cleveland:M.SCARDELLETTI Kansai:T.KASHIWA Norway:Y.THODESEN SouthAustralia:C.O.FUMEAUX WashingtonDC/NorthernVirginia: Columbus:C.CAGLAYAN Kingston:C.E.SAAVEDRA/ OrangeCounty:H.J.DELOSSANTOS SouthBrazil:C.KRETLY R.R.BENOIT Connecticut:C.BLAIR Y.M.ANTAR Oregon:K.MAYS SoutheasternMichigan:A.GRBIC WesternSaudiArabia:A.SHAMIM Croatia:S.HRABAR Kitchener-Waterloo:R.R.MANSOUR Orlando:M.SHIRAZI Spain:M.FERNANDEZBARCIELA Winnipeg:P.MOJABI Czech/Slovakia:V.ZAVODNY Kolkata/India:D.GUHA Ottawa:Q.ZENG SriLanka:A.W.GUNAWARDENA Xian:X.SHI Dallas:R.PANDEY Lebanon:E.HUIJER Philadelphia:A.S.DARYOUSH St.Louis:D.BARBOUR Dayton:A.J.TERZUOLI Lithuania:I.NAIDIONOVA Peru:G.RAFAEL-VALDIVIA Sweden:M.GUSTAFSSON EditorialBoardof IEEETRANSACTIONSONMICROWAVETHEORYANDTECHNIQUES Editors-In-Chief EditorialAssistants JOSECARLOSPEDRO LUCAPERREGRINI LAURAGRIMOLDI ANARIBEIRO UniversidadedeAveiro Univ.ofPavia Italy Portugal Aveiro,Portugal Pavia,Italy AssociateEditors STE´PHANEBILA ALESSANDRACOSTANZO ANDREAFERRERO KAMRANGHORBANI ARUNNATARAJAN CHRISTOPHERSILVA XLIM Univ.ofBologna KeysightTechnol. RMITUniv. OregonStateUniv. TheAerospaceCorporation Limoges,France Bologna,Italy SantaRosa,CA,USA Melbourne,Vic.,Australia Corvallis,OR,USA ElSegundo,CA,USA XUDONGCHEN CHISTIANDAMM JOSEANGELGARCIA JUSEOPLEE HENDRIKROGIER MARTINVOSSIEK Nat.Univ.ofSingapore Univ.Ulm UniversidaddeCantabria KoreaUniv. Univ.ofGhent Univ.ofErlangen-Nu¨rnberg Singapore Ulm,Germany Santander,Spain Seoul,SouthKorea Ghent,Belgium Erlangen,Germany TA-SHUNCHU PATRICKFAY XUNGONG TZYH-GHUANGMA MIGUELANGELSANCHEZSORIANO JOHNWOOD NationalTsingHuaUniversity Univ.ofNotreDame UniversityofCentralFlorida NTUST UniversityofAlicante ObsidianMicrowave,LLC Hsinchu,Taiwan NotreDame,IN,USA Orlando,FL,USA Taipei,Taiwan Alicante,Spain Raleigh-Durham,NC,USA IEEEOfficers JAMESA.JEFFERIES, President WITOLDM.KINSNER, VicePresident,EducationalActivities JOSE´M.F.MOURA, President-Elect SAMIRM.EL-GHAZALY, VicePresident,PublicationServicesandProducts WILLIAMP.WALSH, Secretary MARTINJ.BASTIAANS, VicePresident,MemberandGeographicActivities JOSEPHV.LILLIE, Treasurer FORRESTD.“DON”WRIGHT, President,StandardsAssociation KARENBARTLESON, PastPresident SUSAN“KATHY”LAND, VicePresident,TechnicalActivities SANDRA“CANDY”ROBINSON, President,IEEE-USA JENNIFERT.BERNHARD,Director,DivisionIV—ElectromagneticsandRadiation IEEEExecutive Staff STEPHENP.WELBY, ExecutiveDirector&ChiefOperatingOfficer THOMASSIEGERT, BusinessAdministration CHERIFAMIRAT, InformationTechnology JULIEEVECOZIN, CorporateGovernance KARENHAWKINS, Marketing DONNAHOURICAN, CorporateStrategy CECELIAJANKOWSKI, MemberandGeographicActivities JAMIEMOESCH, EducationalActivities MICHAELFORSTER, Publications JACKBAILEY, GeneralCounsel&ChiefComplianceOfficer KONSTANTINOSKARACHALIOS, StandardsAssociation VACANT, HumanResources MARYWARD-CALLAN, TechnicalActivities CHRISBRANTLEY, IEEE-USA IEEEPeriodicals Transactions/JournalsDepartment SeniorDirector,PublishingOperations: DAWNMELLEY Director,EditorialServices:KEVINLISANKIE Director,ProductionServices: PETERM.TUOHY AssociateDirector,EditorialServices:JEFFREYE.CICHOCKI AssociateDirector,InformationConversionandEditorialSupport: NEELAMKHINVASARA ManagingEditor: CHRISTOPHERPERRY JournalsCoordinator: CHRISTINAM.REZES IEEETRANSACTIONSONMICROWAVETHEORYANDTECHNIQUES(ISSN0018-9480)ispublishedmonthlybytheInstituteofElectricalandElectronicsEngineers,Inc.Responsibilityforthecontentsrestsupon theauthorsandnotupontheIEEE,theSociety/Council,oritsmembers.IEEECorporateOffice:3ParkAvenue,17thFloor,NewYork,NY10016-5997.IEEEOperationsCenter:445HoesLane,Piscataway, NJ08854-4141.NJTelephone:+17329810060.Price/PublicationInformation:Individualcopies:IEEEMembers$20.00(firstcopyonly),nonmember$217.00percopy.(Note:Postageandhandlingcharge notincluded.)Memberandnonmembersubscriptionpricesavailableuponrequest.CopyrightandReprintPermissions:Abstractingispermittedwithcredittothesource.Librariesarepermittedtophotocopyfor privateuseofpatrons,providedtheper-copyfeeof$31.00ispaidthroughtheCopyrightClearanceCenter,222RosewoodDrive,Danvers,MA01923.Forallothercopying,reprint,orrepublicationpermission, writetoCopyrightsandPermissionsDepartment,IEEEPublicationsAdministration,445HoesLane,Piscataway,NJ08854-4141.Copyright(cid:2)c 2018byTheInstituteofElectricalandElectronicsEngineers,Inc. Allrightsreserved.PeriodicalsPostagePaidatNewYork,NYandatadditionalmailingoffices.Postmaster:SendaddresschangestoIEEETRANSACTIONSONMICROWAVETHEORYANDTECHNIQUES,IEEE, 445HoesLane,Piscataway,NJ08854-4141.GSTRegistrationNo.125634188.CPCSalesAgreement#40013087.ReturnundeliverableCanadaaddressesto:PitneyBowesIMEX,P.O.Box4332,StantonRd., Toronto,ONM5W3J4,Canada.IEEEprohibitsdiscrimination,harassmentandbullying.Formoreinformationvisithttp://www.ieee.org/nondiscrimination.PrintedinU.S.A. Digita l Objec tIdentifier 10.1109/TMTT.2018.2873875 This article has been accepted for inclusion in a future issue of this journal. Content is final as presented, with the exception of pagination. IEEETRANSACTIONSONMICROWAVETHEORYANDTECHNIQUES 1 Application of Belevitch Theorem for Pole-Zero Analysis of Microwave Filters With Transmission Lines and Lumped Elements Eng Leong Tan , Senior Member, IEEE, and Ding Yu Heh , Member, IEEE Abstract—This paper presents the application of Belevitch change considerably.The numbers of poles and zeros as well theorem for pole-zero analysis of microwave filters synthesized as their locations for a synthesized filter may be different with transmission linesand lumpedelements. Thescattering (S) from the originally specified or designed ones. For instance, matrix determinant ((cid:2)) based on the Belevitch theorem, aptly the poles of a synthesized Chebyshev filter may be deviated called Belevitch determinant,comprises polesandzeros that are separated in different half-plane regions. Using the Belevitch from the reference Chebyshev ellipse after realization. It is determinant, the poles and zeros of filter transfer functions to be shown later that the number of poles of a synthesized can be determined separately with certainty, e.g., by applying classical coupled-line filter [1] turns out to be more. Hence, the contour integration method based on argument principle. the number of poles and zeros should be ascertained after Note that the contour integration can be evaluated numerically synthesis,whiletheirlocationsshouldbeverifiedandchecked without requiring complicated overall analytical expressions. The proposed method is able to solve the poles and zeros for for further tuning wherever needed, especially for filters with filterssynthesizedwithnoncommensuratetransmissionlinesand transmission lines transformed from lumped element proto- lumpedelements,wherethetransformmethodandtheeigenvalue types. Since poles and zeros control the amplitude response, approach are inapplicable. Several applications are discussed phase, and group delay of a filter directly, they need to be to demonstrate the use of Belevitch theorem and the contour analyzedforbetterunderstandingandmanipulationofthefilter integration method to determine the poles and zeros of various microwave filters on the complex plane. response. Unfortunately, for microwave filters synthesized from IndexTerms—Argumentprinciple,Belevitchtheorem,contour transmission-line structures, it is often difficult to determine integration,microwavefilters,numericalmethod,polesandzeros, S matrix determinant. the poles and zeros since the overall analytical expressions maybe verycomplicateddueto the involvementof(nonpoly- nomial) transcendental functions. The difficulties in solving I. INTRODUCTION poles and zeros are further exacerbated when the filters con- MICROWAVE filters have been the subject of much sist of both transmission lines and lumped elements. Their research due to their wide applications in wireless resultant scattering (S) parameters’ expressions, if one could communication.Variousfilter synthesismethods,designs,and ever derive, would comprise combinations of various integer transformations can be found in numerous classic microwave power terms and trigonometric power terms. This makes the textbooks [1]–[4]. Over the years, many filter structures existing methods, such as Richard, Euler, or digital trans- have also been analyzed and synthesized involving (lossless) form method [1], [9], [10], and coupling matrix eigenvalue lumped elements and transmission-line structures, such as approach [11]–[15] inapplicable for the pole-zero analysis steppedlines,stubs,andcoupledlines [5]–[8].However,most of the filters. For blind modeling methods such as vector analyses have not considered in detail the order and locations fittingtechnique [16],theyalsoinvolvemanyuncertaintiesand of poles and zeros of synthesized filters. In particular, most inconsistencies of poles and zeros depending on the chosen filter analyses of microwave textbooks are subjected only for number of poles and fitting bandwidth. low-pass filter prototype and typically in lumped elements. In this paper, we present the application of Belevitch theo- When replaced by transmission-line structures after various remforthepole-zeroanalysisofmicrowavefilterssynthesized transformations, the transfer functions of these filters may with transmission lines and lumped elements. In Section II, the challenges in solving poles and zeros for such filters Manuscript received October 15, 2017; revised February 28, 2018 and will be exemplified in detail, demonstrating the difficulties June 25, 2018; accepted July 29, 2018. This work was supported in part byresearchprojects throughDSOunderGrantDSOCL12016andinpartby and deficiencies of existing methods, including transform theSingaporeMinistryofEducationTertiaryEducationResearchFundunder method, eigenvalue approach, and vector fitting technique. In Grant2015-1-TR15. (Correspondingauthor:DingYuHeh.) Section III, the Belevitch theorem will be discussed. Using The authors are with the School of Electrical and Electronic Engi- neering, Nanyang Technological University, Singapore 639798 (e-mail: the Belevitch theorem,it will be shownthat the scattering (S) [email protected]; [email protected]). matrix determinant (), called Belevitch determinant, com- Color versions of one or more of the figures in this paper are available prisespolesandzerosthatareseparatedindifferenthalf-plane onlineathttp://ieeexplore.ieee.org. Digital ObjectIdentifier 10.1109/TMTT.2018.2865928 regions. Once the poles and zeros are separated completely, 0018-9480©2018IEEE.Personaluseispermitted, butrepublication/redistribution requires IEEEpermission. Seehttp://www.ieee.org/publications_standards/publications/rights/index.html formoreinformation. This article has been accepted for inclusion in a future issue of this journal. Content is final as presented, with the exception of pagination. 2 IEEETRANSACTIONSONMICROWAVETHEORYANDTECHNIQUES the contour integration method based on argument principle will be applied, which allows them to be determined with certainty.Several applicationswill be discussed in Section IV to demonstrate the use of the Belevitch theorem and the contour integration method to determine the poles and zeros of various microwave filters on the complex plane. II. POLE-ZERO ANALYSIS OFMICROWAVE FILTERS WITH Fig.1. Schematic ofaparallel transmission-line filter. TRANSMISSION LINES AND LUMPEDELEMENTS roots may be solved directly [9], [10]. Following [9], we set For microwave filters involving transmission lines, it is Z =50, Z = Z =2Z ,andtheelectricallengthsofthetwo often challenging to solve for the poles and zeros of transfer 0 1 2 0 parallel transmission lines are θ = π/2, θ = 5π/2, which functions due to the involvement of transcendental functions 1 2 give a = 5. Substituting this a into (6) and upon changing such as the exponentialand trigonometricfunctions.Consider the variable z =u1/2, we obtainan equationthat is consistent the example from [9], where the filter section consists of two parallel transmission-line sections with characteristics with [9, eq. (3)]by notingu =ejπ(f/fd) and the same ak, bk impedances and electrical lengths Z , Z and θ , θ , respec- in [9, Table I]. 1 2 1 2 Thus far, the transform method has been applied for the tively [see Fig. 1(a)]. The expression of S for such a filter 21 case of integera to solve for polynomialroots. However,it is section is given by no longer able to simplify root solving for noncommensurate S21(ω)= NS21(ω)/DS21(ω) (1) lines or when a is not an integer. To demonstrate this, we let θ =3.6514π instead of 5π/2 in Fig. 1(a). Now, substituting where 2 a = 7.3028 into the denominator of (6) gives an equation, NS21(ω) which consists of fractional powers that can no longer be = −2jZ Z Z (Z sin(θ ω/ω )+Z sin(θ ω/ω )) (2) solved easily as polynomial roots. In fact, when a is not 0 1 2 1 1 0 2 2 0 D (ω) integer, it is also mentioned in [9] that the filter cannot be S21 analyzed within the digital-inspired (z or u) framework (e.g., = Z02Z12+ Z02Z22+Z12Z22 (sin(θ1ω/ω0)sin(θ2ω/ω0)) only approximate when θ =0.497π). 1 +Z02Z1Z2(cos2(θ1ω/ω0)+cos2(θ2ω/ω0)+sin2(θ1ω/ω0) The difficulties in solving poles and zeros of transfer + sin2(θ ω/ω )−2cos(θ ω/ω )cos(θ ω/ω )) functions are further exacerbated when the filters consist of 2 0 1 0 2 0 −2j Z Z2Z cos(θ ω/ω )sin(θ ω/ωa ) both transmission lines and lumped elements. To illustrate 0 1 2 2 0 1 0 this, the same filter section in Fig. 1(a) is added with two + Z Z Z2 cos(θ ω/ω )sin(θ ω/ω ) (3) 0 1 2 1 0 2 0 capacitors of capacitance C1 and C2, as shown in Fig. 1(b). ω0 is the center angular frequency. To convert (1)–(3) into The denominator of S21 can be derived as complex frequency s domain, where s = σ + jω, ω is D (ω) substitutedwith−j·s =−j·(σ+jω).Itisworthpointingout S21 = −105C C ω2 p4−2p2p2− j8p2p q +2p2q2 that the S21 denominator of the filter section is already rather 1 2 1 1 2 1 2 2 1 1 complicated even for one section only. As exemplified from +2p2q2− j8p p2q +24p p q q + j8p q q2+ p4 1 2 1 2 1 1 2 1 2 1 1 2 2 (1)–(3),therootscannotbedirectlysolvedforthesenonlinear +2p2q2+2p2q2+ j8p q2q +q4−2q2q2+q4 2 1 2 2 2 1 2 1 1 2 2 functions. Furthermore, as the trigonometric functions are +200ω(2p p + jp q + jp q )(2C p q − jC p p periodic, there are infinitely many poles on the whole s 1 2 1 2 2 1 1 2 1 2 1 2 − jC p p +2C p q + jC q q + jC q q ) domain. It is therefore, difficult to ascertain and solve for the 1 1 2 2 1 2 1 1 2 2 1 2 poles within certain frequency band of interest. −(2p p + jp q + jp q )2 (7) 1 2 1 2 2 1 For commensurate lines whose electrical lengths are the where same (ormultiples)forallsections,onepossible wayto solve for the poles is by utilizing the Richard transform or Euler’s θ θ p =cos 1ω/ω , p =cos 2ω/ω (8) identities [1]. In particular, using 1 2 0 2 2 0 z+1/z z−1/z θ θ cos(−jsθ /ω ) = , sin(−jsθ /ω )= (4) q =sin 1ω/ω , q =sin 2ω/ω . (9) 1 0 1 0 1 0 2 0 2 2j 2 2 za +1/za za −1/za cos(−jsθ /ω ) = , sin(−jsθ /ω )= It is evident that due to the combinations of various inte- 2 0 2 0 2 2j ger power terms (ω2, ω) and trigonometric power terms (5) [sin((θ/2)ω/ω ),cos((θ/2)ω/ω )],(7)cannotbetransformed 0 0 where a =θ /θ , S of (1) is transformed into z domain as intoanothersimplerformtosolvefortherootsdirectly.Hence, 2 1 21 thetransformmethodisinapplicableforthepole-zeroanalysis 2z2a +2za+1−2za−1−2 S (z)= . (6) of filters with transmission lines and lumped elements. 21 4z2a+1−z2a−1−2za −z For certain microwave filter synthesis [11]–[15], its cou- If a is an integer, the numerator and denominator of (6) plingmatrixmayalsobeusedtofindtheS-parameterpolesvia contain only polynomials(with integer powers) for which the eigenvalue approach. Considering, for example, an n-coupled This article has been accepted for inclusion in a future issue of this journal. Content is final as presented, with the exception of pagination. TANANDHEH:APPLICATIONOFBELEVITCHTHEOREMFORPOLE-ZEROANALYSISOFMICROWAVEFILTERS 3 lumpedelementresonatorfilter,its S and S canbederived network can be represented as [17] 21 11 using the coupling matrix as [3] S S 1 u v S = √ 2 A−1 (10a) S= S1211 S1222 = g v∗ −u∗ (13) 21 q q n1 e1 en where g is a strictly Hurwitz polynomial and g, u, and v are 2 −1 S = ± 1− A (10b) related by 11 q 11 e1 uu∗+vv∗ = gg∗. (14) where A = q+ pI − jm (11) The notation g∗ (with subscript “∗”) denotes g∗ =[g(−s∗)]∗, j ω ω wherethesuperscript“*”indicatescomplexconjugation.Note p = − 0 (12) that the effect of g∗ is to reflect all complex roots of g FBW ω ω 0 symmetricallyabouttheimaginaryaxis.Therefore,therootsof qe1 andqen arethescaledexternalqualityfactorsofinputand g and g∗ are symmetricalabout the imaginary axis. Applying outputresonators, FBW is the fractionalbandwidth,ω is the theBelevitchtheorem(13)and(14),itcanbeshownthattheS 0 center angular frequency, I is the n×n identity matrix, q is matrix determinant, aptly called Belevitch determinant, takes an n × n matrix with all entries zeros except q = 1/q the form 11 e1 and qnn = 1/qen, and m is the n ×n coupling matrix. If a = S S −S S = −g∗. (15) narrowband filter can be represented in the coupling matrix 11 22 21 12 g form of (11), the poles can be found by first solving p from Since g is strictly Hurtwiz, (15) shows that all poles of are the eigenvalues of matrix comprising m in (11). Thereafter, located in the left half-plane (LHP) of complex s plane, with the poles are obtained via (12). However, a broadbandfilter can no longer be representable thesamenumberofzeros(rootsofg∗)locatedexactlyopposite of jω axis in the right half-plane (RHP). Note that there will in the coupling matrix form of (11) should it comprise com- be no zeros in the LHP and no poles in the RHP for . As binations of noncommensuratetransmission lines and lumped such, the poles and zeros have been effectively separated into elements. Indeed, the impedance or admittance matrix or A in (11) for such a filter should have various powers of ω, different half-plane regions through (15). This is crucial to (1/ω), as well as sin and cos terms, in much the same allow them to be determined separately with certainty in the subsequent analysis. Furthermore, the Belevitch determinant manner like (7). In other words, the p term of (12) not only involves (ω/ω ) and (ω /ω), while the m matrix of (11) not in (15) is an irreducible expression due to g being strictly 0 0 Hurtwiz. Consequently, the degree of a lossless 2-port is the only contains constant coupling coefficients, but also various nonlinear functions of ω. This makes it difficult or impos- degree of the polynomial g of the canonic form (13) [17]. Byretrievingthepolesof,wehaveessentiallyretrievedthe sible to extract the poles using the eigenvalue approach same poles for S sharing the common denominator g. directly.Hence,theeigenvalueapproachisinapplicableforthe 21 Oncethe polesandzerosare separatedcompletely,one can pole-zeroanalysisoffilterswithtransmissionlinesandlumped proceed to solve for them using several possible ways. To elements. that end, the contour integration method based on argument The blind system modeling methods such as the vector principleisonepossiblewaythatmaybeappliedconveniently. fitting technique, Loewner matrix, or Cauchy method are The argument principle is given by [18], [19] sometimes used to locate the system poles and zeros through fitting of measurement data points. We shall comment briefly 1 1 f(s) on the vector fitting technique [16] that is one of the most Z − P = 2πCargf(s)= 2πj f(s)ds (16) C popularmodelingmethods.Whileitisablindfitting/modeling method, the number of poles to estimate is often uncertain where Cargf(s) is the change in argument of f(s) along and the poles obtained are very much dependent on the closed contour C in the counterclockwise direction, f(s) is a fittingbandwidth.Forafilter withnocomplexzeros,spurious meromorphicfunction, Z and P are the number of zeros and complex zeros may even appear that are not symmetrical poles of f(s), including their multiplicities within the closed about the jω axis, thus violating the S-parameters’ unitary contour, and f(s) is the derivative of f(s). The contour C conditions [11]. Henceforth, this technique and other blind alsoshouldnotpassthroughanypolesandzeros.Thecomplex fitting/modeling methods are beyond the scope of this paper. function f(s) in (16) is to be replaced by the Belevitch determinant (s) of a microwave filter on complex s plane. Since has no zeros in the LHP (Z = 0), we are able III. APPLICATION OFBELEVITCH THEOREM to evaluate (16) to find the number of poles in the LHP. It FORPOLE-ZERO ANALYSIS should be emphasized here that without the use of Belevitch To overcome the difficulties and deficiencies mentioned determinant in (15), the argument principle could not be above, we shall present the application of Belevitch theorem applied directly if f(s) contains both poles and zeros in the for the poles and zeros analysis of microwave filters syn- LHP. Most argument principle-based methods, such as [20] thesized with transmission lines and lumped elements. The and [21], are able to solve only for the zeros of f(s) when it Belevitch theoremstates that the S matrix of a lossless 2-port is analytic within C (P =0). This article has been accepted for inclusion in a future issue of this journal. Content is final as presented, with the exception of pagination. 4 IEEETRANSACTIONSONMICROWAVETHEORYANDTECHNIQUES The contour integration of f ≡ can be evaluated numerically via 1 (s) 1 (s+δs/2)−(s−δs/2) ds = . 2πj (s) πj (s+δs/2)+(s−δs/2) C s (17) Here, centralfinite differenceis used for the derivative(s), while the averaging scheme is used for (s). This reduces the number of stencils from three to two for every point on the contour. δs is the spatial step size chosen along the path to ensure convergence.It should be emphasized that to evalu- ate(17),theoverallanalyticalexpressionofisnotrequired, while only the numerical values of along the contour path are needed. In practice, for microwave filter circuits, the net- work parametersfor individualsections (such as transmission lines,lumpedelements,etc.)canbecomputedandmanipulated readily to obtain the values of the overall circuits. Fig. 2. S-parameters of a coupled-line filter with Chebyshev response Tolocatethepoles,thecontourcanbesuccessivelydivided (N =3). Inset: layout of N+1sections coupled-line filter along witheven andoddcharacteristic impedances oftheindividual section. into smaller sections until they are tightly enclosed. The contour method could even be combined with other efficient iterative rootsearching algorithms, such as Newton’s method, can then be obtained subsequently by reflecting all the Muller’s method, etc., which often require rather good ini- zeros of f(s) symmetrically along jω axis. This will bypass tial guess. Having divided the contour into small sections, the computation of 1/ inversion for greater simplicity and the contour centroid may serve as good initial guess to robustness.Apartfrom the sole determinantin (15), the Bele- these iterativealgorithms.Afterthe poleshavebeenretrieved, vitch theorem can also be applied for other S-parameter we can further determine (if any) the complex zeros of a expressions or combinations. For instance, one may consider microwave filter. This can be done by repeating (16) for f(s) = S / to solve directly for the transmission zeros 21 f ≡ S21 and dividing the contour into smaller sections (zeros of S21) or consider f(s) = S11/ or f(s) = S22/ successively to find Z with known P. to solve directly for the reflection zeros of a microwave filter. Thepolesandzerosmayalso be solvedusingthe approach All in all, many opportunities arise from the application of in [22], employing Newton’s identities. However, one would the Belevitch theorem, which could be further explored for need to estimate the number of distinct zeros and poles, the pole-zero analysis of microwave filters. which is done through forming successive system matrices and checking if they are singular. This is often difficult in IV. APPLICATIONS FOR MICROWAVE FILTERS practice to numericallyascertain whether a matrix is singular. A. Applications I: Classical Filters Furthermore, a large number of zeros and poles often result We shall demonstrate the applications of the Belevitch in high polynomial order and ill-conditioned problem. On theoremandthe contourintegrationmethodtosolve forpoles the other hand, alternative searching technique, such as [23], and zeros of several microwave filters. First, consider the involvesexhaustive triangulationof 2-D domain.Still, it does classical coupled-line filter in Fig. 2 [1]. The specifications not guarantee that all the roots can be found and the risk of of the bandpass filter are Chebyshev response with N = 3, missingrootishigherduetotheirregulardiscretizationofthe passband ripple L =0.5 dB, center frequency f =2 GHz, domain.Overall,thereare stilluncertaintiesanddifficultiesin Ar 0 andfractionalbandwidthFBW= 0.1.Thefilterisrealizedby solving for poles and zeros using these methods. Unlike the four (N +1 = 4) coupled-line sections. The S-parameters of approachabove,by using the Belevitch determinant (15), one thefilterareshowninFig.2.FromtheS plot,threereflection will have the poles and zeros separated completely, with only 11 zeros can be observed, which may lead one to deduce that poles and no zeros in the LHP. There is no need to estimate there are three poles within its passband in line with the the number of poles or zeros, thus obviating the forming of third-orderspecification.Usingthecontourintegrationmethod successive system matrices and checking if they are singular. fortheBelevitchdeterminantwithcontourpathsenclosingthe The poles can then be solved as the roots of polynomial. passband around f = 2 GHz, the poles are determined as Hence, there is no need to divide the contour excessively 0 (normalized by 1 GHz) until the poles are tightly enclosed. In practice, one would still consider dividing the contour sufficiently to reduce the s = 2π(−0.03124+ j1.89849) p1 number of enclosed poles in each contour. This should lower s = 2π(−0.06259+ j2.00000) p2 the polynomial order to avoid ill-conditioned problem. s = 2π(−0.03124+ j2.10151) The Belevitch determinant in (15) further allows us to p3 evaluate (16) on the contours in the RHP with f(s) = sp4 = 2π(−1.88540+ j2.00000) [instead of in the LHP with f(s) = 1/]. The poles of s = 2π(−1.88549+ j2.00000). (18) p5 This article has been accepted for inclusion in a future issue of this journal. Content is final as presented, with the exception of pagination. TANANDHEH:APPLICATIONOFBELEVITCHTHEOREMFORPOLE-ZEROANALYSISOFMICROWAVEFILTERS 5 passband ripple L = 0.01 dB. f is still maintained at Ar 0 2 GHz. The S-parameters are also shown in Fig. 2. Three reflection zerosare still visible fromthe S plot. Meanwhile, 11 the numberof polesdeterminedfromthe determinantcontour integrationmethodisstillfivearoundthepassband.Thepoles are subsequently solved as s = 2π(−0.15601+ j1.67367) p1 s = 2π(−0.71330+ j2.00000) p2 s = 2π(−0.58349+ j2.00000) p3 s = 2π(−0.44182+ j2.00000) p4 s = 2π(−0.15601+ j2.32633). (21) p5 For verification, the transform method is again utilized to obtain the roots in the z domain as z = 0.22427+ j0.85578 p1 z = j0.57108 p2 aFnigd.(3d.)im(aa)gRineaarlyapnadrt(sbo)fimS2a1g.inary parts of S21.Absolute errorsin(c)real zp3 = j0.63237 z = j0.70680 p4 It can be seen that there are five poles within the contour.For z = −0.22427+ j0.85578. (22) p5 verification,we also solvethe polesbyutilizingthe transform methodandfindthepolynomialroots(withpositiveimaginary Using(20),thepoles(ins domain)areinagreementwith(21). part only) in the z domain as It can be seen that all five poles are now near to each other and to jω axis. In addition, note that the two extra poles z = 0.07771+ j0.97266 p1 in (18) and (21) do not lie on/near the Chebyshev ellipse zp2 = j0.95203 on the complex s plane. While the other three poles that z = −0.07771+ j0.97266 are around the Chebyshev ellipse may contribute most to the p3 z = j0.22746 third-order Chebyshev response, the presence of two extra p4 poles would still affect the filter response to a certain extent z = j0.22744. (19) p5 (see Fig. 3). From these two examples, one should appreciate the importance of poles analysis, whereby the number and Thepolesins domaincanthenbededucedfromtheroots(19) locations of the poles of a synthesized filter may be different through the inverse of Richard transform or Euler’s identities from those originally specified. s =ω /θ ·lnz. (20) 0 0 It is found that they agree well with (18) obtained via B. Applications II: Advanced Filters our Belevitch determinant contour integration method, thus Complex zeros may sometimes be utilized to optimize validating our proposed method. the filter phase response for better group delay equalization. Interestingly, from (18), it can be seen that the last two We now apply the Belevitch theoremand the contourintegra- poles sp4 and sp5 are located further away from the jω tionmethodtoanalyzealinearphasefilterwithcomplexzeros axis. To investigate further the consequence of retaining from [3].Thefour-polelinearphasefilterhasapassbandfrom all five poles or only first three poles sp1, sp2, and sp3, 920to975MHzandcanbesynthesizedbycouplingmatrixas ⎡ ⎤ Fig. 3(a) and (b) shows the real and imaginary parts of S 21 0 0.9371 0 0.1953 cthoenystraurectecdomuspianrgedaltlofivSe21poobletsaianneddodnirleycfitlrystfrthormeespimoluelsa,tiaonnd. m =⎢⎢⎣0.90371 0.60196 0.60196 0.90371⎥⎥⎦. (23) Therealzeroslocatedonthe jωaxishavealsobeenincluded, 0.1953 0 0.9371 0 consideringaswellthenegativefrequenciesandperiodicityof this coupled-line filter. Fig. 3(c) and (d) shows the absolute The poles and zeros are solved using the contour integration errors in the real and imaginary parts of S for the results method with Belevitch determinant. The poles and zeros are 21 constructed using five and three poles compared to the simu- plotted on the complex s plane in Fig. 4, marked by “x” and lation.Itcanbeobservedthatifoneexcludesthelasttwopoles “o,”respectively.Itcanbeobservedthatfourpolesarevisible in (18), the errors are larger in both real and imaginary parts withinthepassband.Moreover,acomplexzeroisalsopresent. of S . These errors may at times be tolerable, which explain This further shows the effectiveness of our proposed method 21 whymostliteraturestill regardthe (N =3)coupled-linefilter in retrieving both poles and zeros on the complex plane. For as “third” order effectively with three poles only. verification,the polesare also calculated using the eigenvalue Let us now modify slightly the specifications of the approach.Tothatend, pcanbesolvedfromtheeigenvaluesof coupled-line filter with increased FBW = 0.2 and reduced matrix comprising (23) in (11). Then, the poles are obtained This article has been accepted for inclusion in a future issue of this journal. Content is final as presented, with the exception of pagination. 6 IEEETRANSACTIONSONMICROWAVETHEORYANDTECHNIQUES Fig.4. Polesandzeros plotofalinearphasefilter. Fig. 5. Poles and zeros plot of a linear phase filter cascaded with parallel transmissionlines andlumpedelements. via (12) and it is found that they are consistent with those TABLEI shown in Fig. 4. ROOTSOBTAINEDUSINGMULLER’SMETHODWITHDIFFERENT Consider nextthe linear phase filter of (23)being cascaded INITIALGUESSES(NORMALIZEDTO1GHz) with the parallel transmission lines and lumped elements in Fig. 1(b). The capacitances are set as C = 1 pF and 1 C = 2 pF, respectively, while θ = π/2 and θ = 3.6514π. 2 1 2 This cascaded filter couldno longer be analyzedby the trans- formmethodortheeigenvalueapproach,asitnowconsistsof transmission lines and lumped elements. Using the Belevitch theorem and the contour integration method, we are still able to determine the poles and zeros of the cascaded filter, as shown in Fig. 5. It is observed that there are additional poles and complex zeros introduced due to the cascaded parallel transmissionlinesandcapacitors.Thisexamplehasshownthe capabilities of the proposed method in solving the poles and zeros of filters with transmission lines and lumped elements. The contour integration method is useful to identify the regionthatiscertaintocontainoneormorezeros/poles.Once responses that meet different specifications of each band. In identified, one may use various efficient methods, such as our case, it is a Butterworth response in the first passband Muller’smethod,to locate the rootsforzeros/poles.Although (lowerband)andaChebyshevresponseinthesecondpassband being more efficient, the roots may at times go out of (higher band). Fig. 6 (inset) shows the photograph of the range or may not even converge as exemplified for some fabricated filter, realized by noncommensurate transmission cases in Table I. To demonstrate this, the zeros of filters lines and stubs. The simulation and measurement of the in Figs. 4 and 5aresearchedwithinthecomplexsplaneusing S-parameters have been performed and found to be in good Muller’smethodwith differentinitialguesses.Fromthetable, agreement [24, Fig. 7].Actually,ourmainconcernhereisnot we can see that if the initial guesses are not sufficiently close only the S-parameters but also the poles and zeros, as well to the zeros, the roots may go out of range, e.g., converging as whether they have been conforming to our specifications to the RHP zeros or other harmonic band zeros, or may even (e.g., different passband Butterworth/Chebyshev responses). failtoconverge.Forsufficientlycloseinitialguesses,theroots Withthepresenceofnoncommensuratetransmissionlinesand eventually converge within range and they agree well to the stubs, it is difficult to analyze such a filter for its poles using zeros in Figs. 4 and 5. Hence, this underlines the importance existingmethods.UsingtheBelevitchtheoremandthecontour of identifying the correct zeros/poles region for sufficiently integration method, the poles can be solved and plotted in closeinitialguess,whichcouldbeprovidedforbythecontour Fig. 6, showing three poles present around each passband integration method. for the third-order responses. To further verify each passband We next proceed to analyze a realized dual-band filter that response, the reference Butterworth circle (first passband) we have designed, fabricated,and measured as in [24]. In the and Chebyshev ellipse (second passband) are also drawn on design therein, the poles of each passband are properly dis- the complex s plane based on the specification. It can be tributed on the complex s plane to provide differentpassband seen that the poles of each passband lie quite closely on