眼图matlab仿真

N=2^10; %采样点数
L = 2^6; %每码元的采样点数
M = N/L; 
Alpha = 0.25; %升余弦滚降系数
Na = 3; %在示波器上显示3个元素
Rb = 2; %码元速率
Ts = 1/Rb; %比特间隔
dt = Ts/L; 
df = 1/（N*dt）;
T = N*dt; 
BW = N*df/2; %系统带宽
fs = N/T; 
f = -BW+df/2:df:BW; %频域采样点
t = -T/2+dt/2:dt:T/2; %时域采样点

h = sin（pi*t/Ts）.*cos（pi*t*Alpha/Ts）./（pi*t/Ts.*（1-4*t.^2*Alpha^2/Ts^2））;
Hf = t2f（hfs）;  % 傅里叶变换
EP = zeros（size（f））+eps;

for ii = 1:40
    b = 1-2*sign（randn（1M））; %双极性码
    %产生的PAM信号
    s =