Efficient Filter Design to Compensate Fabrication Imperfections in 6G Communication Systems

In this paper, we present a consistent methodology for the reliable design of 6G-oriented filters with enhanced endurance to construction imperfections. The systematic formulation does not depend on the filter’s operating frequency and employs a robust strategy for obtaining new roots and poles of the filtering function. Essentially, it requires that all the local maxima of the filtering function do not fluctuate beyond the design attenuation levels for a set of predefined roots/poles distortions. To this purpose, two novel algorithms for the derivation of the appropriate filtering functions are developed, in the prior basis, together with a versatile optimization criterion and a heuristic comparison approach that guarantee optimal outcomes. Specifically, the principal idea of the first technique is to accurately extract the roots of the new polynomial from a system of equations on condition that the maximum local peaks of the distorted (due to imperfections) initial polynomial are below a prefixed threshold, such as the unit. Conversely, the second method develops an alternative polynomial, compressed in the amplitude and frequency range, so that a similar prerequisite regarding the maximum local peaks, is satisfied. It is stressed that both methods are fully generalized and may be applied to any polynomial combination, without increasing the overall complexity. The proposed framework is successfully verified in terms of theoretical examples and the numerical simulation of realistic waveguide and mictrostrip line filters, operating at frequencies from 2GHz to 65GHz, which unveil its superiority over existing schemes and implementations.


Introduction
Millimeter wave filters, such those designed for modern 6G systems [1][2][3], require meticulous standardization and high manufacturing accuracy.This is because the sensitivity of a distributed circuit increases with frequency and thus fabrication imperfections may ruin the behavior of the resulting filters in the spectrum of wireless communications [4][5][6][7].Since the latter use printed-circuit multi-layer inserts as discontinuities in order to achieve the required performance, all the parameters and features of the substrate, such as the dielectric permittivity, the thickness, or the losses, can impact the behavior of the filter in the case of wrongly selected values [8].Also, the inaccurate operation of the machining techniques, the milling processes, the chemical etching, the laser cutting, or the metallization methods can create several significant geometric or structural imperfections [9].Amid them, one may discern the fluctuating dimensions, the incorrect digging depth in the substrate, the unwanted curvatures, and the wrongly positioned inserts [10].These imperfections can seriously affect the efficiency of 6G filters by irregularly modifying their frequency response, i.e., they actually change the their transfer function.So, the bandwidth, the amplitude response, and the central frequency of the filter are likely to deviate from the initially designed ones, thus completely degrading the performance of the device.In essence, the contribution of filters in 6G systems is decisive, since their role is to separate useful signals from noise and manage interference in the frequency spectrum, thus aiding the correct data transmission [11].Hence, any deviation in the behavior of the initially designed filtering function will immediately degrade the operation of the entire arrangement.
To evaluate the efficiency of a filter and determine potential imperfections, several key performance criteria should be satisfied.The most important are the variation of the filtering (filter transfer) function, the frequency and amplitude response of the filter, the localization of its central (operating) frequency, and its bandwidth [12].For this purpose, current research focuses on the following areas: (a) the optimization of existing filter fabrication processes (e.g., surface roughness, metallization, and selection of proper materials), (b) the redesign of typical filter architectures according to modern standards (e.g., smooth surfaces, adaptive geometries), and (c) the development of new design strategies.The latter are deemed the most promising and are divided into two classes.The techniques of the first class assess fabrication imperfections via mathematical methods (like the one proposed herein), based on the improvement of the filtering functions, the polynomial chaos, or the enhancement of stochastic techniques [13].The second class launches prototype manufacturing processes, such as the 3D printing, the additive manufacturing, the digital manufacturing, and the free-form fabrication [14].These techniques build filters by selectively adding material (layer by layer), instead of subtracting, as in usual procedures.
Bearing in mind these notions and apart from the meticulous inspection of any conventional fabrication technique, modern design filter strategies try to establish standardized steps to avoid imperfections [15].The first is the generation of the proper 3D computeraided design (CAD ) model of the filter via the respective computer software or a reverseengineering procedure [16].Then, this CAD model is converted to the appropriate stereolithography format, which describes the filter as a discretized mesh, and placed with extreme accuracy in the building device to avoid staircase approximations [17].Next, the fully automated manufacturing process starts, while some minor faults are removed via thermal treatment, sandpapering, shot-peening, surface-coating, or infiltration [18].
Recently, there has been a remarkable interest for the incorporation of emerging technologies into the compensation of filter fabrication imperfections [19,20].The most important is the development of advanced computational packages with artificial intelligence modules to estimate the impact of potential faults or even indicate means for their correction, prior to any fabrication [21].However, the use of thermoplastic and biogradable materials with unique properties is, also, a topic of tense research.Among such media, we can mention the acrylonitrile butadiene styrene, the PolyLactic acid, the PolyEtherEtherKetone, and various advanced ceramics [22].Not to mention that most 3D printing techniques are revisited in order to replace inkjet schemes with material jetting, i.e., the deposition of droplets of material through moveable printer heads [23].In addition, powder bed fusion technologies are considered to selectively fuse particles through a laser or an electron beam at specific filter areas [24].Another promising technology for 6G filters is the 3D screenprinting process, which can design net-shaped components with complex shapes, with negligible imperfections [25].Finally, photopolymerization platforms consisting of liquid resin have become fairly popular due the faultless filter surface finishing they attain [26].Nonetheless, it is emphasized that the above technologies should overcome diverse challenges, before they will become widely popular.In particular, they must properly address: (a) the high cost of the fabrication equipment, (b) the constantly increasing needs for dimensional accuracy and repeatability, (c) the requirement for new materials with advanced properties, and (d) the absence of fully established fabrication standards.
A crucial feature of a filter design technique is to be fully independent of frequency.This is exactly one of the novelties of the technique proposed in this paper and achieved by directly modifying the filtering function after the systematic analysis of the detected imperfections.On the other hand, traditional approaches attempt to adapt to different frequencies by suggesting various formulations for different spectra or approximations that could offer adequate (yet not optimal) results.It becomes apparent that low sensitivity filtering functions must be employed to compensate for the frequently encountered manufacturing tolerances [27][28][29].Considering that these functions are created by a polynomial, the Chebyshev filters (derived by the Chebyshev polynomials) exhibit the best (steepest) behavior in the out-of-band zone, yet at the expense of the highest sensitivity to mechanical defects.On the other hand, the Butterworth filters (derived by the x n function) are the less efficient (smoothest behavior) in the out-of-band zone, but, also, have the highest resistance to mechanical imperfections among all the filters created by polynomials of the same degree.Thus, the need for enhanced and flexible polynomials is escalating, particularly in the rapidly evolving area of 6G communications [30][31][32].
Two state-of-the-art approaches can be found in the relevant literature that, however, do not clearly associate the filter sensitivity to construction tolerances with the choice of the polynomial.Moreover, they do not promptly clarify which solution is the best before its post simulation.Nonetheless, they have the advantage of simplicity and provide the polynomial roots in a detailed way (a feature of the herein-proposed techniques as well).The first approach is the chained-function filter formulation, where the filtering function is the product of various Chebyshev and Butterworth polynomial combinations [33][34][35][36].According to this scheme, the most efficient combination-i.e., the one that produces the steepest plot in the out-of-band zone for a given variation of the peaks in the S 11 -parameter-is obtained via an appropriate criterion and a specific optimization algorithm [35].This variation is verified through a Monte Carlo method, where the input data is the variation of the characteristic resistances of the printed circuit components.It is noteworthy to mention that the chained functions, also, consider the case of rational functions.Conversely, in the second approach, the filtering function imposes a specific pattern on the S 11 -parameter, with low peaks at the edges of a slightly wider passband zone [37].The method is validated, again, by a Monte Carlo implementation, where the input data is the variation of the coupling coefficients of the filter circuit.Principally, both of the above techniques lead to low-sensitivity filtering functions, which, as an upshot, exhibit a lower performance; an issue that requires the design of higher-order filters to retain the same performance [34].
This work addresses the impact of fabrication imperfections on the behavior of millimeter wave filters for 6G communication systems, by, again, lowering the performance of the filtering function.Nevertheless, the key difference from the aforementioned literature is that, now, the new filtering function is systematically and fully related to these manufacturing tolerances.The proposed concept stems from the roots/poles variation of the auxiliary filtering function, caused by a structural variation, which can enforce its peaks in the bandpass zone to surpass a permissible value, namely the unit in the prototype function.Then, the roots and poles are precisely relocated on the real axis, so that, in the enhanced function, the peaks decrease in magnitude and for, the same roots/poles variation, they are not larger than (the ideal case is to be equal to) unit.Such a variation should be computed in advance and this is, herein, conducted by means of the COMSOL Multiphysics ® simulation software [38].Note that, throughout our analysis, the filtering function is presumed to be a real function of the frequency and, equivalently, the prototype filtering function is deemed a real function of the real variable x.
Based on these notions, two novel techniques for the accurate derivation of enhanced and consistent filtering functions are presented in this paper.Both algorithms are independent of the filter's operating frequency and can be used for polynomial and ratio of polynomials (elliptic-like) filtering functions, which can lead to the straightforward and fast design of narrowband low-order millimeter filters.To ensure that the derived filtering functions are optimum (i.e., they have the steepest possible out-of-band behavior, under certain optimization constraints), an instructive criterion is introduced.Furthermore, a fully unified way of comparing the various filtering functions found in the literature with the proposed ones is described.In this manner, and although the differences among the diverse polynomials could be small, the optimum filtering function is promptly obtained.
In summary, the key novelties of the paper are: (a) the development of a generalized design methodology that does not depend on the filter's central frequency and can effectively handle any fabrication imperfection, (b) the introduction of two precise schemes that can be directly applied to any filtering function and lead to significant enhancements, even for demanding imperfections, (c) the proposed formulation is much simpler, as it involves only the magnitude of the S-parameters and not the phase, and (d) the straightforward and fast extraction of the compensated filter design parameters.The new schemes are extensively validated through several real-world waveguide and microstrip line filters, with operating frequencies in the wide range between 2 GHz and 65 GHz, pertinent for 6G systems.Numerical results prove the advantages of the featured method, which exhibits a much better performance than existing filtering functions.

