3-2.m

来自「经典数字信号处理滤波器的源代码 重点是利用巴特沃斯模拟滤波器转而设计其它数字滤」· M 代码 · 共 16 行

%例程3-2  计算随机信号的相关函数

x=[1 3 2 6 2 1 -2 0 1 5 3 2 -3 0 1 2 0 3 1];
y=[0 0 0 1 3 2 6 2 1 -2 0 1 5 3 2 -3 0 1 2 0 3 1];
k=length(y);
e=randn(1,k);
y=y+e;
xk=fft(x,2*k);
yk=fft(y,2*k);
rm=real(ifft(conj(xk).*yk));
rm=[rm(k+2:2*k) rm(1:k)];
m=(-k+1):(k-1);
stem(m,rm) 
xlabel('m');
ylabel('幅度');