diff options
Diffstat (limited to 'macros/freqs.sci')
| -rw-r--r-- | macros/freqs.sci | 29 |
1 files changed, 29 insertions, 0 deletions
diff --git a/macros/freqs.sci b/macros/freqs.sci new file mode 100644 index 0000000..88773fe --- /dev/null +++ b/macros/freqs.sci @@ -0,0 +1,29 @@ +function h = freqs (b, a, w) +//Compute the s-plane frequency response of the IIR filter. +//Calling Sequence +//h = freqs (b, a, w) +//Parameters +//b: vector containing the coefficients of the numerator of the filter. +//a: vector containing the coefficients of the denominator of the filter. +//w: vector containing frequencies +//Description +//This is an Octave function. +//It computes the s-plane frequency response of the IIR filter B(s)/A(s) as H = polyval(B,j*W)./polyval(A,j*W). +//If called with no output argument, a plot of magnitude and phase are displayed. +//Examples +//B = [1 2]; +//A = [1 1]; +//w = linspace(0,4,128); +//freqs(B,A,w); + +funcprot(0); +rhs = argn(2) +if(rhs<3 | rhs>3) +error("Wrong number of input arguments.") +end + +select (rhs) + case 3 then + y = callOctave("freqs",b, a, w) + end +endfunction |