homew1.m

来自「西交大刘树棠写的《数字信号处理》那本书的所有源代码」· M 代码 · 共 10 行

Mlag=50;
n=1:Mlag;%100个采样点2M+1
Ts=0.1;
x=exp(-Ts*n);
[y,ny]=conv_m(x,n,x,n);
rx=xcorr(x,Mlag,'biased');%求自相关
subplot(211);stem(ny,y);%图1为用conv_m计算得到的x与自身的卷积
title('x与自身的卷积');xlabel('ny');ylabel('y');
subplot(212);stem((-Mlag:Mlag),rx);%图2为用xcorr计算得到的x的自相关
title('x的自相关');xlabel('m');ylabel('rx');