STFT 变换求解

zfmeng1993 2018-04-16 07:12:28

STFT 变换相关参数求解，采样频率和目标频率之间是什么关系？

例如：采样频率=1000.0Hz, 窗大小2s, 想要得到4-8Hz的频谱，参数应该如何设置？是直接设置STFT变换后Y轴的坐标吗？

(time, data) = np.loadtxt(input_file,unpack=True, delimiter=",",usecols = (0,1))

fs = 1000.0 # 采样频率
f,t,Sxx = signal.spectrogram(data, fs, nperseg=512)

...全文

788 1 打赏收藏转发到动态举报

写回复

1 条回复

切换为时间正序

请发表友善的回复…

发表回复

半吊子Py全栈工程师 2018-04-17

打赏
举报

回复

短时傅里叶变换~~对这个没有研究~~

求解短时傅里叶变换，使用matlab自带库，以及信号类型识别

针对目前阶段对预应力锚索的无损检测数据分析不够系统、分析方法仍有不足的现状，设计并实现了预应力锚索无损检测数据分析系统。 (1)在理论、算法研究及创新方面: 首先，通过快速傅里叶变换FFT对采集的数据进行滤波；通过短时傅里叶变换STFT对信号进行时频特征分析；在已有的信号处理方法的基础上提出“STFT主频幅值比法”求解预应力锚索自由段、锚固段长度；最后，利用频差法求解预应力大小。 (2)在软件系统实现方面: 根据对锚索无损检测数据分析的任务需求及模块化编程思想将软件分为四个层面（数据访问及管理层、数据计算层、绘图及操作层、结果计算及展示层）；同时基于MVC框架思想、DirectX设计软件核心内容：模型计算、绘图及控制框架。最后，通过Delphi编程设计并实现软件系统。

[本科大学生创新项目] 结构健康监测SHM：预应力锚索无损检测分析(技术：DSP、FFT、STFT、Delphi) 针对目前阶段对预应力锚索的无损检测数据分析不够系统、分析方法仍有不足的现状，设计并实现了预应力锚索无损检测数据分析系统。 (1)在理论、算法研究及创新方面: 首先，通过快速傅里叶变换FFT对采集的数据进行滤波；通过短时傅里叶变换STFT对信号进行时频特征分析；在已有的信号处理方法的基础上提出“STFT主频幅值比法”求解预应力锚索自由段、锚固段长度；最后，利用频差法求解预应力大小。 (2)在软件系统实现方面: 根据对锚索无损检测数据分析的任务需求及模块化编程思想将软件分为四个层面（数据访问及管理层、数据计算层、绘图及操作层、结果计算及展示层）；同时基于MVC框架思想、DirectX设计软件核心内容：模型计算、绘图及控制框架。最后，通过Delphi编程设计并实现软件系统。

文章目录一、傅里叶变换1、傅里叶变换基础2、傅里叶变换形式规律3、一些关键概念3、python 实现傅里叶变换二、短时傅里叶变换（STFT）一、傅里叶变换 1、傅里叶变换基础参考：深入浅出的讲解傅里叶变换傅里叶变换：一段信号可以由若干频率不同的正弦信号叠加构成，DFT将将信号从时域变换到频域用途：滤波：从一个时域图里剔除组成它的某个频域成分很不容易；但是从频域图里剔除这个成分很容易求解微分方程故障诊断：健康的设备的时域曲线图是相似的，异常设备的是不同的，但这一点在时域图种不易被发现，转换到频

librosa是一个应用广泛的音频处理python库。在librosa中有一个方法叫做stft，功能是求音频的短时傅里叶变换。音频短时傅里叶变换后，在对音频取幅值，可以得到音频的线性谱。对线性谱进行mel刻度的加权求和，可以得到语音识别和语音合成中常用的mel谱。短时傅里叶变换的过程是先对音频分帧，再分别对每一帧傅里叶变换。在应用stft方法求解短时傅里叶变换时，发现求出的特征帧的数目有点反常。比如我有一个长度是400个点的音频，如果帧长是100，那么我自然而然的想到，最后应当得到4帧。但实际不是这个

37,720

社区成员

34,239

社区内容

发帖

与我相关

我的任务

社区管理员

加入社区

近7日
近30日
至今

加载中

查看更多榜单

试试用AI创作助手写篇文章吧

+ 用AI写文章