huxiangguan.zip_互相关 matlab_互相关MATLAB_时延估计_时延估计 matlab_相关时延估计

%%原始信号 dt=0.001;N=100; n=0:N-1;t=dt*n; f1=100;f2=200; s1=5*sin(2*pi*f1*t)+10*sin(2*pi*f2*t); s2=5*sin(2*pi*f1*t+0.5*pi)+10*sin(2*pi*f2*t+0.5*pi); figure(1); subplot(2,1,1);plot(t,s1);title('s1'); xlabel('t/ms');ylabel('Amplitude (V)'); subplot(2,1,2);plot(t,s2);title('s2'); xlabel('t/ms');ylabel('Amplitude (V)'); %%互相关分析 x=zeros(N+2*(N-1),1); x(N:2*N-1)=s1; y=s2; c=zeros(2*N-1,1); T=(-(N-1):1:(N-1)) for i=-(N-1):1:N-1 c(i+N)=(dot(x(i+N:i+2*N-1),y))/(norm(x(i+N:i+2*N-1))*norm(y)) end figure(2); plot(T,c);title('相关分析'); xlabel('时延/采样间隔');ylabel('相关系数'); %%确定最大相关系数位置 N=length(c) pmax=c(2) Delay=0 for i=2:1:N-1 if abs(c(i+1))>abs(pmax) pmax=c(i+1) Delay=i+1 else pmax=pmax Delay=Delay end end