王志鵬，王紅娟，蔣華龍，秦 怡

（南陽師范學(xué)院 物理與電子工程學(xué)院，河南 南陽473061）

引言

在如今社會上，信息安全問題越來越重要。由于具有高并行性和高加密維度等優(yōu)點(diǎn)，基于光學(xué)方法的加密技術(shù)得到了越來越多的關(guān)注和研究。近年來，國內(nèi)外專家學(xué)者提出了很多種不同的光學(xué)加密技術(shù)［1－13］。其中，最具代表性的光學(xué)加密技術(shù)就是由Refregier和Javidi于1995年提出的雙隨機(jī)相位光學(xué)加密技術(shù)（double random phase encoding，DRPE）［1］。雙隨機(jī)相位光學(xué)加密技術(shù)利用置于光學(xué)4-f系統(tǒng)的輸入面和頻譜面的2個隨機(jī)相位掩膜實現(xiàn)光學(xué)加密，在輸出平面上得到的復(fù)振幅滿足平穩(wěn)隨機(jī)白噪聲分布。加密結(jié)果包含振幅和相位2個分量。解密過程是加密過程的逆過程，利用加密得到的復(fù)振幅分布以及解密密鑰可以恢復(fù)原始圖像。該技術(shù)要求光學(xué)器件精確對準(zhǔn)，加密結(jié)果為復(fù)振幅分布，不便于記錄和保存，并且在解密時需要利用加密密鑰的復(fù)共軛。雖然DRPE系統(tǒng)能夠很好地抵御暴力攻擊，但是由于其本質(zhì)上是線性系統(tǒng)，該系統(tǒng)易受到“選擇密文 攻 擊 ”［14］、“已 知 明 文 攻 擊 ”［15］、“選 擇 明 文 攻擊”［16］和“唯密文攻擊”［17］等幾種形式的攻擊，安全性存在很大的隱患。為了改進(jìn)傳統(tǒng)的DRPE系統(tǒng)，李榕等人提出了一種基于雙隨機(jī)相位的圖像實值編碼方法［18］。該方法取加密結(jié)果的實部作為編碼圖像，但是在頻域需要將編碼范圍擴(kuò)大4倍。最近，Yishi Shi等人提出了一種基于ptychography的光學(xué)加密技術(shù)［10］。該方法將ptychography技術(shù)和光學(xué)加密技術(shù)相結(jié)合，加密結(jié)果是原始圖像的子區(qū)域經(jīng)過DRPE加密系統(tǒng)以后，在輸出平面上得到的多幅強(qiáng)度分布圖像，在解密時利用迭代算法可以恢復(fù)出原始圖像。這種新的光學(xué)加密方法只需要記錄強(qiáng)度分布，簡化了光學(xué)系統(tǒng)結(jié)構(gòu)，并且可以實現(xiàn)復(fù)值圖像的加密，但是在加密時需要多次曝光，記錄多個強(qiáng)度分布。

相位恢復(fù)算法（phase retrieval algorithm，PRA）［19－20］是一種由光場的強(qiáng)度分布通過迭代或搜索找到最符合真實場強(qiáng)數(shù)據(jù)的相位分布的方法，在天文學(xué)、衍射光學(xué)、全息成像、電子顯微鏡學(xué)以及逆源問題等方面都有著重要的應(yīng)用。已經(jīng)有不 少文獻(xiàn)［6－7，21－26］研究了 相 位 恢 復(fù) 算 法 在 光 學(xué) 加 密中的應(yīng)用。在這些文獻(xiàn)中，都是在加密端利用相位恢復(fù)算法生成一個或多個隨機(jī)相位板（phase only mask，POM）。因此這些方法的加密過程一般都采用數(shù)字方式。在解密端，利用生成的隨機(jī)相位板和密文，可以采用光學(xué)方法解密原始圖像。

