已知连续时间信号为ft=e-0.5t次方u(t),对该信号进行取样,取样间隔为0.1s,用matlab绘制取样信号的三个周期的频谱图
时间: 2024-03-26 18:35:12 浏览: 70
首先,按照上一个问题中的方法生成采样信号,可以使用Matlab中的`linspace`函数生成时间轴上的采样点,然后通过计算得到采样信号:
```
t = linspace(0, 15, 151); % 生成0到15,共151个采样点的时间轴
ft = exp(-0.5*t) .* (t>=0); % 计算采样信号
```
然后,可以使用Matlab中的`fft`函数对采样信号进行傅里叶变换得到频谱,再使用`abs`函数计算频谱的幅值,最后使用`plot`函数绘制频谱图:
```
fs = 1 / (t(2) - t(1)); % 计算采样频率
f = linspace(0, fs/2, length(t)/2+1); % 生成频率轴
F = fft(ft) / length(t); % 对采样信号做傅里叶变换
F = F(1:length(t)/2+1); % 取一半的频率范围
Famp = abs(F); % 计算频谱的幅值
```
接下来,可以使用`repmat`函数复制采样信号,生成三个周期的采样信号,并对每个周期的采样信号都做上述的频谱计算和绘制:
```
N = 3; % 生成3个周期的采样信号
ft_N = repmat(ft, 1, N); % 复制N个周期的采样信号
for i = 1:N
ft_i = ft_N((i-1)*length(t)+1:i*length(t)); % 选择第i个周期的采样信号
Famp_i = abs(fft(ft_i) / length(t)); % 对第i个周期的采样信号做傅里叶变换并计算幅值
Famp_i = Famp_i(1:length(t)/2+1); % 取一半的频率范围
plot(f, Famp_i); % 绘制第i个周期的频谱图
hold on;
end
hold off;
xlabel('Frequency (Hz)');
ylabel('Amplitude');
```
运行上述代码,即可得到采样信号的三个周期的频谱图。
阅读全文