Firstly, the digital angular frequency of wp and ws is converted into analog angular frequency. The conversion method is (here the analog angular frequency is represented by Wp and Ws): WP = 2/t * tan (WP/2), WS = 2/t * tan (WS/2). This formula is available on the internet and becomes the prewarp formula of bilinear transformation.

What you want is an oval shape, which I have never used. I got Chebyshev type, but the principle should be similar, that is, look up its modulus square equation and use WP, WS, DP and DS above to solve the parameters of the modulus square function. Chebyshev's parameters are n and e.

Find out the form of elliptic H(s) expression equation, bring in the parameters just now, and write the specific H(s) in pole expression.

After H(s) is obtained, it is transformed into H(z) by bilinear transformation. The transformation method is s=2/T*( 1-z minus 1 times) /( 1+z minus 1 times), and this formula can be found on the Internet.

After replacement, H(z) is obtained.

Then write H(z) in the form of difference equation, that is, the expression of y(n), and you can draw a picture.

Every word is handwritten, I hope the landlord will adopt it.

Although I didn't give you a paper or simulation results, I have given the operation process completely, step by step.