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
|
function y = czt(x, varargin)
//Chirp Z Transform
//Calling Sequence
//czt (x)
//czt (x, m)
//czt (x, m, w)
//czt (x, m, w, a)
//Parameters
//x: Input scalar or vector
//m: Total Number of steps
//w: ratio between points in each step
//a: point in the complex plane
//Description
//This is an Octave function.
//Chirp z-transform. Compute the frequency response starting at a and stepping by w for m steps. a is a point in the complex plane, and w is the ratio between points in each step (i.e., radius increases exponentially, and angle increases linearly).
//Examples
// m = 32; ## number of points desired
// w = exp(-j*2*pi*(f2-f1)/((m-1)*Fs)); ## freq. step of f2-f1/m
// a = exp(j*2*pi*f1/Fs); ## starting at frequency f1
// y = czt(x, m, w, a);
funcprot(0);
lhs= argn(1);
rhs= argn(2);
if(rhs<1 | rhs > 4)
error("Wrong number of input arguments")
end
select (rhs)
case 1 then
y= callOctave("czt", x);
case 2 then
y= callOctave("czt", x, varargin(1));
case 3 then
y= callOctave("czt", x, varargin(1), varargin(2));
case 4 then
y= callOctave("czt", x, varargin(1), varargin(2), varargin(3));
end
endfunction
|
