correlation1.rar_信号互相关

% 本程序通过两个正弦信号验证互相关系数 n=0:100; x=sin((pi/25)*n); %产生一个周期为50的远端正弦信号 subplot(2,3,1); plot(n,x);grid; title('参考信号'); xjin=sin((pi/18)*n); %产生一近端端正弦信号 subplot(2,3,1); plot(n,xjin);grid; title('参考信号'); y=0.5*x; %将正弦信号幅度降低以此模仿回声信号 subplot(2,3,2); plot(n,y);grid; title('纯的回声信号'); z=0.3*sin((pi/15)*n)+0.2*sin((pi/46)*n); %生成另外一个不同频率，不同幅值的正弦波 subplot(2,3,3); plot(n,z);grid; title('加噪声的回声信号'); % xn=x'; % dn=y'; % M=20; % mu=0.02; % [yn,W,en]=LMS_Algorithm(xn,dn,M,mu); %调用自适应滤波器函数 % subplot(2,3,5); % plot(n,yn);grid; % title('纯回声信号拟合输出'); % % %%%%%%%%%%%%%以下程序求参考信号与滤波器输出纯回声信号的互相关系数%%%%%%%%%%%%%%%%%%% % for i=20:90 % x1=x(i:i+10); % y1=y(i:i+10); % a=corrcoef(x1',y1'); % r(i)=a(1,2); % end % j=1:length(r); % subplot(2,3,4); % plot(j,r); xn=x'; dn=z'; M=30; %mu=0.01; [yn,W,en]=functionLMS(xn,dn,M,mu); %调用自适应滤波器函数 subplot(2,3,6); plot(n,yn);grid; title('实际回声信号拟合输出'); %%%%%%%%%%%%%以下程序求参考信号与滤波器输出纯回声信号的互相关系数%%%%%%%%%%%%%%%%%%% for i=20:90 x1=xjin(i:i+10); y1=yn(i:i+10); a=corrcoef(x1',y1); r(i)=a(1,2); end j=1:length(r); subplot(2,3,4); plot(j,r); title('互相关系数'); %%%%%以下程序用两向量夹角法%%%% t1=norm(y*y'); t2=sqrt(norm(y*y'+(y-yn)^2+xjin*xjin'))*sqr(norm(y*y')); p=t1/t2