ip_06_11.m
上传用户:loeagle
上传日期:2013-03-02
资源大小:1236k
文件大小:1k
- % MATLAB script for Illustrative Problem 6.11.
- echo on
- T=1;
- for n=-2:2,
- for k=-2:2,
- temp=0;
- for i=-2:2, temp=temp+(1/(1+(n-i)^2))*(1/(1+(k-i)^2)); end;
- X(k+3,n+3)=temp;
- echo off ;
- end;
- end;
- echo on;
- N0=0.01; % assuming that N0=0.01
- Ry=X+(N0/2)*eye(5);
- Riy=[1/5 1/2 1 1/2 1/5].';
- c_opt=inv(Ry)*Riy; % optimal tap coefficients
- % find the equalized pulse...
- t=-3:1/2:3;
- x=1./(1+(2*t/T).^2); % sampled pulse
- equalized_pulse=conv(x,c_opt);
- % Decimate the pulse to get the samples at the symbol rate.
- decimated_equalized_pulse=equalized_pulse(1:2:length(equalized_pulse));
- % Plotting command follows.