DescrizioneElliptic Filter Response (4th Order, Detail).svg	English: The response of a 4-th order elliptic filter (detail of the transistion band). Also shown is the maximum extent of the passband ripples, where ε is the ripple factor and Ln is the discrimination factor.
Data	29 gennaio 2009
Fonte	Opera propria
Autore	Inductiveload
Licenza (Riusare questo file)	Public ___domainPublic ___domainfalsefalse Io, detentore del copyright su quest'opera, la rilascio nel pubblico dominio. Questa norma si applica in tutto il mondo. In alcuni paesi questo potrebbe non essere legalmente possibile. In tal caso: Garantisco a chiunque il diritto di utilizzare quest'opera per qualsiasi scopo, senza alcuna condizione, a meno che tali condizioni siano richieste dalla legge.

xp2[zeta_] :=
  Module[{g, num, den},
   g = Sqrt[4*zeta^2 + (4*zeta^2*(zeta^2 - 1))^(2/3)];
   num = 2*zeta^2*Sqrt[g];
   den = Sqrt[8*zeta^2*(zeta^2 + 1) + 12*g*zeta^2 - g^3] - Sqrt[g^3];
   num/den
   ];
xz2[zeta_] := zeta^2/xp2[zeta];

t[zeta_] := Sqrt[1 - 1/zeta^2];

r1[zeta_, x_] := x;
r2[zeta_, x_] := ((t[zeta] + 1)*x^2 - 1)/((t[zeta] - 1)*x^2 + 1);
r3[zeta_, x_] := 
  x*((1 - xp2[zeta])*(x^2 - xz2[zeta]))/((1 - xz2[zeta])*(x^2 - 
        xp2[zeta]));
r4[zeta_, x_] :=
  Module[{num, den},
   num = (1 + t[zeta]) (1 + Sqrt[t[zeta]])^2*x^4 - 
     2 (1 + t[zeta]) (1 + Sqrt[t[zeta]])*x^2 + 1;
   den = (1 + t[zeta]) (1 - Sqrt[t[zeta]])^2*x^4 - 
     2 (1 + t[zeta]) (1 - Sqrt[t[zeta]])*x^2 + 1;
   num/den
   ];

ellgain[xi_, w_, w0_, ep_] := 1/Sqrt[1 + ep^2*r4[xi, w/w0]^2];

Plot[
 ellgain[1.05, w, 1, 0.5],
 {w, 0, 3}]