zixiangguan2.m

来自「用MATLAB实现的自相关的算法源程序」· M 代码 · 共 40 行

%auto correlation about book by huguangshu(sin function with Gassian
%disturbance

N=500;p1=1;p2=0.1;f=1/8;Mlag=60;
u=randn(1,N),u2=u*sqrt(p2);n=[0:N-1];
u_mean=mean(u);
power_u=var(u);
figure(3);
%rxu=xcorr(u,Mlag,'biased');
rxu=xcorr(u,Mlag);
plot(rxu);
s=sin(2*pi*f*n);
figure(1);
subplot(3,1,1)
stem(n,s);
axis([0 40 -1 +1]);
x1=u(1:N)+s;
rx1=xcorr(x1,Mlag,'biased');
%rx1=xcorr(x1);
subplot(3,1,2);plot(rx1);
axis([0 240 -2 +2]);
u2=u*sqrt(p2);x2=u2(1:N)+s;
rx2=xcorr(x2,Mlag,'biased');
%rx2=xcorr(x2)
subplot(3,1,3),plot(rx2);
axis([0 240 -2 +2])

figure(2);
subplot(3,1,1);
%stem(n,s);
plot(s);
x1=u(1:N)+s;
rx1=xcorr(x1,Mlag);
%rx1=xcorr(x1,Mlag,'biased');
subplot(3,1,2);plot(rx1);
%axis([0 40 -1 +1]);
u2=u*sqrt(p2);x2=u2(1:N)+s;
rx2=xcorr(x2,Mlag,'biased')
subplot(3,1,3),plot(rx2);
axis([0 40 -1 +1])