本文在傳統(tǒng)的雙隨機(jī)相位光學(xué)加密系統(tǒng)的基礎(chǔ)上，提出了一種新的光學(xué)加密技術(shù)。在加密端，將原始圖像置于DRPE系統(tǒng)的輸入平面上，同時在輸出平面上利用CCD等感光器件記錄光強(qiáng)分布，作為密文輸出。在解密端，利用相位恢復(fù)迭代算法就可以由密文和解密密鑰恢復(fù)出所加密的原始信息。由于加密時僅需要記錄光強(qiáng)分布，簡化了系統(tǒng)結(jié)構(gòu)，易于光學(xué)實現(xiàn)。經(jīng)過分析也可以發(fā)現(xiàn)，該光學(xué)加密技術(shù)可以很好地抵抗各種攻擊，安全性較高。需要說明的是，雖然在本文中是以DRPE為基礎(chǔ)提出的一種單強(qiáng)度記錄光學(xué)加密技術(shù)，但是實際上該技術(shù)同樣可以應(yīng)用于“基于分?jǐn)?shù)傅里葉變換的加密技術(shù)”和“基于菲涅爾域的加密技術(shù)”等光學(xué)加密系統(tǒng)中，從而可以進(jìn)一步提高系統(tǒng)的安全性。

1 算法描述

1.1 加密過程

該光學(xué)加密過程采用4-f光學(xué)系統(tǒng)。將原始圖像置于4-f系統(tǒng)的輸入平面上進(jìn)行雙隨機(jī)相位加密，在輸出平面上記錄光強(qiáng)分布作為加密輸出。加密系統(tǒng)原理圖如圖1所示。

圖1 光學(xué)加密系統(tǒng)原理圖Fig.1 Schematics of optical encryption system

圖中f（x，y）是待加密的原始圖像，f（x，y）與相位掩膜exp［i2πn（x，y）］相乘后，經(jīng)過透鏡L1的傅里葉變換，然后在頻率域經(jīng)過另一相位掩膜exp［i2πb（u，v）］的濾波，再經(jīng)過透鏡L2的傅里葉逆變換，在輸出平面上利用CCD等強(qiáng)度感應(yīng)器件可以接收到光強(qiáng)分布I（x，y），I（x，y）也就是加密得到的密文。它滿足平穩(wěn)隨機(jī)白噪聲分布。其中，n（x，y）和b（u，v）是2個相互獨(dú)立的分布在0～1之間的隨機(jī)函數(shù)。加密過程可以用數(shù)學(xué)公式表示為

傳統(tǒng)的DPRE系統(tǒng)的密文是復(fù)數(shù)場，記錄時候需要使用干涉裝置，而干涉裝置對于系統(tǒng)的穩(wěn)定性很高，因而光學(xué)實現(xiàn)比較困難。相比于傳統(tǒng)的DRPE系統(tǒng)，該光學(xué)加密過程僅需記錄光強(qiáng)分布作為密文，便于光學(xué)實現(xiàn)，可以直接利用CCD等感光器件進(jìn)行記錄。同時由于CCD等感光器件輸出的是數(shù)字信號，因而也可以非常容易地進(jìn)行存儲和傳輸。

1.2 解密過程

傳統(tǒng)的DRPE系統(tǒng)解密過程是加密過程的逆過程，利用與加密密鑰互為共軛的解密密鑰就可以根據(jù)密文的復(fù)振幅分布恢復(fù)原始圖像。但是，在我們所提出的改進(jìn)的雙隨機(jī)相位加密系統(tǒng)中，在加密時僅保留強(qiáng)度分布作為加密密文，不包含相位分布，密文與明文之間不滿足線性關(guān)系。如果直接用只包含強(qiáng)度分布的密文通過與加密過程相反的逆過程進(jìn)行解密，是無法得到原始圖像的，得到的仍然是隨機(jī)噪聲，這一點(diǎn)在后面的仿真分析中得到了驗證。

由于改進(jìn)系統(tǒng)的加密密文是4-f系統(tǒng)輸出面上的強(qiáng)度分布，并且已知原始圖像是實值圖像，因此解密過程實際上是已知一個強(qiáng)度約束的相位恢復(fù)問題［19］。因此，在解密時可以利用相位恢復(fù)算法來恢復(fù)原始圖像。如果已知密文和解密密鑰（也就是在光學(xué)加密過程中使用的2個隨機(jī)相位掩膜exp［i2πn（x，y）］和exp［i2πb（u，v）］），就可以利用相位恢復(fù)算法通過迭代的方法恢復(fù)原始圖像。解密過程采用數(shù)字方式，具體步驟如下：

