clc,clear,close all;
dt=0.01;
Ts=1;
t=-5*Ts:dt:5*Ts;
N=Ts/dt;
eye_num=5;%眼图数量
N_data=1000;%码元数
% 产生双极性数字信号，+1，-1 序列
d=sign(randn(1,N_data)); %randn 随机生成数字,sign 函数进行一个正负判断，正数为+1，负数为-1
dd=sigexpand(d,N); %将输入的序列扩成间隔为 N-1 个 0 的序列
% 部分响应系统冲击响应
ht=sinc(t./Ts);
ht1=sinc((t-Ts)./Ts);
gt=ht+ht1;%另一种产生 gt 的方法，即 ht 与 ht 延迟一个周期 Ts 后的叠加
st=conv(dd,gt);
tt=-5*Ts:dt:(N_data+5)*N*dt-dt; %设置采样时间
subplot(211);
plot(tt,st);
hold on;
axis([0 20 -3 3]);
xlabel('t');
ylabel('部分响应基带信号');
title('第一类部分响应基带信号');
grid on;
% 画眼图
subplot(212);
ss=zeros(1,eye_num*N);
ttt=0:dt:eye_num*N*dt-dt;
for k=5:50
ss=st(k*N+1:(k+eye_num)*N);
drawnow;
plot(ttt,ss);
hold on;
end;
xlabel('t');
ylabel('部分响应信号眼图');
title('第一类部分响应信号眼图');
grid on;
gtext('图 3-1 第一类部分响应系统的基带信号和其眼图波形')
dt=0.01;
Ts=1;
t=-5*Ts:dt:5*Ts;
N=Ts/dt;
eye_num=5;%眼图数量
N_data=1000;%码元数
% 产生双极性数字信号，+1，-1 序列
d=sign(randn(1,N_data)); %randn 随机生成数字,sign 函数进行一个正负判断，正数为+1，负数为-1
dd=sigexpand(d,N); %将输入的序列扩成间隔为 N-1 个 0 的序列
% 部分响应系统冲击响应
ht=sinc(t./Ts);
ht1=sinc((t-Ts)./Ts);
gt=ht+ht1;%另一种产生 gt 的方法，即 ht 与 ht 延迟一个周期 Ts 后的叠加
st=conv(dd,gt);
tt=-5*Ts:dt:(N_data+5)*N*dt-dt; %设置采样时间
% 比较观察 t 对 h(t)和 g(t)影响
plot(t,ht)
hold on
plot(t,gt,'r--')
hold on
grid;
legend('h(t)','g(t)');
gtext('图 3-2 t 对 g(t)影响')
