xiangguan_qxy.m

来自「相关分析例程」· M 代码 · 共 32 行

% 信号来源文件：rough乌鲁木齐首站出站a；rough鄯善站进站a
clear
close all
%x1:乌鲁木齐首站出站； y1:鄯善站进站
load xinpu_lanzhou_11172005
x1=x1-mean(x1);
y1=y1-mean(y1);

 N=length(x1);
n=0:N-1;
 Fs=50;
 t=n/Fs;
 Lag=100;
 %x1=[x]';
 
%[c,lags]=xcorr(x1,x1,Lag,'unbiased');
[c,lags]=xcorr(x1,y1,'unbiased');
 %c=xcorr(x1,y1);
 subplot(121)
 plot(t,x1,'r',t,y1,'b');title('信号序列起始时间：xinpu_lanzhou_11172005');% 注意：　文件中数据起始时间：6年11月8日31时36分41秒
 xlabel('相对时间/秒');
 %legend('x1:西靖出站','y1:新堡进站');
 axis tight;grid on 
 subplot(122);
 plot(lags/Fs,c,'r');axis tight,title('x1与y1的互相关函数'),xlabel('相对时间/秒');
 grid
 [maxc,lo]=max(c);
 lo=(lo-N)/50
 ds=lo*1628
 vs=148.8*1000/lo