求教 相位差的问题

landzenki

有一组SASW的数据，表面波分析， 2个接收器接收到的同一信号。

分别是时间域信号 x1 y1

做fft  得出频率域信号 X1 Y1
X1=fft(x1)
Y1=fft(y1)

用angle函数求出相位
phase1=angle(X1）
phase2=angle(Y1)

相位差是不是 p=phase2-phase1 呢？

参考文献中利用相速度和波长的关系分析结构内部情况 如下图



如果用相位差求出波从第一个接收器到第二个接收器的时间延迟
t=p*f / 360

然后2个接收器剪距离已知，通过时间延迟求出相速度，
Vr=D/t （D是2个信号接收器之间的距离）

最后求出不同频率的波长
W=Vr / f

但是结果波长有负值，检查发现相位差有负值

自己算出的图是这样的

  
求助~问题到底出在哪里了？
matlab code 和数据文件 附件中附上
code 应该如何修改呢？
怎样去掉下载附件消耗的体力呢，请大家帮忙，当然要0消耗下载了……本人新手……