通信原理實驗報告-LABVIEW2.1要點

1、通信原理實驗報告實驗項目：實驗1掌握用Labview 產(chǎn)主隨機(jī)數(shù)的方法實驗2 統(tǒng)計隨機(jī)數(shù)的概率分布密度函數(shù)及相關(guān)函數(shù)特性實驗 3 產(chǎn)生 m 序列信號源，驗證m 序列的偽隨機(jī)性以及偽隨機(jī)序列的自相關(guān)函數(shù)的雙值特性。實驗 4 模擬產(chǎn)生 AWGN 及 ISI 信道，添加到數(shù)字通信仿真系統(tǒng)中實驗1隨機(jī)數(shù)產(chǎn)生及直方圖統(tǒng)計實驗1隨機(jī)數(shù)產(chǎn)生及直方圖統(tǒng)計一、實驗?zāi)康模?）掌握在一般微型計算機(jī)上產(chǎn)主隨機(jī)數(shù)的方法。（2）統(tǒng)計隨機(jī)數(shù)的概率分布密度函數(shù)。二、實驗內(nèi)容1 .用計算機(jī)產(chǎn)生0, 1均勻分布的（偽）隨機(jī)數(shù)。2 .由0, 1均勻分布隨機(jī)數(shù)產(chǎn)生其它分布的隨機(jī)數(shù)，例：正態(tài) N （0, 1）分布的隨機(jī)數(shù)。3 .用直方

2、圖統(tǒng)計隨機(jī)數(shù)的分布密度。三、實驗設(shè)備微型計算機(jī)及其高級程序語言編譯環(huán)境，例C+、FORTRAN、PASCAL等，也可以應(yīng)用工程計算工具軟件如 MATLAB等。四、實驗原理1 .計算機(jī)產(chǎn)生均勻分布隨機(jī)數(shù)在計算機(jī)算法中，為實現(xiàn)方便，通常使用偽隨機(jī)數(shù)（序列）來代替（真）隨機(jī)數(shù)。偽隨 機(jī)序列是有周期性的數(shù)值序列，當(dāng)其周期N相對很大時，統(tǒng)計特性一定程度上逼近隨機(jī)序列，故效果與（真）隨機(jī)數(shù)相近。2 .高斯分布隨機(jī)數(shù)的獲得實際研究當(dāng)中，高斯（正態(tài)）分布是經(jīng)常被使用到的數(shù)學(xué)模型，可以近似描述很多隨機(jī) 事件的統(tǒng)計特性。，我們可以采用非線性變換法，對比較容易產(chǎn)生的均勻分布隨機(jī)序列進(jìn)行變 換，（近似）得到高斯分布隨

3、機(jī)序列。_1/2_ _Xc =（-2ln Ri） cos2二R2公式中，若Ri和R2是0, 1區(qū)間兩個均勻分布隨機(jī)變量，理論上可以證明 Xc是標(biāo)準(zhǔn)正態(tài)分布（均值為0,方差為1的高斯分布）的隨機(jī)變量。3 .直方圖對于一個隨機(jī)變量，假如我們知道它是正態(tài)的（或其它分布形式）我們可以從隨機(jī)變量 的抽樣估計它的均值和方差，從而得到它的分布密度函數(shù)。預(yù)先對一個隨機(jī)變量分布一無所知，要估計它的分布密度函數(shù)可借助于直方圖統(tǒng)計方法：設(shè)有圖1所示密度函數(shù)fx(x)把隨機(jī)變量X的取值量化，量化階為 2e ,例如對于以 x=2 為中心的量化階內(nèi)，如果 e足夠小。有2 fx(2) : P2 一 ；：二 X :二 2 ;

