生成一个基于MATLAB的利用正弦波验证频域采样定理的代码

好的，以下是一个基于MATLAB的利用正弦波验证频域采样定理的代码：

% 设置采样频率和信号频率
fs = 1000; % 采样频率
f = 100; % 信号频率

% 生成采样信号
t = 0:1/fs:1; % 采样时间
x = sin(2*pi*f*t); % 信号

% 进行离散傅里叶变换
X = fft(x);

% 计算频率轴
f_axis = linspace(-fs/2, fs/2, length(X));

% 绘制原始信号和频域信号
subplot(2, 1, 1);
plot(t, x);
xlabel('时间');
ylabel('幅值');
title('原始信号');

subplot(2, 1, 2);
plot(f_axis, abs(fftshift(X))/length(X));
xlabel('频率');
ylabel('幅值');
title('频域信号');

% 根据采样频率计算最大可还原频率
f_max = fs/2;

% 绘制还原信号
figure;
% 生成还原信号
t_recon = 0:1/f_max:1; % 还原时间
x_recon = sin(2*pi*f*t_recon); % 还原信号
plot(t_recon, x_recon, 'r-');
hold on;
stem(t, x, 'b-');
xlabel('时间');
ylabel('幅值');
title('还原信号与采样信号');
legend('还原信号', '采样信号');

这个代码会生成一个基于正弦波的信号，并进行采样和离散傅里叶变换。然后，它会绘制原始信号和频域信号。接下来，它会计算最大可还原频率，并生成还原信号。最后，它会绘制还原信号和采样信号，以验证频域采样定理。