《電子技術應用》
您所在的位置:首頁 > 可編程邏輯 > 設計應用 > 基于Logistic映射PN序列的FPGA實現
基于Logistic映射PN序列的FPGA實現
摘要: 偽噪聲序列(PN序列)可應用于擴頻通信、信息加密、計算機仿真等領域。PN序列發生器需要一個隨機信號源和一系列的離散、量化算法及其硬件實現技術。確定性的混沌可以復制,具有長期不可預測性,且很難區分一個信號是來自于非確定性系統還是混沌系統。因此,混沌滿足密碼系統設計的基本原則,利用混沌系統作為PN序列的信號源已引起了國內外學者的廣泛關注與研究。
關鍵詞: FPGA Logistic 混沌系統
Abstract:
Key words :

  0 引 言

  偽噪聲序列(PN序列)可應用于擴頻通信、信息加密、計算機仿真等領域。PN序列發生器需要一個隨機信號源和一系列的離散、量化算法及其硬件實現技術。確定性的混沌可以復制,具有長期不可預測性,且很難區分一個信號是來自于非確定性系統還是混沌系統。因此,混沌滿足密碼系統設計的基本原則,利用混沌系統作為PN序列的信號源已引起了國內外學者的廣泛關注與研究。

  基于Logistic映射產生PN序列已有不少研究。文獻[7]利用模擬電路研究了其實現技術,由于混沌對初值和參數的敏感性,兩個模擬電路實現的PN序列,其電路參數不可能完全匹配,且受環境條件的變化而出現失配現象,很難保證保密通信或信息加密中兩個PN序列的完全同步。因此,近年來開展了基于數字電路(如FPGA)的PN序列的實現研究,但現有的研究中大多進行仿真實驗,無法保證硬件實際輸出PN序列的產生及其性能的測試。為此,基于FPGA技術,利用Logistic混沌映射作為隨機信號源對實現PN序列的硬件進行了深入研究,提取Logistic數值序列中二進制數的某一位bi(bi∈{0,1})為PN數字序列,實驗中獲得了硬件輸出序列,并通過一個串口通信電路對硬件輸出序列進行取樣,對其進行序列統計性能分析。分析結果表明,基于Logistic的PN序列性能良好,滿足PN序列測試的標準。

  1 Logistic映射

  Logistic方程如下:

公式

  式中:Xn∈(0,1);μ∈(0,4),當μ取值[3.571 448,4]時,Logistic映射進入混沌態,并表現出復雜的動力學特性。在此取μ=4,令Xn的初值Xo=O.312 5,則xn的時間序列和吸引子相圖如圖1所示。

xn的時間序列和吸引子相圖

  2 Logistic映射的電路設計及時間序列的量化

  DSP Builder。將Matlab和Simlalink系統級設計工具的算法開發、仿真和驗證功能與VHDL綜合、仿真和開發工具整合在一起,實現了這些工具的集成,將系統級設計的實現與DSP算法的開發相鏈接,涵蓋了算法和存儲功能等基本操作。可利用DSP Builder模塊迅速生成算法硬件電路,并轉成VHDL語言,大大縮短了FPGA的設計周期。圖2為Logistic的DSP Builder電路實現模型。

Logistic的DSP Builder電路實現模型

  圖2中為消除毛刺,在各輸出端加入了一個延時器。其中,xout為Xn+1的時間序列;yout為其延時一個時鐘周期后的輸出,即Xo。xout與yout都是以二進制數表示的一系列小數(xn∈(0,1)),可表示為:

公式

  因為時間序列Xn的數值是用二進制數表示的,隨著迭代的不斷進行,Xn將隨之變化,此時其二進制表示中的某一位(O或1)也隨之變化。因此可以提取Xn時間序列中某一位二進制位作為量化值,當Xn的時間序列值隨時間變化時即生成一PN序列。圖2中的XoutBit為PN序列的輸出端,它取自Xn二進制數表示中的第9位。后面的分析表明,這種量化方法同樣具有很好的混沌性能和隨機性。

  為了驗證圖2電路模塊的正確性,對Logistic的模塊電路進行仿真,其時間序列和吸引子如圖3所示。比較圖3和圖1發現,利用DSP Builder設計的電路有效地實現了Logistic映射的功能。

