untitled2.m

来自「对数字图象进行分步处理,第一步滤波, 第二步去噪,第三步显示」· M 代码 · 共 13 行

%1.4(2)
p=0.1;N=500;Mlag=60;
u=rand(1,N);u=u-mean(u);
a=sqrt(12*p);u=u*a;
n=[0:N-1];f=1/16;
x1=sqrt(2)*sin(f*pi*n);
x=x1+u(1:N);
subplot(211);
plot(x(1:100));grid on;
rx=xcorr(x,Mlag,'biased');
subplot(212)
plot((-Mlag:Mlag),rx);grid on;