正弦波与其他信号自相关函数的matlab实现

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('正弦和三角的互相关函数');