ip_06_07.m
上传用户:loeagle
上传日期:2013-03-02
资源大小:1236k
文件大小:1k
源码类别:

通讯编程文档

开发平台:

Matlab

  1. % MATLAB script for Illustrative Problem 6.7.
  2. echo on
  3. N=31;
  4. T=1;
  5. alpha=1/4;
  6. n=-(N-1)/2:(N-1)/2; % the indices for g_T
  7. % The expression for g_T is obtained next.
  8. for i=1:length(n),
  9.   g_T(i)=0;
  10.   for m=-(N-1)/2:(N-1)/2,
  11.     g_T(i)=g_T(i)+sqrt(xrc(4*m/(N*T),alpha,T))*exp(j*2*pi*m*n(i)/N);
  12.     echo off ;
  13.   end;
  14. end;
  15. echo on ;
  16. g_T=real(g_T) ; % The imaginary part is due to the finite machine precision
  17. % Derive g_T(n-(N-1)/2).
  18. n2=0:N-1;
  19. % Get the frequency response characteristics.
  20. [G_T,W]=freqz(g_T,1);
  21. % normalized magnitude response
  22. magG_T_in_dB=20*log10(abs(G_T)/max(abs(G_T)));  
  23. % impulse response of the cascade of the transmitter and the receiver filters
  24. g_R=g_T;
  25. imp_resp_of_cascade=conv(g_R,g_T);
  26. % Plotting commands follow.