利用多光柵傅里葉頻譜的卷積實(shí)現(xiàn)數(shù)的表達(dá)和加密

胡夢(mèng)瑩， 李佳悅， 潘永華

(南京大學(xué) 物理學(xué)院，江蘇 南京210093)

0 引 言

在阿貝成像實(shí)驗(yàn)中，通過對(duì)物面上光柵進(jìn)行調(diào)制，并研究光柵傅里葉頻譜的卷積圖像，可以延伸出大量有趣而有應(yīng)用價(jià)值的拓展［1-6］。本實(shí)驗(yàn)搭建一個(gè)簡(jiǎn)單的光路系統(tǒng)，如圖1 所示。在物面上放上光柵，激光器發(fā)出的激光經(jīng)過透鏡組擴(kuò)束后垂直照在光柵上，首先從兩光柵頻譜的卷積定理出發(fā)，研究了多光柵傅里葉頻譜的卷積圖像，然后利用這種圖像特征，通過對(duì)各個(gè)光柵與水平方向夾角的調(diào)制，賦予其頻譜面上各分立譜圖像及其特定角度以某種定義，實(shí)現(xiàn)了數(shù)的表達(dá)以及加密的功能，最后探討了這種全新方案的優(yōu)越性。

圖1 實(shí)驗(yàn)的光路系統(tǒng)圖

1 多光柵傅里葉頻譜的卷積

1.1 阿貝成像原理［7，8］

在相干光照射下，顯微鏡物鏡成像過程可以看作二次衍射過程。第一次衍射發(fā)生在物平面到物鏡后焦面，通過物的衍射光在物鏡后焦面上形成夫瑯禾費(fèi)衍射圖;第二次衍射發(fā)生在物鏡后焦面到像平面，將物鏡后焦面上夫瑯禾費(fèi)衍射圖上各點(diǎn)發(fā)出的球面波在像平面上相干疊加復(fù)合為中間像，可以通過目鏡觀察到。這兩個(gè)過程本質(zhì)就是兩次傅里葉變換。其物光精確的傅里葉變換關(guān)系［9］為:

其中:fx，fy是x 和y 方向的空間頻率，g(x，y)是物平面的光場(chǎng)復(fù)振幅分布;G(fx，fy)是g(x，y)的空間頻譜。

1.2 卷積定理［10］

設(shè)有兩個(gè)二維函數(shù)g(x，y)和h(x，y)，它們的傅里葉頻譜分別為G(fx，fy)和H(fx，fy)，即

其中，F(xiàn)［］表示傅里葉變換。

根據(jù)卷積的定義(以* 表示):

則有

由此可知，若物面上復(fù)振幅分布為g(x，y)與h(x，y)之乘積，則由卷積定理可知在頻譜面上將觀察到的復(fù)振幅分布將是他們各自頻譜的卷積［11］，同樣若在物面上是g(x，y)和h(x，y)的卷積，則在頻譜面上觀察到的將是他們各自頻譜的乘積。

實(shí)驗(yàn)中分別將兩個(gè)具有不同空間頻率的一維光柵放入光路中作為物，設(shè)g(x，y)為高頻300 條/mm 一維光柵，h(x，y)為低頻100 條/mm 一維光柵，分別觀察其頻譜G(fx，fy)、H(fx，fy)，如圖2 所示。然后將兩光柵重疊在一起作為物，則總的振幅透過率為兩光柵各自振幅透過率之乘積g(x，y)·h(x，y)，在頻譜面上則可觀察到兩光柵各自傅里葉頻譜的卷積G(fx，fy)*H(fx，fy)，如圖3 所示。

圖2 兩個(gè)一維光柵各自的頻譜圖

圖3 兩個(gè)一維光柵頻譜的卷積圖像

1.3 卷積定理的推廣——多光柵傅里葉頻譜的卷積