4、當(dāng)X抽樣Xi ,X2,Xn ,是獨(dú)立且同分布時，則X1,X2，,XN中落在(2-%2 +后)中的個數(shù)2 ；N當(dāng)N較大能用頻率代替概率。所以可以得到概率密度的估計值fx(2):Xk在(2 -8,2+勤中的個數(shù)2 Nk=1, - - ,N5顯然x= 2并不特殊，可以重復(fù)進(jìn)行其它 x值上的密度估計。這樣我們就得到密度函數(shù) fx(x) 的離散估計值。這就是直方圖方法。這種直方圖估計的準(zhǔn)確度可以進(jìn)行專門分析。分析它的方差可以知道，N和£的選擇有一定的要求。顯然，數(shù)據(jù)樣點數(shù) N越大越好，而區(qū)間密度 £選擇與N值相適應(yīng)為女?。有表 1 可供參考。表1樣點數(shù)與直方單元數(shù)N500100020

5、00500010000K2230395674N-數(shù)據(jù)樣點數(shù)；K-隨機(jī)變量取值范圍量化的單元數(shù)。實驗中要求用計算機(jī)對所產(chǎn)生的隨機(jī)數(shù)用直方圖估計其概率密度函數(shù)，井打印出分布函 數(shù)的圖形五、實驗步驟1 .弄懂實驗原理，設(shè)計結(jié)構(gòu)框圖:2 .利用LabVIEW編寫程序:產(chǎn)生隨匯股列長度N隨機(jī)變量宜值定圍量化的聿元激KIdblB六、實驗結(jié)果0.493410.7370 0,71O6c0.06021 0.D015(0.47375直方互幅直實驗2相關(guān)噪聲模型和相關(guān)函數(shù)計算12實驗2相關(guān)噪聲模型和相關(guān)函數(shù)計算一、實驗?zāi)康氖煜は嚓P(guān)噪聲模型，掌握相關(guān)函數(shù)的計算方法。二、實驗內(nèi)容1 .計算一下白噪聲模型的相關(guān)函數(shù)。2

6、.計算相關(guān)噪聲的相關(guān)函數(shù)，它應(yīng)是指數(shù)型的。3 .計算正弦信號加噪聲的相關(guān)函數(shù)，期望得到指數(shù)與余弦迭加的波形。三、實驗設(shè)備微型計算機(jī)及其它外圍設(shè)備。四、實驗原理各態(tài)遍歷的平年I隨機(jī)過程x（t）自相關(guān)函數(shù)可表示為1 Tr（）寸可0 x（t）x（t ）出而對于周期性信號，自相關(guān)函數(shù)可表示為R()1 T00 x(t)x(t )dtT0其中T0為周期。自相關(guān)函數(shù)的離散數(shù)值計算公式1 N _rR（S =Z xkxk +r=0,1,2,，m m < NN -r 口如果我們有N個數(shù)據(jù)記錄，上式就可以在數(shù)字計算機(jī)上進(jìn)行相關(guān)函數(shù)的估計。實驗內(nèi)容中的三個模型：1 .白噪聲模型：產(chǎn)生的0,1均勻分布隨機(jī)序列，都

7、強(qiáng)調(diào)要求它的分布均勻性要好，獨(dú)立性要強(qiáng)。這獨(dú)立性11要求指的是前后相鄰數(shù)據(jù)不相關(guān)。實際上這就是一種白噪聲模型，如果把它變換為-彳，2 均勻分布，則就是一種零均值的白噪聲模型。我們記理想的白噪聲序列為 w©,則其自相關(guān)函數(shù)可表示為2、Rw (i, j) = E Wi, Wj = ； w、j計算機(jī)產(chǎn)生的均勻分布隨機(jī)數(shù)具有白噪聲性質(zhì)，就應(yīng)該有尖峰自相關(guān)函數(shù)。反過來說, 這相關(guān)函數(shù)的尖峰形狀能用來檢驗隨機(jī)數(shù)獨(dú)立性的好壞。2 .相關(guān)噪聲模型Wi f +nf niaa 單位延時對于相應(yīng)的數(shù)據(jù)序列來說，可由白噪聲模型變換出相關(guān)噪聲的模型，白噪聲Wi,經(jīng)圖示裝置處理，產(chǎn)生的輸出為 ni =wi +a

8、ni（0<a<i）我們看到6與其前面幾個值（ni-i,ni-2,，）有依賴關(guān)系。這相關(guān)性的強(qiáng)弱取決于系數(shù)a的取值，這樣得到的數(shù)據(jù)序列 Q就是一個相關(guān)噪聲模型。相關(guān)噪聲的相關(guān)函數(shù)應(yīng)是指數(shù)型的。3 .正弦信號加噪聲這是加噪聲模型，用一個正弦信號與相關(guān)噪聲迭加，得到xi =ni bcos（2：i/T）其中正弦信號初相為零，當(dāng)然也可隨意指定一個初相，這隨機(jī)相位正弦信號的數(shù)據(jù)樣值 只要利用函數(shù)賦值就可得到。正弦信號加噪聲的相關(guān)函數(shù)應(yīng)該是指數(shù)與余弦迭加的波形。五、實驗步驟1 .弄懂實驗原理，設(shè)計結(jié)構(gòu)框圖:2 .利用LabVIEW編寫程序（由于程序框圖太大，這里分成部分貼出）（a）白噪聲及其自相

9、關(guān)函數(shù)京.花后三E數(shù)N(b)相關(guān)噪聲及正弦信號與相關(guān)噪聲迭加數(shù)據(jù)樣本總數(shù)nEJ正弦信號與柜關(guān)噪聲送比伙。(c)相關(guān)噪聲的自相關(guān)函數(shù)及正弦信號與相關(guān)噪聲迭加的自相關(guān)函數(shù)六、實驗結(jié)果數(shù)據(jù)N力叫:什川口川川川干川?，端樣本（自曝聲前3Kj556la2233t4SgR525| 0,0783!白然聲言相關(guān)函數(shù)的離散數(shù)值序列。-0-0189 0.006點 0,0545 -0,1621 M3S4C曲蛙0白噪聲自相關(guān)國教系數(shù)3 (0.8<a<0.9i幅關(guān)應(yīng)言模型數(shù)據(jù)序?如0.044610.0604'，相關(guān)噪聲主相關(guān)函數(shù)0350J-0.25-hiiii 02=01-0.05-曲線口 0510

10、152025303540時間通信號周期T 110正信號與相美賓聲送枕（蝴時間實驗3常量信號檢測的計算機(jī)模擬（新版）13S111111w rIoIoIo1IoIoNfflL J-C.261CC.1678C-C.36981-0,17010.1937-0.4155尸-C.261CC.1678C-C.3698C.8298C0.1937；-0.4165Pevmin-vmaxFC.32910.3239；O.31S3；10,3129C.3C72：O.3O13E最佳學(xué)快門限電平vt |0.63Pmvmin-vmaxPe Min.1233；F0.6501C.6S86EC.627CEC.615C.6C35O.59

11、12tFHfvmin-vmax)F0.0087：。，頤 IfD.0096fC.CICot0.0109；10,0115Pdvmin-vmax；F0.34990.3613；0.3729；0.3345t0.39&4E0.40&7rFFFXV圖P(H1)。，力均值looooei均方差 均方根 d0.1997104469-2.2361曲線口時間實驗4偽隨機(jī)序列產(chǎn)生及其特性研究實驗3偽隨機(jī)序列產(chǎn)生及其特性研究一、實驗?zāi)康? . 了解偽隨機(jī)序列產(chǎn)生的方法，觀察其變換的不同碼型。2 .研究m序列本原多項式與線性移位寄存器的反饋方式間的關(guān)系。3 .驗證m序列的偽隨機(jī)性。4 .驗證偽隨機(jī)序列的自相

