ip_06_10.m

这是通信系统仿真的第六章的例题的仿真程序代码！挺好的

% MATLAB script for Illustrative Problem 6.10.
close all
clear all
echo on
T=1;
Fs=2/T;
Ts=1/Fs;
c_opt=[-2.2 4.9 -3 4.9 -2.2];
t=-5*T:T/2:5*T;			
x=1./(1+((2/T)*t).^2);		        % sampled pulse
t1=-5*T:T:5*T;	
xx=1./(1+((2/T)*t1).^2)
equalized_x=filter(c_opt,1,[x 0 0]);	% since there will be a delay of two samples at the output
% to take care of the delay
equalized_x=equalized_x(3:length(equalized_x));  
% Now, let us downsample the equalizer output.

for i=1:2:length(equalized_x),
  downsampled_equalizer_output((i+1)/2)=equalized_x(i);
  echo off;
end;
echo on ;
figure(1)
subplot(2,1,1)
stem([1:length(x)],x)
subplot(2,1,2)
stem([1:length(equalized_x)],equalized_x)

figure(2)
subplot(2,1,1)
stem([1:length(xx)],xx)
subplot(2,1,2)
stem([1:length(equalized_x)/2+1],downsampled_equalizer_output)

y=[0 0 1 0 0]
tao=T/2
m=[-2:2]

%for m=-2:2
m=-2
    n=[-2:2]
       xx(m+3,n+3)=1./(1+((2/T)*(m*T-n*tao)).^2);	
   %end
 %xx=reshape(xx,5,5)
%c_opt1=inv(xx)*y        %x*c_opt1=[0 0 1 0 0]
% Plotting commands follow.