Development of the Polynomial Filtering Function Methodology
In this section, we formulate two versatile algorithms for the extraction of the necessary filtering function polynomial, significantly less vulnerable to fabrication tolerances, when the variation of its roots has, already, been estimated.According to the first method, we calculate the roots of the new polynomial from a system of equations, requiring that the maximum peaks of the distorted polynomial (i.e., its roots are varied) do not exceed a predetermined value, e.g., the unit.The resulting polynomial is shown to be optimal if it satisfies a specific optimization theorem.Conversely, the second method creates a Chebyshev polynomial, compressed in the amplitude and frequency range, so that, as above, the maximum peaks of the distorted polynomial do not exceed the unit.

The System-Based Optimization Method (SbOM)
Let us start from the transfer function, S 12 , of a filter where ε is the ripple factor and F n (x) a Chebyshev polynomial for the case of a narrowband microwave filter.This simply implies that Equation (1) resembles an inverted Chebyshev polynomial.Bear in mind that Chebyshev polynomials have been, initially, defined as those with the smallest possible deviation from the horizontal axis among all monomial polynomials (with the unit factor at the x n term) of the same degree n.This minimum deviation is achieved by forcing all local maxima of the polynomial to be equal to each other in absolute value.In essence, Chebyshev polynomials became popular since they can produce filtering functions that are maximized outside the interval of their roots (in the prototype problem for |x| ≥ 1), while retaining their peaks below 1, for |x| ≤ 1 [39,40].Therefore, they provide the steepest possible plot in the out-of-band zone from all other polynomials of the same degree.Apparently, this feature renders Chebyshev filters sensitive to fabrication imperfections, because they affect the roots of the polynomial.As Chebyshev polynomials exhibit the smallest possible deviation from the horizontal axis, a change of their roots causes a definite increase of, at least, one local maximum (peak) beyond 1.This means that the filtering function does not satisfy the typical design standards, e.g., it leads to |S 12 | < −0.043 dB or |S 11 | > −20 dB, for ε = 0.1.On the other hand, the roots in Butterworth filters and thus the local maxima are all at 0; so, a change of the roots does create peaks, yet there is enough space for them to remain below the maximum level of variation.
In the present work, we consider that the seeking polynomial, F n (x), has real roots in [−1, 1] (i.e., arranged as −1 ≤ x 1 ≤ • • • ≤ x n ≤ 1), its peaks in [−1, 1] do not exceed 1, and F n (x) crosses the (1, 1) point.Hence, F n (x) is given by or F n (x) crosses the (1, 1) point.For even symmetry, F n (x) can be written as Owing to mechanical imperfections, the roots of the polynomial become x i ± d i , with d i ≥ 0 the maximum variance of the ith root.For simplicity, we assume that the increment and decrement of the ith root are equal, while, below, this constraint can be relaxed.Obviously, the total number of x i ± d i combinations of the roots produces a large set of distorted polynomials whose local peaks vary in size.A local peak receives its maximum value for a specific combination.Hence, to find the proper x i roots, we require that all these maxima should not to exceed 1 in the prototype problem.This yields a system of n − 1 equations, half of which are used to set the maximum peaks to 1 and half to acquire their abscissas.It is emphasized that in the case of even symmetry, the number of equations for the maximum peaks is n/2 and for the abscissas (n/2) − 1.
Conversely, in the case of odd symmetry, these numbers are (n − 1)/2 for both the peaks and the abscissas.Furthermore, it will be numerically shown in the following paragraphs that, via an optimization criterion, the closest each local maximum peak is located to 1, the steeper the plot of the resulting polynomial is produced outside Prior to formulating the new SbOM, we focus on some important characteristics of F n (x), considering that the arrangement of its −1 ≤ x 1 ≤ • • • ≤ x n ≤ 1 roots is maintained.Specifically, we provide the subsequent propositions: Proposition 1. Denoting as i the peak between the x i and x i+1 root, it is stated that this peak increases if x i decreases (similarly, if the x i+1 root increases).Indeed, if the abscissa point of the ith peak is denoted as q i , then the distorted polynomial, G n (x, d i ), at q i , can be defined as and as x i < q i , it holds that Since the prior inequality is valid for q i , it is, also, valid for the new ξ i points that are the abscissas of the G n (x, d i ) peaks; a fact which proves the proposition.Proposition 2. For a d i small relative to the x i root, the maximum of the ith peak occurs when the x 1 . . .x i roots depart from the x i+1 . . .x n ones.The reasoning is the same as the one used to compare the ith maximum of the F n (x) with that of the G n (q i , d i ) above.It is simply, now, considered that the roots move consecutively (i.e., one at a time).Hence, the maximum of the ith peak occurs for the Proposition 3. The absolute value of every F n (x) peak is smaller than 1, since its maximum is not larger than 1.