利用DSP Builder設計的電路有效地實現了Logistic映射的功能

  3 基于FPGA的PN序列實現

  圖2中,xout與yout具有相同的k及l值,k=1,l=9,即k+l=10。因為所使用的FPGA開發板自帶THS5651DA轉換器接受的是10引腳數據,k+l=10可實現匹配。盡管xout在最后輸出前轉成了10位長度,但在整個內部運算環節,是50位長度的,有限字長明顯大于10,這樣可以最大限度地保持精度,減少有限字長效應帶來的誤差,有效實現混沌的非周期特性。

  理論上混沌序列是非周期的,但因為在使用有限字長表示混沌狀態的情況下,加上運算結果的近似化,混沌序列最終會演化為一個周期序列。但是可以通過加長有限字長來延長周期。從實用的角度看,只要設計得好,映射和分叉參數選擇合適,所產生的序列周期將足夠長,序列的相關特性將足夠好,可用的序列數量將足夠多,可認為位轉換就已夠用,因為通信的時間不可能無限長。研究還發現,混沌映射產生的有限字長周期序列具有正的最大Lyapllnov指數。

  將圖2所示的圖形化模塊通過“Signal Corepiler”直接轉化為硬件描述語言,再利用QuartusⅡ進行編譯,其仿真時序圖見圖4。其中,Output為抽取的混沌數字PN序列,即圖2中的XoutBit;Outputl和Output2分別為Xn+1和Xn。

仿真時序圖

  把工程文件在CycloneⅡ的EP2C35F484C8進行編譯和引腳分配并下載到該芯片進行硬件實驗,從數字示波器上觀察到的PN序列如圖5(a)所示。為了驗證硬件實現的精度,通過THS5651芯片對得到的混沌數字信號進行D/A轉化,再通過示波器觀察其時間序列和吸引子,如圖5(b),圖5(c)所示。

 

