PriskriboElliptic 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.
Dato	la 29-a de januaro 2009
Fonto	Propra verko
Aŭtoro	Inductiveload
Permeso (Reuzo de la dosiero)	Public domainPublic domainfalsefalse Mi, la posedanto de la aŭtorrajto por ĉi tiu verko, liberigas ĝin tutmonde kiel publikan havaĵon. En iuj landoj tio povas esti laŭleĝe neebla; en tiu okazo: Mi rajtigas ĉiun ajn uzi ĉi tiun verkon por ajna celo, sen iuj ajn kondiĉoj, krom se tiuj kondiĉoj estas postulataj per leĝo.

Mathematica Code

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}]