（1）任意賦予原始圖像一個初始估計值f0（x，y），然后開始迭代過程；

（2）在第k次迭代過程中，輸入圖像fk（x，y）在輸入平面和傅里葉平面分別被相位掩膜exp［i2πn（x，y）］和exp［i2πb（u，v）］調(diào)制，在輸出平面上得到光波分布函數(shù)gk（x，y）：

（3）保留gk（x，y）的相位，振幅用密文I（x，y）的平方根代替，得到修正了的光波分布函數(shù)gkm（x，y）：

（5）更新輸入圖像的取值，得到fk＋1（x，y）：

（6）計算輸出平面上得到的分布函數(shù)gk（x，y）的振幅與密文的平方根［I（x，y）］1／2的相關(guān)系數(shù)［27］。相關(guān)系數(shù)代表著恢復(fù)圖像與明文圖像的符合程度，相關(guān)系數(shù)越接近于1，就表示解密圖像與原始圖像越相似。令［I（x，y）］1／2＝f，｜gk（x，y）｜＝fm，則相關(guān)系數(shù)CC可表示為

（7）重復(fù)步驟（2）～（6），直到相關(guān)系數(shù)達(dá)到所設(shè)定的閾值，迭代過程結(jié)束。

經(jīng)過以上的迭代過程，就可以恢復(fù)出原始圖像。但是由于在加密過程中引入了隨機(jī)相位掩膜，恢復(fù)的原始圖像存在著一定的散斑噪聲干擾。為了提高解密結(jié)果的圖像質(zhì)量，還需要利用散斑噪聲抑制算法來抑制散斑噪聲。由于該解密算法采用純數(shù)字方法，因此可以非常方便地引入各種數(shù)字圖像處理技術(shù)來抑制噪聲，改善解密的效果。在這里，我們采用非常簡單的中值濾波的方法抑制散斑噪聲。

2 安全性分析

我們所提出的方法相對于傳統(tǒng)的DRPE加密系統(tǒng)具有更高的安全性。這主要是基于以下2個原因：

一方面，我們所提出的方法在解密時需要的密鑰更多。傳統(tǒng)的DRPE加密系統(tǒng)在解密時只需要利用位于傅里葉平面的一個解密密鑰就可以得到正確的解密輸出。而我們所提出的方法中，要想正確恢復(fù)原始圖像，必需同時知道2個隨機(jī)相位板的相位分布。因此，我們所提出的方法需要2個解密密鑰，密鑰空間更大，安全性更高。同時，還可以將2個解密密鑰分配給不同的用戶持有，只有2個密鑰持有人同時提供各自的密鑰，解密系統(tǒng)才可以正確恢復(fù)原始圖像，這樣就可以進(jìn)一步提高加密系統(tǒng)的安全性。

另一方面，我們所提出的方法的明文和密文之間不滿足線性關(guān)系，可以更好地抵抗各種攻擊。傳統(tǒng)的DRPE加密系統(tǒng)的密鑰空間大，能夠很好地抵抗暴力攻擊。但是，由于其本質(zhì)上是一種線性系統(tǒng)，因而它存在著巨大的安全隱患，容易受到“選擇密文攻擊”、“已知明文攻擊”、“選擇明文攻擊”和“唯密文攻擊”等幾種形式的攻擊。相對來說，我們所提出的這種改進(jìn)的單強(qiáng)度記錄光學(xué)加密系統(tǒng)就可以很好地抵抗這幾種形式的攻擊。因為這幾種攻擊一般都需要將密文經(jīng)過傅里葉變換得到頻譜面上的密文［14－17］，我們所提的這種光學(xué)加密系統(tǒng)只保留輸出面上的光強(qiáng)分布作為密文，沒有記錄相位信息，破壞了密文和明文之間的線性關(guān)系，因此無法利用同樣的方法來有效地破解加密系統(tǒng)。所以說，該改進(jìn)的單強(qiáng)度記錄光學(xué)加密系統(tǒng)具有更高的安全性。后面的計算機(jī)仿真也驗證了這里的分析。

