请问由功率谱密度怎么得到时域信号

ialan

有MATLAB程序吗？

[ 本帖最后由 zhlong 于 2007-6-20 22:12 编辑 ]

zhlong

这个问题论坛有讨论过的，你找找看。

VibrationMaster

由功率谱只能得到自相关，无法得到原始的“时域信号”

zhlong

功率谱不是可以通过自相关做FFT 和 直接对信号做FFT（周期图法）两种方法得到吗？

zhangnan3509

同问

wanyeqing2003

论坛上有相关的帖子介绍这些内容。

可以用傅里叶变换处理，不过需要加入相位因素。

[ 本帖最后由 zhangnan3509 于 2007-7-21 16:58 编辑 ]

w89986581

原帖由 wanyeqing2003 于 2007-6-20 20:44 发表
论坛上有相关的帖子介绍这些内容。

可以用逆富里叶变换处理，不过需要加入相位因素。

没有相位信息是无法回复的,否则严重畸变.

zhlong

那么这个呢？
t=1:64;
signal=sin(t);
s=fft(signal,64);
X=real(ifft(s,64));
subplot(211),plot(X,'b'),grid on;title('IFFT得到的信号')
subplot(212),plot(signal,'b'),grid on;title('原始信号')

w89986581

实际测量时,一般是得不到谱密度的相位信息.
s=fft(signal,64);已经包含相位信息拉.
你尝试从ss*恢复到原来的信号,呵呵.

VibrationMaster

功率谱是描述平稳过程的一个概念，其逆傅立叶变换是自相关。
而通常时域信号是一个具体函数，从自相关中恢复不出来。
如果认为对象是一个确定性信号+弱噪声，且噪声为各态历经，那么按照一定时间对齐的傅立叶谱平均（与幅值谱有关）可以消除噪声影响，从而由傅立叶谱可恢复时域信号。

VibrationMaster

你演示的是傅立叶谱，不是功率谱

zhlong

其实和这位网友的问题一样
http://forum.vibunion.com/forum/viewthread.php?tid=44361&extra=page%3D1%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D186

[ 本帖最后由 zhlong 于 2007-6-20 22:37 编辑 ]

zhlong

呵呵，我自己没仔细想，s=fft(signal,64); 是一个复数自然包含了相位信息。
而求功率谱首先对s求了模abs(s)，然后通过式子ss=((abs(s)*2/N).^2)/2求得。所以通过功率谱ss最多也只能求得abs(s)，而无法求得s。

汗自己一个:loveliness:

[ 本帖最后由 zhlong 于 2007-6-20 22:37 编辑 ]

wanyeqing2003

我做过用功率谱生成时程曲线。

wanyeqing2003

提供的相位不同，时程曲线也不同。也就是说同一个功率谱曲线，可以生成无数个时程曲线。这是因为功率谱曲线不包含相位信息。