12、關(guān)函數(shù)的雙值特性。二、實驗內(nèi)容1 .選擇合適的m序列本原多項式，設(shè)計n級(n=38)線性反饋移位寄存器， 產(chǎn)生m序列。2 .分析記錄m序列的周期P與級數(shù)n之間的關(guān)系是否符合 P = 2n-1。3 .討論m序列的性質(zhì)和相關(guān)函數(shù)特性。三、實驗設(shè)備1 .直流穩(wěn)壓電源2 .示波器3 .單片計算機(jī)實驗電路裝置四、實驗原理1. m序列在通信系統(tǒng)中，為了研究隨機(jī)噪聲對系統(tǒng)的影響，往往需要人為生成隨機(jī)噪聲。而在20世紀(jì)60年代，人們發(fā)明了 “偽隨機(jī)噪聲”，才真正滿足了研究的需要。偽隨機(jī)噪聲具有類似于隨機(jī)噪聲的某些統(tǒng)計特性，同時又能夠重復(fù)產(chǎn)生，避免了隨機(jī)噪聲不可重現(xiàn)的缺點，因而 獲得了廣泛的應(yīng)用。而偽隨機(jī)噪聲是

13、由周期性數(shù)字序列經(jīng)濾波等處理后得到的，這種周期性 數(shù)字序列就是“偽隨機(jī)序列”，有時也被稱作偽隨機(jī)信號或偽隨機(jī)碼。至今，最廣為人知的二進(jìn)制偽隨機(jī)序列是“最長線性反饋移位寄存器序列”，簡稱m序列。m序列因其隨機(jī)特性和預(yù)先可確定性及可重復(fù)實現(xiàn)的優(yōu)點，在實際領(lǐng)域中得到廣泛應(yīng)用。2. m序列產(chǎn)生原理m序列是由帶線性反饋的移存器產(chǎn)生的周期最長的序列。圖4-1為一般的線性反饋移存器產(chǎn)生m序列的原理飯框圖。圖4-1線性反饋移位寄存器原理方框圖由于帶有線性反饋，在移位脈沖作用下，線性移位寄存器各級的狀態(tài)將不斷變化，通常移位寄存器的最后一級做輸出，當(dāng)移位寄存器的級數(shù)及時鐘一定時，輸出序列就由初始狀態(tài)和反饋邏輯完全

