clear;close?all;clc;
M=16;
L=512;
P=4;
ini_phase=0;
roll_off=0.35;
a=2*randint（1Lsqrt（M））-（sqrt（M）-1）;
b=2*randint（1Lsqrt（M））-（sqrt（M）-1）;
x=a+j*b;N=L*P;
y=zeros（1N）;
for?n=1:N
y（n）=0;
for?k=1:L
t=（n-1）/P-（k-1）;
y（n）=y（n）+x（k）*（sin（pi*t+eps）/（pi*t+eps））*...
（cos（roll_off*pi*t+eps）/（（1-（2*roll_off*t）^2）+eps））;
end
end
sff