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

通讯编程文档

开发平台:

Matlab

  1. % MATLAB script for Illustrated Problem 7.1. 
  2. echo on
  3. T=1;
  4. delta_T=T/200;           % sampling interval
  5. alpha=0.5;           % rolloff factor
  6. fc=40/T;           % carrier frequency
  7. A_m=1;      % amplitude
  8. t=-5*T+delta_T:delta_T:5*T;       % time axis
  9. N=length(t);
  10. for i=1:N,
  11.   if (abs(t(i))~=T/(2*alpha)),  
  12.     g_T(i) = sinc(t(i)/T)*(cos(pi*alpha*t(i)/T)/(1-4*alpha^2*t(i)^2/T^2));
  13.   else
  14.     g_T(i) = 0;   % The value of g_T is 0 at t=T/(2*alpha)
  15.   end;   % and at t=-T/(2*alpha).
  16.   echo off ;
  17. end;
  18. echo on;
  19. G_T=abs(fft(g_T));                % spectrum of g_T
  20. u_m=A_m*g_T.*cos(2*pi*fc*t);   % the modulated signal
  21. U_m=abs(fft(u_m));   % spectrum of the modulated signal
  22. % actual frequency scale
  23. f=-0.5/delta_T:1/(delta_T*(N-1)):0.5/delta_T;
  24. % Plotting commands follow.
  25. figure(1);
  26. plot(f,fftshift(G_T));  
  27. axis([-1/T 1/T 0 max(G_T)]);
  28. figure(2);
  29. plot(f,fftshift(U_m));