3 計算機(jī)仿真

受實驗條件所限，我們暫時沒有能夠進(jìn)行相應(yīng)的光學(xué)實驗，因此我們主要在計算機(jī)上使用MATLAB 7.1進(jìn)行了仿真實驗。一系列仿真實驗驗證了該單強(qiáng)度記錄光學(xué)加密系統(tǒng)的有效性和安全性。

在有效性仿真中，我們分別選取了二值圖像和灰度圖像作為原始圖像，利用所提出的單強(qiáng)度記錄光學(xué)加密系統(tǒng)進(jìn)行光學(xué)加密和解密。該單強(qiáng)度記錄光學(xué)加密系統(tǒng)對二值圖像（256pixel×256pixel）進(jìn)行加密和解密的仿真結(jié)果如圖2所示。其中，圖2（a）是原始二值圖像，圖2（b）是經(jīng)過單強(qiáng)度記錄光學(xué)加密系統(tǒng)加密后得到的密文，該密文滿足隨機(jī)白噪聲分布。圖2（c）是利用迭代解密算法得到的解密輸出?？梢钥吹?，利用密文和解密密鑰，經(jīng)過相位恢復(fù)算法的迭代解密過程，是可以恢復(fù)原始圖像的。但是解密輸出相比于原始圖像，疊加了一定的散斑噪聲，也就是說該單強(qiáng)度記錄光學(xué)加密系統(tǒng)是有損的。圖2（d）是當(dāng)解密密鑰為錯誤密鑰時的解密結(jié)果。很顯然，這時的解密結(jié)果是隨機(jī)噪聲。也就是說，如果不知道解密密鑰，是無法正確恢復(fù)原始圖像的。圖2（e）是由密文直接利用雙隨機(jī)相位加密的逆過程得到的仿真結(jié)果，該仿真結(jié)果為白噪聲分布，說明僅利用4-f系統(tǒng)輸出平面上的光強(qiáng)分布直接利用DRPE的逆過程是無法恢復(fù)原始圖像的。

圖2 對二值圖像加密和解密的仿真結(jié)果Fig.2 Simulation results of encryption and decryption for binary image

利用我們所提出的單強(qiáng)度記錄光學(xué)加密系統(tǒng)對一幅灰度圖像（Lena圖，256pixel×256pixel）進(jìn)行加密和解密的仿真結(jié)果如圖3所示。其中，圖3（a）是原始灰度圖像，圖3（b）是經(jīng)過單強(qiáng)度記錄光學(xué)加密系統(tǒng)加密后得到的密文，圖3（c）是利用迭代解密算法得到的解密輸出?？梢钥吹?，對于灰度圖像，利用密文和解密密鑰，經(jīng)過相位恢復(fù)算法的迭代解密過程同樣是可以恢復(fù)原始圖像的。圖3（d）是當(dāng)解密密鑰為錯誤密鑰時的解密結(jié)果，圖3（e）是由密文直接利用雙隨機(jī)相位加密的逆過程得到的仿真結(jié)果?？梢钥吹?，圖2和圖3的仿真結(jié)果是相似的，也就說明，我們所提出的單強(qiáng)度記錄光學(xué)加密系統(tǒng)既可以實現(xiàn)對二值圖像的加密，也可以實現(xiàn)對灰度圖像的加密。

圖3 對灰度圖像加密和解密的仿真結(jié)果Fig.3 Simulation results of encryption and decryption for grayscale image

