summaryrefslogtreecommitdiff
path: root/modules/signal_processing/macros/ell1mag.sci
blob: 15cbfc7c7c9daa71d9f2e73f674f7a42c1e3c107 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) INRIA
//
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution.  The terms
// are also available at
// http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.txt

function [v]=ell1mag(eps,m1,z)
    //Function used for squared magnitude of an elliptic filter
    //Usually m1=eps*eps/(a*a-1);
    //  eps     :Passband ripple=1/(1+eps**2)
    //  m1      :Stopband ripple=1/(1+(eps**2)/m1)
    //  z       :Sample vector of values in the complex plane
    //  v       :Elliptic filter values at sample points
    //
    //!
    s=%sn(z,m1);un=ones(z);
    v=real(un./(un+eps*eps*s.*s))


endfunction