Primero convierta las frecuencias angulares digitales de wp y ws en frecuencias angulares analógicas. El método de conversión es (aquí las frecuencias angulares simuladas están representadas por Wp y Ws): WP = 2/t * tan (WP/2), WS = 2/t * tan (WS/2). Esta fórmula está disponible en línea y se convierte en la fórmula previa a la deformación de la transformación bilineal.
Lo que quieres es una forma ovalada, yo nunca la he usado. Obtuve el tipo Chebyshev, pero el principio debería ser similar, es decir, verificar su ecuación cuadrada modular y usar los WP, WS, DP y DS anteriores para resolver los parámetros de la función cuadrada modular. Los parámetros de Chebyshev son n y e.
Descubra la forma de la ecuación de expresión elíptica H(s), introduzca los parámetros ahora mismo y escriba las H(s) específicas en ella usando la expresión polar.
Después de obtener H(s), conviértalo en H(z) mediante transformación bilineal. El método de transformación es s=2/T*(1-z menos 1 vez)/(1+z menos 1 vez). Esta fórmula se puede encontrar en línea.
Después de la sustitución, se obtiene H(z).
Luego escribe H(z) en forma de ecuación diferencial, que es la expresión de y(n), y podrás dibujar una gráfica.
Cada palabra está escrita a mano, espero que el cartel la adopte.
Aunque no te he entregado un trabajo ni resultados de simulación, te he dado un proceso de operación completo, solo síguelo paso a paso.