关于时间域转频率域并画出相位差
本帖最后由 wdhd 于 2016-9-19 11:02 编辑请问我有2笔资料分别是aa.txt跟bb.txt,这2笔数据是由信号仪所得到的同一时间,放置于不同地点的2个接收器的时间域信号,请问小弟的程序这样编写对吗?要做2个接收器的相位差是不是把aa.txt与bb.txt的相角相减即可?烦请各位大大不吝指教~十分感激!
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 编辑 ]
页:
[1]