关于时间域转频率域并画出相位差

bonju13

　　请问我有2笔资料分别是aa.txt跟bb.txt,这2笔数据是由信号仪所得到的同一时间,放置于不同地点的2个接收器的时间域信号,请问小弟的程序这样编写对吗?要做2个接收器的相位差是不是把aa.txt与bb.txt的相角相减即可?烦请各位大大不吝指教~十分感激! aa.TXT (31 KB, 下载次数: 5)

2009-7-16 23:31 上传

点击文件名下载附件

bb.TXT (31 KB, 下载次数: 4)

2009-7-16 23:31 上传

点击文件名下载附件

　　load aa.txt -ascii
　　N=1024;
　　plot(aa(:,1),aa(:,2))
　　Y=fft(aa(:,2),N)/length(aa(:,2));
　　Fs=1/(aa(2,1)-aa(1,1));
　　t=(0:N-1)/Fs;
　　x=Fs/2*linspace(0,1,N/2);
　　subplot(2,1,1);
　　plot(x,2*abs(Y(1:N/2))) ;
　　p=unwrap(angle(Y));
　　f = (0:length(Y)-1)'/length(Y)*Fs;
　　subplot(2,1,2);
　　plot(f,p)
　　[ 本帖最后由 bonju13 于 2009-7-16 23:32 编辑 ]