Proposition 4.
The abscissas of all roots will be in [−1, 1], due to the fact that F n (x) crosses the (1, 1) point.Indeed, if there were an F n (x) root outside [−1, 1], then the (1, 1) point would act as a peak between two roots, implying that there would exist a combination of the roots for which the specific peak would be larger than 1; i.e., a complete contradiction to the way the roots are calculated.Proposition 5.The previously defined coefficient K is always smaller than its K Chebyshev counterpart of the Chebyshev polynomial.This is because, in terms of Proposition 3, the absolute values of the F n (x) peaks are smaller than 1 and thus the peaks of the monic F n (x)/K polynomial are below the level of the corresponding monic Chebyshev polynomial.Since the latter has the smallest possible deviation from the horizontal axis (and is unique), the point where |F n (x)/K| crosses x = 1 is higher than the maximum of the monic Chebyshev polynomial, C n (x), and less than 1.This implies that |C n (1 It can be numerically verified that |K| has its maximum value for the optimum polynomial, as explained in the next paragraphs.Similarly, we can derive that the larger the root variation, the more the new filtering function resembles a Butterworth function.Moreover, in the worst-case scenario, it is presumed that the roots fluctuate independently to each other, which may produce less steep filtering functions.However, the relation between fluctuations depends on the filter implementation.

Calculation of Maxima
In general, if d i are large compared with the roots, then we have to take into account many variations of root combinations to retrieve every maximum local peak of the polynomial.Nonetheless, for small d i , we may consider that the arrangement of x i remains and, by means of Proposition 2, obtain the ith maximum of F n (x).On the other hand, if d i are large compared with the distance between the roots or some of the roots coincide, then the auxiliary arrangement of x i can no longer remain and movements of roots along different directions may provide the local maximum peaks.To overcome this hindrance, the algorithm we employ controls the values of the local peaks along every root movement.Actually, if there are doubts for the maxima of a local peak, then other variations of root combinations can be used.Nevertheless, this happens only for a few cases in our analysis.
To extract x i , let us, now, consider that the increment, d + i , and the decrement, d − i , of the ith root differ.According to the preceding notions, for relatively small root variations, the maximum of the local ith peak occurs for the roots and the polynomial For the case of even symmetry, the corresponding root variation is while the distorted polynomial reads (6)   with 1 ≤ i ≤ (n/2) − 1, and for i = 0 (the distorted polynomial with maximum local peak between the ) and K as defined in Equation (3).
In this context, for the required x i , we must solve the next system of n − 1 equations for 1 ≤ i ≤ (n/2) − 1, and for i = 0, and for 1 ≤ i ≤ (n/2) − 1.In Equations ( 8) and ( 10), ξ i is a point where the partial derivative of 8) and ( 9) is the value of the maximum local ith peak of the n (x).Conversely, for the case of odd symmetry (i.e., x 1 = 0), we consider that d − 1 = d + 1 = d 1 and the respective root variation becomes whereas the distorted polynomial is given by for 2 ≤ i ≤ (n + 1)/2 − 1, and ) .Thus, Equations ( 8)- (10) can, now, be written as for 2 ≤ i ≤ (n + 1)/2 − 1, and for i = 1, and for 2 ≤ i ≤ (n + 1)/2 − 1, and for i = 1.Notice that in Equations ( 8), (14), and (15), we have p i ≤ 1.The reason for not imposing p i = 1 is that, generally, the system of Equations ( 8)-( 10)-or, equivalently, the system of Equations ( 14)-(17) − is not solvable for every p i = 1.Hence, in our algorithm, we, initially, set all p i = 1 and, in the end, some roots may be found equal to each other, implying that for some p i (whose calculation is not required) it holds that p i < 1.