在物理學(xué)中，任何一個(gè)符合疊加原理的線性系統(tǒng)都存在卷積［12］。一般來(lái)說，卷積是一種通過兩個(gè)函數(shù)f 和g 生成第三個(gè)函數(shù)的數(shù)學(xué)算子，表征函數(shù)f 與經(jīng)過翻轉(zhuǎn)和平移的函數(shù)g 的重疊部分的面積。如果將參加卷積的一個(gè)函數(shù)看做區(qū)間的指示函數(shù)，則卷積還可以看做是“移動(dòng)平均”的推廣。數(shù)學(xué)上，卷積運(yùn)算滿足交換律和結(jié)合律，也即可以實(shí)現(xiàn)多個(gè)函數(shù)的連續(xù)卷積運(yùn)算，所有參與卷積的函數(shù)地位相同。

在光學(xué)上，為了實(shí)現(xiàn)卷積，應(yīng)用以下傅氏變換公式［13-14］:

式中:A(νX)是a(x)的傅氏變換;F 表示傅氏變換;νX是空間頻率。故我們完全可以將1.2 中2 光柵傅里葉頻譜的卷積推廣到多個(gè)光柵［15］。在實(shí)驗(yàn)條件足夠理想時(shí)，4 個(gè)光柵、5 個(gè)光柵乃至n 個(gè)光柵傅里葉頻譜的卷積圖像都有望實(shí)現(xiàn)。

本文通過實(shí)驗(yàn)觀測(cè)到了3 個(gè)一維光柵(空間頻率分別為300 條/mm、100 條/mm 以及10 條/mm)傅里葉頻譜的卷積圖像，如圖4 所示。

圖4 3 個(gè)光柵傅里葉頻譜的卷積圖像

2 利用多光柵傅里葉頻譜卷積實(shí)現(xiàn)數(shù)的表達(dá)

2.1 基本思想

空間頻率越大的光柵，其衍射分立譜的間距越大。如果在物面上先放置一個(gè)空間頻率最大的一維光柵G1，再放置一個(gè)空間頻率第二大的一維光柵G2，此時(shí)在頻譜面上就可以得到它們各自頻譜的卷積圖像。保持G1 不動(dòng)，旋轉(zhuǎn)G2，可以觀察到G2 的分立譜繞著G1頻譜的每個(gè)點(diǎn)轉(zhuǎn)動(dòng)，即與G1 頻譜的夾角發(fā)生變化。再加入一個(gè)空間頻率第三大的一維光柵G3，此時(shí)頻譜面上即為3 個(gè)光柵傅里葉頻譜卷積圖像。保持G1、G2固定不動(dòng)，旋轉(zhuǎn)G3，則G3 的分立譜也會(huì)繞著G1、G2傅里葉頻譜卷積圖像的每個(gè)頻譜點(diǎn)轉(zhuǎn)動(dòng)。按照同樣的方式，保持先加入的光柵固定不動(dòng)，依次加入并轉(zhuǎn)動(dòng)空間頻率由大到小的一維光柵，便可觀察到同樣變化規(guī)律的現(xiàn)象。

基于上述方式調(diào)制光柵，通過定義一種簡(jiǎn)明的規(guī)則，即可實(shí)現(xiàn)數(shù)的表達(dá)。規(guī)則如下:將各個(gè)光柵的分立譜圖像與十進(jìn)制數(shù)的各個(gè)數(shù)的位權(quán)相對(duì)應(yīng)，通過依次調(diào)制各一維光柵與水平軸之間的夾角得到特定的頻譜分布圖，并認(rèn)為不同夾角表示不同的數(shù)值，再將各個(gè)數(shù)乘以各自的位權(quán)后相加，即可得到一個(gè)特定的數(shù)，從而在多光柵傅里葉頻譜卷積圖像上實(shí)現(xiàn)了數(shù)的表達(dá)。

2.2 實(shí)驗(yàn)實(shí)現(xiàn)

