Section: New Results

Matching problems and their applications

Participants : Laurent Baratchart, Martine Olivi, David Martinez Martinez, Fabien Seyfert.

This is collaborative work with Stéphane Bila (Xlim, Limoges, France), Yohann Sence (Xlim, Limoges, France), Thierry Monediere (Xlim, Limoges, France), Francois Torrès (Xlim, Limoges, France).

Filter synthesis is usually performed under the hypothesis that both ports of the filter are loaded on a constant resistive load (usually 50 Ohm). In complex systems, filters are however cascaded with other devices, and end up being loaded, at least at one port, on a non purely resistive frequency varying load. This is for example the case when synthesizing a multiplexer: each filter is here loaded at one of its ports on a common junction. Thus, the load varies with frequency by construction, and is not purely resistive either. Likewise, in an emitter-receiver, the antenna is followed by a filter. Whereas the antenna can usually be regarded as a resistive load at some frequencies, this is far from being true on the whole pass-band. A mismatch between the antenna and the filter, however, causes irremediable power losses, both in emission and transmission. Our goal is therefore to develop a method for filter synthesis that allows us to match varying loads on specific frequency bands, while enforcing some rejection properties away from the pass-band.

Figure 5 shows a filter with scattering matrix $S$, plugged at its right port on a frequency varying load with reflexion parameter $L_{1,1}$. If the filter is lossless, simple algebraic manipulations show that on the frequency axis the reflexion parameter satisfies:

The matching problem of minimizing $|G_{1,1}|$ amounts therefore to minimize the pseudo-hyperbolic distance $\delta$ between the filter's reflexion parameter $S_{2,2}$ and the load's reflexion $L_{1,1}$, on a given frequency band. On the contrary enforcing a rejection level on a stop band, amounts to maintaining the value of $\delta (L_{1,1},S_{2,2})$ above a certain threshold on this frequency band. For a broad class of filters, namely those that can be modeled by a circuit of $n$ coupled resonators, the scattering matrix $S$ is a rational function of McMillan degree $n$ in the frequency variable. The matching problem thus appears to be a rational approximation problem in the hyperbolic metric.

Approach based on interpolation

When the degree $n$ of the rational function $S_{2,2}$ is fixed, the hyperbolic minimization problem is non-convex and led us to seek methods to derive good initial guesses for classical descent algorithms. To this effect, if $S_{2,2}=p/q$ where $p$, $q$ are polynomials, we considered the following interpolation problem $𝒫$: given $n$ frequency points $w_1\cdots w_n$ and a transmission polynomial $r$, to find a monic polynomial $p$ of degree $n$ such that:

where $q$ is the unique monic Hurwitz polynomial of degree $n$ satisfying the Feldtkeller equation

which accounts for the losslessness of the filter. The frequencies $\left(w_k\right)$ are perfect matching points, as $\delta (S_{2,2}\left(w_k\right),L_{1,1}\left(w_k\right))=0$ holds, while the real zeros $\left(x_k\right)$ of $r$ are perfect rejection points (i.e. $\delta (S_{2,2}\left(x_k\right),L_{1,1}\left(x_k\right))=1$). The interpolation problem is therefore a point-wise version of our original matching-rejection problem. The monic restriction on $p$ and $q$ ensures the realisability of the filter in terms of coupled resonating circuits. If a perfect phase shifter is added in front of the filter, realized for example with a transmission line on a narrow frequency band, these monic restrictions can be dropped and an interpolation point $w_{n+1}$ added, thereby yielding another interpolation problem $\widehat{𝒫}$. Our main result, states that $𝒫$ as well as $\widehat{𝒫}$ admit a unique solution. Moreover the evaluation map defined by $\psi \left(p\right)=\left(p/q\left(x_1\right),\cdots ,p/q\left(x_n\right)\right)$ is a homeomorphism from monic polynomials of degree $n$ onto ${𝔻}^n$ ($𝔻$ the complex open disk), and ${\psi }^{-1}$ is a diffeomorphism on an open, connected, dense set of ${𝔻}^n$. This last property has shown crucial for the design of an effective computational procedure based on continuation techniques. Current implementation of the latter tackles instances of $𝒫$ or $\widehat{𝒫}$ for $n=10$ in less than $0.1sec$, and allows for us a recursive use of this interpolation framework in multiplexer synthesis problems. We presented these techniques at the European Microwave Week 2015 in the workshop dedicated to "Recent Advances in the Synthesis of Microwave Filters and Multiplexers". The detailed mathematical proofs can be found in [21] and will be submitted shortly. On a related topic, namely the de-embedding of filters in multiplexers, our work has been published in [13] .

Uniform matching and global optimality considerations

The previous interpolation procedure provides us with a matching/rejecting filtering characteristics at a discrete set of frequencies. This can serve as a starting point for heavier optimization procedures where the matching and rejection specifications are expressed uniformly over the bandwidth. Although the practical results thus obtained have shown to be quite convincing, we have no proof of their global optimality. This led us to seek alternative approaches able to asses, at least in simple cases, global optimality of the derived response. Following the approach of Fano and Youla, we considered the problem of a designing a $2\times 2$ loss-less frequency response, under the condition that a specified load can be "unchained" from one of its port. This classically amounts to set interpolation conditions on the response at the transmission zeros of the Darlington extension of the load. When the load admits a rational representation of degree 1, and if the transmission zeros of the overall system are fixed, then we were able to show that the uniform matching problem over an interval reduces to a convex minimization problem with convex constraints over the set of non-negative polynomials of given degree. In this case, which is already of some practical interest for antenna matching (antenna usually exhibit a single resonance in their matching band which is reasonably approximated at order 1), it is therefore possible to perform filter synthesis with a guarantee on the global optimality of the obtained characteristics. Procedures to derive the solution are currently being investigated, and lie at the heart of our contribution to the ANR-project Cocoram.