相位谱方法求速度引出的加窗问题

zsure79

首先两个波形信号是在不同位置采集的，有色散，对它有色散的部分加窗求FFT，根据得到的各频率的相位角求速度，c=2*pi*f*L / angel,
其中L是两个不同接受位置的距离，c为速度，但是好像加窗的位置有一点变化，速度就变化很大？？？
如图：其中下面的波形加窗位置3800-6200不变，上面的波形一个加窗位置是7800-10000，另一个是7800-10200. 得到的色散如下面两个图。
波形中其他部分强迫赋值为零。
否则如何很好的加窗？？？？？？

[ 本帖最后由 zsure79 于 2008-6-2 10:02 编辑 ]

zsure79

?????????????????????????????

songzy41

楼主实际上是要求相速度。但是楼主用的方法是“根据得到的各频率的相位角求速度”，但怎么得到相位角angle呢？

zsure79

相位角好像是这样来的：

f1=fft(signal1); f2=fft(signal2);
N=length(f1);
s12=f1(1:N/2).*conj(f2(1:N/2);
ang=angle(s12);
就是两个信号的fft相乘，其中一个是fft的共轭复数。
然后相位角等于这个乘积的相位。
这里的方法是：ultrasonics 34（1996）793-799
好像附件没有办法上传 pdf 文件。懂行的人帮给看一下
另外说点闲话啊：

不是每一个人都闲得发慌，整天都在这里灌水，看看帖子数就知道了，
这个论坛是讨论问题的，我觉得还是很纯净，一些人在这里拿着大棒乱挥，不知所谓！！！
好像就他自己是搞学问的，其他人都是在灌水。给人的印象不好。。。

songzy41

原帖由 zsure79 于 2008-6-19 17:23 发表

                               
登录/注册后可看大图

相位角好像是这样来的：

f1=fft(signal1); f2=fft(signal2);
N=length(f1);
s12=f1(1:N/2).*conj(f2(1:N/2);
ang=angle(s12);
就是两个信号的fft相乘，其中一个是fft的共轭复数。
然后相位角等于这个乘积的相 ...

从s12求出的是互谱的相角，但互谱的相角便是两信号之间在该频率上的相移吗？

zsure79

原帖由 songzy41 于 2008-6-20 09:06 发表

                               
登录/注册后可看大图

从s12求出的是互谱的相角，但互谱的相角便是两信号之间在该频率上的相移吗？

按道理是的。一个信号由于传播时间的差异，其相位会逐渐拉开，看看第一个图（信号图）会理解的清楚写，我不是专门做信号处理的，所以也搞得不是很清楚，我觉得有些信号这样子是可以的，但是有些信号相位会变化的剧烈，就应该不能用这种方法。