14、確定。輸出序列為ak = a0alan，是一個周期序列。經(jīng)一次移位線性反饋，左端新得到的輸入an為nan =Cinn5 c2anN 55 Cna0 =£ Qan-L (模 2) i 1因此，一般說來，對于任意一個輸入ak,有nak = '，Gak_ti 4上式稱為“遞推方程”，它給出移位輸入與移位前各級狀態(tài)的關(guān)系。按照遞推方程，就可 以用軟件產(chǎn)生 m序列。當(dāng)初始狀態(tài)為全零狀態(tài)是，移位寄存器輸出全零序列，因此，因避免 設(shè)置全0電路。在遞推方程中，G的取值決定了序列的結(jié)構(gòu)，所以G是一個很重要的參量。G的取值情況可以用特征多項式(或特征方程)f (x)描述：nf (x) = co

15、c1x . cnxn 二、Gxii=0對于n階m序列，最大能產(chǎn)生的周期為 2n -1o其中能滿足此要求的最簡 f (x)也稱“本源多項式”。其中38階序列的本源多項式如表 4-1所示。表4-1本源多項式表級數(shù)n8進(jìn)制表小f (x)代數(shù)式313x3+x+1423x4+x+1545x5+x2+16103x6+x+17211x7+x3+18435x8+x4+x3+x2+13. m序列的性質(zhì)m序列具有幾個有趣的性質(zhì)和統(tǒng)計特性。1 .對n級m序列周期為P =2n -12 .在一個周期中，“0"出現(xiàn)2n" 1次，“1”出現(xiàn)2nl次3 .在一個周期中，共有 2n.個游程，其中:1人長度為

