python绘制时频图

最近一直在做轴承RUL。在进行xinh信号处理上面遇见了很多的问题，想要用python进行信号时频图的绘制。但是在选择绘制图的函数上也是煞费苦心。

首先来说plt.imshow()函数，虽然他的功能强大，但是他并没有相对应的可以将（时间，频率，CWT处理后的信号）进行联合绘制的参数设定。在使用plt.pcolor（）这个函数时，它运行速度极慢。后来选择使用了plt.pcolormesh()函数进行时频图的绘制。

为了保证所有的图像时频图一致性，需要设置X,Y轴，以及colorbar的大小。在设置colorbar的时候可以直接使用plt.pcolormesh()函数。plt.pcolormesh(t, frequencies, abs(cwtmatr), cmap=cm,vmin=0, vmax=3)这两个标红的参数进行调整。

代码见：https://download.csdn.net/download/x_uhen/10614778

简单的绘制时间，频率和数据之间的图代码如下：

# -*- coding: utf-8 -*-
import numpy as np
import pylab as plsampling_rate = 8000
fft_size = 512
t = np.arange(0, 1.0, 1.0/sampling_rate)
x = np.sin(2*np.pi*156.25*t)  + 2*np.sin(2*np.pi*234.375*t)
xs = x[:fft_size]
xf = np.fft.rfft(xs)/fft_size
freqs = np.linspace(0, sampling_rate/2, fft_size/2+1)
xfp = 20*np.log10(np.clip(np.abs(xf), 1e-20, 1e100))
pl.figure(figsize=(8,4))
pl.subplot(211)
pl.plot(t[:fft_size], xs)
pl.xlabel(u"时间(秒)")
pl.title(u"156.25Hz和234.375Hz的波形和频谱")
pl.subplot(212)
pl.plot(freqs, xfp)
pl.xlabel(u"频率(Hz)")
pl.subplots_adjust(hspace=0.4)
pl.show()

python绘制时频图相关推荐

python 绘制时频图 plt.specgram
时频图以横轴为时间,纵轴为频率,用颜色表示幅值.在一幅图中表示信号的频率.幅度随时间的变化 matplotlib.pyplot.specgram(x, NFFT=None, Fs=None, Fc=N ...
MATLAB短时傅里叶变换并绘制时频图
data=load('dataset/rawdata/eegdata4.mat'); Au=data.eegdata(:,1); Fs = 1200; [B, F, T, P] = spectrogr ...
Python实现“EMD\EEMD\VMD+Hilbert时频图”与“CWT小波时频图”
Python实现"EMD\EEMD\VMD+Hilbert时频图"与"CWT小波时频图" 信号处理中常需要分析时域统计量.频率成分,但不平稳信号的时域波形往 ...
融合DE 端和FE端数据，利用小波变换生成时频图，再分别利用DCNN、KNN和DNN进行对比实验(python代码)
1.数据集介绍: 试验台如图所示,试验台左侧有电动机,中间有扭矩收集器,右侧有动力测试仪,控制电子设备在图中没有显示.SKF6203轴承使用16通道数据采集卡采集轴承的振动数据,并在驱动端部分(DE) ...
基于小波时频图和2D-CNN的滚动轴承故障检测
目录一.研究思路 1.基于小波时频图和CNN的滚轴故障诊断方法的研究思路如下: 二.数据集介绍与数据处理 1.数据集介绍 2.数据集分割与合并 3.数据集分析三.小波时频图导出四.CNN网络的构 ...
信号时频图的横纵坐标对调
博客1:如何翻转信号的时频图论文有时候需要绘制信号的时频图.Matlab有专用的函数spectrogram可以用来帮助我们绘制信号的时频图.默认的图例是横轴为频率,纵轴为时间.如果希望将横轴和纵轴对 ...
python绘制语谱图（详细注释）
用python 绘制语谱图 1.步骤: 1)导入相关模块 2)读入音频并获取音频参数 3)将音频转化为可处理形式(注意读入的是字符串格式,需要转换成int或short型) 代码如下: import ...
python 股票图表_k线图分析法_【趣味案例】用Python绘制K线图，一眼看清股市状况...
本文介绍关于[趣味案例]用Python绘制K线图,一眼看清股市状况及神一般的裸k交易法,精髓就这三步,一目了然!精髓就这三步,一目了然!(附图解析)的相关内容. [趣味案例]用Python绘制K线图, ...
Python绘制桑基图Sankey,Pyecharts不显示html页面,桑基图只显示标题,原因总结
Python绘制桑基图Sankey,Pyecharts不显示html页面,桑基图只显示标题,原因总结说说使用上的三点注意: 1.nodes中的"name"不要重新命名,否则会不识 ...

python绘制时频图

python绘制时频图相关推荐

最新文章

热门文章