spectrogram（频谱图函数）matlab使用STFT参数选择以及滑动窗的问题

最新推荐文章于 2024-03-02 21:01:06 发布
最新推荐文章于 2024-03-02 21:01:06 发布
阅读量2.3w 收藏 199
点赞数 44
分类专栏： matlab
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_43158059/article/details/115032195
版权

spectrogram-翻译为频谱图
matlab原文档中的介绍如下：
该频谱图用于短时傅里叶画频谱图。
[S,F,T,P]=spectrogram(x,window,noverlap,nfft,fs)
[S,F,T,P]=spectrogram(x,window,noverlap,F,fs)

spectrogram(x)
s = spectrogram(x)
s = spectrogram(x, window)
s = spectrogram(x, window, noverlap)
s = spectrogram(x, window, noverlap, nfft)
s = spectrogram(x, window, noverlap, nfft, fs)
[S,F,T,P] = spectrogram(x, window, noverlap, nfft, fs)
[S,F,T,P] = spectrogram(x, window, noverlap, f, fs)
[S,F,T,P] = spectrogram(x, window, noverlap, f, fs)

参数：

  • x表示输入信号；

  • window表示窗函数，如果window的值是一个整数，那么被分段的x的每一段的长度都等于window，并采用默认的Hamming窗；如果window是一个向量，那么被分段后每一段的长度都等于length(window)，且输入的向量即为所要加的窗函数；

  • noverlap表示两段之间的重合点数，window-overlap=步长。overlap的值必须要小于窗长，如果没有指定overlap，默认是窗长的一半，即50%的overlap；

  • nfft表示fft的点数，fft的点数跟窗长可以是不同的，当没有指定该参数时，Matlab会取max(256,2^(ceil(log2(length(window)))))，即当窗长小于256时，fft的点数是256；当窗长大于256时，fft的点数取大于窗长的最小的2的整数次幂；这个数值应该是double或者single，不是向量。区分函数中nfft（离散点大小）和f（频率区间）的使用不同

  • f 表示显示的频谱范围，f是一个向量 (eg:f=1:30;%1-30hz的stft)，长度跟s的行数相同；

    • 当x是实信号且nfft为偶数时，s的行数为(nfft/2+1)

    • 当x是实信号且nfft为奇数时，s的行数为(nfft+1)/2

    • 当x是复信号时，s的行数为nfft

    • 当在输入的参数列表中指定f后，函数会在f指定的频率处计算频谱图，返回的f跟输入的f是相同的；

  • fs表示采样率，用来归一化显示使用；

%输出变量

  • S—输入信号x的短时傅里叶变换。它的每一列包含一个短期局部时间的频率成分估计，时间沿列增加，频率沿行增加。如果x是长度为Nx的复信号，则S为nfft行k列的复矩阵，其中k取决于window， 如果window为一个标量，则k = fix((Nx-noverlap)/(window-noverlap)).如果window为向量，则k = fix((Nx-noverlap)/(length(window)-noverlap))

  • F—在输入变量中使用F频率向量，函数会使用Goertzel方法计算在F指定的频率处计算频谱图。
    对于实信号x，如果nfft为偶数，则S的行数为(nfft/2+1)，如果nfft为奇数， 则行数为(nfft+1)/2，列数同上。
    指定的频率被四舍五入到与信号分辨率相关的最近的DFT容器(bin)中。而在其他的使用nfft语法中，短时傅里叶变换方法将被使用。对于返回值中的F向量，为四舍五入的频率，其长度 等于S的行数。

  • T—频谱图计算的时刻点，其长度等于上面定义的k，值为所分各段的中点。

  • P—能量谱密度PSD(Power Spectral
    Density)，对于实信号，P是各段PSD的单边周期估计；对于复信号，当指定F频率向量时，P为双边PSD。

