『壹』 Matlab中短时傅里叶变换中一些函数变换步骤不明,求解释!
(1)matlab是按列储存的,当然可以使1*256的了,不过需要自己设。傅里叶变化是对称变化,一般都回要求数据个数是答2的幂次。
(2)tfr对谁做fft变换,fft就和谁是同样大小的矩阵。
(3)上面说了,fft变换后的波谱是对称的,如果tfr是256的话,只需取256/2=128个就行。
(4)随机数组或原数据体中数据之间无关的话,那么变换后的振幅图和相位图都是上下波动的,密集线就是波动情况,密集程度与点距和线长有关。
『贰』 短时傅里叶变换窗函数窗口宽度的选择
首先,我的答案是针对于matlab时频分析工具箱的,在这个工具箱里,段师傅里叶变换的函数为tfrstft
。下面介绍一下用法
格式:
[tfr,
t,
f]
=
tfrstft(x)
%
计算时间序列x的短时傅里叶变换,参数tfr为短时傅里叶变换系数,t为系数tfr对应的时刻,f为归一化频率向量
[tfr,
t,
f]
=
tfrstft(x,
t)
%
计算对应时刻t的短时傅里叶变换
[tfr,
t,
f]
=
tfrstft(x,
t,
n)
%
计算n点对应时刻t的短时傅里叶变换
[tfr,
t,
f]
=
tfrstft(x,
t,
n,
h)
%
参数h为归一化频率平滑窗
[tfr,
t,
f]
=
tfrstft(x
t,
n,
h,
trace)
%
trace显示算法进程
说明:
x--信号
t--时间(缺省值为1:length(x))
n--频率数(缺省值为length(x))
h--频率滑窗,h归一化为单位能量(缺省值为hamming(n/4))
trace--如果非零,显示算法的进程(缺省值为0)
tfr--时频分解(为复值),频率轴观察范围为-0.5~0.5
也就是说,如果你想改变窗函数,你需要修改h这个变量,希望对你有帮助