观察离散时间和连续时间信号频谱关系的matlab程序

clear; T = 1; fs = 1000; t = -T :1/(3*fs): T; Ts = 0.2; %脉冲宽度 x = pulse(t, Ts); f = -500 : 0.1 : 500; w = 2*pi*f; x_jw = trapz(t, x'.*exp(-1i*t'*w)); figure; subplot(2, 1, 1); plot(t, x); title('连续时间脉冲函数'); subplot(2, 1, 2); plot(f, abs(x_jw)); xlabel('Hz'); title('连续时间信号幅频特性'); %进行采样得到离散时间信号 n = -T : 1/fs : T; x_n = pulse(n, Ts); w = 0 : 1/(100*pi) : 4*pi; x_ejw = x_n * exp(-1i*fs*n'*w); figure; subplot(2, 1, 1); stem(n*fs, x_n); title('离散时间脉冲函数'); subplot(2, 1, 2); plot(w/pi, fftshift(abs(x_ejw))); xlabel('\pi'); title('离散时间信号幅频特性');