n=0:499;

x=sin(0.1*n*pi);%产生正弦波

s=-2+4*rand(1,length(x));%产生随机噪声

X=x+s;

[b,a]=xcorr(x);

[d,e]=xcorr(x,s,'unbiased');

[m,l]=xcorr(x,X,'unbiased');

subplot(621);plot(n,x);xlabel('正弦信号');

subplot(622);plot(n,s);xlabel('噪声信号');

subplot(623);plot(n,X);xlabel('正弦和噪声合成的信号');

subplot(624);plot(a,b);xlabel('正弦信号的自相关函数');

subplot(625);plot(e,d);xlabel('正弦和噪声的互相关函数');

subplot(626);plot(l,m);xlabel('正弦波与自身加噪声的互相关函数');

subplot(627);plot(Lags,c);xlabel('合成信号的自相关函数');

fb=square(0.1*pi*n);%产生方波

subplot(6,2,9);plot(n,fb);xlabel('方波信号');

[o,p]=xcorr(x,fb,'unbiased');

subplot(6,2,10);plot(p,o);xlabel('正弦和方波的互相关函数');

sgb=sawtooth(0.1*pi*n,0.5);%产生三角波

subplot(6,2,11);plot(n,sgb);xlabel('三角波信号');

[q,r]=xcorr(x,sgb,'unbiased');

subplot(6,2,12);plot(r,q);xlabel('正弦和三角的互相关函数');