matlab 绘制原始信号的谐波

转载于 2019-07-21 18:05:00 发布 · 3k 阅读

3 ·

本内容遵循CC 4.0 BY-SA版权协议

原文链接：http://www.cnblogs.com/darkchii/p/11222049.html

本文通过Matlab代码演示了锯齿信号的生成及其快速傅立叶变换(FFT)的频谱分析过程。采用分段采样策略，展示了信号在时间和频率域上的变化，为理解信号处理中的频谱特性提供了直观的视觉效果。

　　这里以锯齿信号为例：

clear;clc;
figure

time = 0:1/20:1000;
wave = sawtooth(time);
subplot(3, 1, 1);
plot(time, wave);

samples = length(wave);
sampling = 100;
ampl = zeros(sampling, floor(samples/sampling));
frame = 1;
for i=0:sampling:samples
    first = max(1,i+1);
    last = min(samples,i+sampling);
    ampl(1:sampling, frame) = abs(fft(wave(first:last)));
    frame = frame + 1;
end

subplot(3, 1, 2);
for i=1:sampling
    plot3(0:1:samples/sampling, ones(1, floor(samples/sampling) + 1) * i, ampl(i,:));
    hold on
end

subplot(3, 1, 3);
for i=1:samples/sampling
    plot3(0:1:sampling - 1, ones(1, sampling) * i, ampl(:,i));
    hold on
    %drawnow
end

　　运行结果：