Research Article  Open Access
Design of Microwave Multibandpass Filters with Quasilumped Resonators
Abstract
Design of RF and microwave filters has always been the challenging engineering field. Modern filter design techniques involve the use of the threedimensional electromagnetic (3D EM) solvers for predicting filter behavior, yielding the most accurate filter characteristics. However, the 3D EM simulations are time consuming. In this paper, we propose electriccircuit models, instead of 3D EM models, suitable for design of RF and microwave filters with quasilumped coupled resonators. Using the diakoptic approach, the 3D filter structure is decomposed into domains that are modeled by electric networks. The coupling between these domains is modeled by capacitors and coupled inductors. Furthermore, we relate the circuitelement values to the physical dimensions of the 3D filter structure. We propose the filter design procedure that is based on the circuit models and fast circuitlevel simulations, yielding the element values from which the physical dimensions can be obtained. The obtained dimensions should be slightly refined for achieving the desired filter characteristics. The mathematical problems encountered in the procedure are solved by numerical and symbolic computations. The procedure is exemplified by designing a triplebandpass filter and validated by measurements on the fabricated filter. The simulation and experimental results are in good agreement.
1. Introduction
Traditional approach in filter design has several distinct stages: specification, approximation, synthesis, implementation, and study of imperfections. Main filter characteristics, such as filter type, central frequency, bandwidth, insertion loss, and return loss, are given in the specification stage. Mathematical expressions of physically realizable frequency characteristics are formulated in the approximation stage. Development of an electric circuit with ideal elements, from the approximation stage, is carried out in the synthesis stage. The electric circuit is transformed into the physical device in the implementation stage for a desired technology. Study of imperfections is the final filter design stage when the physical structure is tuned or optimized to compensate for various parasitic effects [1, 2].
Particular implementation of the filter elements depends on the chosen technology. If the planar technology is used there are various approaches for filter implementation. These planar filters have been implemented as microstrip filters with openloop resonators [3], microstrip interdigital filters [4], microstrip SIR (stepped impedance resonator) filters [5], duallayer filters [6], dualmetalplane patchviaspiral filters [7], multilayer filters with apertures on the ground plane [8], and multilayer stacked LC resonator filters [9]. For the above mentioned filter design examples, a common approach is to find relations between geometrical characteristics of the element implementations and the corresponding element values, such as inductances, capacitances, resistances, coupling coefficients, and quality factors (factor). These relations are typically generated by extensive 3D EM simulations but can be obtained experimentally as well. Practically, these relations are presented graphically as a set of curves describing, for example, quality factor as a function of coupling gaps [5], coupling as a function of interresonator spacing and quality factor as a function of tappedline feed position [9], resonant frequency as a function of resonator dimensions [7], and electric, magnetic, and mixed coupling as a function of apertures in the ground plane [10]. Using these data, initial dimensions and mutual positions of element implementations are chosen.
The initial implementation guess, obtained from the set of curves, is refined iteratively in order to meet the filter specifications. The analysis of the actual filter structure and its parts is performed by 3D EM simulators, such as IE3D [11], yielding the most accurate filter characteristics. However, the 3D EM simulations are time consuming and require extensive computer resources, such as CPU speed and the operating memory size. Moreover, the filter parameter refinement often requires repeated 3D EM simulations which further extend the simulation time.
In this paper, we focus on electriccircuit models, instead of 3D EM models, suitable for design of RF and microwave multibandpass filters with quasilumped coupled resonators. The circuit models are composed of linear timeinvariant passive lumped elements: capacitors, inductors, coupled inductors, and resistors. Using the diakoptic approach, the 3D filter structure is decomposed into domains that are modeled by linear timeinvariant passive electric networks with finite number of elements. Electromagnetic coupling between the domains is modeled by capacitors (for electrical coupling) and coupled inductors (for magnetic coupling). Compared to the 3D EM simulators, the circuitlevel simulators analyze the filter almost instantaneously. Therefore, repeated circuit simulations can be efficiently used to identify the key filter parameters which predominantly affect the filter characteristics. Consequently, these key parameters can be effectively used to finetune the filter performance and meet the specifications [12, 13].
In recent papers, electriccircuit models for filter design have been often used, as in [5, 7, 14]. They have been used to model filters and to analyze the influence of various filter elements on the filter characteristics. Prior to deploying a circuit model, 3D models of the filter have already been made and elements of the circuit model have been extracted from that 3D model. Circuit model is actually based on the 3D model that is to be implemented.
In this paper, the concept of circuit model is used as a primary tool to analyze and design a filter. A structure of the quasilumped coupledresonator bandpass filter is proposed, for which the electriccircuit model is derived. The filter structure is decomposed and every filter 3D component is modeled by lumped electric elements. These elements are analytically described and relations between the electricelement values and the geometrical dimensions of the filter structure are established in software [15]. Using circuit models, three single bandpass filters are designed. After fine tuning, these filters are analyzed by 3D EM solvers. Finally, a triplebandpass filter is designed using the three single bandpass filters. The filter is simulated using the 3D EM model and validated by measurements on the fabricated filter. The simulated results and experimental results are in good agreement.
2. ElectricCircuit Model of Bandpass Filter
In this section, we propose a novel microwave bandpass filter design based on coupled resonators. The resonators are implemented with quasilumped elements, using the planar multilayer technology. This approach has been inspired by filter design techniques reported in [7, 9, 16].
We have made a comprehensive analysis of one type of the multilayer quasilumped resonator and the corresponding feeders [17, 18]. From the design viewpoint, the filter implementation using these multilayer resonators is very flexible, but the 3D EM simulation is time consuming. In this work, one of our objectives is to reduce the simulation time as much as possible.
The multilayer quasilumped resonator is shown in Figure 1. The resonator is realized using doublesided microstrip. The mid layer is the common ground plane.
The resonator contains three elements: a spiral inductor, a viahole, and a parallelplate capacitor. The via connects the inner end of the spiral inductor at the top plane and the capacitor at the bottom plane. The outer end of the spiral inductor is open line. The via passes through the structure without electrical connection to the common ground plane.
The points 1 to 3 are marked on the resonator structure (Figure 1) in order to decompose the resonator into domains that are modeled by electrical networks. The points of the filter structure correspond to electriccircuit nodes. An electriccircuit model of the resonator is shown in Figure 2.
In this circuit model, the squarespiral inductor is modeled with electric network between nodes 2 and 3, and the element parameters are inductance (), resistance (), and parasitic capacitances ( and ) with their resistances ( and ). The interspiral capacitance () is capacitance between the turns of the spiral inductor.
The parameters of the cylindrical viahole are inductance () and resistance () between nodes 1 and 3.
The parameters of the capacitor are capacitance () and resistance () between nodes 1 and 0.
Neglecting the via inductance and all circuit losses modeled by resistors, the parallel resonant frequency of the unloaded resonator (Figure 2) is which is determined by the inductance and capacitances , , , and .
This resonator is used as a main building block for the bandpass filter design. In this research, the secondorder bandpass filter is designed, using coupled resonators, as shown in Figure 3. The nodes are marked on the filter resonators (1 to 4 for the first and 1′ to 4′ for the second resonator), as shown in Figure 3.
(a)
(b)
There are two possibilities of the interresonator coupling. In the first case, we have inductortoinductor coupling on one side of substrate and capacitortocapacitor coupling on the other side. This type of coupling is shown in Figure 3(a) and it will be denoted as coupling (meaning inductor to inductor and capacitor to capacitor ). In the second case, we have inductortocapacitor coupling on both sides of the layers. This type of coupling is shown in Figure 3(b) and will be denoted as coupling. These two possibilities require different models of the interresonator coupling. Two electric circuits are shown in Figure 4, corresponding to these two cases of interresonator coupling.
(a)
(b)
The coupling is modeled by capacitor (with capacitance ) between two patch capacitors (nodes 1 and 1′). The capacitance between the end of two inductors (nodes 2 and 2′) is modeled by capacitor . The mutual inductance () is modeled with inductive coupling between the inductors. The sign of mutual inductance () can be positive or negative, depending on the mutual inductor orientation. Symmetric orientation of the inductors is also depicted in Figure 3(a). If we adopt the dot marking convention as in Figure 4(a), then is positive for the symmetric inductor orientation and negative for the asymmetric inductor orientation.
For the interresonator coupling, there is dominantly capacitive coupling between inductor and capacitor, and it is modeled by capacitance (between nodes 1 and 2′, 1′ and 2), as in Figure 4(b).
The feeder is realized as a microstrip line whose width corresponds to the characteristic impedance of 50 Ω (). The feeder is capacitively coupled to the resonator and it is dominantly coupled to the outer end of the resonator inductor. This is modeled by capacitance between nodes 2 and 4 (2′ and 4′).
3. Analytic ClosedForm Expressions for Quasilumped Element Values
In order to analyze and design filter characteristics, using the electriccircuit model, it is necessary to make a correlation between the circuitelement values and the physical dimensions of the 3D filter structure. For that purpose, we use closedform expressions for inductance, capacitance, and resistance and a software package for finding matrix parameters of the transmission lines [15]. In this research, the substrate used for filter design is RT/Duroid 5880 (, , substrate thickness mm, metallization thickness μm).
The capacitance of a rectangular parallelplate structure must be calculated considering the effect of fringing field. This is important when the size of the gap separating the electrodes (i.e., substrate thickness) is comparable to the dimensions of the capacitor plates. In that case, the capacitance cannot be calculated as for a simple model of parallelplate devices. The effective capacitance is greater in this case and should be calculated as [19] where and are outer capacitor dimensions (in our case ) and substrate thickness . The capacitance of the patch capacitor with respect to the capacitor size is show in Figure 5.
The inductance of cylindrical via hole in microstrip can be calculated as a modified inductance of a cylindrical conductor and it will be calculated using the following formula [20]: For the chosen substrate and intended via radius (), the obtained value is nH (length mm and radius mm). The via resistance is neglected.
The relation between the squarespiral inductance and the squarespiral physical dimensions is given by the following formula [12]: The parameter denotes the number of inductor turns, is the outside diameter and is the inside diameter of the inductor, is the fill ratio, and is a correction factor for the presence of a ground plane. is calculated as follows: Using formula (4) one can find the spiral line width (), the turn spacing (), and the outside diameter () for the required squarespiral inductance. The approximate maximum number of the inductor turns (for given ) is , in which case . Using this formula and setting in (4), we can derive formula for the maximum theoretically achievable inductance for given , , and : This formula can be useful for the first estimation of the inductor size. The squarespiral inductance with respect to the number of turns is shown in Figure 6, without considering physical limitations during fabrication.
The parasitic capacitance of the inductor toward ground plane may be calculated using (2) with . Although layout of the inductor is different from the layout of the capacitor, our research has shown that (2) can be used with great accuracy in the case of greater inductor fill ratio. This is certainly correct in the case of our filter design objectives, where one of the goals is miniaturization. In our model, the parasitic capacitance of the inductor toward ground plane is equally distributed to and ; that is, .
The interspiral capacitance causes parallel resonances of the spiral inductor and it has influence on effective inductance at higher frequencies. The satisfactory estimation of the capacitance value is obtained with [21] where is the effective relative permeability for the spiral line of width and the given substrate. The parameter is the average spiral length.
The inductor resistance () is calculated by multiplying the spiral length with per unit length resistance of microstrip transmission line of width . Per unit length resistance is frequency dependent, and it is obtained using a software package for calculating matrix parameters of the transmission lines [15]. This dependence is shown in Figure 7.
The mutual inductance between two squarespiral inductors can be calculated using approach presented in [12]. The method is based on formula for the mutual inductance for two adjacent sections of lengths and , of the same width and mutual spacing (Figure 8(a)): where all dimensions are in microns. The parameter is the effective length of the two sections between which the mutual inductance is being calculated. As an approximation, can be considered as an average length of these two sections, .
(a)
(b)
The method for calculating the total mutual inductance is shown in Figure 8(a). The total mutual inductance is calculated as a sum of mutual inductances of each pair of parallel sections, positioned opposite to each other. There is mutual inductance between each segment of the left spiral inductor (denoted from to ) and each opposite segment of the right spiral inductor (denoted from to ). In this case, we assume that is the number of complete turns. Mutual inductances of nonopposite segments are neglected. For the adopted current direction, positive and negative mutual inductances are marked. The total mutual inductance can be calculated as follows: where is given by (8) and for to , to . In the case of the asymmetrical inductor orientation, the total mutual inductance has the opposite algebraic sign. Formula (9) is primarily used to obtain the maximum achievable value of the mutual inductance, because a greater mutual inductance value can cause appearance of transmission zeros. As an illustration, the coupling coefficient () is shown in Figure 8(b), with respect to the distance between spiral inductors (), for several inductor sizes () and .
The capacitances between the resonators and capacitance between the feeder and the resonator can be found using a software package for calculating matrix parameters of the coupled transmission lines [15]. This approach can be explained by the details shown in Figures 9 and 3.
The capacitance is equal to the mutual capacitance (sometimes referred to as the coupling capacitance) of a short section (length ) of coupled lines of the same width (), mutually separated by the distance (Figure 3(a)). The capacitance is equal to the mutual capacitance of a short section (length ) of coupled lines of the same width (), mutually separated by the distance (Figure 3(a)).
The capacitance is equal to the mutual capacitance of a short section (length ) of coupled lines of widths and , mutually separated by the distance (Figure 9, detail 1).
The capacitance is the sum of the capacitance between the central part of the feeder and the resonator () and the capacitance between the side part of the feeder and the resonator (). The first one is equal to the mutual capacitance of a short section (length ) of coupled lines of widths and , mutually separated by the distance (Figure 9, detail 2). The second one is equal to the mutual capacitance of a short section (length ) of coupled lines of the same width (), mutually separated by the distance (Figure 9, detail 3).
The values of these capacitances are presented in Figures 10 and 11, for dimensions of interest.
We are able to estimate the values of each element in the circuit model and the corresponding physical dimensions using formulas and graphs presented in this section.
4. Influence of CircuitElement Values on Filter Frequency Response
Using the electriccircuit model, it is now possible to make fast analysis of the influence of every element on the filter frequency response. It has been shown that the variation of the spiral inductor parameters has the most significant influence on the resonant frequency. The resonant frequency is mostly determined by the spiral inductance and capacitances , , and . The value of capacitance also influences resonant frequency, but it is less significant. The influence of the via becomes more pronounced at higher frequencies. The interspiral capacitance also has significant influence on filter selectivity. Greater increases filter selectivity but also decreases resonant frequency. A possible way to increase this value is to use greater number of spiral turns and smaller turn spacing . As expected, , , , , and determine filter bandwidth, insertion loss, and return loss. The above observations imply that the main resonator element for adjusting filter characteristics is the spiral inductor. Variation of the inductor parameters changes resonant frequency and selectivity simultaneously.
By analytically solving [22] the circuit model for the resonator (Figure 4(a) with and , or Figure 4(b) with ), the following formula for the resonant frequency is obtained:
The most important phenomenon is that filters using two different types of interresonator coupling have different overall characteristics. interresonator coupling always generates transmission zeros, while interresonator coupling does not. This effect for and interresonator coupling is shown in Figure 12. Using interresonator coupling, there is a negative electric coupling between the capacitors (through ) and positive or negative magnetic coupling between the inductors (through ). If these two couplings have the opposite signs and has sufficiently high value, transmission zeros can emerge. The opposite signs of and can be provided only with the symmetric inductor orientation in the resonators. Sufficiently high positive value of can be achieved using a different patchviaspiral approach [7]. However, the resonator realization in [7] has grounded inductors which results in a greater magnetic coupling . Using the equivalent circuit proposed here, it can be found that, in the case of interresonator coupling, the value of would be needed in order to have transmission zeros. Such a great value cannot be implemented using the proposed resonator layout.
In the case of interresonator coupling, the existence of zeros is caused by another phenomenon—the resonance of the inductor and certain capacitors. The resonator equivalent circuit proposed here has a complex inductor model. Therefore, the frequencies at which the transmission zeros occur can be found using the electriccircuit model and modern mathematical tools for expression simplification [22]. This yields to the following formulas for interresonator coupling (Figure 4(b)):
According to formulas (11) it is evident that the resonance of the spiral inductor and the surrounding capacitors is causing the appearance of the transmission zeros.
5. Filter Design Using ElectricCircuit Model
In this section, we present the novel procedure for the bandpass filter design using the electriccircuit model (from Figure 4(b)) and relations between lumpedelement values and their physical dimensions (Section 3). The procedure steps are as follows.(1)In the first step, lumpedelement values of the electriccircuit model are adjusted to meet the designed central frequency and the bandwidth. The most efficient way is to adopt initial dimensions of the resonator in terms of and () and the dimensions of the spiral inductor and . In our case, with the goal of miniaturization and available milling technique [23], we will adopt dimensions of the squarespiral inductor mm and mm. This will give us approximate values of , , , and . The only elements now left for adjusting are for the central frequency adjustment; for return loss; , , , and for insertion loss and bandwidth adjustment. We obtain values of , , , , and after adjusting frequency response for the circuit model. Although this overall model gives the value of , this inductance is notably affected by the influence of feeder capacitance and interresonator capacitances. The observed value of is actually the effective spiral inductance. Therefore, we propose to estimate the actual value of by observing unloaded resonator. The actual value of the inductance will give us actual value of the spiral inductor size. This approach means removing interresonator coupling by setting , , and reducing the influence of feeder by decreasing value of , while observing parameter of the circuit model. should be decreased until it has no influence on resonant frequency (observing the parameter). Now, setting the value of to adjust the desired resonator central frequency will give us the actual value of the spiral inductance. Tuned this way, the spiral inductance is found to be the most accurate.(2)In the second step, feasibility of the obtained lumpedelement values is checked using formulas and design curves. If necessary, the first step is repeated and values are modified in order to find feasible solution. In a few iterations, the first variant should be easily achieved.(3)The third step is realization of 3D EM filter model, performing full 3D EM simulation and comparing the obtained and specified frequency responses.(4)In the final step, necessary modifications are performed on the 3D EM model and the final simulation and validation has to be done.
In the rest of the text, interresonator coupling will be considered only, because it provides realizable transmission zeros.
Using the presented approach, we have designed three bandpass filters, whose frequency responses are shown in Figure 13. These filters are designed for central frequencies of 1.6 GHz (BandI), 3.5 GHz (BandII), and 5.8 GHz (BandIII), with fractional bandwidths (FBW) of 4% (BandI), 4.5% (BandII), and 7.5% (BandIII). These bands are used for GPS navigation, WiMAX access, and WiFi, respectively.
With adopted resonator dimensions mm, mm, and mm (superscripts I, II, and III will denote parameters for 1.6 GHz, 3.5 GHz, and 5.8 GHz bands, resp.) corresponding lumpedelement values are given in Table 1. The via inductance and all circuit losses are neglected, except the squarespiral inductor resistance ().

