递推子空间方法辨识结构的时变固有频率——非平稳随机信号处理

比丘尼萧竹

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

yugang2010

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

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

比丘尼萧竹

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

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

ccc5458

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

比丘尼萧竹

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

yghit08

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

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

比丘尼萧竹

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

yghit08

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

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

比丘尼萧竹

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

看庞世伟的博士论文，发现N4SID方法更新时对广义HANKEL矩阵的定义有错，按照他的方法编程循环的时候出现了HANKEL矩阵的行数需要变化才能循环下去，后来找了他的另一篇期刊论文，发现了问题所在。可能我之前也没讲清楚，还是感谢您的解答