首先水平放置空間頻率最大的一維光柵G1(300條/mm)，頻譜如圖5 所示。再水平放置空間頻率第二大的一維光柵G2(100 條/mm)，它們頻譜的卷積圖像如圖6 所示，規(guī)定該頻譜圖像表示0 ×10°。

圖5 300 條/mm 一維光柵頻譜圖

圖6 0 ×10°的頻譜圖像

再旋轉(zhuǎn)G2，規(guī)定其分立譜與水平軸(即最大空間頻率光柵G1 分立譜位置)夾角α 為π/10、2π/10、3π/10、4π/10、5π/10、6π/10、7π/10、8π/10、9π/10 依次表示1 ×10°、2 ×10°、3 ×10°、4 ×10°、5 ×10°、6 ×10°、7 ×10°、8 ×10°、9 ×10°。此時(shí)G2 的分立頻譜所代表的位權(quán)為10°。通過兩光柵傅里葉頻譜卷積圖像，即可表示0 ～9 的所有整數(shù)。頻譜圖像如圖7 所示。

圖7 表示0 ～9 的傅里葉頻譜卷積圖像

利用多光柵傅里葉頻譜卷積圖像可以實(shí)現(xiàn)更高位的數(shù)的表達(dá)。加上空間頻率第三大的一維光柵G3(10條/mm)，在保持G1、G2 固定不動(dòng)的基礎(chǔ)上旋轉(zhuǎn)G3，規(guī)定其分立譜與水平軸夾角β 為0、π/10、2π/10、3π/10、4π/10、5π/10、6π/10、7π/10、8π/10、9π/10 依次表示0、1 ×101、2 ×101、3 ×101、4 ×101、5 ×101、6 ×101、7 ×101、8 ×101、9 ×101。此時(shí)G3 分立頻譜所代表的位權(quán)為101，3 個(gè)光柵傅里葉頻譜卷積圖像所表示的數(shù)即為β×10/π×101+α×10/π×100。所以，利用三個(gè)光柵傅里葉頻譜卷積圖像，即可表示0 到99 間所有整數(shù)。例如，如圖8 所示，α =5π/10，β =8π/10，所表示的數(shù)即為8 ×101+5 ×100=85;如圖9 所示，α =7π/10，β=2π/10，所表示的數(shù)即為2 ×101+7 ×100=27;如圖10 所示，α=7π/10，β =5π/10，所表示的數(shù)即為5 ×101+7 ×100=57。

圖8 表示85 的傅里葉頻譜卷積圖像

圖9 表示27 的傅里葉頻譜卷積圖像

圖10 表示57 的傅里葉頻譜卷積圖像

依照上述規(guī)定，依次加入并旋轉(zhuǎn)空間頻率由大到小的各個(gè)光柵，將每個(gè)角度所表示的數(shù)乘以各自的位權(quán)后相加，即可得到更高位的數(shù)。理論上，可利用n 個(gè)光柵傅里葉頻譜卷積圖像來(lái)表示0 ～10n－1－1 之間的所有整數(shù)，即利用多光柵傅里葉頻譜卷積實(shí)現(xiàn)了數(shù)的表達(dá)。

3 利用多光柵傅里葉頻譜卷積實(shí)現(xiàn)數(shù)的加密

在某種規(guī)則下，多光柵傅里葉頻譜的卷積圖像可以表示數(shù)。由于這種方式是用圖來(lái)間接地表示數(shù)，并且若不熟悉規(guī)則就不可能知道表達(dá)的數(shù)是什么，所以完全可以將其作為一種數(shù)的加密手段，而多光柵傅里葉頻譜卷積圖像就是藏有密碼信息的“加密圖”。

3.1 基本思路

加密的關(guān)鍵在于規(guī)則的制定與應(yīng)用。為使得這種基于光柵調(diào)制的加密方式具有較強(qiáng)的安全性與靈活性，制定了以下規(guī)則:

