用MATLAB实现连续时间信号的采样及重构 以f(t)=Sa(t)=(sin t)/t为例，分别以过采样、等采样、欠采样三种情况, 绘出原信号、采样信号、重构信号的时域波形图。

以下是一些可能的MATLAB代码，用于实现f(t) = Sa(t) = (sin t)/t 连续时间信号的采样和重构：

% 定义采样率和采样时长
fs = 100; % 采样率
T = 2*pi; % 采样时长
ts = 1/fs; % 采样间隔

% 定义原始信号
t = -T/2:0.001:T/2;
x = sinc(t/T);

% 过采样，采样率为 400 Hz
fs1 = 400; % 过采样率
ts1 = 1/fs1; % 过采样间隔
n1 = ceil(T/ts1); % 过采样点数
t1 = 0:ts1:(n1-1)*ts1; % 过采样时间轴
x1 = sinc(t1/T); % 过采样信号

% 等采样，采样率为 100 Hz
t2 = -T/2:ts:T/2; % 等采样时间轴
x2 = sinc(t2/T); % 等采样信号

% 欠采样，采样率为 50 Hz
fs3 = 50; % 欠采样率
ts3 = 1/fs3; % 欠采样间隔
n3 = ceil(T/ts3); % 欠采样点数
t3 = 0:ts3:(n3-1)*ts3; % 欠采样时间轴
x3 = sinc(t3/T); % 欠采样信号

% 重构信号
x_r1 = sinc((t1 - n1*ts1/2)/T)*fs1*ts1; % 过采样信号重构
x_r2 = sinc((t2 - T/2)/T)*fs*ts; % 等采样信号重构
x_r3 = sinc((t3 - n3*ts3/2)/T)*fs3*ts3; % 欠采样信号重构

% 绘制图形
figure;
subplot(3,3,1); plot(t, x); title('Original Signal');
subplot(3,3,4); stem(t1, x1); title('Oversampled Signal, fs1 = 400 Hz');
subplot(3,3,7); plot(t, x_r1); title('Reconstructed Signal, fs1 = 400 Hz');
subplot(3,3,2); plot(t, x); title('Original Signal');
subplot(3,3,5); stem(t2, x2); title('Uniformly Sampled Signal, fs = 100 Hz');
subplot(3,3,8); plot(t, x_r2); title('Reconstructed Signal, fs = 100 Hz');
subplot(3,3,3); plot(t, x); title('Original Signal');
subplot(3,3,6); stem(t3, x3); title('Undersampled Signal, fs3 = 50 Hz');
subplot(3,3,9); plot(t, x_r3); title('Reconstructed Signal, fs3 = 50 Hz');

运行上述代码后，将会生成一个包含 9 个子图的图形窗口，其中每行对应一种采样情况，包括原始信号、采样信号以及重构信号的时域波形图。