概述

卷積是一種線性運算,其本質(zhì)是滑動平均思想,廣泛應(yīng)用于圖像濾波。而隨著人工智能及深度學(xué)習(xí)的發(fā)展，卷積也在神經(jīng)網(wǎng)絡(luò)中發(fā)揮重要的作用，如卷積神經(jīng)網(wǎng)絡(luò)。本參考設(shè)計主要介紹如何基于INTEL 硬浮點的DSP Block實現(xiàn)32位單精度浮點的卷積運算，而針對定點及低精度的浮點運算，則需要對硬浮點DSP Block進(jìn)行相應(yīng)的替換即可。

原理分析

設(shè):f(x), g(x)是兩個可積函數(shù)，作積分：

隨著x的不同取值，該積分定義了一個新的函數(shù)h(x)，稱為函數(shù)f(x)與g(x)的卷積，記為h(x)=f(x)*g(x)。

如果卷積的變量是序列x(n)和h(n)，則卷積的結(jié)果為

其中*表示卷積。因此兩個序列的卷積，實際上就是多項式的乘法，用個例子說明其工作原理。a = [7,5,4]; b = [6,7,9];則實現(xiàn)a和b的卷積，就是把a和b作為一個多項式的系數(shù)，按多項式的升冪或降冪排列，即為：

因此得到a*b=[42,79,122,73,36]；與Matlab運算結(jié)果一致。而二維卷積可以采用通用多項式乘積方法實現(xiàn)卷積運算。

基于INTEL FPGA的實現(xiàn)分析

如上我們確定了兩個序列的卷積等同于兩個多項式的乘法，因此當(dāng)我們需要計算序列[a0,a1,a2, …,an-1]與[b0,b1,b2, …,bn-1]的卷積結(jié)果時，可以成立a,b兩個n階多項式，如下所示：

則[a0,a1,a2, …,an-1]與[b0,b1,b2, …,bn-1]的卷積結(jié)果即為由a*b得到的多項式的各項系數(shù)所組成的序列。令c=a*b，得到

則由多項式c的各階系數(shù)所組成的新的序列[c0,c1,c2, …,c2n-1]即為[a0,a1,a2, …,an-1]與[b0,b1,b2, …,bn-1]的卷積結(jié)果。則按照高階多項式計算展開可得到：

┆┆

┆┆

因此卷積的運算可以轉(zhuǎn)化為行向量與列向量相乘的結(jié)果，即乘累加的運算結(jié)構(gòu)。

Intel FPGA在Arria10DSP Block中首次支持了單精度硬浮點DSP block,是行業(yè)內(nèi)第一個支持單精度DSP block，硬浮點DSP block架構(gòu)如圖1所示：

圖1 硬浮點DSPblock架構(gòu)

硬浮點DSP Block包含硬浮點乘法器，硬浮點加法器，支持乘累加運算，因此采用硬浮點DSPblock實現(xiàn)行列向量相乘是非常好的方式。下面我們針對一個實際的卷積運算，介紹如何基于INTEL硬浮點DSP block實現(xiàn)。假設(shè)我們需要求隨機數(shù)組a=[4,8,9,11]與b=[10,5,7,13]的卷積運算結(jié)果，則根據(jù)上面的分析，保持?jǐn)?shù)組a順序不變，而數(shù)組b需根據(jù)上述分析結(jié)果，針對每一個卷積結(jié)果產(chǎn)生新的序列。所以整個實現(xiàn)包括數(shù)列重組模塊和硬浮點乘法器模塊及輸出處理。下面是實現(xiàn)框圖及仿真結(jié)果。

圖2 實現(xiàn)框圖

圖3 Modelsim仿真結(jié)果

仿真結(jié)果與Matlab實現(xiàn)結(jié)果一致，并且該設(shè)計中充分考慮了FPGA并行擴展特性，對于低速率要求的設(shè)計可采用DSP Block復(fù)用的方式節(jié)約DSP block數(shù)量。