Optimal Formulation
A critical aspect in the development of the proposed SbOM is to investigate whether the solution of Equations ( 8)-( 10) is optimal, namely if it provides the maximum possible value of |F n (x)|, for |x| ≥ 1.This is the case if the solution satisfies a specific optimization theorem.Actually, there are some proofs in the literature concerning the behavior of Chebyshev polynomials outside [−1, 1].Some of them are based on the alternation of the pronoun of the difference between the Chebyshev polynomial and another polynomial function that intersects it [39].Others rely on the representation of the Chebyshev polynomial as a Lagrange interpolation [40].Such approaches, however, can not be applied to the proposed technique, since our objective is to draw a conclusion about the behavior of F n (x) by observing for some random |x 0 | ≥ 1, i.e., a linear relationship between the roots [41].In this framework, our optimization problem opts for the minimization for some random |x 0 | ≥ 1, as, also, alternatively described in [42].
Focusing on our approach, we firstly obtain a solution from ( 8)- (10) and then examine if this is optimal, according to whether it satisfies a specific optimization theorem.Explicitly, we want to optimize the polynomial function F n (x 0 ), for x 0 > 1, whose variables are its roots In essence, the formulation, introduced herein, can be deemed as an appropriately tailored version of the optimization theorem of a constrained function [43], which uses the Lagrange λ multipliers.
Our aim is to prove that λ i > 0, when the x i are such that | max i F (i) n (x)| = p i are valid.Differently speaking, F n (x 0 ) will be maximized, when | max i F (i) n (x)| = p i apply and therefore F n (x 0 ) will be larger than the value of any other polynomial at the same point x 0 , whose local maxima are not larger than p i and at least one is smaller.Hence, we launch the representation of and the system of equations ∂L(x 1 , x 2 , . . ., x n/2 )/∂x i = 0, for 1 ≤ i ≤ n/2.Then, owing to Equations ( 8)-( 10), we derive the following system in matrix form where for j > i Consequently, when the above theorem is satisfied, namely when λ i > 0, we can reliably draw the subsequent significant deductions, whose validity does not depend (by any means) on the filter's operating frequency:

•
The polynomial in the out-of-band zone, i.e., for |x| ≥ 1, is maximized when the maxima of the local peaks receive the values of p i and is not smaller than them.

•
The polynomial in the out-of-band zone, i.e., for |x| ≥ 1, is better than any other polynomial, whose local maxima are not larger than p i and at least one of them is smaller.• Since p i can be defined arbitrarily, we can set all of them equal to 1, on the condition that Equations ( 8)-( 10) are solvable for this selection.This is the maximum normalized value and provides the maximum possible polynomial value for |x| ≥ 1, namely the steepest possible plot in the out-of-band zone.• When all the maximum peaks p i of a polynomial are one-to-one smaller than or equal to the maximum peaks p i of another polynomial and at least one p i is less than one p i , then the former polynomial is less steep than the latter in the out-of-band zone.
As presented in the next sections, we successfully apply the prior theorem to various solutions of the Equations ( 8)-( 10) system, for different d i combinations during the design of sixth-order filters.Results reveal that, in all cases, the λ i s are found to be positive and therefore our solutions are, indeed, the optimal ones.Moreover, since x i and ξ i are known analytically, it is possible that the above theorem can provide an alternative interpretation to the behavior of Chebyshev polynomials outside [−1, 1].It is stressed that the specific procedure may be used to demonstrate the advantages of the proposed filtering function over existing options, namely it can precisely produce the steepest possible out-of-band plot.Finally, we can deduce that the SbOM leads to the optimum polynomial when the d i s are none-zero, whereas the Chebyshev polynomial is the optimum one when the d i s vanish.

