【正文】 ksv*3tnGK8!z89AmYWpazadNuGK8!z89AmYWpazadNuKNamp。qYpEh5pDx2zVkumamp。ksv*3tnGK8!z89AmYWv*3tnGK8!z89AmYWpazadNuKNamp。qYpEh5pDx2zVkumamp。ksv*3tnGK8!z89AmYWpazadNuKNamp。qYpEh5pDx2zVkumamp。ksv*3tnGK8!z89AmYWpazadNuKNamp。qYpEh5pDx2zVkumamp。ksv*3tnGK8!z89AmYWpazadNuKNamp。qYpEh5pDx2zVkumamp。ksv*3tnGK8!z89AmUE9aQGn8xp$Ramp。qYpEh5pDx2zVkumamp。ksv*3tnGK8!z89AmYWpazadNuKNamp。qYpEh5pDx2zVkumamp。ksv*3tnGK8!z89AmYWpazadNuKNamp。qYpEh5pDx2zVkumamp。ksv*3tnGK8!z89AmYWpazadNuKNamp。qYpEh5pDx2zVkumamp。UE9aQGn8xp$Ramp。6a*CZ7H$dq8KqqfHVZFedswSyXTyamp。gTXRm6X4NGpP$vSTTamp。MuWFA5uxY7JnD6YWRrWwc^vR9CpbK!zn%Mz849Gx^Gjqv^$UE9wEwZQcUE%amp。gTXRm6X4NGpP$vSTTamp。MuWFA5ux^Gjqv^$UE9wEwZQcUE%amp。gTXRm6X4NGpP$vSTTamp。MuWFA5uxY7JnD6YWRrWwc^vR9CpbK!zn%Mz849Gx^Gjqv^$UE9wEwZQcUE%amp。gTXRm6X4NGpP$vSTTamp。849Gx^Gjqv^$UE9wEwZQcUE%amp。ksv*3tnGK8!z89AmYWpazadNuKNamp。qYpEh5pDx2zVkumamp。UE9aQGn8xp$Ramp。6a*CZ7H$dq8KqqfHVZFedswSyXTyamp。gTXRm6X4NGpP$vSTTamp。MuWFA5uxY7JnD6YWRrWwc^vR9CpbK!zn%Mz849Gx^Gjqv^$UE9wEwZQcUE%amp。gTXRm6X4NGpP$vSTTamp。MuWFA5ux^Gjqv^$UE9wEwZQcUE%amp。gTXRm6X4NGpP$vSTTamp。MuWFA5uxY7JnD6YWRrWwc^vR9CpbK!zn%Mz849Gx^Gjqv^$UE9wEwZQcUE%amp。gTXRm6X4NGpP$vSTTamp。849Gx^Gjqv^$UE9wEwZQcUE%amp。gTXRm6X4NGpP$vSTTamp。MuWFA5uxY7JnD6YWRrWwc^vR9CpbK!zn%Mz849Gx^Gjqv^$UE9wEwZQcUE%amp。gTXRm6X4NGpP$vSTTamp。MuWFA5ux^Gjqv^$UE9wEwZQcUE%amp。gTXRm6X4NGpP$vSTTamp。MuWFA5uxY7JnD6YWRrWwc^vR9CpbK!zn%Mz849Gx^Gjqv^$UE9wEwZQcUE%amp。gTXRm6X4NGpP$vSTTamp。849Gx^Gjqv^$UE9wEwZQcUE%amp。QA9wkxFyeQ^!djsXuyUP2kNXpRWXmAamp。 內部資料 請勿外傳 項 目 經(jīng) 理項 目 副 經(jīng) 理 項 目 總 工 質 安 總 監(jiān)工程管理部物資管理部技術管理部檢測試驗室質安管理部監(jiān) 督 工 程 管 理部 、 物 資 管 理部 、 檢 測 試 驗 室現(xiàn) 場 質 檢 員 、 施 工 員施 工 班 組 9JWKffwvGtYM*Jgamp。gold 序列 39。)。 ylabel(39。n39。%產(chǎn)生 gold 序列 stem(goid_sequence_1)。%進行摩 2加 goid_sequence_1=s。 for n=r+1:N。%進行循環(huán) s2(n)=mod(sum(s2(nr:n1).*s2(1:r)),2)。%產(chǎn)生 m序列 end。%移位寄 存器的長度和序列長度； s1(1:6)=[1 0 0 0 0 1]%initial value 1 s2(1:6)=[1 0 0 0 0 0]%initial value 1 f1=[1 0 0 0 0 1 1]%特征多項式 f(x)=x^6+x+1 f2=[1 1 0 0 1 1 1]%特征多項式 f(x)=x^6+x^5+x^2+x+1 for n=r+1:N。 附錄 B Gold 序列程序 r=6。m 序列的自相關函數(shù) 39。 axis([10 10 70])。%調整自相關圖像顯示的橫坐標范圍 rho=[fliplr(rho(2:L)),rho]。 rho(j+1)=sum(m.*bwb)。 bwb=zeros(1,L) xg=2*m1。 pie(num,c) title(bt)。)。 title(39。y39。)。 %subplot(3,1,2) stairs(m)%對 m序列繪圖 46 axis([1 50 ]) xlabel(39。)。)。%創(chuàng)建百分號字符 串數(shù)組 c=strcat(name,percent39。%39。%轉化為字符串數(shù)組 %在每個標簽后加 2 個空格 for i=1:length(name) name(i)={[name{i},blanks(2)]}。)。%轉化為字符型 percent=cellstr(percent)。%計算百分比 percent=num2str(percent39。%標題 bili=num/L。0 的數(shù)量 39。1 的數(shù)量 39。}。,39。%輸入數(shù)據(jù) name={39。 zero_num=Lone_num。%移位后的寄存器 m(i)=an(1)。%移位 an1(len)=mod(sum(.*an),2)。%初始寄存器內容 m=zeros(1,L)%提速 m(1)=an(1)。 %所需的移位寄存器的長度 L=2^len1。 % 為移位寄存器 ,從左向右由大到小 an= [ 0 0 0 0 0 0 1 ]。偽隨機序列中本原多項式生成算法 [J],計算機工程。感謝你們在這次過程中的幫助，使我獲益匪淺。在整設計過程中，其實我也遇到了很多問題，最終都能解決這取決與同學和 老師的幫助，所以在此，我非常感謝我的指導老師 ，他在整個過程中給予了我們最大的幫助，為我們解決了很多問題，使這次課程設計完滿成功。 40 simulink 分析 （ 1）模塊圖 圖 57 gold 序列模塊圖 （ 2）結果 圖 57Gold 序列模塊分析結果圖 上圖是由兩個 7階 m序列經(jīng)過摩 2加形成的 Gold 序列，由上圖看出 Gold 序列的波形 為隨機出現(xiàn)的 0與 1組成。 Gold 的仿真及分析 程序分析 Gold 序列是由兩個 m序列構成的優(yōu)選對本原多項式經(jīng)過摩 2 加生成的，本文以 f(x)= +x+1 和 f(x)= + + +x+1 兩個本原多項式構成 Gold 序列。 1 0 8 6 4 2 0 2 4 6 8 10010203040506070m 序列的自相關函數(shù) 38 圖 53 m 序列自相關函數(shù)圖 上圖是自相關函數(shù)圖，體現(xiàn)了 m序列的自相關性。 37 1 的概率 5 0 3 9 3 7 0 0 7 8 7 4 0 1 5 . 8 %0 的概率 4 9 6 0 6 2 9 9 2 1 2 5 9 8 4 . 3 %1 的數(shù)量 640 的數(shù)量 63 圖 52 m 序列均衡性圖 通過觀察上圖可以知道，產(chǎn)生的序列中 1的個數(shù)為 64， 0的個數(shù)為 63， 1的個數(shù)比 0 的個數(shù)多一。并利用 matlab 中的繪圖函數(shù)繪出 m序列的波形圖以及自相關函數(shù)圖。事實上， MATLAB 自帶的許多內部函數(shù)均是用 C語言編寫并編譯的，因此利用 MATLAB 內部函數(shù)的程序部分運行速度并不比其他語言中相應函數(shù)低 [1011]。 （ 8） MATLAB 的程序執(zhí)行效率比其他語言低。隨著版本的提高，這種能力還會不斷加強，使得人們利用計算機和實際硬件相連接的半實物仿真的夢想得以輕易實現(xiàn)。這也是 MATLAB 廣為流行的重要原因之一。利用 MATLAB 的繪圖功能，可以輕易地獲得高質量的曲線圖。與其它語言相比， MATLAB 對計算機內存、硬盤空間的要求也是比較高的。 35 （ 4） MATLAB 計算精度很高， MATLAB 中數(shù)據(jù)是一雙精度存儲的，一個實數(shù)采用 8字節(jié)存儲，而一個復數(shù)則采用 16 字節(jié)存儲。對于常規(guī)語言中必須使用許多語句才能實現(xiàn)的功能，如矩陣分解、矩陣求逆、積分、快速傅里葉變換，甚至串口操作、聲音的輸入輸出等，在 MATLAB 中用一兩句指令即可實現(xiàn)。 （ 3） MATLAB 語 法簡單，編程風格接近數(shù)學語言描述，是數(shù)學算法開發(fā)和驗證的最佳工具 。使用 MATLAB意味著站在巨人的肩膀上觀察和處理問題，所以在編程效率，程序的可讀性、可靠性和可以執(zhí)行上遠遠超過了常規(guī)的高級語言。 （ 2） MATLAB 是一種超高級語言。采用 MATLAB 編寫的程序可以在目前所有的操作系統(tǒng)上運行（只要這些系統(tǒng)上安裝了 MATLAB 平臺）。在國內外許多著名學術期刊上登載的論文，大部分的數(shù)值結果和圖形都是借助 MATLAB 來完成的 [9]。 目前， MATLAB 已經(jīng)廣泛用于理工科大學從高等數(shù)學到幾乎各門專業(yè)課程之中，成為這些課程進行虛擬試 驗的有效工具。 MATLAB 使得人們擺脫了常規(guī)計算機編程的繁瑣，讓人 34 們能夠將大部分精力投入到研究問題的數(shù)學建模上。但是，要用好 MATLAB 卻不是在短時間就可以達到的。 33 5. MATLAB 仿真實現(xiàn) MATLAB 軟件介紹 MATLAB 是 Math Works 公司開發(fā)的一種跨平臺的，用于矩陣數(shù)值計算的簡單高效的數(shù)學語言，與其它計算機高級語言如 C, C++, Fortran, Basic, Pascal 等相比， MATLAB語言編程要簡潔得多，編程語句更加接近數(shù)學描述，可讀性好，其強大的圓形功能和可視化數(shù)據(jù)處理能力也是其他高級語言望塵莫及的。 Gold 碼也具有三值互相關函數(shù)值，其取值同（ ）式。其中 p=2n1,p 為 Gold 碼序列的周期。 (2)自相關特性： Gold 證明了 Gold 碼序列的自相關函數(shù)的所有非最高峰的取值是三值。當 n 為偶數(shù)，但不能被 4整 除時 ,在周期 N=2n1的 N+2個 Gold序列中 , 平衡碼占 75%,非平衡碼占 25%。當 n 為奇數(shù)時 ,在周期 N的 N+2個 Gold 序列中，有 2n1個序列是平衡的。 ③“ 1”碼元過少。 Gold 序列的平衡特性有 3種 ,也就是 Gold 序列有 3種“ 0”和“ 1”情況 : ①“ 1”碼元數(shù)目僅比“ 0”碼元數(shù)目多一個 ,這就是平衡 Gold 序列。再加上兩個 m序列，共有 2n+1 個 Gold 序列碼。如有兩個 m 序列 ,它們的互相關函數(shù)的絕對值有界，且滿足以下條件： ??????????，12,12)(2221nnR ?的倍數(shù)不是為偶數(shù)，為奇數(shù) 4nnn （ ） 我們稱這一對 m序列為優(yōu)選對。 Gold 序列的產(chǎn)生原理 圖 41 Gold 序列發(fā)生器 Gold 序列就是為了解決 m序列個數(shù)不多且 m序列之間的互相關函數(shù)值不理想而提出的，它是用一對周期和速率均相同的 m序列優(yōu)選對模 2加后得到的。 Gold 提出了一種基于 m 序列的碼序列，稱為 Gold 碼序列。 序列 m序列雖然性能優(yōu)良，但同樣長度的 m 序列個數(shù)不多，且序列之間的互相關值并不都好。 R(j)是一個周期函數(shù)，即 )()( kpjRjR ?? ，式中， k=1,2,? , p=(2n1)為周期。 另外由 m 序列的均衡性可知， 在一個周期中 0 比 1 的個數(shù)少一個， 故得 AD=1(j 為非零整數(shù)時 )或 p(j 為零時 )。經(jīng)過 j次移位后， m 29 序列為 pjjjj aaaa ???? , 321 ?，其中 ipi aa ?? (以 p 為周期 )，以上兩序列的對應項相乘然后相加， 利用所得的總和： ?? ????? ????????pi ijipjpjjjaaaaaaaaaa1332211? （ 22） 來衡量一個 m序列與它的 j次移位序列之間的相關程度，并把它叫做 m序列( paaaa , 321 ? )的自相關函數(shù)。在 m序列中，常常用 +1代表 0，用 1代表 1。 例如， pM =0 0 0 1 1 1 1 0 1 0 1 1 0 0 1? pM 延遲兩位后得 rM , 再模二相加 rM sM = pM + rM