(1)用各個(gè)光柵的分立譜圖像與水平軸之間的特定夾角表示特定數(shù)值。這種規(guī)則使得加密方式靈活多變:例如，加密者可以用夾角為0π/10、1π/10、2π/10、3π/10、4π/10、5π/10、6π/10、7π/10、8π/10、9π/10 的光柵分立譜圖像依次表示數(shù)0、1、2、3、4、5、6、7、8、9，也可用來(lái)表示9、8、7、6、5、4、3、2、1、0，或用來(lái)表示0、2、4、6、8、1、3、5、7、9，等等。另外還可以用其他夾角來(lái)表示數(shù)字:如完全可以用1π/8 來(lái)表示1，等等?？傊瑠A角與數(shù)值的對(duì)應(yīng)關(guān)系完全由加密者所確定。

(2)將各個(gè)光柵的分立譜圖像與各密碼數(shù)的排列位置相對(duì)應(yīng)。在(1)中已明確了各個(gè)光柵的分立譜圖像與水平軸之間的特定夾角所表示的數(shù)，加密者可再制定規(guī)則來(lái)確定每個(gè)光柵的分立譜圖像所對(duì)應(yīng)各個(gè)數(shù)的排列順序。例如，可令空間頻率第二大的一維光柵的分立譜代表密碼的第一個(gè)數(shù)，空間頻率第三大的一維光柵的分立譜代表密碼的第二個(gè)數(shù)……，以此類推，空間頻率最小的光柵分立譜代表密碼的最后一個(gè)數(shù)。加密者便可通過第2 節(jié)中方式依次調(diào)制光柵，使得到的多光柵傅里葉頻譜卷積圖像即為某個(gè)特定密碼的“加密圖”。同樣，確定數(shù)的排列順序的規(guī)則也是靈活多變的:完全可以用空間頻率第二大到最小的光柵的分立譜倒序表示密碼，這樣同一幅圖所存儲(chǔ)的密碼的順序相反。如果再采用另外的規(guī)則來(lái)確定各密碼數(shù)的排列順序，則同一幅圖所存儲(chǔ)的密碼又會(huì)改變。由于規(guī)則的制定方式由加密者確定，這大大增加了解密難度。

(3)若密碼數(shù)較多，還可用多幅多光柵傅里葉頻譜卷積圖像實(shí)現(xiàn)加密。由于每幅卷積圖像所能存儲(chǔ)的數(shù)字信息量是有限的，可以讓每幅圖存儲(chǔ)密碼的某一部分，而完整的密碼內(nèi)容就是每幅圖像存有密碼的某種組合。加密者可自行規(guī)定每幅圖存儲(chǔ)的是密碼的哪一部分，這個(gè)不確定因素使得加密程度進(jìn)一步加深，解密難度也進(jìn)一步加大。

3.2 實(shí)驗(yàn)實(shí)現(xiàn)

(1)利用一幅3 個(gè)光柵傅里葉頻譜的卷積圖像實(shí)現(xiàn)對(duì)2 位密碼的加密

加密者采用如下規(guī)則:夾角為0π/10、1π/10、2π/10、3π/10、4π/10、5π/10、6π/10、7π/10、8π/10、9π/10的光柵分立譜圖像依次表示數(shù)0、1、2、3、4、5、6、7、8、9，并規(guī)定空間頻率第二大的一維光柵的分立譜表示密碼的第一個(gè)數(shù)，空間頻率第三大的一維光柵的分立譜表示密碼的第二個(gè)數(shù)。則若要對(duì)密碼75 實(shí)現(xiàn)加密，可用第2 節(jié)中的方式依次調(diào)制光柵G1、G2、G3，使得α =7π/10，β=5π/10，如圖11 所示，該3 光柵傅里葉頻譜卷積圖像即為存有密碼75 的“加密圖”。

圖11 存有密碼75 的“加密圖”