圖4 是迭代過程中恢復(fù)圖像與明文圖像之間的相關(guān)系數(shù)CC的變化曲線，其中相關(guān)系數(shù)的計算如公式（6）所示。從圖4可以看到，針對二值圖像和灰度圖像進(jìn)行迭代解密時相關(guān)系數(shù)的兩條變化曲線都是單調(diào)遞增的，說明該迭代解密過程是收斂的。對比兩條曲線也可以發(fā)現(xiàn)，經(jīng)過400次迭代之后，恢復(fù)的灰度圖像與原始圖像的相關(guān)系數(shù)可以達(dá)到0.937 2，而恢復(fù)的二值圖像與原始圖像的相關(guān)系數(shù)只有0.724 1。這主要是由于二值圖像的圖案相對更加簡單，每個像素的取值只有0和1，因此同樣的散斑噪聲對于二值圖像的影響就相對更大些，導(dǎo)致恢復(fù)的二值圖像與原始圖像的相關(guān)系數(shù)比較小。不過，根據(jù)圖2的仿真結(jié)果也可以看到，雖然疊加的散斑噪聲使得通過迭代解密算法恢復(fù)的二值圖像與原始圖像的相關(guān)系數(shù)并不是很大，但是仍然是可以由恢復(fù)的二值圖像獲得準(zhǔn)確的原始信息的。

圖4 相關(guān)系數(shù)CC變化曲線Fig.4 Curve of correlation coefficient

由前面的仿真結(jié)果可以看到，經(jīng)過迭代解密算法恢復(fù)的圖像疊加有一定的散斑噪聲，影響了解密輸出的質(zhì)量。由于我們所提出的解密方法是采用純數(shù)字方式進(jìn)行迭代計算，因此可以在解密算法中非常方便地引入相應(yīng)的數(shù)字圖像處理技術(shù)來減小散斑噪聲。這里我們采用比較簡單的中值濾波方法來減小散斑噪聲。經(jīng)過迭代解密算法恢復(fù)的圖像經(jīng)過窗口尺寸為4pixel×4pixel的中值濾波器之后的輸出圖像如圖5所示。其中，圖5（a）是經(jīng)過中值濾波之后的二值圖像，它與原始二值圖像的相關(guān)系數(shù)為0.915 6；圖5（b）是經(jīng)過中值濾波之后的灰度圖像，它與原始灰度圖像的相關(guān)系數(shù)為0.937 4。由圖5可以看到，經(jīng)過中值濾波之后，輸出圖像相比于迭代解密結(jié)果（圖2（c）和圖3（c）），疊加的散斑噪聲在一定程度上減小了，輸出圖像和原始圖像的相關(guān)系數(shù)在一定程度上增大了，但代價是圖像變得比較模糊。可以預(yù)見的是，如果應(yīng)用更加復(fù)雜的數(shù)字圖像處理技術(shù)是可以得到更好的恢復(fù)效果的。

圖5 經(jīng)過中值濾波之后的輸出圖像Fig.5 Output images after median filtering

上面的仿真實驗驗證了我們所提出的雙隨機(jī)相位加密方法的有效性，接下來我們通過與傳統(tǒng)DRPE系統(tǒng)的對比實驗來驗證我們所提方法的安全性。在這里主要利用在文獻(xiàn)［15］和［16］中提出的“已知明文攻擊”（known plaintext attack，KPA）和“選擇明文攻擊”（chosen plaintext attack，CPA）兩種攻擊方法分別對傳統(tǒng)DRPE方法和我們所提的方法進(jìn)行破解。篇幅所限，在這里僅列出對灰度圖像進(jìn)行破解的仿真結(jié)果。利用破解密鑰進(jìn)行解密的仿真結(jié)果如圖6所示。其中，圖6（a）和圖6（b）是采用“已知明文攻擊”分別對傳統(tǒng)DRPE和我們所提出的加密方法進(jìn)行破解，然后利用破解密鑰解密的仿真結(jié)果。圖6（c）和圖6（d）是采用“選擇明文攻擊”分別對傳統(tǒng)DRPE和我們所提出的加密方法進(jìn)行破解，然后利用破解密鑰解密的仿真結(jié)果。

圖6 安全性分析的仿真結(jié)果Fig.6 Simulation results of security analysis

