ms2r4.m

《数字信号处理实验(MATLAB版)》这

f=200;Um=1;nt=2;            %输入信号频率、振幅和显示周期个数
Fs=4000;N=Fs/f;                 %采样频率,求采样点数N
T=1/f;                          %T为信号的周期
dt=T/N;                         %采样时间间隔
n=0:nt*N-1;                     %建立离散信号的时间序列
tn=n*dt;                        %确定时间序列样点在时间轴上的位置
x=Um*square(2*f*pi*tn);         %产生时域信号
subplot(2,1,1),plot(n,x);   %连续信号
axis([0,nt*N,1.1*min(x),1.1*max(x)]);
subplot(2,1,2),stem(n,x);   %离散序列
axis([0,nt*N,1.1*min(x),1.1*max(x)]);