由于規(guī)則為加密者靈活制定，所以如果加密者變換加密規(guī)則，那么圖11 解密后將是另外一個(gè)密碼。例如，加密者用夾角為0π/10、1π/10、2π/10、3π/10、4π/10、5π/10、6π/10、7π/10、8π/10、9π/10 的光柵分立譜圖像依次表示數(shù)9、8、7、6、5、4、3、2、1、0，則圖11就變?yōu)榇嬗忻艽a24 的“加密圖”;或者加密者若不改變上段中各夾角的光柵分立譜圖像所表示的數(shù)，而是改變各個(gè)光柵的分立譜圖像與各密碼數(shù)的排列位置的對(duì)應(yīng)關(guān)系，如用空間頻率第二大的一維光柵的分立譜表示密碼的第二個(gè)數(shù)，空間頻率第三大的一維光柵的分立譜表示密碼的第一個(gè)數(shù)，則圖11 就成為存有密碼57的“加密圖”。所以，一幅多光柵傅里葉頻譜卷積圖像可實(shí)現(xiàn)對(duì)不同密碼的加密。同樣，同個(gè)密碼也可以用不同的多光柵傅里葉頻譜卷積圖像實(shí)現(xiàn)加密。由于真正存儲(chǔ)的密碼取決于加密者制定的規(guī)則，而這規(guī)則千變?nèi)f化，故成功解密的難度非常大，密碼的安全性很高。

(2)對(duì)多位密碼進(jìn)行加密。若加密者要對(duì)多位密碼進(jìn)行加密，可利用多幅多光柵傅里葉頻譜卷積圖像。

例如，加密者要對(duì)6 位密碼進(jìn)行加密，可利用3 幅3 個(gè)光柵傅里葉頻譜卷積圖像，每幅圖存有其中2 位密碼。如果對(duì)于每個(gè)2 位密碼加密方式都采用如同2.2 節(jié)中密碼75 的“加密圖”的加密規(guī)則，并令從左往右排列的3 幅圖分別存有6 位密碼的第一第二位、第三第四位和第五第六位，則若要對(duì)密碼725855 實(shí)現(xiàn)加密，可用第2 節(jié)中方式依次調(diào)制光柵，得到3 幅3 光柵傅里葉頻譜卷積圖像并依次排列，如圖12 所示，這3幅從左往右排列的三光柵傅里葉頻譜卷積圖像即為密碼725855 的“加密圖”。

圖12 存有密碼725855 的“加密圖”

事實(shí)上，加密者還可以規(guī)定從左往右排列的3 幅圖分別存有密碼的第一第六位、第二第五位、以及第三第四位，則圖12 實(shí)現(xiàn)了密碼755582 的加密。或改變讀圖順序，令從右往左排列的三幅圖分別存有密碼的第一第二位、第三第四位和第五第六位，則圖12 又成為密碼555872 的“加密圖”。僅僅是改變每幅圖所存有的密碼的排位或讀圖的順序，3 幅3 個(gè)光柵傅里葉頻譜卷積圖像就可實(shí)現(xiàn)多個(gè)不同的六位密碼的加密，若再改變每一幅圖的加密規(guī)則，則利用這3 幅3 個(gè)光柵傅里葉卷積圖像可實(shí)現(xiàn)非常多種不同的六位密碼的加密。同時(shí)，同一個(gè)密碼還可以用不同的多幅多光柵傅里葉頻譜卷積圖像實(shí)現(xiàn)加密。所以，利用多幅多光柵傅里葉頻譜的卷積圖像實(shí)現(xiàn)加密的安全性比利用一幅要高得多，破解密碼難度極大。

3.3 用此種方式實(shí)現(xiàn)加密功能的優(yōu)越性

本實(shí)驗(yàn)利用多光柵傅里葉頻譜卷積圖像實(shí)現(xiàn)加密，其優(yōu)越性體現(xiàn)在以下幾點(diǎn):

(1)通過各光柵頻譜圖與水平方向夾角來(lái)存儲(chǔ)數(shù)字，在實(shí)驗(yàn)上易于實(shí)現(xiàn)，且現(xiàn)象清楚明顯，既易于調(diào)制，也易于讀取。