16、k的游程有個，1<k<n2；2k長度為n1的有1個“0”游程；長度為n的有1個“1”游程。4 .歸一化自相關(guān)函數(shù)為1（j=°）-1/P （1 < j| < P 1）且Wj）的周期為P五、設(shè)計要求1 .利用單片機(jī)編程實現(xiàn)，設(shè)計n級（n=38） m序列的發(fā)生器，每級可選擇實現(xiàn)1種序列碼型。2 .序列碼元速率：（推薦）100Baud 。3 .為便于示波器觀測 m序列波形，給出序列周期同步信號作為示波器觸發(fā)源，如圖 4-2 示意。以n=3為例：m序列一m m序列輸出11 1 1 0 0日01 1 1 1 0 0目0發(fā)生器序列周期同步脈沖 nn圖4-2 m序列發(fā)生和序列

17、周期同步信號六、實驗內(nèi)容1 .m序列生成本次試驗由MCS51系列的AT89S52單片機(jī)完成 m序列機(jī)器周期同步脈沖的生成。其主 要部分的程序設(shè)計如下： a） m序列的移位寄存器生成法m序列由移位寄存器法產(chǎn)生（詳細(xì)原理見上面“實驗原理”部分） ，反饋環(huán)路的設(shè)置采 用了表4-1中的本源多項式的取值方法。 每當(dāng)定時器T0中斷輸出m序列的一個二進(jìn)制位， 便把全局變量v置零。在主程序中，每當(dāng)檢測到v的值為零，則根據(jù)當(dāng)前要求的階數(shù)移位產(chǎn)生新的位，將 v置1,并等待定時器中斷輸出。為保證在38階都可生成，并避免移位寄存器出現(xiàn)全零狀態(tài)，移存器（實際上是一個數(shù)組）的初值設(shè)置為1,0,1,0,1,0,1,0, n

18、階m序列則取其前n個作為m序列的初值。（詳見“程序清單”主程序中無限循環(huán)部分） b） 定時中斷輸出m序列我們選擇m序列的波特率為 400Baud,而單片機(jī)晶振為12.0592MHz ,所以定時器的計 數(shù)量為 pwm_time=11059200/12/400=2304。我們采用定時器 T0作為m序列輸出定時器，工作于模式 1 （即16位定時器）。中斷服 務(wù)程序如下：timer0（） interrupt 1 using 1/T0 中斷，發(fā)送 m 序列 EA = 0;TH0 = vth0;TL0 = vtl0;TR0 = 1;pwm_out=pwm_status;sync=sync_sta;v=0;

19、EA = 1;2 .實驗結(jié)果及m序列性質(zhì)驗證a）實驗結(jié)果n = 3級m序列 n = 4級m序列 n = 5級m序列 n = 6級m序列 n = 7級m序列 n = 8級m序列b）序列性質(zhì)1.M序列的周期實驗觀察到的m序列為:1010011101011001000111101011001000111101011001000111101010110011011101101001001110-101010100110011101110100101100-101010100110010111010001110110-由表一可得生成的M序列的周期p=2n -1 ,如n=3時，p=7; n=4時，p=15

20、2. M序列的游程4 1由上表所得的M序列所得，當(dāng)n=4時，共有2=8個游程，分別為1111, 0, 1, 0, 11,00, 1, 000;其中長度為4的游程有1個（1111）,長度為3的游程有1個（000）,長度為2 的游程有2個（11, 00）,長度為1的游程有4個（0, 1, 0, 1）綜上長度為1 （k=1）的游程占了總游程數(shù)的一半（1/2）,長度為2 （k=2）的游程占了11長度為3 （n-1）的游程是連續(xù)0的游程（000）,長度為4 （n）的游程為連續(xù)1的2* k 4游程（1111）。同理分析n=3, 5, 6, 7, 8時得到結(jié)論：m序列的一個周期（p=2n-1）中，游程總數(shù)為

