File: cheb1mag.sci

package info (click to toggle)
scilab 2.6-4
  • links: PTS
  • area: non-free
  • in suites: woody
  • size: 54,632 kB
  • ctags: 40,267
  • sloc: ansic: 267,851; fortran: 166,549; sh: 10,005; makefile: 4,119; tcl: 1,070; cpp: 233; csh: 143; asm: 135; perl: 130; java: 39
file content (24 lines) | stat: -rw-r--r-- 747 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
 
function [h2]=cheb1mag(n,omegac,epsilon,sample)
//<h2>=cheb1mag(n,omegac,epsilon,sample)
//Square magnitude response of a type 1 Chebyshev filter
//omegac=passband edge
//epsilon such that 1/(1+epsilon**2)=passband ripple
//sample vector of frequencies where the square magnitude
//is desired.
//  n       :Filter order
//  omegac  :Cut-off frequency
//  epsilon :Ripple in pass band
//  sample  :Vector of frequency where cheb1mag is evaluated
//  h2      :Chebyshev I filter values at sample points
//
//!
//Author F.D.
// Copyright INRIA

   [n1,n2]=size(sample);
   un=ones(n1,n2);
   Tn=chepol(n,'x');  //n-th Chebyshev polynomial
   fr=freq(Tn,1,sample/omegac);   //fr=Tn(sample/omegac)
   h2=un./(un+epsilon*epsilon*fr.*fr)   //magnitude