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.