21、2n-1。其中長度為1的游程個數(shù)占游程總數(shù) 的1/2;長度為2的游程個數(shù)占游程總數(shù)的 1/22=1/4;長度為3的游程個數(shù)占游程總 數(shù)白1 1/23=1/8;一般地，長度為k的游程個數(shù)占游程總數(shù)的1/2k=2-k,其中1wkw（n-2）。而且，在長度為k游程中，連1游程與連0游程各占一半，長為（n-1）的游程是連0游程，長為n的游程是連1游程。七、程序清單#include <absacc.h>#include <reg51.h>/*/7/*/常量定義/*/*/*/*定時器初值*/#define V_TH1 0xee /* 5ms , ! notice: crystal

22、11.0592MHz */#define V_TL1 0x00/*數(shù)碼管段驅(qū)動寄存器地址*/#define ADDR_8SEG XBYTE0x2000/*數(shù)碼管位驅(qū)動和指示燈驅(qū)動寄存器地址*/#define ADDR_SEL XBYTE0x4000/*按鍵*/sbit KEY1=P1A0;sbit KEY2=Pil；sbit KEY3=Pi2;sbit KEY4=Pi3;/*/*/*/變量定義/*/*/*/*數(shù)碼管位驅(qū)動和指示燈驅(qū)動信號輸出緩存，定義了一個可位尋址的變量 */ unsigned char bdata output_sel;sbit led_1 = output_se|A5;sb

23、it led_2 = output_se|A6;sbit led_3 = output_se|A7；sbit led_4 = output_se|A4;/*數(shù)碼管掃描驅(qū)動指針，為測試外部存儲器（ U3 6264）,特使用xdata類型*/ unsigned char xdata digi_scaner;/*測試用計數(shù)值十進(jìn)制表示，為測試外部存儲器（ U3 6264）,特使用xdata類型*/ unsigned char xdata digi4=0,0,0,3;unsigned char temp4;/按鍵消抖/*pwm 控制的變量*/unsigned int pwm_time;unsigned

24、 char vth0,vtl0;bit pwm_status;/ 信號bit sync_sta;/同步int power;/2的乘方的結(jié)果，等于 2的n次方unsigned char j;/用于 for 循環(huán)int prd;/period ,由于序列周期的計數(shù)sbit pwm_out = P3A5;/m 序列輸出sbit sync = P3A4;序列周期同步信號 unsigned char v=0;/ 發(fā)送完一位置零/*/*/*/函數(shù)定義/*/*/*/* 7段數(shù)碼顯示譯碼參數(shù)：DATA :需要顯示的數(shù)字或符號；返回彳1:7段譯碼結(jié)果(D7。=PGFEDCBA )*/unsigned char

25、NUMTOSEG7(unsigned char DA TA) unsigned char AA;switch (DA TA) case 0: AA=0xc0;break; /*' 0' */case 1: AA=0xf9;break; /*' 1' */case 2: AA=0xa4;break; /*'2' */case 3: AA=0xb0;break; /*'3'*/case 4: AA=0x99;break; /*'4'*/case 5: AA=0x92;break;/*'5'*/case

26、6: AA=0x82;break;/*'6'*/case 7: AA=0xf8;break;/*7*/case 8: AA=0x80;break; /*'8'*/case 9: AA=0x90;br eak; /*' 9' */case 10: AA=0x88;break; /*' */ 'Acase 11: AA=0x83;break; /*' B'*/case 12: AA=0xc6;break; /*' C'*/case 13: AA=0xa1;break; /*' D'*/c

27、ase 14: AA=0x86;break; /*' E'*/case 15: AA=0x8e;break; /*' F'*/case '-':AA=0xbf;break;/*破折號，此處原誤為 0xdf,系04級王資凱同學(xué)指正*/case '_':AA=0xf7;break; /* 下劃線 */case '':AA=0xff;break;/* 消隱*/default: AA=0xff;return(AA);/* T1時鐘中斷服務(wù)程序在本例中，每5ms被執(zhí)行一次*/time)interrupt 3 using 0