通過仿真結(jié)果可以看到，傳統(tǒng)的DRPE系統(tǒng)容易受到KPA和CPA等幾種形式的攻擊，這主要是由于它在本質(zhì)上是一個線性系統(tǒng)，系統(tǒng)輸出和輸入之間滿足簡單的線性關(guān)系。而我們所提出的雙隨機(jī)相位加密系統(tǒng)在輸出平面上只記錄光強(qiáng)分布，這樣使得系統(tǒng)輸出和輸入不再滿足線性關(guān)系，因而利用同樣的方法破解密鑰進(jìn)行解密，其輸出仍然為隨機(jī)噪聲，無法正確恢復(fù)原始圖像。光學(xué)加密方法的應(yīng)用范圍。因此如何減小散斑噪聲，進(jìn)一步提高解密輸出圖像的質(zhì)量是今后繼續(xù)深入研究的重點(diǎn)。

4 結(jié)論

本文提出了一種基于相位恢復(fù)算法的單強(qiáng)度記錄光學(xué)加密技術(shù)。該光學(xué)加密技術(shù)將原始圖像置于4-f系統(tǒng)的輸入平面上進(jìn)行雙隨機(jī)相位加密，記錄輸出平面上的光強(qiáng)分布作為密文。加密過程可以采用純光學(xué)方法實現(xiàn)。在解密時，利用相位恢復(fù)迭代算法就可以由密文恢復(fù)出所加密的原始圖像。解密過程采用數(shù)字方式，因此可以非常方便地引入中值濾波等數(shù)字圖像處理方法來進(jìn)一步提高解密圖像的質(zhì)量。該光學(xué)加密技術(shù)僅需一次曝光，只需要記錄輸出面上的光強(qiáng)分布作為密文，便于光學(xué)實現(xiàn)。同時加密結(jié)果可以利用CCD等感光器件直接記錄，記錄的密文直接是數(shù)字信號，可以非常容易地進(jìn)行存儲和傳輸。此外，由于該加密系統(tǒng)的輸出與輸入之間不滿足線性關(guān)系，因而不易受到各種攻擊，安全性較高。計算機(jī)仿真實驗結(jié)果驗證了系統(tǒng)的有效性和安全性。需要指出的是，該光學(xué)加密系統(tǒng)的迭代解密結(jié)果是有損的，解密輸出疊加有一定的散斑噪聲，這將會限制該

［1］ Philippe R，Bahram J.Optical image encryption based on input plane and Fourier plane random encoding［J］.Optics Letters，1995，20（7）：767-769.

［2］ Nomura T，Bahram J.Optical encryption using a joint transform correlator architecture ［J］.Optical Engineering，2000，39（8）：2031-2035.

［3］ Unnikrishnan G，Joseph J，Singh K.Optical encryption by double-random phase encoding in the fractional Fourier domain［J］.Optics Letters，2000，25（12）：887-889.

［4］ Situ Guohai，Zhang Jingjuan.Double random-phase encoding in the Fresnel domain［J］.Optics Letters，2004，29（14）：1584-1586.

［5］ Zhang Yan， Wang Bo.Optical image encryption based on interference［J］.Optics Letters，2008，33（21）：2443-2445.

［6］ Li Youzhi，Kreske K，Rosen J.Security and encryption optical systems based on a correlator with significant output images［J］.Applied Optics，2000，39（29）：5295-5301.

［7］ Chang H T.Image encryption using separable amplitude-based virtual image and iteratively retrieved phase information ［J］.Optical Engineering，2001，40（10）：2165-2171.

［8］ Chen Linfei，Zhao Daomu.Optical image encryption based on fractional wavelet transform ［J］.Optics Communications，2005，254（4-6）：361-367.

［9］ Pérez-CabréE，Cho M，Javidi B.Information authentication using photon-counting double-randomphase encrypted images［J］.Optics Letters，2011，36（1）：22-24.

［10］Shi Yishi，Li Tuo，Wang Yali，et al.Optical image encryption via ptychography ［J］.Optics Letters，2013，38（9）：1425-1427.

［11］Qin Yi，Wang Hongjuan，Wang Li，et al.Virtual optical image encryption based on random amplitudephase encoding ［J］.Opto-Electronic Engineering，2012，39（7）：126-131.秦怡，王紅娟，王麗，等.隨機(jī)振幅－相位編碼虛擬光學(xué)加 密 系 統(tǒng) ［J］.光 電 工 程，2012，39（7）：126-131.

