The design of linear phase filters has been intensively discussed in literature. For the two most common error criteria (l2l2 and l∞l∞), optimal solution algorithms exist. The least squares norm filter can be found by solving an overdetermined system of equations, whereas the Chebishev norm filter is easily found by using either the Remez algorithm or linear programming. For many typical applications, linear phase filters are good enough; however, when arbitrary magnitude and phase constraints are required, a more complicated approach must be taken since such design results in a complex approximation problem. By replacing CC in the linear phase algorithm with a complex Fourier kernel matrix, and the real desired frequency vector DD with a complex one, one can use the same algorithm from (Reference) to design complex lplp filters.