28、EA=0; /*關(guān)中斷*/*重新對計數(shù)器賦初值，并啟動定時計數(shù) */TH1=V_TH1;TL1=V_TL1;TR1=1;EA=1;output_sel = 0xf0;/*初值，令數(shù)碼管驅(qū)動位無效，指示燈全滅*/*按鍵4消抖和按鍵4確認(rèn)的處理*/if (KEY4=0) led_4=0;temp4+;if(temp4>20 && KEY4=1)temp4=0;if(+digi3>8) digi3=3;power=1;prd=1;for(j=0;j<digi3;j+) power*=2; */*數(shù)碼管掃描驅(qū)動指針值從1到4重復(fù)變換，每5ms間隔對一個數(shù)碼管進(jìn)行驅(qū)動，

29、20ms 一個輪回if (+digi_scaner>=5) digi_scaner = 1;switch (digi_scaner) case 1: /*驅(qū)動第一個數(shù)碼管*/output_sel |= 0x01;ADDR_8SEG = NUMTOSEG7(pwm_status);/* 輸出到鎖存器 U5 (在電路圖中找)*/ break;case 2: /*驅(qū)動第二個數(shù)碼管*/output_sel |= 0x02;ADDR_8SEG =0xff; /* 輸出到鎖存器 U5*/ break;case 3: /*驅(qū)動第三個數(shù)碼管*/output_sel |= 0x04;ADDR_8SEG =

30、 0xff; /* 輸出到鎖存器 U5*/ break;case 4: /*驅(qū)動第四個數(shù)碼管*/output_sel |= 0x08;ADDR_8SEG = NUMTOSEG7(digi3); /* 輸出到鎖存器 U5*/ break;ADDR_SEL = output_sel;/*輸出到鎖存器 U6 (在電路圖中找)*/timer0() interrupt 1 using 1/T0 中斷，發(fā)送 m 序列 EA = 0;TH0 = vth0;TL0 = vtl0;TR0 = 1;pwm_out=pwm_status;sync=sync_sta;v=0;EA = 1;/* 主程序 */main(

31、)int n;int i;int count=0;int temp;int yd8=1,0,1,0,1,0,1,0;/*單片機(jī)系統(tǒng)初始化*/P1=0XFF;P3=0XFF;/*中斷工作寄存器初始化*/IE=0X80;/EA=1IP=0;PT0=1;PT1=0;/*定時器工作寄存器初始化*/TMOD=0X11;/T0/1 work as 16b timerTCON=0;/*定時器1的寄存器初始化*/TH1=V_TH1;TL1=V_TL1;/*全局變量賦初值*/digi_scaner=1;/ 數(shù)碼管計數(shù)pwm_time=2304;/ 機(jī)器周期數(shù)，對應(yīng) 400Hzvth0=(0xffff-pwm_t

32、ime)>>8;vtl0=(0xffff-pwm_time);prd=1;power=8;TH0 = vth0;TL0 = vtl0;/*開定時器，定時器0中斷允許*/TR0=1;TR1=1;ET0=1;ET1=1;output_sel=0xf0;temp4=0;*/*主循環(huán)，本例中，在 T1中斷服務(wù)程序未被執(zhí)行的空余時間里，處理機(jī)在以下程序中不斷循環(huán) while (1)n=digi3-1; 隨機(jī)碼的級數(shù)-1if(n=2|n=3|n=5|n=6)if(v=0)if(prd=1) sync_sta=1; elsesync_sta=0;prd+=1;if(prd>=power)prd=1;temp =(yd0+ydn)%2; 取最低位與修高位異或運(yùn)算for(i=0; i<n; i+) ydi = ydi+1;其余數(shù)字全部向前移動一位ydn= temp; /新的數(shù)據(jù)加入隊列最末pwm_status=temp