［12］Qin Yi，Gong Qiong，LüXiaodong，et al.Double image encryption system based on principles of interference［J］.Journal of Applied Optics，2013，34（2）：260-266.秦怡，鞏瓊，呂曉東，等.基于干涉原理的雙圖像加密系統(tǒng)［J］.應(yīng)用光學(xué)，2013，34（2）：260-266.

［13］Wang Hongjuan，Zhang Yingying，Huang Yiding，et al.Color image encryption technology based on asymmetric cryptosystem ［J］.Journal of Applied Optics，2013，34（6）：980-984.王紅娟，張穎穎，黃義定，等.基于非對稱密碼系統(tǒng)的彩色圖像加密技術(shù)［J］.應(yīng)用光學(xué)，2013，34（6）：980-984.

［14］Carnicer A，Montes-Usategui M，Arcos S，et al.Vulnerability to chosen-cyphertext attacks of optical encryption schemes based on double random phase keys［J］.Optics Letters，2005，30（13）：1644-1646.

［15］Peng Xiang，Zhang Peng，Wei Hengzheng，et al.Known plaintext attack on double random phase encoding encryption technique［J］.Acta Physica Sinica，2006，55（3）：1130-1136.彭翔，張鵬，位恒政，等.雙隨機(jī)相位加密系統(tǒng)的已知明 文 攻 擊 ［J］. 物 理 學(xué) 報，2006，55（3）：1130-1136.

［16］Wei Hengzheng，Peng Xiang，Zhang Peng，et al.Chosen-plaintext attack on double random phase encoding encryption technique［J］.Acta Optica Sinica，2007，27（5）：824-829.位恒政，彭翔，張鵬，等.雙隨機(jī)相位加密系統(tǒng)的選擇明文攻擊［J］.光學(xué)學(xué)報，2007，27（5）：824-829.

［17］Peng Xiang，Tang Hongqiao，Tian Jindong.Ciphertext-only attack on double random phase encoding optical encryption system ［J］.Acta Physica Sinica，2007，56（5）：2629-2636.彭翔，湯紅喬，田勁東.雙隨機(jī)相位編碼光學(xué)加密系統(tǒng)的唯密文攻擊［J］.物理學(xué)報，2007，56（5）：2629-2636.

［18］Li Rong，Li Ping.Research on the image security in double random phase real-value encryption［J］.Acta Photonica Sinica，2005，34（6）：952-955.李榕，李萍.雙隨機(jī)相位圖像加密的實值編碼研究［J］.光子學(xué)報，2005，34（6）：952-955.

［19］Gerchberg R W，Saxton W O.A practical algorithm for the determination of phase from image and diffraction phase pictures［J］.Optik，1972，35（2）：237-246.

［20］Fienup J R.Phase retrieval algorithms：a comparison［J］.Applied Optics，1982，21（15）：2758-2769.

［21］Wang R K，Watson L A，Chatwin C.Random phase encoding for optical security ［J］.Optical Engineering，1996，35（9），2464-2469

［22］Chang HT，Lu WC，Kuo CJ.Multiple-phase retrieval for optical security systems by use of randomphase encoding［J］.Applied Optics，2002，41（23），4815-4834.

［23］Situ G，Zhang J.A lensless optical security system based on computer-generated phase only masks［J］.Optics Communications，2004，232（1-6），115-122.

［24］Liu Zhengjun，Liu Shutian.Double image encryption based on iterative fractional Fourier transform ［J］.Optics Communications，2007，275（2）：324-329.

［25］Liu Zhengjun，Xu Lie，Guo Qing，et al.Image watermarking by using phase retrieval algorithm in gyrator transform domain［J］.Optics Communications，2010，283（24）：4923-4927.

［26］Liu Wei，Liu Zhengjun，Liu Shutian.Asymmetric cryptosystem using random binary phase modulation based on mixture retrieval type of Yang-Gu algorithm［J］.Optics Letters，2013，38（10）：1651-1653.

［27］Situ Guohai，Zhang Jingjuan. Multiple-image encryption by wavelength multiplexing［J］.Optics Letters，2005，30（11）：1306-1308.