序列卷积与相关运算（含matlab演示代码）

---

0 前言

  近日学习“数字信号处理”课程，上课时讲到序列卷积与相关运算，一时半会想不起来，遂趁周末复习了一下信号与系统中相关的内容，将复习所得整理如下。

1 序列卷积运算

1.1 卷积的定义

  序列的卷积运算揭示了 离散信号作用于离散线性非时变系统的机理 。两个离散信号的卷积定义为:
$$y[k]=x_1[k]\ast x_2[k]=\sum _{n=-\infty }^{+\infty }{x}_1[n]x_2[k-n]\text{\hspace{1em}(1.1)}$$

1.2 卷积的性质

1. 交换律
   $$x_1[k]\ast x_2[k]=x_2[k]\ast x_1[k]$$

这里只对交换律进行证明，其他几条性质的证明是类似的
证：由式$(1.1)$可得
$$y[k]=x_1[k]\ast x_2[k]=\sum _{n=-\infty }^{+\infty }{x}_1[n]x_2[k-n]\text{\hspace{1em}(1)}$$
换元，令$t=k-n$,式$(1)$变为
$$\sum _{t=+\infty }^{-\infty }{x}_1[k-t]x_2[t]\text{\hspace{1em}(2)}$$
对于 累加 运算来说，将$+\infty$与$-\infty$颠倒一下是不影响最终结果的，因此可以将式$(2)$写作
$$\sum _{t=-\infty }^{+\infty }{x}_1[k-t]x_2[t]\text{\hspace{1em}(3)}$$
再使用$n$替代$t$（只是换个符号表示而已，不影响式子结果）
$$\sum _{n=-\infty }^{+\infty }{x}_1[k-n]x_2[n]=x_2[k]\ast x_1[k]\text{\hspace{1em}(4)}$$
证毕。

2. 结合律
3. 分配律
4. 位移特性

2 序列相关运算

  为了方便记忆，一般需要将 相关运算 与 卷积运算 联系起来。

2.1 相关的定义

两个能量有限的实序列$x[k]$与$y[k]$的互相关运算定义为
$$r_{xy}[n]=\sum _{k=-\infty }^{+\infty }x[k]y[k+n]\text{\hspace{1em}(2.1)}$$

一眼看过去是不是和卷积运算$(1.1)$很像？事实上它们确实有联系（数学推导上的），这种相似性可以帮助我们进行记忆。
以下就是通过数学推导来辅助记忆的过程
首先对于$(2.1)$进行换元，令$t=k+n$，于是式$(2.1)$变为
$$r_{xy}[n]=\sum _{t=-\infty }^{+\infty }x[t-n]y[t]\text{\hspace{1em}(2.2)}$$
再使用$n$替代$t$，于是式$(2.2)$变为
$$r_{xy}[n]=\sum _{k=-\infty }^{+\infty }x[k-n]y[k]\text{\hspace{1em}(2.3)}$$
是不是已经与卷积运算非常相似了？对卷积运算比较熟悉的应该已经能看出结论了，但这里还是多讲一步方便理解。
接下来的步骤请务必注意$\sum$的下标是$k$还是$n$！！！
我们将$x[-n]$和$y[n]$视作进行卷积运算的两个序列，根据式$(1.1)$，$y[n]\ast x[-n]$定义为
$$y[n]\ast x[-n]=\sum _{k=-\infty }^{+\infty }y[k]x[-(-k)-n]\text{\hspace{1em}(2.4)}$$
化简得
$$y[n]\ast x[-n]=\sum _{k=-\infty }^{+\infty }y[k]x[k-n]\text{\hspace{1em}(2.5)}$$
可以看到式$(2.3)$与式$(2.5)$是完全一致的，即
$$r_{xy}[n]=x[-n]\ast y[n]$$

2.2 相关运算的作用

  相关运算的用途有很多，这里只展示最常见且最有用的两种

2.2.1 检测序列y[k]中是否含有序列x[k]

  验证该用途所使用的matlab代码如下所示，以Akie秋绘的翻唱歌曲《なんでもないや》为例进行演示（鄙人是老绘星了）。

% 用于演示自相关运算的作用
clc;
clear;
close all;

