Rename the variable in the cheb2mag function. It was overriding a Scilab function...
[scilab.git] / scilab / modules / signal_processing / macros / cheb2mag.sci
1 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
2 // Copyright (C) INRIA - F.D
3 // 
4 // This file must be used under the terms of the CeCILL.
5 // This source file is licensed as described in the file COPYING, which
6 // you should have received as part of this distribution.  The terms
7 // are also available at
8 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
9 function [h2]=cheb2mag(n,omegar,A,samplefreq)
10 //<h2>=cheb2mag(n,omegar,A,samplefreq)
11 //Square magnitude response of a type 1 Chebyshev filter
12 //omegar = stopband edge
13 //samplefreq = vector of frequencies where the square magnitude
14 //h2 is desired.
15 //  n       :Filter order
16 //  omegar  :Cut-off frequency
17 //  A       :Attenuation in stop band
18 //  samplefreq  :Vector of frequency where cheb2mag is evaluated
19 //  h2      :Chebyshev II filter values at sample points
20 //
21 //!
22
23
24    [n1,n2]=size(samplefreq);
25    un=ones(n1,n2);
26    Tn=chepol(n,'x');             //n-th Chebyshev polynomial
27    frd=freq(Tn,1,omegar*un./samplefreq);   //frd=Tn(omegar/samplefreq)
28    h2=un./(un+(A*A-1)*un./real(frd.*frd))
29 endfunction
30