Solution of the Equation System
Depending on the symmetry type, we solve sequentially the previously derived system of Equations: (a) ( 8)-( 10) for x 1 , ξ 2 , x 2 , ξ 3 , x 3 , . . ., ξ n/2 , x n/2 , x 1 , . . .(even symmetry) or (b) ( 14)-( 17) for x 2 , ξ 2 , x 3 , . . ., ξ (n−1)/2 , x (n+1)/2 , ξ 1 , x 2 , . . .(odd symmetry), until the full convergence of the roots.Recall that ξ i are the abscissas of the maximum local peaks, while our guess values are the roots of the Chebyshev polynomial of the same degree, each time taking into account that p i = 1.To this goal, our algorithm entails that if a x i is calculated larger than the x i+1 , then x i receives the value of x i+1 (double root), whereas if a x i is calculated smaller than the x i−1 , then x i−1 receives the value of x i (double root).The former case results in p i < 1 and the latter in p i−1 < 1, although we do not have to compute them.Observe that for the system of Equations ( 8)- (10), there is no solution when, finally, x 1 < 0 (even symmetry) or for the system of Equations ( 14)-( 17), when x 2 < 0 (odd symmetry).In the case that the final solution has double roots, it should be checked whether it satisfies other variations of root combinations.Attention must be, also, drawn to the fact that when not all roots are single, the solution is not unique and thus one can not possibly deduce which of them is the optimum, unless their plots are a posteriori examined.
As a first example to comprehend the prior scheme, consider a sixth-order filter with Applying the SbOM, we find that x 1 = 0.268, x 2 = 0.715, x 3 = 0.896, while for the corresponding Chebyshev polynomial the roots are x 1 = 0.2588, x 2 = 0.7071, x 3 = 0.9659.Let us, now, more elaborately outline the initial steps of our algorithm.Focusing on the even symmetry case, we employ Equation (9) to acquire x 1 (while x 2 and x 3 come from the Chebyshev polynomial), namely (x 1 + 0.1) 2 (0.7071 + 0.1) 2 (0.9659 + 0.1) 2 (1 where we selected the positive value.Then, this x 1 is plugged into Equation (10) to compute ξ 1 (again, x 2 and x 3 are obtained from the Chebyshev polynomial) as where we chose the solution between x 1 and x 2 .Next, this ξ 1 is substituted into Equation ( 8) to extract x 2 (now, x 1 is acquired from the previous step and x 3 from the Chebyshev polynomial), i.e., (0.4919 + 0.1) 2 −0.1115 2 0.4919 2 −(x 2 +0.1) 2 0.4919 2 −(0.9659 which is, again, the positive value.Similarly, this x 2 is plugged into Equation (10) to obtain ξ 2 (with x 1 obtained from the previous step and x 3 from the Chebyshev polynomial) and so on until the convergence of the solution.The plot of the new polynomial along with its typical Chebyshev counterpart are illustrated in Figure 1a.Additionally, Figure 1b presents three polynomial implementations, extracted through the SbOM as this is distorted by diverse combinations of the roots, where none of the maximum local peaks is larger than 1.In order to further elaborate with our analysis, Figure 2a shows the magnitude of the S-parameters for our sixth-order filter, derived though the Chebyshev, Butterworth, and SbOM filtering functions of Figure 1a.Observe that the SbOM results lie between the Chebyshev and Butterworth ones, while the |S 11 | obtained via the SbOM is at least 10 dB lower than its Chebyshev counterpart.On the other hand, Figure 2b compares the magnitude of the S-parameters for the sixth-order filter, retrieved through the Chebyshev and SbOM filtering functions that are distorted by the ±( Secondly, we examine a sixth-order filter with d 1 = 0.9, The SbOM leads to x 1 = x 2 = x 3 = 0.26.This is not an acceptable outcome because the third maximum peak (abscissa closest to 1) is larger than 1, for the distortion described by the ±(x 1 + 0.9), ±x 2 − 0.2, ±x 3 − 0.7 combination of the roots (i.e., K x 2 − (x 1 + 0.9 ).Note that this combination does not satisfy the rule of Proposition 2. In contrast, a feasible solution is the set of x 1 = 0.128, x 2 = x 3 = 0.268.Therefore, the third maximum peak is equal to 1, for the distortions denoted by the ±(x 1 + 0.9), ±x 2 − 0.2, ±x 3 − 0.7 and ±x 1 − 0.9, ±x 2 − 0.2, ±(x 3 + 0.7) combinations of the roots (namely, K x 2 − (x 1 + 0.9 as well as K (x + 0.9 ) and the maxima of the other peaks are smaller than 1.Nevertheless, this solution is not unique, since the x 1 = x 2 = x 3 = 0.135 set is, also, acceptable.Thus, the third maximum peak is equal to 1, for the distortion described by the ±(x 1 + 0.9), ±x 2 − 0.2, ±x 3 − 0.7 combination of the roots and the maxima of the other peaks are smaller than 1 as well as smaller than those of the first solution.This reveals that the p i of the latter solution are smaller than those of the former one, which is, finally, preferred.Our selection is, also, verified by means of Figure 3, which shows that the former solution (i.e., the x 1 = 0.128, x 2 = x 3 = 0.268 roots) provides the steepest plot.

The Compressed Chebyshev Polynomial Method (CoCPM)
The key concept of the novel method is that the Chebyshev polynomial (or any other polynomial) is compressed in the amplitude and frequency range to: (a) satisfy the criterion which requires that its maximum local peaks (in absolute value) do not exceed 1 and (b) cross the (1, 1) point.In particular, via the corresponding plot, we find which combination of the Chebyshev polynomial root variations leads to its maximum absolute value, i.e., the maximum of all the maximum local peaks.Then, we compress the polynomial both in amplitude (so that the above maximum absolute value becomes 1) and in frequency (so that the polynomial crosses the (1, 1) point).We designate the final outcome as the compressed Chebyshev polynomial (CoCP) and describe its derivation in detail below.
Let us suppose a nth-order Chebyshev polynomial, C n (x), expressed, in the case of even symmetry (with a similar analysis for odd symmetry), as Next, we introduce factor k and coefficient a < 1, such that the maximum deviation of |kC n (x/a)| from the horizontal axis (due to the variation of its roots) does not exceed 1 and C n (x/a) crosses the (1, 1) point.This, in turn, means that kC n (1/a) = 1 and k = 1/C n (1/a), respectively.Actually, a is used to retain the same bandpass zone.If the maximum deviation of kC n (x/a) from the horizontal axis occurs at the ith peak, we define the following distorted polynomial so that its maximum value, i.e., 1, is at its ith peak, for x = ξ i .Moreover, the unknown a and ξ i are determined by the system of and In this framework, the desired CoCP is denoted as , it remains a typical Chebyshev polynomial in [−a, a], where its roots are the ax i and its maximum value is Notice that, in (a, 1], Equation ( 23) increases until it crosses the (1, 1) point.The solution of the proposed polynomial is always very close to the optimal one, accomplishing, also, a much simpler filtering process and filter implementation, which can be readily utilized for the improvement of existing techniques, as shown in Section 3. Due to the theorem of [43], the roots of Equation ( 23) are calculated from the distorted Hence, the optimal solution is obtained from p i ≤ 1, which yields a definitely steeper plot.
As an example, lets us consider the case of a sixth-order filter, with d 1 = d 2 = d 3 = 0.1.The resulting Chebyshev polynomial obtains its maximum value for the x 1 − d 1 , x 2 − d 2 , x 3 + d 3 root variation, i.e., at the third peak.Then, by employing Equations ( 21) and ( 22), we obtain ξ 3 = 0.8485 and a = 0.9266, so that the roots of the CoCP are ax 1 = 0.2398, ax 2 = 0.6552, ax 3 = 0.895, where x i are the roots of the Chebyshev polynomial.It is must be stated that the optimal solution is acquired with all p i = 1, which is better than the CoCPM, where only p 3 = 1.Lastly, and through Equation ( 24), the maximum value of the new polynomial, in [−0.9266, 0.9266], is found to be 0.1847.The prior outcomes are shown in Figure 4, while Figure 5

Extension to Rational Filtering Functions
A noteworthy asset of the SbOM is its straightforward application to rational polynomial functions.For this aim, let us presume the prototype rational polynomial function . ., y n + d n combinations of the poles, i.e., the zeros approach the ith valley.In fact, due to the close proximity of zeros to poles and the steep transit zone, the aforementioned scheme is efficient for small changes, namely the rational function is more sensitive to zeros and poles variations than the polynomial.However, the principal notion for deriving the F n (x), such that its peaks and valleys are not over predetermined values, is still valid.Observe that the optimization theorem of [43] is applicable to the |∂F n (x)/∂x| x=1 derivative, owing to the steep plot of the transit zone.
Next, we design a sixth-order elliptic filter, whose polynomial function has the the x 1 = 0.4017, x 2 = 0.868, x 3 = 0.9898 zeros and the y 1 = 1.0691, y 2 = 1.2192, y 3 = 2.6344 poles ([44] see pp. 33, 34 and replace m with m 1/2 in (2.5.19) and (2.5.20)).Assume, for instance, that all d i s are equal to 0.1, i.e., the roots and poles change by ±0.1.This leads to a sixth-order prototype rational filtering function with x 1 = 0.379, x 2 = 0.827, x 3 = 0.888 zeros and y 1 = 1.353, y 2 = 1.448, y 3 = 2.704 poles, as shown in Figure 5a.Its solution is obtained by solving Equations ( 8)-( 10) sequentially for x 1 , ξ 2 , x 2 , ξ 3 , x 3 , ξ 4 , y 1 , ξ 5 , y 2 .Then, y 3 is calculated by imposing the K = (1 3 )] = 100 coefficient via the lim x→∞ F n (x) = K constraint.Moreover, ξ 4 and ξ 5 are the abscissas of the lowest valleys between y 1 , y 2 and y 2 , y 3 poles, respectively.These interesting findings are presented in Figures 6 and 7, which, also, include the case when the prior filter is distorted by certain combinations of the roots.As detected, the featured extension to rational polynomial functions provides a promising treatment for several complicated distortions, thus guaranteeing a reliable realization process for effective and robust filters.

Realistic Applications and Numerical Verification
For its comprehensive validation, the novel methodology is applied to the design of several waveguide and microstrip line filters, which constitute indispensable elements of contemporary 6G communication systems.The selected filters operate at fairly different frequencies, covering the popular 6G spectrum between 2 GHz and 65 GHz, in order to prove that our technique does not depend on frequency.In this context, the required polynomials are derived through the proposed SbOM and CoCPM of Section 2, while all numerical simulations are performed via the COMSOL Multiphysics ® simulation software [38].

Design of a Sixth-Order Waveguide Filter
To derive the sixth-order polynomial for the specific waveguide filter, we, firstly, estimate the variations of the roots from the relevant literature.Then, we compare our SbOM and CoCPM solutions with the chained function filter scheme [45] and the filter presented in [37].It is should be emphasized that this comparison is conducted according to a unified way, proposed herein, in order to draw trustworthy conclusions, namely:

•
The maximum amplitude of each prototype polynomial, as distorted owing to the variation of its roots, should be 1.
The required estimation on the variation of the roots is accomplished via the plots in [45], where a distorted sixth-order Chebyshev polynomial is depicted.Note that the position of the roots refers to the abscissas of the minima of the S 11 -parameter magnitude, for which we derive that d 1 ∈ [−0.2588, 0.02], d 2 ∈ [−0.07, 0.07], and d 3 ∈ [−0.14, 0.1].Moreover, to consider the worst case scenario, we assume that d 1 , d 2 , d 3 can move independently to each other, and hence the diverse elevations of the peaks may be larger than those in [45].Subsequently, by means of the SbOM and Equations ( 8)- (10), the roots of the optimal polynomial F 6 (ω) are calculated as x 1 = 0.3993, x 2 = 0.6875, x 3 = 0.8832.Moreover, through the CoCPM and Equations ( 21) and ( 22), the respective roots are x 1 = 0.2363, x 2 = 0.6456, x 3 = 0.8819.For the latter solution, the maximum amplitude occurs for the third peak that should be equal to 1, i.e., we solve Equations ( 21) and ( 22) for the x 1 − 0.2588, x 2 − 0.07, x 3 + 0.1 root variation.Indeed, through the proper plots, it can be promptly detected that all the other peaks have lower maximum values due to the fluctuations of the roots.Furthermore, the maximum of the prior CoCP in [−a, a], with a = 0.913, is found to be 0.1479 from Equation (24).So, the waveguide filter is equivalent to a Chebyshev filter with ε = 0.1479ε = 0.01486 in [−a, a].This means that we lower the plot of the S 11 -parameter magnitude by an additional 16.5 dB (from the initial −20 dB), whereas the (1, 1) point continues to correspond to the level of −20 dB.
For our comparisons, we employ the chained function of [45], which combines the first-, second-, and third-order Chebyshev polynomials.Concerning the specific sixth-order waveguide filter, this results in 8x 6 − 10x 4 + 3x 2 = 8x 2 (x 2 − 0.70712)(x 2 − 0.8662) that crosses the (1, 1) point.Since the maximum value of its second peak is larger than 1, the chained function is compressed in frequency and amplitude, by means of Equations ( 21) and (22), thus leading to the 0, 0.6972, and 0.8539 roots.Next, we consider the polynomial presented in [37], whose roots are 0.3036, 0.7683, 0.9659 and its largest peak, due to root variation, is the second one.However, the specific function does not cross the (1, 1) point; therefore, we compress it in frequency and amplitude, via Equations ( 21) and ( 22), to acquire the 0.2731, 0.6911, and 0.8688 roots.Summarizing all the optimal F 6 (ω), Figure 8a illustrates their initial plots, which practically coincide for |x| > 1. Conversely, Figure 8b shows their behavior owing to the variation of their roots, as discussed in the previous paragraph.Note that none of the F 6 (ω) maximum amplitudes surpass 1.In addition, Figure 9a,b present the magnitude of the S 12 -parameter (ε = 0.1) and the S 11 -parameter (ε = 0.1), respectively, both computed in terms of the corresponding polynomials of Figure 8a.It can be detected that the outcomes of the proposed CoCPM and the compressed polynomial [37] are practically identical for |x| > 1. Lastly, the impact of the various F 6 (ω) root variations of Figure 8b on the |S 11 | (ε = 0.1) is examined in Figure 9c.Therefore, we prove that none of the plots exceeds the level of −20 dB, unlike typical formulations that lack to offer adequate results.In this manner, one can comprehend the merits of our method which, together with the compression concept for the significant improvement of existing approaches, can provide flexible 6G filters.

Design of a Fourth-Order Chebyshev Microstrip Line Filter at 2.4 GHz
The verification of our theoretical framework, herein, focuses on the analysis of a microstrip line filter from the Chebyshev family, which comprises a frequent selection in modern communication systems.Therefore, after studying the filter in its initial form, we deform the microstrip line to emulate potential fabrication tolerances and then compute the new filtering function, via our methodology, which leads to the new microstrip filter.Finally, the resulting structure is deformed, again, in order to systematically investigate its overall behavior and substantiate the effectiveness of the featured design procedure.
Based on these aspects, we study the fourth-order Chebyshev microstrip line filter of Figure 10a with a central frequency of 2.4 GHz, described in [46].In its initial form, the dimensions of the filter are L 1 = 11.4 mm, L 2 = L 3 = 11.5 mm, w = 0.4 mm, s 1 = 0.7 mm, s 2 = 1.6 mm, and s 3 = 1.8 mm.Also, tapers are used at the input-output ends to reduce the fringing fields, due to the transition from the 50 Ω input microstrip lines to the thinner coupled lines of the structure.Our basic design restriction is |S 11 | < −9.636 dB.As previously described, this implies that the filtering function is given by 8(x 2 − 0.383 2 )(x 2 − 0.924 2 ), with ε = 0.35 in the |S 12 | formula.Therefore, and after the appropriate numerical simulations, Figure 10b presents the magnitude of the S-parameters, which, as promptly observed, are not symmetrical, with the first |S 11 | peak computed at −7.98 dB.Subsequently, the prior filter is deformed both in terms of the width and length of its microstrip lines.For example, let us consider Figure 11, which presents the top right-hand part of the filter shown in Figure 10a.In particular, width deformations are created by adding ±ε i = 30 µm (for i = 1, 3, 4, 5, 7, 8), while length deformations by adding ±ε i = 30 µm (for i = 2, 6, 9).Moreover, regarding the inner microstrip lines, there exist six possible ways of deformation, i.e., six ε i s: one for each of their short edges and two for each of their long ones.On the other hand, concerning the outer microstrip lines, there are three possible ways of deformation, i.e., three ε i s: one for each of their short edges and one for each of their long ones.Hence, the total number of deformations is 30.Furthermore, to avoid excessive simulations and motivated by the sensitivity function for parasitic elements of [47], we presume a Taylor series dependence of any filter characteristic, ch, on the ε i deformation, as where α i abd γ i are unknown weights.Note that, throughout our analysis, we have selected ch to stand for the local maximum of the S 11 -parameter magnitude, although equivalent outcomes could be derived if, for example, ch represented a root of the filtering function.
Considering these notions, we perform two numerical simulations for every ε i (i.e., for the +ε i and −ε i value), setting the rest of the ε i s to zero, and record each local maximum.Actually, we keep the ε i for which the local peak is maximized, namely the +ε i , −ε i , or ε i = 0 case.Thus, a total of 60 simulations are conducted and three sequences of ε i s are obtained to maximize each local peak.Finally, we simulate the deformed filter for these sequences and derive the local peaks of the |S 11 | to be maximized, as illustrated in Figure 12.
The last stage of the validation process involves the design of the new fourth-order Chebyshev microstrip line filter of Figure 15, consistent with the filtering function derived above via the CoCPM and shown in Figure 13a.Basically, this filtering function leads to a filter with a bandwidth of 9.387% and a ripple of 0.105 dB.To calculate the symmetric odd-mode, Z oi , and even-mode, Z ei , impedances (for i = 1, 2, 3) of the device, we employ the well-known formulae from ( [12] i.e., (8.108) and (8.121), respectively).Thus, we obtain Z o1 = 38.4426Ω, Z e1 = 74.7126Ω, Z o2 = 44.6487Ω, Z e2 = 56.8367Ω, and Z o3 = 45.6354Ω, Z e3 = 55.2984Ω.Moreover, the dimensions of the structure in Figure 15 are calculated by means of the relevant well-known expressions, found in [12,48].So, after some mathematical manipulations, one acquires L 1 = 11.3593mm, L 2 = 11.0938mm, L 3 = 11.0615mm, w 1 = 0.801 mm, s 1 = 0.415 mm, w 2 = 1.074 , s 2 = 1.494 mm, w 3 = 1.089 mm, and s 3 = 1.819 mm.The promising performance of the numerically simulated filter is verified through the magnitude of its S-parameters in Figure 16a, which indicate the significantly reduced |S 11 | local peaks and the considerable bandwidth.Similar deductions can be drawn from the correspondingly enhanced phase of the filter's S-parameters in Figure 16b.Finally, the new filter is deformed, again, via the same sequence of microstrip line deformation, already applied to the initial structure.Figure 17 presents the magnitude and the phase of the S-parameters.Hence, the maximum value (and the corresponding phase) of the |S 11 | first peak (Figure 17a,d) is −8.05 dB, second peak (Figure 17b,e) is −10.4 dB, and third peak (Figure 17c,f) is −10.92 dB.As observed, apart from the slight elevation of the first peak, attributed, chiefly, to design inaccuracies, the rest of the peaks do not exceed the limitation of −9.696 dB.Furthermore, all the peaks are definitely below the levels of their initial filter counterparts.These observations verify the robustness of the featured filter (even if deformed) and the efficiency of our technique to treat construction imperfections in demanding design scenarios.Having successfully determined the efficiency of our technique at the low part of the 6G spectrum, we extend our investigation to much higher frequencies.In this manner, we will be able to prove that its overall performance is completely independent of the filter's operating frequency and certify its reliable use at any frequency range.As already described in the previous examples, after studying the filter in its initial form, the microstrip line is deformed to imitate potential fabrication tolerances.Next, we compute the new filtering function, which leads to the new (and far more resilient to fabrication imperfections) microstrip filter.The resulting device is then deformed, again, in order to thoroughly assess its overall behavior and affirm the advantages of the design process.
Hence, let us examine the fourth-order microstrip parallel coupled-line bandpass filter, depicted in its initial form in Figure 18a, with a central frequency at 62.5 GHz, as presented in [49].Actually, such devices, operating in the vicinity of 60 GHz, have triggered a significant research interest, due to anticipated impressive usage of these frequencies in a multitude of future scenarios [1,2,19].The dimensions of the filter are L 1 = 0.995934 mm, L 2 = 0.924559 mm, L 3 = 0.984250 mm, w 1 = 0.643128 mm, w 2 = 0.604012 mm, w 3 = 0.642112 mm, s 1 = 0.079248 mm, s 2 = 0.062992 mm, and s 3 = 0.078740 mm.Its numerical simulation leads to the outcomes of Figure 18b, which illustrates the magnitude of the S-parameters.Observe that all |S 11 | peaks are below the desired threshold of −10 dB.To apply our method, the prior filter is, next, deformed by adding ±ε i = 5 µm (for i = 1, 2, . . ., 30) to each of its microstrip lines.As a consequence, a total of 60 simulations are performed and three sequences of ε i 's are acquired in order to maximize each local peak.Then, the deformed filter is computationally simulated for these sequences to extract the local peaks of the |S 11 | to be maximized, which are shown in Figure 19.Explicitly, the maximum value of the |S 11 | first peak (Figure 19a) is −10.65 dB, second peak (Figure 19b) is −6.84 dB, and third peak (Figure 19c) is −11.02 dB. a system of equations with the prerequisite that the local maxima of the deformed initial polynomial do not exceed a certain limit, without the need of considering their phase.On the other hand, the CoCPM develops an equivalent polynomial through a compression process in the amplitude and frequency regime along with a corresponding threshold concerning the local maxima.For the optimum profile of the featured filtering functions, a versatile criterion has been formulated along with a useful procedure for objectively comparing them to their various counterparts from existing approaches.Both techniques are completely universal in their applicability (with regard to frequency spectrum or device geometry) and can handle any filtering function.Furthermore, they are straightforward to implement, efficient in the compensation of fabrication imperfections, fast to deliver accurate results, and frugal in computational resources.Actually, these aspects constitute the key differences of the novel theoretical analysis with those of the related literature and pave the way for future research, such as the consideration of substrate losses.The proposed methodology has been validated via the design of real-world waveguide and microstrip filters, covering the broad range between 2 GHz and 65 GHz, which indicate its potential to serve as a trustworthy tool for the construction of high-end 6G filters.
combination of the roots.It can be, clearly, detected that the |S 11 | acquired via the distorted SbOM never exceeds the level of −20 dB (complying fully with our design requirements), unlike the distorted Chebyshev outcome.The latter deductions are, also, verified by Figure2cwhich indicates a definite lag in the phase of the S 12 -parameters, acquired from the distorted Chebyshev and SbOM filtering functions.

Figure 1 .Figure 2 .Figure 3 .
Figure 1.Diverse sixth-order polynomial realizations for the design of the corresponding filter (absolute values are depicted).(a) The Chebyshev polynomial and the polynomial derived via the proposed SbOM, with d 1 = d 2 = d 3 = 0.1.(b) Three different polynomials, derived via the proposed SbOM (with d 1 = d 2 = d 3 = 0.1) as this is distorted by various combinations of the roots.The first, second, and third maximum peak, of each case, respectively, do not exceed 1.

Figure 4 .Figure 5 .
Figure 4. Diverse sixth-order polynomial realizations for the design of the corresponding filter (absolute values are depicted).(a) The Chebyshev polynomial and the polynomial derived via the proposed CoCPM, with d 1 = d 2 = d 3 = 0.1.(b) Two different polynomials, derived via the proposed CoCPM, with d 1 = d 2 = d 3 = 0.1 and the ±x 1 − d 1 , ±x 2 − d 2 , ±(x 3 + d 3 ) combination of the roots, respectively.The third maximum peak does not exceed 1.

Figure 6 .Figure 7 .
Figure 6.Design of a sixth-order rational filter.(a) The sixth-order (sixth-order polynomial for the numerator and the denominator) prototype rational filtering function (absolute values depicted) derived by considering that the zeros and poles vary by ±0.1 (all d i s are equal to 0.1).(b) Magnitude of the S-parameters (ε = 0.1), extracted through the sixth-order elliptic polynomial and the rational filtering function of Figure 6a.

Figure 8 .Figure 9 .
Figure8.Optimal sixth-order polynomial, F 6 (ω), realizations for the design of the corresponding waveguide filter, derived via the proposed SbOM and CoCPM, the compressed polynomial[37], and the compressed chained function[45] (absolute values are depicted).(a) The initial form of the different F 6 (ω).All plots are almost identical for |x| > 1.(b) Behavior of the diverse F 6 (ω) owing to the variation of their roots.The maximum amplitudes do not exceed 1.

Figure 10 .
Figure 10.The initial form of the fourth-order Chebyshev microstrip line filter terminated by means of a taper at both sides for the necessary input-output signal routing.(a) Geometry and (b) magnitude of the S-parameters.The maximum value of the |S 11 | first peak is −7.98 dB.

Figure 11 .
Figure 11.Geometry of the deformed, due to fabrication imperfections, coupled microstrip lines of the fourth-order Chebyshev filter of Figure 10a.The sketch illustrates the top right-hand part of the filter and the dashed lines depict the deformations at the microstrip lines.

Figure 13 .Figure 14 .
Figure 13.(a) The new fourth-order filtering function derived via the SbOM and the CoCPM and (b) magnitude of the respective S-parameters.

3 Figure 15 .Figure 16 .Figure 17 .
Figure 15.Geometry of the new fourth-order Chebyshev microstrip line filter designed via the CoCPM-extracted filtering function of Figure 13a.

3. 3 .
Design of a Fourth-Order Microstrip Parallel Coupled-Line Bandpass Filter at 62.5 GHz

FrequencyFigure 18 .
Figure 18.The initial form of the fourth-order microstrip parallel coupled-line bandpass filter terminated by means of a taper at both sides for the appropriate input-output signal routing.(a) Geometry and (b) magnitude of the S-parameters.The maximum value of the |S 11 | first peak is −10.59 dB.

Figure 21 .
Figure 21.The S-parametrs of the new fourth-order microstrip parallel coupled-line bandpass filter.(a) Magnitude and (b) phase.The bandwidth of the filter is 41.331% and its ripple 0.1651 dB.
Magnitude of the S-parameters of the deformed fourth-order Chebyshev microstrip line filter.The deformation occurs according to a sequence of ε i s, such that the |S 11 | first, second, or third peak is maximized.Explicitly, the maximum value of the |S 11 | (a) first peak is −4.27 dB, (b) second peak is −6.16 dB, and (c) third peak is −7.88 dB.