matlab例子：

   
    x=[7.40000000000000e-05;6.50000000000000e-05;8.20000000000000e-05;0.000104000000000000;9.00000000000000e-05;8.70000000000000e-05;6.30000000000000e-05;4.80000000000000e-05;4.70000000000000e-05;2.60000000000000e-05;1.40000000000000e-05;2.20000000000000e-05;4.30000000000000e-05;2.00000000000000e-05;2.60000000000000e-05;2.40000000000000e-05;8.00000000000000e-06;4.00000000000000e-06;1.40000000000000e-05;1.70000000000000e-05;1.60000000000000e-05;2.30000000000000e-05;3.90000000000000e-05;4.10000000000000e-05;2.90000000000000e-05;2.10000000000000e-05;1.70000000000000e-05;1.20000000000000e-05;9.00000000000000e-06;1.80000000000000e-05;2.10000000000000e-05;2.10000000000000e-05;2.40000000000000e-05;1.40000000000000e-05;3.00000000000000e-05;5.30000000000000e-05;7.00000000000000e-05;4.40000000000000e-05;5.40000000000000e-05;9.10000000000000e-05;3.00000000000000e-06;-1.50000000000000e-05;-7.00000000000000e-06;-4.00000000000000e-05;-6.10000000000000e-05;-5.90000000000000e-05;-8.60000000000000e-05;-0.000108000000000000;-0.000102000000000000;-9.10000000000000e-05;-9.10000000000000e-05;-8.20000000000000e-05;-6.40000000000000e-05;-5.90000000000000e-05;-5.40000000000000e-05;-6.10000000000000e-05;-5.50000000000000e-05;-5.80000000000000e-05;-6.20000000000000e-05;-6.40000000000000e-05;-5.60000000000000e-05;-6.00000000000000e-05;-6.00000000000000e-05;-6.50000000000000e-05;-7.10000000000000e-05;-7.30000000000000e-05;-7.10000000000000e-05;-6.30000000000000e-05;-5.20000000000000e-05;-5.80000000000000e-05;-5.10000000000000e-05;-4.80000000000000e-05;-4.30000000000000e-05;-6.40000000000000e-05;-9.20000000000000e-05;-0.000106000000000000;-9.90000000000000e-05;-0.000112000000000000;-0.000117000000000000;-0.000112000000000000;-0.000106000000000000;-8.10000000000000e-05;-8.30000000000000e-05;-9.40000000000000e-05;-0.000101000000000000;-9.30000000000000e-05;-9.20000000000000e-05;-8.70000000000000e-05;-0.000114000000000000;-0.000117000000000000;-7.10000000000000e-05;-3.50000000000000e-05;-1.10000000000000e-05;1.10000000000000e-05;2.40000000000000e-05;1.50000000000000e-05;-2.60000000000000e-05;-0.000103000000000000;-0.000152000000000000;-7.60000000000000e-05;-4.90000000000000e-05;-8.10000000000000e-05;-4.70000000000000e-05;-1.10000000000000e-05;-7.00000000000000e-06;-1.20000000000000e-05;-1.00000000000000e-05;-2.30000000000000e-05;-3.90000000000000e-05;-2.60000000000000e-05;-1.20000000000000e-05;-2.50000000000000e-05;-3.00000000000000e-05;-1.10000000000000e-05;-3.00000000000000e-06;-2.00000000000000e-06;-8.00000000000000e-06;-3.50000000000000e-05;-3.70000000000000e-05;-2.50000000000000e-05;-2.00000000000000e-06;-1.10000000000000e-05;-1.40000000000000e-05;-2.70000000000000e-05;-3.10000000000000e-05;-1.10000000000000e-05;2.00000000000000e-06;1.00000000000000e-05;2.10000000000000e-05;4.20000000000000e-05;3.10000000000000e-05;3.10000000000000e-05;4.40000000000000e-05;3.10000000000000e-05;2.90000000000000e-05;3.20000000000000e-05;3.20000000000000e-05;4.50000000000000e-05;4.00000000000000e-05;4.40000000000000e-05;5.30000000000000e-05;3.60000000000000e-05;4.00000000000000e-05;4.20000000000000e-05;4.30000000000000e-05;4.10000000000000e-05;5.40000000000000e-05;9.70000000000000e-05;8.10000000000000e-05;6.10000000000000e-05;5.60000000000000e-05;7.00000000000000e-05;7.40000000000000e-05;6.10000000000000e-05;7.80000000000000e-05;9.10000000000000e-05;7.30000000000000e-05;6.20000000000000e-05;4.80000000000000e-05;4.50000000000000e-05;5.30000000000000e-05;5.20000000000000e-05;2.50000000000000e-05;1.90000000000000e-05;1.10000000000000e-05;2.00000000000000e-05;4.40000000000000e-05;3.00000000000000e-05;1.70000000000000e-05;1.80000000000000e-05;3.00000000000000e-05;3.00000000000000e-05;3.30000000000000e-05;4.50000000000000e-05;1.80000000000000e-05;8.00000000000000e-06;2.30000000000000e-05;3.00000000000000e-05;1.50000000000000e-05;2.80000000000000e-05;2.10000000000000e-05;4.10000000000000e-05;4.30000000000000e-05;4.60000000000000e-05;5.90000000000000e-05;3.70000000000000e-05;2.60000000000000e-05;2.60000000000000e-05;9.00000000000000e-06;1.30000000000000e-05;2.10000000000000e-05;2.90000000000000e-05;3.70000000000000e-05;5.30000000000000e-05;5.00000000000000e-05;3.90000000000000e-05;2.50000000000000e-05;2.60000000000000e-05;3.00000000000000e-05;4.40000000000000e-05;5.80000000000000e-05;5.20000000000000e-05;5.40000000000000e-05;8.00000000000000e-05;9.20000000000000e-05;9.20000000000000e-05;7.50000000000000e-05;4.20000000000000e-05;3.20000000000000e-05;3.60000000000000e-05;3.20000000000000e-05;4.20000000000000e-05;3.30000000000000e-05;3.90000000000000e-05;6.10000000000000e-05;7.30000000000000e-05;7.70000000000000e-05;5.60000000000000e-05;7.30000000000000e-05;6.40000000000000e-05;4.50000000000000e-05;5.30000000000000e-05;4.10000000000000e-05;5.30000000000000e-05;4.20000000000000e-05;1.00000000000000e-05;2.20000000000000e-05;2.70000000000000e-05;3.10000000000000e-05;4.50000000000000e-05;5.50000000000000e-05;6.10000000000000e-05;6.20000000000000e-05;4.70000000000000e-05;4.80000000000000e-05;4.80000000000000e-05;3.50000000000000e-05;2.90000000000000e-05;3.70000000000000e-05;3.90000000000000e-05;2.00000000000000e-05;-2.20000000000000e-05;-4.20000000000000e-05;-3.70000000000000e-05;5.00000000000000e-06;1.50000000000000e-05;1.90000000000000e-05;3.50000000000000e-05;4.60000000000000e-05;1.20000000000000e-05;-3.60000000000000e-05;-3.90000000000000e-05;-4.80000000000000e-05;-5.00000000000000e-05;-4.50000000000000e-05;-4.80000000000000e-05;-3.30000000000000e-05;-3.00000000000000e-05;-4.40000000000000e-05;-2.70000000000000e-05;-2.50000000000000e-05;-1.00000000000000e-05;-1.20000000000000e-05;-1.80000000000000e-05;-3.10000000000000e-05;-1.70000000000000e-05;3.70000000000000e-05;2.40000000000000e-05;0;7.00000000000000e-06;3.40000000000000e-05;1.10000000000000e-05;-7.00000000000000e-06;-1.40000000000000e-05;-2.30000000000000e-05;-7.00000000000000e-06;1.70000000000000e-05;3.40000000000000e-05;3.50000000000000e-05;1.80000000000000e-05;1.80000000000000e-05;4.10000000000000e-05;4.00000000000000e-05;4.60000000000000e-05;4.40000000000000e-05;4.60000000000000e-05;3.50000000000000e-05;2.60000000000000e-05;2.20000000000000e-05;2.40000000000000e-05;9.00000000000000e-06;-3.00000000000000e-06;-4.00000000000000e-06;7.00000000000000e-06;-1.30000000000000e-05;-3.50000000000000e-05;-3.20000000000000e-05;-3.40000000000000e-05;-4.00000000000000e-06;-1.40000000000000e-05;-3.90000000000000e-05;-9.00000000000000e-06;2.00000000000000e-06;-1.00000000000000e-05;1.80000000000000e-05;2.10000000000000e-05;2.00000000000000e-06;2.20000000000000e-05;3.30000000000000e-05;4.10000000000000e-05;3.90000000000000e-05;2.90000000000000e-05;1.70000000000000e-05;2.10000000000000e-05;2.20000000000000e-05;2.30000000000000e-05;2.10000000000000e-05;2.20000000000000e-05;1.50000000000000e-05;1.40000000000000e-05;3.80000000000000e-05;2.80000000000000e-05;4.10000000000000e-05;4.10000000000000e-05;4.30000000000000e-05;4.20000000000000e-05;3.90000000000000e-05;3.00000000000000e-05;3.40000000000000e-05;3.40000000000000e-05;5.70000000000000e-05;4.10000000000000e-05;2.10000000000000e-05;3.00000000000000e-05;2.10000000000000e-05;5.00000000000000e-05;4.80000000000000e-05;4.50000000000000e-05;7.00000000000000e-05;3.50000000000000e-05;2.70000000000000e-05;4.20000000000000e-05;5.00000000000000e-05;2.40000000000000e-05;1.80000000000000e-05;7.00000000000000e-06;2.00000000000000e-05;3.20000000000000e-05;2.40000000000000e-05;3.30000000000000e-05;1.00000000000000e-05;7.00000000000000e-06;1.90000000000000e-05;1.70000000000000e-05;6.00000000000000e-06;-5.00000000000000e-06;1.00000000000000e-05;2.00000000000000e-05;3.10000000000000e-05;2.00000000000000e-05;2.10000000000000e-05;-1.00000000000000e-06;-1.80000000000000e-05;-2.20000000000000e-05;-4.00000000000000e-06;9.00000000000000e-06;3.00000000000000e-06;-1.00000000000000e-06;-1.60000000000000e-05;-2.60000000000000e-05;-2.10000000000000e-05;1.00000000000000e-06;1.20000000000000e-05;2.00000000000000e-05;2.00000000000000e-06;2.00000000000000e-06;2.50000000000000e-05;5.00000000000000e-06;4.00000000000000e-06;3.00000000000000e-06;1.20000000000000e-05;-2.00000000000000e-06;7.00000000000000e-06;1.60000000000000e-05;1.00000000000000e-05;-5.00000000000000e-06;-1.90000000000000e-05;-3.50000000000000e-05;-3.70000000000000e-05;-2.30000000000000e-05;-3.30000000000000e-05;-4.20000000000000e-05;-4.50000000000000e-05;-3.50000000000000e-05;-2.90000000000000e-05;-5.00000000000000e-06;1.00000000000000e-05;2.70000000000000e-05;4.70000000000000e-05;7.40000000000000e-05;6.70000000000000e-05;2.00000000000000e-05;-1.10000000000000e-05;-1.60000000000000e-05;-2.00000000000000e-05;-2.60000000000000e-05;-2.70000000000000e-05;-3.70000000000000e-05;-3.90000000000000e-05;-5.30000000000000e-05;-7.70000000000000e-05;-4.40000000000000e-05;-3.40000000000000e-05;-4.30000000000000e-05;-5.70000000000000e-05;-5.50000000000000e-05;-5.90000000000000e-05;-6.90000000000000e-05;-5.30000000000000e-05;-4.30000000000000e-05;-5.30000000000000e-05;-4.50000000000000e-05;-4.20000000000000e-05;-3.30000000000000e-05;-4.30000000000000e-05;-6.90000000000000e-05;-6.90000000000000e-05;-5.20000000000000e-05;-5.20000000000000e-05;-7.90000000000000e-05;-8.80000000000000e-05;-5.40000000000000e-05;-4.20000000000000e-05;-1.70000000000000e-05;-1.20000000000000e-05;-2.00000000000000e-05;-1.90000000000000e-05;-2.00000000000000e-06;3.00000000000000e-06;-5.00000000000000e-06;-1.60000000000000e-05;-3.00000000000000e-05;-2.00000000000000e-05;-8.00000000000000e-06;-2.30000000000000e-05;-2.10000000000000e-05;-2.60000000000000e-05;-1.80000000000000e-05;8.00000000000000e-06;2.30000000000000e-05;3.50000000000000e-05;2.50000000000000e-05;-2.30000000000000e-05;-9.00000000000000e-06;-5.00000000000000e-06;-1.00000000000000e-05;-4.50000000000000e-05;-6.80000000000000e-05;-7.30000000000000e-05;-6.70000000000000e-05;-5.00000000000000e-05;-3.70000000000000e-05;-5.10000000000000e-05;-5.60000000000000e-05;-7.70000000000000e-05;-9.30000000000000e-05;-9.30000000000000e-05;-0.000106000000000000;-0.000118000000000000;-0.000130000000000000;-9.10000000000000e-05;-8.40000000000000e-05;-0.000107000000000000;-9.10000000000000e-05;-9.00000000000000e-05;-0.000113000000000000;-0.000124000000000000;-0.000147000000000000;-0.000139000000000000;-0.000142000000000000;-0.000137000000000000;-0.000152000000000000;-0.000148000000000000;-0.000143000000000000;-0.000152000000000000;-0.000139000000000000;-0.000133000000000000;-0.000152000000000000;-0.000159000000000000;-0.000145000000000000;-0.000144000000000000;-0.000145000000000000;-0.000101000000000000;-0.000130000000000000;-0.000153000000000000;-0.000146000000000000;-0.000177000000000000;-0.000218000000000000;-0.000223000000000000;-0.000223000000000000];%x是一个500*1的输入向量
    ;
    %x是一个500*1的输入向量
    %%%
    window=100;%滑动窗大小
    noverlap=99;%重叠样本数为100重叠99平移样本,每隔1个单位画一个框计算
    %%%
    f=1:30;%1-30hz变换
    Fs=250;%采样频率
	[B,F, T, P] = spectrogram(x',window,noverlap,f,Fs);  
    % B是F大小行T大小列的频率峰值，P是对应的能量谱密度
	figure
	imagesc(T,F,P);
	set(gca,'YDir','normal')
	colorbar;
	xlabel('时间 t/s');
	ylabel('频率 f/Hz');
	title('短时傅里叶时频图');

结果：在这里插入图片描述

那么有人要问了。为什么我选取信号时间长度和频谱图中时间长度不同。
原因是短时傅里叶变化选取的框的大小会把时间点长度给覆盖，不满足框大小的时间段。
所以上面的程序中。
原始为500个采样点（250hz，2秒），matlab使用了间隔为1个步长的100大小的框。那么（500-100）/1=400计算区间。转换到时间就是400/250hz=1.6秒，
matlab原始时间是0-2s，现在成了0.2-1.8s。实际上就是1.6s，matlab将时间中点作为舍取中心。于是左右就合并了50个采样点也就是说50/250=0.2s的单位。

如果希望解决这个问题，那么避免时间框选择太大，（选择小框缺点是计算损失一些低频信息，如1hz的信号分辨精度下降）
我修改窗函数大小的结果，
window=60；
noverlap=59;%滑动间隔依旧保持1个步长

  x=[7.40000000000000e-05;6.50000000000000e-05;8.20000000000000e-05;0.000104000000000000;9.00000000000000e-05;8.70000000000000e-05;6.30000000000000e-05;4.80000000000000e-05;4.70000000000000e-05;2.60000000000000e-05;1.40000000000000e-05;2.20000000000000e-05;4.30000000000000e-05;2.00000000000000e-05;2.60000000000000e-05;2.40000000000000e-05;8.00000000000000e-06;4.00000000000000e-06;1.40000000000000e-05;1.70000000000000e-05;1.60000000000000e-05;2.30000000000000e-05;3.90000000000000e-05;4.10000000000000e-05;2.90000000000000e-05;2.10000000000000e-05;1.70000000000000e-05;1.20000000000000e-05;9.00000000000000e-06;1.80000000000000e-05;2.10000000000000e-05;2.10000000000000e-05;2.40000000000000e-05;1.40000000000000e-05;3.00000000000000e-05;5.30000000000000e-05;7.00000000000000e-05;4.40000000000000e-05;5.40000000000000e-05;9.10000000000000e-05;3.00000000000000e-06;-1.50000000000000e-05;-7.00000000000000e-06;-4.00000000000000e-05;-6.10000000000000e-05;-5.90000000000000e-05;-8.60000000000000e-05;-0.000108000000000000;-0.000102000000000000;-9.10000000000000e-05;-9.10000000000000e-05;-8.20000000000000e-05;-6.40000000000000e-05;-5.90000000000000e-05;-5.40000000000000e-05;-6.10000000000000e-05;-5.50000000000000e-05;-5.80000000000000e-05;-6.20000000000000e-05;-6.40000000000000e-05;-5.60000000000000e-05;-6.00000000000000e-05;-6.00000000000000e-05;-6.50000000000000e-05;-7.10000000000000e-05;-7.30000000000000e-05;-7.10000000000000e-05;-6.30000000000000e-05;-5.20000000000000e-05;-5.80000000000000e-05;-5.10000000000000e-05;-4.80000000000000e-05;-4.30000000000000e-05;-6.40000000000000e-05;-9.20000000000000e-05;-0.000106000000000000;-9.90000000000000e-05;-0.000112000000000000;-0.000117000000000000;-0.000112000000000000;-0.000106000000000000;-8.10000000000000e-05;-8.30000000000000e-05;-9.40000000000000e-05;-0.000101000000000000;-9.30000000000000e-05;-9.20000000000000e-05;-8.70000000000000e-05;-0.000114000000000000;-0.000117000000000000;-7.10000000000000e-05;-3.50000000000000e-05;-1.10000000000000e-05;1.10000000000000e-05;2.40000000000000e-05;1.50000000000000e-05;-2.60000000000000e-05;-0.000103000000000000;-0.000152000000000000;-7.60000000000000e-05;-4.90000000000000e-05;-8.10000000000000e-05;-4.70000000000000e-05;-1.10000000000000e-05;-7.00000000000000e-06;-1.20000000000000e-05;-1.00000000000000e-05;-2.30000000000000e-05;-3.90000000000000e-05;-2.60000000000000e-05;-1.20000000000000e-05;-2.50000000000000e-05;-3.00000000000000e-05;-1.10000000000000e-05;-3.00000000000000e-06;-2.00000000000000e-06;-8.00000000000000e-06;-3.50000000000000e-05;-3.70000000000000e-05;-2.50000000000000e-05;-2.00000000000000e-06;-1.10000000000000e-05;-1.40000000000000e-05;-2.70000000000000e-05;-3.10000000000000e-05;-1.10000000000000e-05;2.00000000000000e-06;1.00000000000000e-05;2.10000000000000e-05;4.20000000000000e-05;3.10000000000000e-05;3.10000000000000e-05;4.40000000000000e-05;3.10000000000000e-05;2.90000000000000e-05;3.20000000000000e-05;3.20000000000000e-05;4.50000000000000e-05;4.00000000000000e-05;4.40000000000000e-05;5.30000000000000e-05;3.60000000000000e-05;4.00000000000000e-05;4.20000000000000e-05;4.30000000000000e-05;4.10000000000000e-05;5.40000000000000e-05;9.70000000000000e-05;8.10000000000000e-05;6.10000000000000e-05;5.60000000000000e-05;7.00000000000000e-05;7.40000000000000e-05;6.10000000000000e-05;7.80000000000000e-05;9.10000000000000e-05;7.30000000000000e-05;6.20000000000000e-05;4.80000000000000e-05;4.50000000000000e-05;5.30000000000000e-05;5.20000000000000e-05;2.50000000000000e-05;1.90000000000000e-05;1.10000000000000e-05;2.00000000000000e-05;4.40000000000000e-05;3.00000000000000e-05;1.70000000000000e-05;1.80000000000000e-05;3.00000000000000e-05;3.00000000000000e-05;3.30000000000000e-05;4.50000000000000e-05;1.80000000000000e-05;8.00000000000000e-06;2.30000000000000e-05;3.00000000000000e-05;1.50000000000000e-05;2.80000000000000e-05;2.10000000000000e-05;4.10000000000000e-05;4.30000000000000e-05;4.60000000000000e-05;5.90000000000000e-05;3.70000000000000e-05;2.60000000000000e-05;2.60000000000000e-05;9.00000000000000e-06;1.30000000000000e-05;2.10000000000000e-05;2.90000000000000e-05;3.70000000000000e-05;5.30000000000000e-05;5.00000000000000e-05;3.90000000000000e-05;2.50000000000000e-05;2.60000000000000e-05;3.00000000000000e-05;4.40000000000000e-05;5.80000000000000e-05;5.20000000000000e-05;5.40000000000000e-05;8.00000000000000e-05;9.20000000000000e-05;9.20000000000000e-05;7.50000000000000e-05;4.20000000000000e-05;3.20000000000000e-05;3.60000000000000e-05;3.20000000000000e-05;4.20000000000000e-05;3.30000000000000e-05;3.90000000000000e-05;6.10000000000000e-05;7.30000000000000e-05;7.70000000000000e-05;5.60000000000000e-05;7.30000000000000e-05;6.40000000000000e-05;4.50000000000000e-05;5.30000000000000e-05;4.10000000000000e-05;5.30000000000000e-05;4.20000000000000e-05;1.00000000000000e-05;2.20000000000000e-05;2.70000000000000e-05;3.10000000000000e-05;4.50000000000000e-05;5.50000000000000e-05;6.10000000000000e-05;6.20000000000000e-05;4.70000000000000e-05;4.80000000000000e-05;4.80000000000000e-05;3.50000000000000e-05;2.90000000000000e-05;3.70000000000000e-05;3.90000000000000e-05;2.00000000000000e-05;-2.20000000000000e-05;-4.20000000000000e-05;-3.70000000000000e-05;5.00000000000000e-06;1.50000000000000e-05;1.90000000000000e-05;3.50000000000000e-05;4.60000000000000e-05;1.20000000000000e-05;-3.60000000000000e-05;-3.90000000000000e-05;-4.80000000000000e-05;-5.00000000000000e-05;-4.50000000000000e-05;-4.80000000000000e-05;-3.30000000000000e-05;-3.00000000000000e-05;-4.40000000000000e-05;-2.70000000000000e-05;-2.50000000000000e-05;-1.00000000000000e-05;-1.20000000000000e-05;-1.80000000000000e-05;-3.10000000000000e-05;-1.70000000000000e-05;3.70000000000000e-05;2.40000000000000e-05;0;7.00000000000000e-06;3.40000000000000e-05;1.10000000000000e-05;-7.00000000000000e-06;-1.40000000000000e-05;-2.30000000000000e-05;-7.00000000000000e-06;1.70000000000000e-05;3.40000000000000e-05;3.50000000000000e-05;1.80000000000000e-05;1.80000000000000e-05;4.10000000000000e-05;4.00000000000000e-05;4.60000000000000e-05;4.40000000000000e-05;4.60000000000000e-05;3.50000000000000e-05;2.60000000000000e-05;2.20000000000000e-05;2.40000000000000e-05;9.00000000000000e-06;-3.00000000000000e-06;-4.00000000000000e-06;7.00000000000000e-06;-1.30000000000000e-05;-3.50000000000000e-05;-3.20000000000000e-05;-3.40000000000000e-05;-4.00000000000000e-06;-1.40000000000000e-05;-3.90000000000000e-05;-9.00000000000000e-06;2.00000000000000e-06;-1.00000000000000e-05;1.80000000000000e-05;2.10000000000000e-05;2.00000000000000e-06;2.20000000000000e-05;3.30000000000000e-05;4.10000000000000e-05;3.90000000000000e-05;2.90000000000000e-05;1.70000000000000e-05;2.10000000000000e-05;2.20000000000000e-05;2.30000000000000e-05;2.10000000000000e-05;2.20000000000000e-05;1.50000000000000e-05;1.40000000000000e-05;3.80000000000000e-05;2.80000000000000e-05;4.10000000000000e-05;4.10000000000000e-05;4.30000000000000e-05;4.20000000000000e-05;3.90000000000000e-05;3.00000000000000e-05;3.40000000000000e-05;3.40000000000000e-05;5.70000000000000e-05;4.10000000000000e-05;2.10000000000000e-05;3.00000000000000e-05;2.10000000000000e-05;5.00000000000000e-05;4.80000000000000e-05;4.50000000000000e-05;7.00000000000000e-05;3.50000000000000e-05;2.70000000000000e-05;4.20000000000000e-05;5.00000000000000e-05;2.40000000000000e-05;1.80000000000000e-05;7.00000000000000e-06;2.00000000000000e-05;3.20000000000000e-05;2.40000000000000e-05;3.30000000000000e-05;1.00000000000000e-05;7.00000000000000e-06;1.90000000000000e-05;1.70000000000000e-05;6.00000000000000e-06;-5.00000000000000e-06;1.00000000000000e-05;2.00000000000000e-05;3.10000000000000e-05;2.00000000000000e-05;2.10000000000000e-05;-1.00000000000000e-06;-1.80000000000000e-05;-2.20000000000000e-05;-4.00000000000000e-06;9.00000000000000e-06;3.00000000000000e-06;-1.00000000000000e-06;-1.60000000000000e-05;-2.60000000000000e-05;-2.10000000000000e-05;1.00000000000000e-06;1.20000000000000e-05;2.00000000000000e-05;2.00000000000000e-06;2.00000000000000e-06;2.50000000000000e-05;5.00000000000000e-06;4.00000000000000e-06;3.00000000000000e-06;1.20000000000000e-05;-2.00000000000000e-06;7.00000000000000e-06;1.60000000000000e-05;1.00000000000000e-05;-5.00000000000000e-06;-1.90000000000000e-05;-3.50000000000000e-05;-3.70000000000000e-05;-2.30000000000000e-05;-3.30000000000000e-05;-4.20000000000000e-05;-4.50000000000000e-05;-3.50000000000000e-05;-2.90000000000000e-05;-5.00000000000000e-06;1.00000000000000e-05;2.70000000000000e-05;4.70000000000000e-05;7.40000000000000e-05;6.70000000000000e-05;2.00000000000000e-05;-1.10000000000000e-05;-1.60000000000000e-05;-2.00000000000000e-05;-2.60000000000000e-05;-2.70000000000000e-05;-3.70000000000000e-05;-3.90000000000000e-05;-5.30000000000000e-05;-7.70000000000000e-05;-4.40000000000000e-05;-3.40000000000000e-05;-4.30000000000000e-05;-5.70000000000000e-05;-5.50000000000000e-05;-5.90000000000000e-05;-6.90000000000000e-05;-5.30000000000000e-05;-4.30000000000000e-05;-5.30000000000000e-05;-4.50000000000000e-05;-4.20000000000000e-05;-3.30000000000000e-05;-4.30000000000000e-05;-6.90000000000000e-05;-6.90000000000000e-05;-5.20000000000000e-05;-5.20000000000000e-05;-7.90000000000000e-05;-8.80000000000000e-05;-5.40000000000000e-05;-4.20000000000000e-05;-1.70000000000000e-05;-1.20000000000000e-05;-2.00000000000000e-05;-1.90000000000000e-05;-2.00000000000000e-06;3.00000000000000e-06;-5.00000000000000e-06;-1.60000000000000e-05;-3.00000000000000e-05;-2.00000000000000e-05;-8.00000000000000e-06;-2.30000000000000e-05;-2.10000000000000e-05;-2.60000000000000e-05;-1.80000000000000e-05;8.00000000000000e-06;2.30000000000000e-05;3.50000000000000e-05;2.50000000000000e-05;-2.30000000000000e-05;-9.00000000000000e-06;-5.00000000000000e-06;-1.00000000000000e-05;-4.50000000000000e-05;-6.80000000000000e-05;-7.30000000000000e-05;-6.70000000000000e-05;-5.00000000000000e-05;-3.70000000000000e-05;-5.10000000000000e-05;-5.60000000000000e-05;-7.70000000000000e-05;-9.30000000000000e-05;-9.30000000000000e-05;-0.000106000000000000;-0.000118000000000000;-0.000130000000000000;-9.10000000000000e-05;-8.40000000000000e-05;-0.000107000000000000;-9.10000000000000e-05;-9.00000000000000e-05;-0.000113000000000000;-0.000124000000000000;-0.000147000000000000;-0.000139000000000000;-0.000142000000000000;-0.000137000000000000;-0.000152000000000000;-0.000148000000000000;-0.000143000000000000;-0.000152000000000000;-0.000139000000000000;-0.000133000000000000;-0.000152000000000000;-0.000159000000000000;-0.000145000000000000;-0.000144000000000000;-0.000145000000000000;-0.000101000000000000;-0.000130000000000000;-0.000153000000000000;-0.000146000000000000;-0.000177000000000000;-0.000218000000000000;-0.000223000000000000;-0.000223000000000000];%x是一个500*1的输入向量
    ;%x是一个500*1的输入向量
    %%%
    window=60;%滑动窗大小
    noverlap=59;%重叠样本数为250重叠249平移样本,每隔1个单位画一个框计算
    %%%
    f=1:30;%1-30hz变换
    Fs=250;%采样频率
	[B,F, T, P] = spectrogram(x',window,noverlap,f,Fs);  
    % B是F大小行T大小列的频率峰值，P是对应的能量谱密度
	figure
	imagesc(T,F,P);
	set(gca,'YDir','normal')
	colorbar;
	xlabel('时间 t/s');
	ylabel('频率 f/Hz');
	title('短时傅里叶时频图');

采样损失好了很多。
在这里插入图片描述
时间上恢复接近了2s，但是stft的弊端凸显出来了。和上图相比
短时傅里叶的缺点
1窗函数宽度固定不变，时间和频率的局部化不能同时满足
2窗函数越大，越接近实际的信号傅里叶变化处理，但是时间长度缩短。
3窗函数越小，时间长度基本保持连续，但是频谱的低频信息（幅值大小）会损失。

所以要根据实际信号分析的特点需求选取合适大小的滑动窗是关键。
另外一位博主讲解matlab中2019直接引入了STFT函数的使用
可以看一下。
[https://blog.csdn.net/zhanghaijun2013/article/details/106331270]

高山仰止景
关注
  • 44
    点赞
  • 踩
  • 199
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 15
    评论
专栏目录
短时傅里叶变换（stft）Matlab源码
06-07
短时傅里叶变换（stft）的Matlab源码，已封装为Function函数，可直接根据源码注释调用
使用 Matlab 进行短时傅立叶变换 (STFT)：通过短时傅立叶变换 (STFT) 进行时频分析。-matlab开发
05-29
本代码是一个 Matlab 函数，它提供给定信号 x[n] 的短时傅立叶变换 (STFT)。 该函数是 Matlab 命令“spectrogram”的替代方法。 函数的输出是： 1) 具有复数 STFT 系数的矩阵，列上有时间，行上有频率； 2) 频率向量； 3）时间向量。 为了阐明该功能的使用，给出了一个例子。 为了方便起见，在函数的开头给出了输入和输出参数。 该代码基于以下内容中描述的理论： [1] H.日沃米罗夫。 关于 STFT 分析和 ISTFT 合成例程的开发及其实际实现。 TEM 期刊，ISSN：2217-8309，DOI：10.18421/TEM81-07，卷。 8，第 1 期，第 56-64 页，2019 年 2 月。（ http://www.temjournal.com/content/81/TEMJournalFebruary2019_56_64.pdf ）
Matlab中短时傅里叶变换 spectrogram和stft的用法
zhanghaijun2013的博客
05-25 1万+
  在Matlab中，做短时傅里叶变换需要使用函数spectrogram，而在Matlab2019中，引入了一个新的函数stft，下面我们就来看下这两个函数都如何使用。   短时傅里叶变换的基本原理就是将数据分段加窗，做fft，在分段时会有overlap，因此一个向量的短时傅里叶变换结果是一个矩阵。了解了这点，下面的函数及参数就更加容易理解了。 spectrogram 参数列表   先来看spectrogram函数，在更早期的版本中，这个函数的名字是specgram，几种常用的用法如下： spectrogr
使用 Matlab 进行短时傅里叶变换 (STFT)（matlab实现）
最新发布
你的Matlab小助手
03-02 836
短时傅里叶变换（STFT）是一种在信号处理领域广泛应用的技术，用于分析信号在时频域上的特性。相比于对整个信号进行傅里叶变换，STFT将信号分成短时段，然后对每个时段进行傅里叶变换，以获取信号在不同时间和频率上的频谱信息。MATLAB作为一种强大的数学工具和编程环境，提供了丰富的函数库，方便了短时傅里叶变换的实现。通过使用MATLAB，工程师和研究人员能够轻松地对复杂信号进行时频分析，并可视化这些分析结果。文章中一些内容引自网络，会注明出处或引用为参考文献，难免有未尽之处，如有不妥，请随时联系删除。
语音.zip_spectrogram_spectrogram函数_自编语谱图函数_语谱图
07-14
自编语谱图函数，可实现基本的spectrogram功能
spectrogram函数
LYLYC_3的博客
04-02 6739
360百科：https://baike.so.com/doc/7161179-7385189.html spectrogram是一个MATLAB函数，使用短时傅里叶变换得到信号的频谱图。当使用时无输出参数，会自动绘制频谱图;有输出参数，则会返回输入信号的短时傅里叶变换。 中文名称：波谱图; 频谱图; 光谱图; 外文名称：spectrogram 功能：short-time...
Matlab中stft函数的使用与原理
飞行空间
07-15 2万+
1 引入 频谱分析的方法是建立在如下假设之上： 连续时间信号是频谱上的带限信号，其采样后的离散值可以提供原信号频谱的信息。带限信号对应的时域信号定义域往往是 。因此，往往采用抗混叠滤波后再采样的方式。 由于实际的信号往往是信号参数时变的信号，即信号的频率，幅度，相位随时间变化。使用傅里叶分析就不能体现出时变的特性。 Questions 1、如果需要分析工程实际中的信号，会有什么样的具体需求？...
解析MATLAB短时傅里叶变换函数spectrogram()
平淡多好
03-15 3万+
BB： 最近做脑电信号的时频分析，考虑采用该函数，但在参数设置和输出结果上有些疑惑，参考了众多资料，现对出现的困惑进行一些总结。 如何理解spectrogram()函数的输入输出参数？ Abstract: 我想这篇博文可以帮你弄清楚这几个问题，spectrogram()函数有无返回值调用时图像如何转换，函数输入输出参数的的含义，以及其对谱图分辨力的影响。 Note: spectrogram()函数...
matlab时频分析之短时傅里叶变换 spectrogram
热门推荐
hyhhyh21的博客
03-22 6万+
matlab时频分析之短时傅里叶变换 spectrogram 短时傅里叶变换常用于缓慢时变信号的频谱分析，可以观察沿时间变化的频谱信号。 其优点如下图所示，弥补了频谱分析中不能观察时间的缺点，也弥补了时域分析不能获取频率的缺点。 1 STFT的基本原理 基本原理可以理解为对一段长信号，截取每一段时间的短信号做fft，将得到的频谱图时间沿时间轴排列，及可得到时频的云图。 2 matlab中实现 ...
matlab spectrogram 函数用法
cocofisher的博客
08-05 1万+
项目说明 介绍 matlab 中 spectrogram 函数的用法，其功能为实现短时傅里叶变换并得到信号的频谱图。 语法： [S,F,T,P]=spectrogram(x,window,noverlap,nfft,fs) 短时傅里叶变换 短时傅里叶变换的基本原理如下： 使用一个窗函数在时域信号上（有重叠地）不断平移，如下图所示： 其中，window 指窗的长度，noverlap 指窗的上一个位置和下一个位置之间重叠的长度（实际上noverlap>window/2noverlap>wind
MATLAB 短时傅里叶变换函数运用 - spectrogram
wanjiac的博客
07-09 3821
函数形式： spectrogram support documents
语音特征提取: 看懂梅尔语谱图(Mel-spectrogram)、梅尔倒频系数(MFCCs)的原理
weixin_50547200的博客
05-27 4万+
本文内容主要来自于:Speech Processing for Machine Learning: Filter banks, Mel-Frequency Cepstral Coefficients (MFCCs) and What’s In-Between | Haytham Fayek 1. 什么是梅尔语谱图和梅尔倒频系数？ 机器学习的第一步都是要提取出相应的特征(feature)，如果输入数据是图片，例如28*28的图片，那么只需要把每个像素(pixel)作为特征，对应的像素值大小(代表颜色的强.
spectrogram_fft.rar_FFT音频_Spectrogram matlab_绘制语谱图_语音频谱图_音频spect
07-15
matlab用来绘制wav音频文件的语谱图的程序
全面剖析STFT(短时傅里叶变换)spectrogram函数用法
06-21
短时傅里叶变换（STFT）可以有效地处理非平稳数字信号，给出瞬时频率的信息。STFT的本质就是将一段信号截成多段，对每一段进行单独傅里叶变换（FFT），最后在时域上接起来，就形成了STFT的结果。Matlab中关于STFT的函数是spectrogram，有界面化与代码两种方式，本视频对这两种方式进行了深度解剖，一看就会。下载本视频的朋友，请加QQ714362241免费索取完整版的视频（完整版视频85分钟） Part1：spectrogram函数是啥，有何作用？ Part2：为什么要用到spectrogram函数？ Part3：spectrogram函数的过程 是怎样的？ Part4：界面化的spectrogram。 Part5：如何准确无误地使用代码spectrogram。 Part6：酷炫地展示spectrogram结果。
spectrogram.rar_matlab语谱图_spectrogram stft_傅里叶_时谱图_短时傅
07-14
基于MATLAB的短时傅里叶变换（STFT）的语谱图的画法，MATLAB环境下
spectrogram:显示语音的频谱图-matlab开发
05-30
用颜色图 (jet) 显示语音的频谱图。
基于matlab使用脉冲波形分析仪应用程序绘制频谱图
03-22
例如，使用该应用程序显示连续FMCW波形的频谱图。 1. 将波形设置为FMCW 2. 将扫描间隔设置为Symmetric 3. 将扫描次数设置为4 4. 在“分析”选项卡中，选择Spectrogram 然后，您将看到与此类似的信号频谱图。
空气质量指数计算公式
高山
02-28 1万+
空气质量指数是通过计算各个污染物空气质量分指数得到的 计算每一个污染物的空气质量分指数在对所有的分指数取最大值就可得到： 其中IAQI_p为污染物项目p的空气质量分指数；Cp为污染物项目p的质量浓度值，BPhi为污染物浓度限高位值；BPlo为污染物浓度限低位值；IAQIhi为BPhi对应的空气质量分指数；IAQIlo为BPlo对应的空气质量分指数。 这些数值在表格中有相应对应的值 matlab代码如下： function o = caliaqi(i) %UNTITLED %%根据限度值和浓度值计算iaq
matlab STFT 高斯窗
09-06
高斯窗在matlab中的STFT函数中可以使用。在matlab中，使用高斯窗可以通过指定窗长来选择适当的滑动窗口大小。根据实际信号分析的需求，选择合适大小的滑动窗口是关键。在MATLAB 2019版本中，引入了STFT函数的使用，可以通过该函数来进行频谱分析。 在使用STFT函数时，可以设置nfft参数来表示FFT的点数。当没有指定nfft时，MATLAB会根据窗长自动选择一个合适的FFT点数，取值为max(256,2^(ceil(log2(length(window)))))。当窗长小于256时，FFT点数为256；当窗长大于256时，FFT点数为大于窗长的最小的2的整数次幂。需要注意的是，nfft参数应该是double或single类型的，而不是一个向量。 在MATLAB的spectrogram函数中，可以使用高斯窗来进行短时傅里叶变换，生成频谱图。spectrogram函数的参数包括x（输入信号）、window（窗函数）、noverlap（重叠部分的样本数）、nfft（FFT点数）和fs（采样率）。可以根据需求选择相应的参数来进行频谱分析。 总结起来，MATLAB中可以使用高斯窗进行STFT分析，可以通过设置窗长和FFT点数来选择适当的滑动窗口大小和频谱分辨率。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *2* *3* [spectrogram（频谱图函数）matlab使用STFT参数选择以及滑动窗的问题](https://blog.csdn.net/qq_43158059/article/details/115032195)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_1"}}] [.reference_item style="max-width: 100%"] [ .reference_list ]

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 15
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

高山仰止景

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值