求大神帮忙

Ts=1;
N=15;
eye_num=6;
a=1;
N_data=1000;
dt=Ts/N;
t=-3*Ts:dt:3*Ts;
% 产生双极性数字信号
d=sign(randn(1,N_data));
dd=sigexpand(d,N); %
% 基带系统冲击响应（升余弦）
ht=sinc(t/Ts).*(cos(a*pi*t/Ts))./(1-4*a^2*t.^2/Ts^2+eps);
st=conv(dd,ht);
tt=-3*Ts:dt:(N_data+3)*N*dt-dt;
subplot(211)
plot(tt,st);
axis([0 20 -1.2 1.2]);
xlabel('t/Ts');
ylabel('基带信号');
subplot(212)
% 画眼图
ss=zeros(1,eye_num*N);
ttt=0:dt:eye_num*N*dt-dt;
for k=3:50
ss=st(k*N+1:(k+eye_num)*N);
drawnow;
plot(ttt,ss);
hold on;
end;
xlabel('t/Ts');
ylabel('基带信号眼图');
% 将输入的序列扩成间隔为N-1个0的序列
运行错误：??? At运tempt to execute SCRIPT sigexpand as a function.
Error in ==> <a href="error:D:\MATLABwu\work\Untitled3.m,10,1">Untitled3 at 10</a>
dd=sigexpand(d,N); %
运行的时候老是出错，百度的方法也不能解决!哪位大神帮我解决下呀！