matlab绘制取样信号频谱图

明日7时前随时采纳，需要代码，过程能详细一点就更好了，非常感谢

基于new bing的编写，有帮助记得采纳一下哦！：
【运行截图：】

【代码：】

% 定义信号 f(t)
syms t;
f = exp(-0.5*t)*heaviside(t);

% 取样间隔 Ts 和采样时间范围 t_range
Ts = 0.1;
t_range = 0:Ts:2;

% 进行采样
f_sampled = subs(f, t, t_range);

% 绘制取样信号的频谱图
N = length(f_sampled); % 采样点数
f_fft = fft(double(f_sampled)); % 进行 FFT 变换并转换为双精度数组
f_fft_shifted = fftshift(f_fft); % 进行 FFT 移频
freq_range = [-N/2:N/2-1]/(N*Ts); % 频率范围

figure;
plot(freq_range, abs(f_fft_shifted));
xlim([-3, 3]); % 限制频率范围
xlabel('Frequency (Hz)');
ylabel('Amplitude');
title('Sampled Signal Spectrum');

前言
创作不易,如果对你有帮助,给个采纳哦!

效果截图

代码及逻辑分析

根据采样定理，采样频率应该大于等于信号最高频率的两倍，即

f_s >= 2*f_max

其中f_max为信号的最高频率，对于连续时间信号f(t) = $$e^{-0.5t}u(t)$$ ，其最高频率为0.5Hz。因此，可以选择采样频率f_s = 1Hz，对原信号进行等间隔采样。

Matlab代码如下：

% 设置参数
fs = 1; % 采样频率
Ts = 1 / fs; % 采样周期
t = 0:Ts:(3*Ts); % 时间序列，绘制三个周期即可
f = exp(-0.5*t) .* (t >= 0); % 原信号
x = f .* cos(2*pi*5*t); % 等间隔采样信号

% 绘制时域图形
subplot(2, 1, 1);
stem(t, x, 'Linewidth', 2);
xlabel('Time (s)');
ylabel('Amplitude');
title('Sampled Signal');

% 绘制频谱图
subplot(2, 1, 2);
Fs = length(x) / t(end); % 实际采样频率
N = length(x); % 采样点数
f = (-N/2:N/2-1) * Fs / N; % 频率序列
X = fft(x) / N; % FFT变换
X = fftshift(X); % 将0频率移动到中心
plot(f, abs(X), 'Linewidth', 2);
xlim([-3, 3]);
xlabel('Frequency (Hz)');
ylabel('Magnitude');
title('Frequency Spectrum');

运行代码后，将会得到取样信号的时域图形和频谱图。其中时域图形为一个包含多个抽样点的离散信号，频谱图显示了取样信号的频谱，能够看出取样信号的频率内容。