% 部分响应信号眼图示意，pres.m
clear all;
close all;
Ts=1;
N_sample=16;
eye_num=11;
N_data=1000;
dt=Ts/N_sample;
t=-5*ts:dt:5*Ts;
%产生双极性数字信号
d=sign(randn(1,N_data));
dd=sigexpand(d,N_sample);
%部分响应系统冲激响应
ht=sinc((t+eps)/Ts)./(1-(t+eps)./Ts);
ht(6*N_sample+1)=1;
st=conv(dd,ht);
tt=-5*Ts:dt:(N_data+5)*N_sample*dt-dt;
figure(1);
subplot(211);
plot(tt,st);
axis([0 20 -33]);xlabel('t/Ts');ylabel('部分响应基带信号');
subplot(212);
%画眼图
ss=zeros(1,eye_num*N_sample);
ttt=0:dt:eye_num*N_sample*dt-dt;
for k=5:50
    ss=st(k*N_sample+1:(k+eye_num)*N_sample);
    drawnow;
    plot(ttt,ss);hold on;
end
%plot(ttt,ss);
xlabel('t/Ts');ylabel('部分响应信号眼图');