二进制基带脉冲成形及眼图的绘制

clear all;

fs=107;%采样频率
Ts=1;%基带信号码元周期
N_sample=fs;%码元采样点数
N=1000;%码元数目

% ------ 时域参数 ------ %
dt=1/N_sample;
t=0:dt:（N*N_sample-1）*dt;

% ------ 绘制升余弦滚降滤波器的冲激响应函数 ------ %
a=1;
t0=-3:dt:3;
ht=sin（pi*t0/Ts）./（pi*t0./Ts）.*cos（a*pi*t0/Ts）./（1-4*a.^2*t0.^2/Ts.^2）;

ht（（length（t0）+1）/2）=1;

figure（1）;
plot（t0ht‘LineWidth‘1.5）;
axis（[-3 3 -0.2 1.2]）;
title（‘升余弦滤波器的冲击响应函数h（t）‘）;
xlabel（‘时间（s）‘）;
ylabel（‘电压值（V）‘）;
grid on;

% ------ 初始化二进制数字信号码元 ------ %
bc0=[-1zeros（1106）];
bc1=[1zeros（1106）];

% ------ 生成二进制随机基带信号 ------ %
rand_code=randi（[01]1N）;

% ------ 初始化二进制基带脉冲信号 ------ %
mt=zeros（1N*N_sample-1）;

% ------ 计算二进制基带脉冲信号 ------ %
for i=1:N
    if rand_code（i）==1
        mt（（i-1）*N_sample+1:i*N_sample）=bc1;
    else
         mt（（i-1）*N_sa