1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
|
function h = cl2bp (m, w1, w2, up, lo, gridsize)
//Constrained L2 bandpass FIR filter design.
//Calling Sequence
//h = cl2bp (m, w1, w2, up, lo, gridsize)
//h = cl2bp (m, w1, w2, up, lo)
//Parameters
//m: degree of cosine polynomial, i.e. the number of output coefficients will be m*2+1
//w1 and w2: bandpass filter cutoffs in the range 0 <= w1 < w2 <= pi, where pi is the Nyquist frequency
//up: vector of 3 upper bounds for [stopband1, passband, stopband2]
//lo: vector of 3 lower bounds for [stopband1, passband, stopband2]
//gridsize: search grid size; larger values may improve accuracy, but greatly increase calculation time.
//Description
//This is an Octave function.
//Constrained L2 bandpass FIR filter design. Compared to remez, it offers implicit specification of transition bands, a higher likelihood of convergence, and an error criterion combining features of both L2 and Chebyshev approaches.
//Examples
//h = cl2bp(5, 0.3*pi, 0.6*pi, [0.02, 1.02, 0.02], [-0.02, 0.98, -0.02], 2^11)
//h =
//
// 0.038311
// 0.082289
// -0.086163
// -0.226006
// 0.047851
// 0.307434
// 0.047851
// -0.226006
// -0.086163
// 0.082289
// 0.038311
funcprot(0);
lhs = argn(1)
rhs = argn(2)
if (rhs < 5 | rhs > 6)
error("Wrong number of input arguments.")
end
select(rhs)
case 5 then
res = callOctave("cl2bp", m, w1, w2, up, lo)
case 6 then
res = callOctave("cl2bp", m, w1, w2, up, lo, gridsize)
end
endfunction
|
