gaussf.m
上传用户:m_sun_001
上传日期:2014-07-30
资源大小:1115k
文件大小:1k
源码类别:

matlab例程

开发平台:

Matlab

  1. % Program 3-13
  2. % gaussf.m
  3. %
  4. % Function to form Gaussaian filter
  5. %
  6. % programmed by H.Harada
  7. function [xh] = gaussf(B,irfn,ipoint,sr,ncc)
  8. %**************************************************************** 
  9. % irfn : Number of symbols to use filtering
  10. % ipoint : Number of samples in one symbol
  11. % sr : symbol rate
  12. %  B :  filter coeficiense
  13. % ncc; : 1 -- transmitting filter  0 -- receiving filter
  14. %**************************************************************** 
  15. point = ipoint;
  16. tr = sr ;  
  17. n = ipoint .* irfn;
  18. mid = ( n ./ 2 ) + 1;
  19. fo=B/sqrt(2*log(2));
  20. for i = 1 : n 
  21.   icon = i - mid;
  22.   ym = icon;
  23.  xt=1/2*(erf(-sqrt(2/log(2))*pi*B*(ym/ipoint-1/2)/tr)+erf(sqrt(2/log(2))*pi*B*(ym/ipoint+1/2)/tr));
  24.       
  25.     
  26.     
  27.   if ncc == 0                         % in the case of receiver
  28.      xh( i ) = xt ;
  29.   elseif ncc == 1                       % in the case of transmitter
  30.      xh( i ) = xt;
  31.   else
  32.     error('ncc error');
  33.   end   
  34.   end  
  35. %******************** end of file ***************************