比丘尼萧竹 发表于 2014-6-3 10:14

那个问题解决了,还真心是里面一个函数的编写出了问题。函数是自己下载的工具箱里面的,别人编的也没细看。多谢了

yugang2010 发表于 2014-6-3 21:49

比丘尼萧竹 发表于 2014-5-28 13:00
这个分别是FFT和STFT的结果,两者之间差两倍了,可能是窗函数的原因,但是还是没有查出来

fft变换时用的是所有的数据点。而STFT用的是窗函数做的,如果STFT函数是用的数据点的一半做窗长度的话,则在频率轴上,所得到的数据点要等于fft方法得到的数据点的一半。这时候需要重新设置频率坐标轴,我估计你应该是这种情况。
无论所用的分析点数怎么变,你做出来的信号的频率只有你的采样频率有关。不可能出现不一致的情况。
分析点数的多少只影响频率分辨率的好坏。

比丘尼萧竹 发表于 2014-6-5 15:27

yugang2010 发表于 2014-6-3 21:49
fft变换时用的是所有的数据点。而STFT用的是窗函数做的,如果STFT函数是用的数据点的一半做窗长度的话, ...

恩,多谢了!程序里面确实用到fftshift来重新排列频率坐标轴,函数的返回值多乘了一个0.5,导致结果变成了fft的一半,现在已经解决了

ccc5458 发表于 2014-7-21 15:26

您所在的用户组无法下载或查看附件{:{44}:}

比丘尼萧竹 发表于 2014-9-16 10:06

不好意思,又有问题要问了,看递推子空间的算法,基于投影估计子空间跟踪(PAST)的递推子空间不是需要首先根据输入和输出构造HANkel矩阵(输入hankel阵和输出hankel阵)吗?单单做实验不是只能获取响应信息吗

yghit08 发表于 2014-9-23 17:52

比丘尼萧竹 发表于 2014-9-16 10:06
不好意思,又有问题要问了,看递推子空间的算法,基于投影估计子空间跟踪(PAST)的递推子空间不是需要首先 ...

确实需要分别构造输入输出的hankel矩阵,实验中有采集输入信号的时候,对于小型构件或者容易激励的件通常用锤击,这时候默认是初始激励条件,不用构造输入的hankel矩阵。

比丘尼萧竹 发表于 2014-12-16 20:56

你好,在编程的时候又遇到问题了,想请教下你。对于输入输出的HANKEL矩阵怎么形式的呢,里面有个行数是未知的呀,循环的时候数据长度是变的,此时HANKEL的行数怎么确定呢?

yghit08 发表于 2014-12-17 12:57

比丘尼萧竹 发表于 2014-12-16 20:56
你好,在编程的时候又遇到问题了,想请教下你。对于输入输出的HANKEL矩阵怎么形式的呢,里面有个行数是未 ...

没明白你说的问题,hankel矩阵既然是构造的,那么它的维数是确定的,至于如何确定这里又涉及到系统的阶的问题,有不少文献可以参考,但是多只是涉及hankel矩阵的列,因为构造的需要是列满秩矩阵(可能错了,我的认识有限),对行的要求不具体,这一点我也太清楚,多是在使用过程中试凑(可能也是我臆断了)

比丘尼萧竹 发表于 2014-12-18 10:29

yghit08 发表于 2014-12-17 12:57
没明白你说的问题,hankel矩阵既然是构造的,那么它的维数是确定的,至于如何确定这里又涉及到系统的阶的 ...

看庞世伟的博士论文,发现N4SID方法更新时对广义HANKEL矩阵的定义有错,按照他的方法编程循环的时候出现了HANKEL矩阵的行数需要变化才能循环下去,后来找了他的另一篇期刊论文,发现了问题所在。可能我之前也没讲清楚,还是感谢您的解答
页: 1 [2]
查看完整版本: 递推子空间方法辨识结构的时变固有频率——非平稳随机信号处理