% 1st_判断输出y[k]中是否含有x[k]
fs = 44100;  % 采样频率
t = 5;  % 采样时间
dt = 0:1/fs:t-1/fs;  % 采样时刻
samples = [1, t*fs];
[src,fs] = audioread('music_akie.mp3', samples);
x = src(:, 1);  % 取单声道数据
% sound(x, fs);
subplot(321);
plot(dt, x);
title('(a)$x[k]$', 'interpreter', 'latex');


% 计算有用信号自相关
[corr_x, lags] = xcorr(x);
subplot(322);
plot(lags, corr_x);
ylim([-2000, 3000]);
title('(b)$r_{x}[n]$', 'interpreter', 'latex');

% 创建噪声信号
noise = normrnd(0, 0.1, [t*fs, 1]);
% sound(noise, fs);
subplot(323);
plot(dt, noise);
title('(c)noise $n$', 'interpreter', 'latex');

% 计算噪声信号自相关
[corr_n, lags] = xcorr(noise);
subplot(324);
plot(lags, corr_n);
title('(d)$r_{n}[n]$', 'interpreter', 'latex');

% 混合音乐信号与噪声信号
y = x + noise;
% sound(y, fs);  % 可以听一下混合后的信号，可以听到歌声，但是从是时域图上看不出来是否还含有歌声信号
subplot(325);
plot(dt, y);
title('(e)$y = x+n$', 'interpreter', 'latex');

% 计算噪声污染后音乐信号自相关
[corr_y, lags] = xcorr(y);
subplot(326);
plot(lags, corr_y);
ylim([-2000, 3000]);
title('(f)$r_{y}[n]$', 'interpreter', 'latex');

(a)音乐信号(b)音乐信号的自相关函数(c)噪声信号(d)噪声信号的自相关函数
(e)受到噪声污染的音乐信号(f)被污染音乐信号的自相关函数

  可以看到，音乐信号本身的自相关函数（图b）与受到污染音乐信号的自相关函数非常相似，其原因在于 噪声信号的自相关函数除了零点处均约等于0 如此一来我们就可以通过相似度计算判断受到污染的信号中是否含有原信号。

  具体的原理在此不再推导，根据“2.1 相关的定义”中的内容按部就班计算即可，如需另外讲解，可以在评论区说一下，鄙人看到会来补充说明的！

  值得一提的是，当我们增大噪声的幅度，就会发现受到污染音乐信号的自相关函数不再像原来的自相关函数，而是更像噪声信号的自相关函数了，换言之，个人认为，相关运算检测的是混合信号中的主要成分 ！！！接下来我以Akie翻唱的另一首歌《lemon》为例演示，结果如下图所示。

(a)音乐信号1(b)音乐信号1的自相关函数(c)音乐信号2(d)音乐信号2的自相关函数(e)混合的音乐信号(f)混合音乐信号的自相关函数

2.2.2 回声消除

  首先需要对回声信号进行建模，根据生活常识，回声信号$echo\_x[k]$就是 原信号$x[k]$加上经过一段时间后返回的有衰减的信号$\lambda x[k-N]$ 。其数学表达式如下
$$echo\_x[k]=x[k]+\lambda x[k-N]$$
  然后我们进行理论的推理，先求回声信号的自相关函数$r_{echo\_x}[n]$
$$\begin{gathered} r_{echo\_x}[n]=\sum _{k=-\infty }^{+\infty }echo\_x[k]echo\_x[k+n] \\ =\sum _{k=-\infty }^{+\infty }(x[k]+\lambda x[k-N])(x[k+n]+\lambda x[k+n-N]) \\ =\sum _{k=-\infty }^{+\infty }(x[k]x[k+n]+\lambda x[k]x[k+n-N]+\lambda x[k-N]x[k+n]+{\lambda }^{2}x[k+n-N]x[k-N]) \\ =(1+{\lambda }^2)r_x[n]+\lambda r_x[n-N]+\lambda r_x[n+N] \end{gathered}$$
  可以看到，回声信号的自相关函数$r_{echo\_x}[n]$其实是原始信号自相关函数$r_x[n]$左右平移得到。