求大神解释这段matlab程序,

%就是对同一源信号想x1作延迟x2并加上噪声,然后用两种方法分析得到x1和x2的%互相关函数!并对加噪声的x1 分析了自相关函数clear alln=99; T=6; t=0:2*pi/(n-1):2*T*pi; N=length(t);fs=1/(n-1);x1=4*sin(t); %源信号 y1=4*sin(t+pi/6); %%源信号延迟pi/6相位 x=awgn(x1,20); %加高斯白噪声,信噪比SNR==20dB y=awgn(y1,0。001); %同上x2=awgn(x1。/4,-3); subplot(311),plot(t,x)subplot(312),plot(t,y)tic %tic和toc一起用于计算程序运算时间[Cc,r]=xcorr(x,y); %xcorr用于求互相关矩阵函数,Cc为互相关矩阵,r为延迟矢量 [Cc1,r1]=xcorr(x1,y1); %同上subplot(313),plot(r*2*pi/(n-1),Cc,'r'),hold on;plot(Cc1,'-');tic S1=fft(x,2*N-1); %求x的2*N-1点fftS2=fft(y,2*N-1);X=S1。* conj(S2); %求x和y的互功率谱R1=real(fftshift(ifft(X)));%利用ifft求反变换得到互相关函数,用fftshift反转toc %对应tic求计算时间figure(2)subplot(211),plot(Cc),title('xcorr')subplot(212),plot(R1);title('ifft')figure(3)plot(xcorr(x2,'unbiased')); %求x2的自相关函数并作图