周期ws=2.4*pi/Ts;%理想低通截止频率n=-100:100;%定义序列的长度是201nTs=n*Ts%抽样点f=sinc(nTs/pi);%抽样信号Dt=0.005;t=-20:Dt:20;fa=f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs'*ones(1,length(t))));%信号重建t1=-20:0.5:20;f1=sinc(t1/pi);subplot(212);plot(t,fa)xlabel('t');ylabel('fa(t)');title('Sa(t)函数');grid;?图7函数(4)单位阶跃信号t=-1:0.01:4;%定义时间样本向量t0=0;%指定信号发生突变的时刻ut=stepfun(t,t0);%产生单位阶跃信号plot(t,ut)%绘制波形axis([-1,4,-0.5,1.5])%设定坐标轴范围title('单位阶跃信号')%设定标题图8单位阶跃信号(5)门函数的波形t=-4:0.01:4;%定义时间样本向量t1=-2;%指定信号发生突变的时刻u1=stepfun(t,t1);%产生左移位的阶跃信号e(t+2)t2=2;%指定信号发生突变的时刻u2=stepfun(t,t2);%产生右移位的阶跃信号e(t-2)g=u1-u2;%表示门函数plot(t,g)%绘制门函数的波形axis([-4,4,-0.5,1.5])%设定坐标轴范围-4<x<4,-0.5<y<1.5title('门函数')%设定标题图9门函数3.2常用离散时间信号的实现(1)单位抽样序列n1=-20:20;x1=[zeros(1,20),1,zeros(1,20)];subplot(3,2,1);stem(n1,x1);title('单位抽样序列');xlabel('n');ylabel('δ(n))');