(2)基于此種光柵調(diào)制方式，可定義各種規(guī)則進(jìn)行加密。規(guī)則的任意性、復(fù)雜性和可靈活變換性，大大地增強(qiáng)了加密強(qiáng)度，也大大提高了密碼信息存儲(chǔ)的安全性。加密者可任意改變規(guī)則，只有在特定時(shí)間知道該時(shí)間段特定規(guī)則的人，才可以成功解密。

4 結(jié) 語(yǔ)

本文從卷積定理的推廣，在實(shí)驗(yàn)上研究并驗(yàn)證了多光柵傅里葉頻譜的卷積圖像。利用多光柵傅里葉頻譜的卷積圖像的特性，通過依次調(diào)制物面上各個(gè)光柵與水平方向的角度，并賦予其頻譜面上各光柵分立譜圖像以及分立譜圖像的特定角度以某種定義，創(chuàng)造出了一種全新的多位數(shù)的表達(dá)方式。同時(shí)還提出，利用這種調(diào)制方式，加上靈活多變的定義規(guī)則，可以實(shí)現(xiàn)數(shù)字密碼的加密，而該種加密方式安全可靠，其優(yōu)越性是顯而易見的。

［1］ 余 欣，潘永華. 光的三進(jìn)制邏輯運(yùn)算與光學(xué)加密［J］. 物理實(shí)驗(yàn)，2008，28(3):1-8.

［2］ 王漢琛，劉 璟，潘永華，周 進(jìn).運(yùn)用阿貝成像原理進(jìn)行平面密鋪結(jié)構(gòu)的頻譜分析［J］.2011，31(10):5-10.

［3］ 李芳菊，董康軍.利用阿貝成像原理制作低頻全息光柵［J］.物理實(shí)驗(yàn)，2008，28(5):1-3.

［4］ 馮少彤.互補(bǔ)分形圖像的光學(xué)傅里葉變換［J］. 激光雜志，2003，24(1):27-28.

［5］ 王 翚，劉香茹，石發(fā)旺. 利用阿貝成像原理制作全息光柵的理論分析［J］.河南科技大學(xué)學(xué)報(bào)(自然科學(xué)版)，2006(2):1-3.

［6］ 袁 霞，王晶晶，金華陽(yáng). 阿貝成像原理和空間濾波實(shí)驗(yàn)的改進(jìn)［J］.物理實(shí)驗(yàn)，2010，30(3):1-3.

［7］ 高文琦等. 光學(xué)(第三版)［M］. 南京:南京大學(xué)出版社，2013:250-251.

［8］ 趙凱華，鐘錫華.光學(xué)(下冊(cè))［M］.北京:北京大學(xué)出版社，1984:75-77.

［9］ 梁昆淼，劉 法，繆國(guó)慶.數(shù)學(xué)物理方法［M］.3 版.北京:高等教育出版社，1998:102.

［10］ 蘇為寧，王思慧，高文莉，等編. 大學(xué)物理實(shí)驗(yàn)(理科)第二冊(cè)［M］.南京:南京大學(xué)出版社，2014:130.

［11］ 潘元?jiǎng)?，馮璧華，于 瑤.大學(xué)物理實(shí)驗(yàn)(第二冊(cè))［M］.南京:南京大學(xué)出版社，1998:102.

［12］ 毛惠良，多重卷積公式及其應(yīng)用［J］. 工科數(shù)學(xué)，1995，11(3):173-176.

［13］ Jack. D. Gaskill. Linear systems，F(xiàn)ourier transforms，and optics［M］. Arizona:John Wiley ＆ Sons，Inc，1978，160-201.

［14］ 馮光彤，金國(guó)鈞，章程軍，等. 光學(xué)分形圖像的產(chǎn)生［J］. 中國(guó)激光，1997，24(8):2.

［15］ 高 峰.傅里葉變換透鏡參數(shù)的研究［J］. 實(shí)驗(yàn)室研究與探索，1997(6):59-60.