從數字示波器上觀察到的PN序列

  4 PN序列隨機性分析

  為了對硬件輸出的PN序列進行性能測試,首先通過一個串口通信模塊電路對FPGA輸出的數字PN序列進行采樣,然后利用計算機對其進行隨機性能分析。這里采用的FPGA開發板其內核為CycloneⅡ的EP2(235F484C8,外圍電路提供了USB口、并行I/O以及串口等一系列可供用戶選擇的端口。抽取Xn二進制表示的小數部分的第9位作為比特序列,并通過串口來實時獲得該位在一定時間段t內的所有值,其串口通信模塊電路如圖6所示。

串口通信模塊電路

  為了確保計算機在接收該序列時不出現比特丟失現象,可預先向計算機發送一組周期序列來保證數據的可靠性。假設串口傳輸波特率為9 600,系統晶振為20 MHz,為了得到整數倍的分頻,可根據公式:

n=(20×106×k)/9600

  令k=3,則可得出n為6 250。其中,比例因子k可通過鎖相環來獲得。上述兩步驟可通過在QuartusⅡ中添加鎖相環代碼與分頻器代碼實現。此時,讓系統時鐘T=1/9 600再經過二分頻后通過串口輸出,如果在接收端收到“01010”的周期序列,則表明FPGA與計算機之間的通信無問題。

  5 PN序列隨機性能分析

  5.1 頻率測試(Frequency Test)

  該測試可用來判斷序列當中0和1的個數是否近似相等,若近似相等,則符合隨機序列的必要條件。令n0,n1分別代表0與1的個數;n代表序列s的長度。

  該統計方程為:

T1=(n0-n1)2/n

  若n≥lO,則T1近似符合1維自由度的x2分布。

  對采集到的PN序列進行自由度為l的x2檢驗,取顯著性水平α=0.05,可以從標準表中查到顯著性水平a=O.05的x2值為3.84。如果T1值小于3.84,則該序列通過測試。n0,n1和T1的測試值如表l所示,由此可知,該序列的T1值遠小于3.84,通過測試。

n0

  5.2 串列測試(Serial Test)

  該測試可用來判斷序列中00,01,10,11子序列的個數是否近似相等。若近似相等,則符合隨機序列的必要條件。令n0,n1分別代表0與1的個數,而n00,n01,n10,n11分別代表s序列中00,01,10,11的個數。由于n00+n01+no+nll=n-l,該統計方程為:

公式

  若n≥21,則T2近似符合2維自由度的x2分布。對序列進行自由度為2的x2檢驗。顯著性水平α置為O.05,查標準表得此時的x2值為5.99。測試結果如表2所示,T2=O.352 5《5.99,通過測試。

測試結果

  5.3 Poker測試(Poker Test)

  令m為符合(n/m)≥5×2m的正整數,并令k=n/m。將序列s分成k個不重復的部分,每部分的長度為m,令ni為第i個長度為m的序列標號,1≤i≤2m。Poker測試可用來確定ni在序列s中,長度為m的各個不同子序列出現的概率是否近似相等。該統計方程為:

公式

  式中:T3近似符合2m-1維自由度的x2分布。值得注意的是,頻率測試其實是Poker測試的一個特例,令m=1,Poker測試即成為了頻率測試。分別取m=1,2,3,4,5,6,對待測序列進行自由度為2m-1的x2檢驗,結果列于表3。由表3可見,對于不同m值的T3值均小于對應的x2值,通過測試。

Poker測試結果

  5.4 游程測試(RLins Test)

  該測試的目的是為了判斷s序列中不同長度的游程是否具有隨機序列的游程特性。理論上,在一個長度為n隨機序列中,長度為i的連續O或連續1的期望值為ei=(n—i+3)/2i+2。令k為當ei≥5時的最大的i;Bi,Gi分別為各個block和gap的個數。游程測試的統計方程為:

公式

  則T4近似符合2k-2維自由度的χ2分布。

  分別取i為1,2,3,4,5,6,7,8,9。對本序列,k值為9,所以可對T4進行自由度為16的χ2檢驗。查標準表知自由度為16的χ2值為26.3。將上述數值代人統計方程,運算得T4=19.883 1<26.3,通過測試。其測試結果如表4所示。

游程測試結果

  5.5 自相關測試(Autocorrelation Test)

  該測試是為了檢測s序列與其(非循環的)移動后的序列的相關性。令d為一個整數,1≤d≤(n/2)。比較s序列和它移動d位之后的序列,兩者不同的位的數目可由公式計算得出。⊕表異或運算。自相關測試的統計方程為:

公式

  當,n-d≥10時,T5近似符合N(O,1)分布。

  為判斷序列是否符合N(O,1)分布,可用余誤差函數計算公式再將erfc(z)與顯著性水平α=O.05比較,若大于它,則通過測試,測試結果如表5所示。由檢驗結果可知,每一個erfc(z)都大于a,通過測試。

 

自相關測試結果

  6 結 語

  利用Logistic混沌映射作為隨機信號源,基于FP-GA技術研究了PN序列的硬件產生。通過對產生的PN序列進行性能分析發現,其具有很好的隨機性,可被應用于擴頻通信和數字信息加密等領域。由于直接以離散混沌作為隨機信號源,相對于將連續混沌離散化后再量化產生PN序列,其精確度更高。由于數字序列是通過抽取二進制數值中某一位產生的,只要迭代值的位數足夠多,序列的可變性就更大,從而增大了密鑰空間和通信的保密性。

此內容為AET網站原創,未經授權禁止轉載。
亚洲一区二区欧美_亚洲丝袜一区_99re亚洲国产精品_日韩亚洲一区二区
亚洲美女精品一区| 99热这里只有精品8| 欧美成人日韩| 久久久亚洲一区| 欧美亚洲综合久久| 午夜精品视频在线观看| 中文国产亚洲喷潮| 9l视频自拍蝌蚪9l视频成人| 亚洲精品1区| 亚洲国产精品视频一区| 欧美中文在线观看国产| 欧美一区二区三区啪啪| 午夜免费电影一区在线观看| 亚洲午夜久久久久久久久电影院| 日韩亚洲精品视频| 亚洲另类自拍| 亚洲精品日产精品乱码不卡| 亚洲国产日韩欧美| 亚洲国产综合视频在线观看| 在线不卡免费欧美| 亚洲第一免费播放区| 亚洲成人在线网站| 亚洲级视频在线观看免费1级| 亚洲成在人线av| 1024成人网色www| 亚洲高清视频在线观看| 一区二区动漫| 日韩午夜高潮| 亚洲性线免费观看视频成熟| 亚洲午夜在线| 翔田千里一区二区| 欧美一区二区视频网站| 久久久777| 久久夜色精品国产噜噜av| 久热精品在线视频| 欧美精品日韩三级| 欧美偷拍另类| 国产午夜精品全部视频播放| 国产一区二区三区高清在线观看| 黑人巨大精品欧美一区二区小视频 | 美女主播精品视频一二三四| 欧美国产日本| 国产精品h在线观看| 国产欧美一区二区三区沐欲| 精品999网站| 日韩视频一区二区三区| 亚洲男人的天堂在线aⅴ视频| 久久国产天堂福利天堂| 日韩午夜在线| 性欧美xxxx视频在线观看| 久久深夜福利免费观看| 欧美全黄视频| 国产区二精品视| 亚洲国产一区二区三区青草影视 | 国产亚洲a∨片在线观看| 在线精品视频一区二区| 99精品国产一区二区青青牛奶 | 亚洲字幕在线观看| 亚洲国产精品悠悠久久琪琪| 一区二区三区高清在线| 久久精品国产v日韩v亚洲| 欧美成人a视频| 国产精品成人一区二区艾草| 国产亚洲欧美日韩在线一区 | aa级大片欧美| 久久国产99| 宅男精品视频| 久久阴道视频| 亚洲久色影视| 一本久久知道综合久久| 欧美亚洲一区二区三区| 欧美成人日韩| 国产酒店精品激情| 亚洲欧洲精品一区二区三区 | 亚洲毛片网站| 久久精品欧洲| 欧美性猛交xxxx乱大交蜜桃| 黄色成人av在线| 亚洲一区二区三区高清| 最新国产乱人伦偷精品免费网站| 亚洲影院免费观看| 免费久久99精品国产自| 国产精品拍天天在线| 最近中文字幕日韩精品| 欧美在线观看视频一区二区三区| 亚洲午夜三级在线| 欧美国产视频在线观看| 国产一区日韩二区欧美三区| 一本色道久久99精品综合| 亚洲激情一区二区| 久久er精品视频| 欧美性大战久久久久| 亚洲人成欧美中文字幕| 亚洲成人在线视频网站| 欧美一级电影久久| 欧美午夜电影一区| 亚洲日本中文字幕免费在线不卡| 久久精品国产第一区二区三区| 亚洲一区在线免费观看| 欧美日本精品| 亚洲黄色成人网| 亚洲第一在线视频| 久久精品国产99精品国产亚洲性色| 国产精品国产a| 日韩一二三在线视频播| 日韩视频在线一区二区| 欧美 日韩 国产精品免费观看| 国产日韩在线一区| 午夜电影亚洲| 欧美伊人久久久久久久久影院| 欧美午夜理伦三级在线观看| 亚洲精品久久久蜜桃| 亚洲伦理一区| 欧美黄色视屏| 亚洲激情成人网| 亚洲精品久久嫩草网站秘色| 欧美/亚洲一区| 亚洲国产精品久久久久久女王| 亚洲激情在线观看| 女女同性精品视频| 亚洲激情av| 99在线热播精品免费99热| 欧美精品日本| 99精品久久久| 午夜视频久久久| 日韩亚洲视频| 欧美高清视频一区二区| 亚洲黄色一区二区三区| 日韩视频在线观看一区二区| 欧美精品久久天天躁| 日韩亚洲欧美在线观看| 亚洲专区免费| 国产欧美精品在线播放| 欧美在线不卡视频| 久热综合在线亚洲精品| 亚洲国产高清aⅴ视频| 亚洲卡通欧美制服中文| 欧美激情精品久久久久久| 91久久精品一区二区别| 99re6这里只有精品| 欧美人交a欧美精品| 一本色道久久综合| 午夜一区不卡| 国产一区二区三区四区在线观看| 久久av二区| 欧美1级日本1级| 亚洲伦理久久| 亚洲欧美一区二区激情| 国产日韩欧美成人| 亚洲高清自拍| 欧美精品久久久久a| 中日韩高清电影网| 久久精品国产96久久久香蕉| 在线看片一区| 亚洲视频大全| 国产午夜久久| 亚洲精品中文字| 欧美午夜电影一区| 午夜精品久久久久99热蜜桃导演| 久久综合狠狠综合久久激情| 91久久精品美女| 亚洲欧美日韩精品久久久| 国产视频在线观看一区二区| 亚洲福利视频一区二区| 欧美日韩亚洲一区三区| 午夜精品免费| 欧美国内亚洲| 亚洲一区二区三区免费在线观看| 久久久精品国产免大香伊 | 一本久久综合亚洲鲁鲁| 国产精品区免费视频| 久久精品水蜜桃av综合天堂| 欧美精品久久久久久久久老牛影院| 亚洲一区二区三区免费在线观看| 久久影音先锋| 一区二区黄色| 久久亚洲午夜电影| 91久久久久久国产精品| 亚洲乱码精品一二三四区日韩在线| 亚洲伊人网站| 亚洲高清自拍| 篠田优中文在线播放第一区| 在线精品在线| 欧美一级淫片aaaaaaa视频| 亚洲国产一区二区三区a毛片| 欧美亚洲一区| 亚洲精品免费观看| 久久久噜噜噜| 亚洲一级影院| 欧美激情综合网| 久久精品国产第一区二区三区| 欧美视频手机在线| 亚洲国产日韩精品| 国产欧美精品一区二区三区介绍| 99re66热这里只有精品4| 国内精品视频在线播放| 亚洲一区免费观看| 亚洲精品免费电影| 毛片精品免费在线观看| 亚洲欧美欧美一区二区三区|