频谱与倒谱分析在Matlab中的实现

阅读量：

一、频谱

原理：傅里叶变换

采样频率需高于原始信号频率的两倍

方法：采用fft和ifft函数（其中fft点数决定了频谱所覆盖的频率范围）

经过fft处理后，频率分布区间呈现为前半部分[0,fs/2]，后半部分[-fs/2,0]
为使频率区间恢复至常规形式，应使用fftshift函数进行调整
若需提取正频率部分，可在执行fft后截取前半段数据，或在使用fftshift后截取后半段数据

在这里插入图片描述

MATLAB程序：

复制代码

    %时域波形
    w=rectwin(61)%矩形窗函数
    subplot(411)
    plot(w,'k')
    xlabel('样本数');ylabel('幅度');title('矩形波时域波形')
    %频域波形
    w1=fft(w,1024)%对窗函数做1024点的fft
    w2=w1/(max(w1))%幅度归一化
    w3=20*log10(abs(w2))%出来的是双边频谱（数学谱）
    w4=fftshif

全部评论 (0)

还没有任何评论哟~