f_getrate.asv

DSP程序 Matlab是一套用于科学工程计算的可视化高性能语言与软件环境。它集数值分析、矩阵运算、信号处理和图形显示于一体

function [x,y,b,fs,user,i_type,i_typeold] = f_getrate (i_type,i_typeold,fs,L,M,c,m,f_type,user,x,y,hc_type);

%F_GETRATE: Compute input and output for GUI module G_MULTIRATE
%
% Usage: [x,y,b,fs,user] = f_getrate (i_type,fs,L,M,c,m,f_type,user,x,y);
%
% Inputs: 
%         i_type  = input type
%
%                   1 = white noise  
%                   2 = damped cosine
%                   3 = amplitude modulated (AM)
%                   4 = frequency modulated (FM)
%                   5 = record sound
%                   6 = user defined
%
%         i_typeold = previous i_type
%         fs        = sampling frequency 
%         L         = interpolation factor (1 to 100)
%         M         = decimation factor (1 to 100)
%         c         = damping factor
%         m         = FIR filter order 
%         f_type    = FIR filter type
%
%                     0 = windowed (rectangular)
%                     1 = windowed (Hanning)
%                     2 = windowed (Hamming)
%                     3 = windowed (Blackman)
%                     4 = frequency-sampled
%                     5 = least-squares
%                     6 = equiripple 
%
%         user      = string containing name of MAT-file that
%                     contains input x and samping frequency fs
%         x         = 1 by N vector containing current input samples
%         y         = 1 by P vector containing current output samples
%         hc_type   = array of handles to input type radio buttons
% Outputs: 
%          x       = 1 by N vector containing input samples
%          y       = 1 by P vector containing output samples where
%                    P = floor(L*N/M);
%          b       = 1 by m+1 coefficient vector for FIR filter
%          fs      = sampling frequency
%          user    = string containing name of MAT-file that
%                    contains input x and samping frequency fs
%          i_type  = input type
%          i_typ

% Initialize

N = length(x);
T = 1/fs;
x_save = x;
fs_save = fs;

% Construct input
 
switch (i_type)
   
case 1,                                          % white noise
   
    clear x
    N = 1000;
    x = f_randu (N,1,-1,1);
   
case 2,                                          % damped cosine
   
   clear x
   N = 500;
   F_0 = fs/10;
   k = 0 : N-1;
   x = c.^k .* cos(2*pi*F_0*k*T);
      
case 3,                                          % amplitude modulated (AM)
   
   clear x
   N = 1500;
   F_0 = fs/180;
   F_1 = 32*F_0;
   k = 0 : N-1;
   x = sin(2*pi*F_0*k*T) .* cos(2*pi*F_1*k*T);
   
case 4,                                          % frequency mopdulated (FM)
   
   clear x
   N = 1500;
   F_0 = fs/80;
   F_1 = 4*F_0;
   k = 0 : N-1;
   f = F_0 + (F_1-F_0)*sin(4*pi*k*T);
   x = cos(2*pi*f.*k*T);
   
case 5,                                          % record x
   
    
   fs = 8000;
   [x,cancel] = f_getsound (x,1.0,fs);
   if cancel
       x = x_save;
       fs = fs_save;
   end
   
case 6,                                          % user defined
   
   caption = 'Select MAT-file containing x,fs';
   [user1,path] = f_getmatfile (caption,user); 
   if user1 ~= 0
       user = user1;
       old_path = pwd;
       cd (path);
       load (user,'x','fs')
       cd (old_path)
   else
       x = x_save;
   end

end

% Check x and compute y using rate conversion

x = f_tocol(x);
[y,b] = f_rateconv (x,fs,L,M,m,f_type);