• 在图像处理领域，Kernel = convolution matrix = mask，它们一般都为一个较小的矩阵；
• 用于：Sharpen，Blur， Edge enhance，Edge detect，Emboss（使凸出；在……上作浮雕图案；装饰）

1. 卷积操作与卷积矩阵的等价性

• （1）创建一维信号

  N = 100;
  
  s = zeros(N, 1);
  
  k = [20, 45, 70];
  
  a = [2, -1, 1];
  
  s(k) = a;

• （2）创建卷积核

  L = 4;
  
  h = ones(L, 1)/L;

• （3）使用 matlab 创建 convolution matrix

  首先考虑，卷积矩阵的 size，卷积操作conv(h, s) 返回的响应其长度为 L+N-1, 因此：

  M = N+L-1;
  
  H = sparse(M, N);
  
  e = ones(N, 1)
  
  for i = 0:L-1,
  
      H = H + spdiags(e*h(i+1), -i, M, N);
  
  end

• （4）验证二者操作的等价性：

  err = H*s - conv(h, s);
  
  max_err = max(abs(err));

2. spy：可视化稀疏矩阵的形式

而不是简单地将稀疏矩阵以一个矩阵的形式显示，显然是因为其规模十分之大，且较为稀疏，不易直观观察。

Kernel (image processing)

8.2. Convolution Matrix