频域采样定理演示matlab

以下是频域采样定理的matlab演示代码：

clc
close all;
clear all;

% 定义采样频率和观测时间
Fs = 2000; % 采样频率为2kHz
Tp = 50/1000; % 观测时间为50ms

% 定义信号
n = 0:1/Fs:Tp-1/Fs;
x = 5*sin(2*pi*100*n) + 3*sin(2*pi*200*n) + 2*sin(2*pi*300*n);

% 绘制原始信号的时域图像
figure;
subplot(2,1,1);
plot(n,x);
xlabel('Time (s)');
ylabel('Amplitude');
title('Original Signal in Time Domain');

% 进行频域采样
M = 64; % 采样点数为64
Xk = fft(x,M); % 对信号进行M点FFT

% 绘制采样信号的频域图像
k = 0:M-1;
fk = k/Tp;
subplot(2,1,2);
stem(fk,abs(Xk));
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('Sampled Signal in Frequency Domain');

该代码演示了如何使用matlab进行频域采样定理的演示。首先定义了采样频率和观测时间，然后定义了一个信号。接着绘制了原始信号的时域图像。然后进行频域采样，采样点数为64，对信号进行M点FFT。最后绘制了采样信号的频域图像。可以看到，采样信号的频域图像中只有三个频率分量，与原始信号相同，说明频域采样定理成立。