* Bug 16290: %sn() set obsolete
[scilab.git] / scilab / modules / signal_processing / macros / ell1mag.sci
1 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
2 // Copyright (C) INRIA
3 //
4 // Copyright (C) 2012 - 2016 - Scilab Enterprises
5 //
6 // This file is hereby licensed under the terms of the GNU GPL v2.0,
7 // pursuant to article 5.3.4 of the CeCILL v.2.1.
8 // This file was originally licensed under the terms of the CeCILL v2.1,
9 // and continues to be available under such terms.
10 // For more information, see the COPYING file which you should have received
11 // along with this program.
12
13 function v = ell1mag(eps, m1, z)
14     //Function used for squared magnitude of an elliptic filter
15     //Usually m1=eps*eps/(a*a-1);
16     //  eps     :Passband ripple=1/(1+eps**2)
17     //  m1      :Stopband ripple=1/(1+(eps**2)/m1)
18     //  z       :Sample vector of values in the complex plane
19     //  v       :Elliptic filter values at sample points
20     //
21     //!
22     s = ellipj(z,m1);
23     un = ones(z);
24     v = real(un ./ (un+eps*eps*s.*s))
25
26
27 endfunction