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

通讯编程文档

开发平台:

Matlab

  1. % MATLAB script for Illustrative Problem 6.11.
  2. echo on
  3. T=1;
  4. for n=-2:2,
  5.   for k=-2:2,
  6.     temp=0;
  7.     for i=-2:2, temp=temp+(1/(1+(n-i)^2))*(1/(1+(k-i)^2)); end;
  8.     X(k+3,n+3)=temp;
  9.     echo off ;
  10.   end;
  11. end;
  12. echo on;
  13. N0=0.01;         % assuming that N0=0.01
  14. Ry=X+(N0/2)*eye(5);
  15. Riy=[1/5 1/2 1 1/2 1/5].';        
  16. c_opt=inv(Ry)*Riy;    % optimal tap coefficients
  17. % find the equalized pulse...
  18. t=-3:1/2:3;
  19. x=1./(1+(2*t/T).^2);    % sampled pulse
  20. equalized_pulse=conv(x,c_opt);    
  21. % Decimate the pulse to get the samples at the symbol rate.
  22. decimated_equalized_pulse=equalized_pulse(1:2:length(equalized_pulse));
  23. % Plotting command follows.