diff options
| -rw-r--r-- | macros/sinetone.sci | 94 |
1 files changed, 64 insertions, 30 deletions
diff --git a/macros/sinetone.sci b/macros/sinetone.sci index 407fa45..f174869 100644 --- a/macros/sinetone.sci +++ b/macros/sinetone.sci @@ -1,32 +1,66 @@ -function y= sinetone(x, varargin) +function y = sinetone(freq, rate, sec, ampl) //Return a sinetone of the input -//Calling Sequence -//y= sinetone(FREQ) -//y= sinetone(FREQ, RATE) -//y= sinetone(FREQ, RATE, SEC) -//y= sinetone(FREQ, RATE, SEC, AMPL) -//Parameters -//FREQ: frequency of sinetone -//RATE: Sampling rate -//SEC: Length in seconds -//AMPL: Amplitude -//Description -//Return a sinetone of frequency FREQ with a length of SEC seconds atsampling rate RATE and with amplitude AMPL.The arguments FREQ and AMPL may be vectors of common size.The defaults are RATE = 8000, SEC = 1, and AMPL = 64. -funcprot(0); -rhs= argn(2); -if(rhs<1 | rhs>4) -error("Wrong number of input arguments") -end - -select(rhs) - case 1 then - y= callOctave("sinetone", x); - case 2 then - y= callOctave("sinetone", x , varargin(1)); - case 3 then - y= callOctave("sinetone", x , varargin(1), varargin(2)); - case 4 then - y= callOctave("sinetone", x , varargin(1), varargin(2), varargin(3)); - -end +//Calling Sequence: +//y= sinetone(freq) +//y= sinetone(freq, rate) +//y= sinetone(freq, rate, sec) +//y= sinetone(freq, rate, sec, ampl) +//Parameters: +//freq: frequency of sinetone +//rate: sampling rate +//sec: length in seconds +//ampl: amplitude +//Description: +//Return a sinetone of frequency 'freq' with a length of 'sec' seconds at sampling rate 'rate' and with amplitude 'ampl'. +//The arguments freq and ampl may be vectors of common size. The defaults are rate = 8000, sec = 1, and ampl = 64. +//Examples: +//sinetone(5, 2, 1, 8) +//ans = [4.89E-15; -9.79E-15] + + funcprot(0); + rhs= argn(2); + + select rhs + case 1 then + rate = 8000; + sec = 1; + ampl = 64; + case 2 then + sec = 1; + ampl = 64; + case 3 then + ampl = 64; + case 4 then + break; + else + error("sinetone: wrong number of input arguments"); + end + + if ( isvector(freq) & isvector(ampl) & ~isequal(size(freq), size(ampl)) ) then + error("sinetone: freq and ampl must be vectors of equal size"); + end + + if ( ~(isscalar(rate) && isscalar(sec)) ) then + error("sinetone: rate and sec must be scalars"); + end + + n = length (freq); + ns = round (rate * sec); + y = zeros (ns, n); + + for k = 1:n + y(:, k) = ampl(k) * sin (2 * %pi * (1:ns) / rate * freq(k))'; + end + endfunction + +//input validation: +//assert_checkerror("sinetone()", "sinetone: wrong number of input arguments"); +//assert_checkerror("sinetone([6, 9, 4, 2], 2, 3, [6, 2, 0])", "sinetone: freq and ampl must be vectors of equal size"); +//assert_checkerror("sinetone(1, [1, 2])", "sinetone: rate and sec must be scalars"); +//assert_checkerror("sinetone(1, 2, [4, 3])", "sinetone: rate and sec must be scalars"); + +//tests: +//assert_checkequal(size(sinetone(18e6, 150e6, 19550/150e6, 1)), [19550, 1]); +//assert_checkequal(sinetone(5), sinetone(5, 8000, 1, 64)); +//assert_checkequal(size(sinetone([1, 2, 3], 4000, 1, [8, 9, 6])), [4000, 3]); |