Using the values from Table 1 and formulas for finding physical dimensions of the resonator, the distance between the feeder and the resonator and the interresonator distance can be found: , , , mm, mm, mm, and mm. Using these dimensions, three single bandpass filters have been realized and simulated. The obtained results, along with the corresponding filter layouts, are shown in Figure 14.
(a)
(b)
(c)
(d)
(e)
(f)
In the next step, it was necessary to make small adjustments with respect to the central frequency, so small modifications of the squarespiral inductor length were made. In addition, adjustments were made to the distance between resonators. It may seem that the interresonator distance has been changed notably during this first modification for BandII and BandIII. Actually, the proposed circuit model implies that both the interresonator coupling and the coupling between the feeder and the resonator have been achieved with the first segments of the spiral inductors. The obtained values of , , , , and correspond to this model. This means that the values for couplings will be more accurate for greater spiral lengths, that is, lower frequency bands. For this modification, we decided to implement the interresonator coupling with the third segment of the spiral, for BandI, and the fourth segment of the spiral, for BandII and BandIII. This decision is motivated by the fact that the coupling with the third and the fourth segment of the spiral allows us to achieve the same coupling with smaller interresonator distance [18, 24]. In addition, this modification reduces the footprint. Coupling with the third segment, instead of the fourth one, for BandI, is used because it provides lower level of parasitic resonance. Frequency responses of the three modified bandpass filters generated by 3D EM and circuitlevel simulations are shown in Figures 15(a), 15(c), and 15(e). 3D EM models of the three bandpass filters are shown in Figures 15(b), 15(d), and 15(f). As can be seen, satisfying results have been achieved after only a few modifications. We must emphasize that each of these modifications as well as simulations has taken little time comparing to time necessary to design a filter without knowing the initial dimensions. Comparison of the responses in Figure 15(a) shows discrepancy with respect to the parasitic resonance that occurs at 4.95 GHz for 3D EM model, which is not present for the electriccircuit model. The reason for the appearance of this resonance is spiral inductor length, which behaves as a transmission line at 1.6 GHz. The spiral length is approximately 36 mm and, at the frequency of 4.95 GHz, its electrical length is approximately . Since the spiral is not grounded at the outer end, the impedance seen at node 3 (Figure 1) is equal to zero at 4.95 GHz. This can be validated by replacing the spiral with the line (at 1.6 GHz) in the circuit model and observing appearance of this parasitic resonance. Because of its simplicity, our model does not include this phenomenon.
(a)
(b)
(c)
(d)
(e)
(f)
6. TripleBandpass Filter Design
By incorporating three single bandpass filters into one triple bandpass filter, we obtained the 3D EM model shown in Figure 16(a). The filter dimensions are shown in Figure 16(b). Due to additional interband influences, small adjustments have been made and the final dimensions have been achieved: mm, mm, mm, mm, mm, mm, and mm.
(a)
(b)
The blockdiagram for the triplebandpass filter is shown in Figure 17(a). Frequency responses of the triplebandpass filter, obtained by 3D EM and circuitlevel simulations, are shown in Figure 17(b).
(a)
(b)
It can be noticed that the parasitic resonance (marked with “P” in Figure 17(b)) of BandI resonators (at 4.95 GHz) has disappeared. This is due to the vicinity of the two transmission zeros and . Upper transmission zero of BandII and lower transmission zero of BandIII have suppressed this effect of unwanted resonance. As it was earlier stated, this resonance is determined by the dimensions of the spiral inductor of BandI. It is possible to adjust spiral dimensions (length and width) to satisfy required central frequency and to simultaneously tune the parasitic resonance to be canceled by the transmission zeros of the other bands.
By comparing the positions of the transmission zeros , , and (Figures 15(a), 15(c), and 15(e)) with their positions for the triplebandpass filter (Figure 17(b)), it is shown that they have been shifted toward the corresponding passbands. This happens for each of the three passbands, when a triplebandpass filter is realized. The zero shifting occurs because each filter, far enough from its passband, behaves predominantly as a capacitor connected between its input and output ports. This phenomenon is also observed by the circuitlevel simulations.
In this final case, the circuit model is used for verification only, not for the design. The values of the circuit lumped elements for the three blockdiagrams (in Figure 17(a)) are tuned to match the frequency response obtained by the 3D EM simulation of each filter. This tuning has provided the actual circuitelement values. The spiral inductance and the interspiral capacitance have been tuned to match transmission zeros and retain resonant frequency. The spiral inductances have been changed to nH, nH, and nH (initial values: nH, nH, and nH) and the interspiral capacitance to pF, pF, and pF (initial values: pF, pF, and pF). Additionally, capacitance (interband capacitance) has been taken into account. This is capacitance between capacitors of the firstband and the secondband resonators and also between the secondband and the thirdband resonators. Its value is approximately pF. Three individual circuit models comprised into one circuit describe the overall filter characteristic. It should be noted that the circuitlevel simulations, using the proposed circuit model, capture all relevant phenomena important for the filter design procedure.
Disagreement between the desired frequency response (obtained by using the circuit models) and the response obtained by 3D EM simulations is greater for higher frequencies, that is, for higher bands. Two reasons explain this behavior. The first one is that our model is based on the frequency independent lumped elements. The second one is that our implementation of feedertoresonator coupling and resonatortoresonator coupling is based on the outer end of the spiral inductor.
7. TripleBandpass Filter Fabrication and Experimental Results
In the previous section, the procedure for the triplebandpass filter design has been presented and satisfactory results have been obtained. Following the frequency response of the 3D EM simulations, the triplebandpass filter is fabricated. The filter is fabricated on RT/Duroid 5880 substrate (, , mm, and μm). Final filter dimensions (see Figure 16(b)) are as follows (all in millimeters): , , , , , , , , , , and . The photograph of the filter is shown in Figure 18.
The fabricated filter has a footprint of , where is the guided wavelength of a 50 Ω microstrip line at the central frequency of the first passband (about 7.2 mm × 8.4 mm, not including the feed lines).
The proposed size is among the best achieved results (smallest footprint) compared to different design approaches. In single layer techniques, various triplebandpass filters are summarized in [25], with the best achieved result of 11.2 mm × 9.3 mm () which confirms the benefits in size reduction of the multilayer approach. Our achieved size is better compared to similar approaches with twolayer substrate 31 mm × 29 mm () [26] and 31.4 mm × 29 mm () [27]. Slightly better results are presented in [28] with triplebandpass filter size of 2.9 mm × 2.4 mm () but with the use of an additional substrate and conductive layer, with total of three substrates () and four conductive layers.
The filter is measured using Agilent N5227A PNA microwave network analyzer. The simulated and measured responses of the proposed filter are in good agreement, as shown in Figure 19.
Transmission zeros are located at the frequencies of 1.5, 1.9, 3.3, 4.1, 5.5, and 7 GHz. Measured 3 dB bandwidths are 4%, 4.6%, and 7.5%, insertion losses are 2.56, 2.68, and 2.85 dB, and return losses are 15, 14, and 12 dB at frequencies of 1.6, 3.5, and 5.6 GHz, respectively. Unwantedfrequency suppression is greater than 30 dB between the bands and greater than 20 dB above the third passband. The most notable difference between the simulated responses and the measured ones is represented by the third passband ripple. This is attributed to the manufacturing process during which vias have been manually soldered to the spiral inductor and patch capacitor. More precise machining process (or tuning the interresonator spacing at the existing model) would certainly lead to the better matching of the obtained results.
8. Conclusions
In this paper, we have proposed the use of electriccircuit models for faster design of the multibandpass filters with quasilumped coupled resonators. In addition, we proposed a novel filter which is realized as a multilayer structure. Using the diakoptic approach, the 3D filter structure is decomposed into domains, and each domain is modeled by linear timeinvariant passive electrical elements. The element characteristics have been analytically described, while relations between their electrical and physical parameters have been established.
The fast circuitlevel simulators have been used to investigate the influence of each element on the filter characteristic. Two types of interresonator coupling have been investigated, and . The limitations of the achievable capacitances, inductances, and mutual inductances have been found. This has yielded the constraints and guidelines in the design of a bandpass filter with the required characteristics. The values of electriccircuit elements have been mapped to the physical dimensions of the filter elements. This procedure has been used to independently design three bandpass filters which could be combined into one tripleband bandpass filter. Using 3D EM simulations, the obtained tripleband filter structure has been analyzed and modified and the final filter model has been established. The triplebandpass filter has been fabricated and measured and the results obtained by simulation and measurement are in good agreement. This triplebandpass filter has also achieved very small size due to the multilayer implementation (7.2 mm × 8.4 mm, i.e., ).
During the filter design, the number of timeconsuming 3D EM simulations has been reduced significantly, compared to the traditional filter design approach. 3D EM simulations have been used at the last stage of the filter design with the element values obtained from the electriccircuit model. Only a few additional simulations have been needed at the end for fine tuning. Along with the very good results regarding the filter size reduction, we have shown that a significant time saving could be achieved, in the filter analysis and design, by means of modeling the 3D EM structure by an accurate electric circuit.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This research is supported by the Ministry of Education, Science and Technological Development of the Republic of Serbia under Grant TR32005.
References
 J.S. Hong, Microstrip Filters for RF/Microwave Applications, John Wiley & Sons, Hoboken, NJ, USA, 2nd edition, 2011.
 M. D. Lutovac, D. V. Tošić, and B. L. Evans, Filter Design for Signal Processing Using MATLAB and Mathematica, Prentice Hall, Upper Saddle River, NJ, USA, 2001, Translated in Chinese, Publishing House of Electronics Industry (PHEI), Beijing, China, 2004.
 Q.X. Chu and H. Wang, “A compact openloop filter with mixed electric and magnetic coupling,” IEEE Transactions on Microwave Theory and Techniques, vol. 56, no. 2, pp. 431–439, 2008. View at: Publisher Site  Google Scholar
 D. G. Swanson Jr., “Narrowband microwave filter design,” IEEE Microwave Magazine, vol. 8, no. 5, pp. 105–114, 2007. View at: Publisher Site  Google Scholar
 R. Dhakal and N.Y. Kim, “A compact symmetric microstrip filter based on a rectangular meanderedline stepped impedance Resonator with a tripleband bandstop response,” The Scientific World Journal, vol. 2013, Article ID 457693, 7 pages, 2013. View at: Publisher Site  Google Scholar
 X. Luo, J.G. Ma, and E.P. Li, “Compact high selectivity and wide stopband bandpass filter using duallayer structure,” in Proceedings of the 41st European Microwave Conference (EuMC '11), pp. 460–463, October 2011. View at: Google Scholar
 S.C. Lin, C.H. Wang, and C. H. Chen, “Novel patchviaspiral resonators for the development of miniaturized bandpass filters with transmission zeros,” IEEE Transactions on Microwave Theory and Techniques, vol. 55, no. 1, pp. 137–146, 2007. View at: Publisher Site  Google Scholar
 J.S. Hong and M. J. Lancaster, “Aperturecoupled microstrip openloop resonators and their applications to the design of novel microstrip bandpass filters,” IEEE Transactions on Microwave Theory and Techniques, vol. 47, no. 9, pp. 1848–1855, 1999. View at: Publisher Site  Google Scholar
 C.H. Chen, C.H. Huang, T.S. Horng et al., “Very compact stacked lc resonatorbased bandpass filters with a novel approach to tune the transmission zeros,” IEEE Microwave and Wireless Components Letters, vol. 19, no. 5, pp. 293–295, 2009. View at: Publisher Site  Google Scholar
 A. Djaiz and T. A. Denidni, “A new compact microstrip twolayer bandpass filter using aperturecoupled SIRhairpin resonators with transmission zeros,” IEEE Transactions on Microwave Theory and Techniques, vol. 54, no. 5, pp. 1929–1936, 2006. View at: Publisher Site  Google Scholar
 IE3DFullWave EM Simulation, Optimization, and Synthesis Package, Zeland Software, Inc., Fremont, Calif, USA, 2007.
 I. J. Bahl, Lumped Elements for RF and Microwave Circuits, Artech House, Norwood, Mass, USA, 2003.
 R. J. Cameron, C. M. Kudsia, and R. R. Mansour, Microwave Filters for Communication Systems: Fundamentals, Design and Applications, John Wiley & Sons, Hoboken, NJ, USA, 2007.
 L.Y. Ren, “Triband bandpass filters based on dualplane microstrip/DGS slot structure,” IEEE Microwave and Wireless Components Letters, vol. 20, no. 8, pp. 429–431, 2010. View at: Publisher Site  Google Scholar
 Microwave Office, Applied Wave Research, El Segundo, Calif, USA, 2009, http://www.appwave.com/.
 M. M. Potrebić and D. V. Tošić, “A novel design of a compact multilayer resonator using doublesided microstrip,” Optoelectronics and Advanced Materials: Rapid Communications, vol. 6, no. 34, pp. 441–445, 2012. View at: Google Scholar
 D. Miljanović, M. Potrebić, D. V. Tošić, and Z. Stamenković, “Feeder realization for quasilumped multilayer resonators with low Qfactor,” in Proceedings of the 11th WSEAS International Conerence on Circuits, Systems, Electronics, Control & Signal Processing (CSECS '12), pp. 63–68, December 2012. View at: Google Scholar
 D. M. Miljanović, M. M. Potrebić, D. V. Tošić, and Z. Stamenković, “Design of miniaturized bandpass filters using quasilumped multilayer resonators,” Journal of Circuits, Systems and Computers, vol. 23, no. 6, Article ID 1450083, 21 pages, 2014. View at: Publisher Site  Google Scholar
 M. Hosseini, G. Zhu, and Y.A. Peter, “A new formulation of fringing capacitance and its application to the control of parallelplate electrostatic micro actuators,” Analog Integrated Circuits and Signal Processing, vol. 53, no. 23, pp. 119–128, 2007. View at: Publisher Site  Google Scholar
 M. E. Golfarb and R. A. Pucel, “Modeling via hole grounds in microstrip,” IEEE Microwave and Guided Wave Letters, vol. 1, no. 6, pp. 135–137, 1991. View at: Publisher Site  Google Scholar
 R. Melati, A. Hamid, L. Thierry, and M. Derkaoui, “Design of a new electrical model of a ferromagnetic planar inductor for its integration in a microconverter,” Mathematical and Computer Modelling, vol. 57, no. 12, pp. 200–227, 2013. View at: Publisher Site  Google Scholar
 Mathematica, Wolfram Research, Inc., Champaign, Ill, USA, 2013, http://www.wolfram.com.
 MITS Electronics, PCB Prototyping Machine MITS FP21TP Precision, MITS Electronics, Tokyo, Japan, 2008, http://www.mitspcb.com.
 J.S. Hong and M. J. Lancaster, “Couplings of microstrip square openloop resonators for crosscoupled planar microwave filters,” IEEE Transactions on Microwave Theory and Techniques, vol. 44, no. 11, pp. 2099–2109, 1996. View at: Publisher Site  Google Scholar
 H.W. Liu, Y. Wang, X.M. Wang et al., “Compact and high selectivity triband bandpass filter using multimode steppedimpedance resonator,” IEEE Microwave and Wireless Components Letters, vol. 23, no. 10, pp. 536–538, 2013. View at: Publisher Site  Google Scholar
 Y.F. Chen, S.H. Huang, H.W. Wu, and H.Y. Lee, “Design of compact triband bandpass filter using multilayer substrate technique,” in Proceedings of the AsiaPacific Microwave Conference (APMC '13), pp. 527–529, Seoul, Republic of Korea, November 2013. View at: Publisher Site  Google Scholar
 P.Y. Hsiao and R.M. Weng, “A triband bandpass filter using 2T4L2Cpattern multilayer structure,” in Proceedings of the Asia Pacific Microwave Conference (APMC '09), pp. 1416–1418, Singapore, December 2009. View at: Publisher Site  Google Scholar
 C.H. Chen, C.H. Huang, T.S. Horng, and S.M. Wu, “Highly miniaturized multiband bandpass filter design based on a stacked spiral resonator structure,” IEEE Transactions on Microwave Theory and Techniques, vol. 60, no. 5, pp. 1278–1286, 2012. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Dejan Miljanović et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.