clc;clear all;
clc;clear all;
load('a0.mat');
load('B.mat');
u=a0;
% [fuzhi,tmax0]=max(u(250000:500000));
t=(1:1000000);
NN=1000000;
t0=1; %t1me 1nterval
% for n=1:11
% u1=A(:,n);
% [fz,tmax]=max(u(250000:500000));
% tk(n)=tmax;
% end
% u=abs(hilbert(u)); % 计算包络线
t1=250000;
t2=500000;
% T=t2-t1;
T=2000;
[m1 tu]=max(u(t1:t2,1));
for n=1
u11=B(:,n);
% u11=abs(hilbert(u11)); % 计算包络线
% [m1 tu]=max(u);
% t1=tu-40; %修改于2016/7/16,t1\t2的取值影响到n(窗口数)的大小
% T=349-t1;
%
% u11=abs(hilbert(u11)); % 计算包络线
for i=1000:5000 %理论值为共平移213个点
for tt=1:NN-i
u1(tt,1)=u11(tt+i,1);
% u1(tt)=u11(tt-i);
end
[m11 tu1]=max(u1(t1:t2,1));
N=length(t1:t1+T);
rr=corrcoef(u(tu-T:tu+T,1),u1(tu1-T:tu1+T,1));
R(i)=rr(1,2);
end
R1(jj,:)=R;
[Rmax,nts]=max(R);
% s1=fft(u,2*N-1);
% s2=fft(u1,2*N-1);
% X=s1.*conj(s2);%conj(s2): 求s2的共轭, X为s1,s2的相关函数【参照复值信号得相关函数定义】
% R=real(fftshift(ifft(X))); %快速变换算法计算随机序列的功率谱密度，再作反变换计算出相关函数。
% [max1,mn]=max(R);
% pnt(n)=mn;
% Rm(n)=max1;
% s(jj-1)=(jj-1)*0.165; %距离
end
评论0