通信原理实验报告材料-LABVIEW2.1

1、通信原理实验报告实验项目：实验1掌握用Labview产主随机数的方法实验2统计随机数的概率分布密度函数及相关函数特性实验3 产生m序列信号源，验证 m序列的伪随机性以及伪随机序列的自相关函数的双 值特性。实验4 模拟产生AWGN及 ISI信道，添加到数字通信仿真系统中实验1随机数产生及直方图统计、实验目的（1）掌握在一般微型计算机上产主随机数的方法。（2）统计随机数的概率分布密度函数。、实验内容1. 用计算机产生0,1均匀分布的（伪）随机数。N （ 0, I ）分布的随机数。2 .由0 , 1均匀分布随机数产生其它分布的随机数，例：正态3 用直方图统计随机数的分布密度。二、实验设备C+、FOR

2、TRAN PASCAL等，也可以应用工程微型计算机及其高级程序语言编译环境，例 计算工具软件如 MATLAB等。四、实验原理1. 计算机产生均匀分布随机数在计算机算法中，为实现方便，通常使用伪随机数（序列）来代替（真）随机数。伪随机序列是有周期性的数值序列，当其周期N相对很大时，统计特性一定程度上逼近随机序列，故效果与（真）随机数相近。2. 高斯分布随机数的获得实际研究当中，高斯（正态）分布是经常被使用到的数学模型，可以近似描述很多随机 事件的统计特性。，我们可以采用非线性变换法，对比较容易产生的均匀分布随机序列进行变 换，（近似）得到高斯分布随机序列。1/2Xc =（-21 n R1） CO

3、S2二R2公式中，若 R和F2是0 , 1区间两个均匀分布随机变量，理论上可以证明Xc是标准正态分布（均值为0，方差为1的高斯分布）的随机变量。3. 直方图对于一个随机变量，假如我们知道它是正态的（或其它分布形式）我们可以从随机变量 的抽样估计它的均值和方差，从而得到它的分布密度函数。预先对一个随机变量分布一无所知，要估计它的分布密度函数可借助于直方图统计方法：设有图1所示密度函数fx(x)把随机变量X的取值量化，量化阶为2 例如对于以x=2为中心的量化阶内，如果&足够小。有2 ;fx(2) ： P2 ； ： X ： 2；fx(2)Xk在(2 - ；,2 ；)中的个数当X抽样Xi，爼,Xn ，

4、是独立且同分布时，则P2 -.22,.人中洛在()中的个数2名N当N较大能用频率代替概率。所以可以得到概率密度的估计值显然x= 2并不特殊，可以重复进行其它 x值上的密度估计。这样我们就得到密度函数 fx(x) 的离散估计值。这就是直方图方法。这种直方图估计的准确度可以进行专门分析。分析它的方差可以知道，N和&的选择有一定的要求。显然，数据样点数 N越大越好，而区间密度 &选择与N值相适应为好。有表 1 可供参考。表1样点数与直方单元数N50010002000500010000K2230395674N-数据样点数；K-随机变量取值范围量化的单元数。实验中要求用计算机对所产生的随机数用直方图估计

5、其概率密度函数，井打印出分布函 数的图形五、实验步骤1.弄懂实验原理,设计结构框图:2.利用LabVIEW编写程序:产生随f咸见民度N产生旳正态N ( 0 J )分布松前匸蕊歹人b UHistogram Graph 2Histogram Graph建机受星靈值矩围量化笳单元数KI蓟鼻六、实验结果； :10.06021 10.0015( 10.4934 1()*6126； 10.4737 10.7370= |0-71D6：直方图直方囹计数齐00.20.40.60.B1实验2相关噪声模型和相关函数计算、实验目的熟悉相关噪声模型，掌握相关函数的计算方法。、实验内容1. 计算一下白噪声模型的相关函数。

6、2. 计算相关噪声的相关函数，它应是指数型的。3. 计算正弦信号加噪声的相关函数，期望得到指数与余弦迭加的波形。二、实验设备微型计算机及其它外围设备。四、实验原理 各态遍历的平稳随机过程x(t)自相关函数可表示为1 TR( HTim- 0 x(t)x(t )dt而对于周期性信号，自相关函数可表示为1 JoR()To 0x(t)x(t x)dt其中T)为周期。自相关函数的离散数值计算公式r =0,1,2,m m N1N_rR( )-XkXkTN -r心如果我们有N个数据记录，上式就可以在数字计算机上进行相关函数的估计。实验内容中的三个模型：1 .白噪声模型：产生的0,1均匀分布随机序列，都强调要

7、求它的分布均匀性要好，独立性要强。这独立1 1 性要求指的是前后相邻数据不相关。实际上这就是一种白噪声模型，如果把它变换为卜？ ,均匀分布，则就是一种零均值的白噪声模型。我们记理想的白噪声序列为w，则其自相关函数可表示为2RwQ, j) = EWi ,Wjj其中Ij.0 j计算机产生的均匀分布随机数具有白噪声性质，就应该有尖峰自相关函数。反过来说, 这相关函数的尖峰形状能用来检验随机数独立性的好坏。2 .相关噪声模型对于相应的数据序列来说，可由白噪声模型变换出相关噪声的模型，白噪声W，经图示装置处理，产生的输出为 n, =w a（ oa43至歆日(0.8a0.9)43正弦信号周期T蛊数b (b

8、 = 0.3 04)11相关秦肓檯型:数垢序niP口xi = ni + b*co5(2*pi*i/n； HEH1(c)相关噪声的自相关函数及正弦信号与相关噪声迭加的自相关函数相黃谍言模甦進序列筒ixJS松-K N丐二忑三三丐二丕芦卧TITDE.L7ama口曲 柜关噪岸至宅芫国取h斋取議直亭列柜关噪岸至宅芫国取h斋取議直亭列正弦信寻杠幌肓茅柜关團議六、实验结果款捱“巳嚴N白嗓車刍咱丟函纹的蕃粧史直喜列3 f0.0189| 00064( 0*0545 4H621 018Mm噪芦吕治关圄毅rnc罢数a (0.8*； a 0,9)口口口0.07910.0785 0.04461 0.06044)1035

9、正弦信号周期T110正信号与牢芫嗓雷動咚“|aii67 |ouw| n i. i H a(0077?10.28021 -0.05641正弦信号控噪声导相黄固数曲纬0価;|oIoIo1IoIoNfflNP(H1)均值均方差 沟方根|Oj00061Jo.1997； 1 (0.4469-C.261C C.1678C -C.3698 -C.17C1 C.1987- -C.165FF| -01261010.167SC | -0.369B 10298? 0.1987E -04165FFFFFFPevmin-vmax103294 103239； 10318J1103129 103072； 1030131最佳

10、!l决门限电平vtPe Min(063ai233；Pmvmin-vmaxC.65C1C.6S86E C.627CE C.615 C.6C35 C.5912tF4vmin-vmaxC.DC87： C.OOSlf C.0&96f C.ClCSt C.C1C9： C.C1150.349910.3613；10.3729；C.385t0.39&4E10.40&7Pdvmin-vmaxX图d12.2361ffigEO对闫实验3伪随机序列产生及其特性研究、实验目的1 了解伪随机序列产生的方法，观察其变换的不同码型。2. 研究m序列本原多项式与线性移位寄存器的反馈方式间的关系。3 .验证m序列的伪随机性。4

11、.验证伪随机序列的自相关函数的双值特性。二、实验内容1 .选择合适的m序列本原多项式，设计 n级（n=38）线性反馈移位寄存器，产生 m序列。2 .分析记录m序列的周期P与级数n之间的关系是否符合 P=2n-1。3 .讨论m序列的性质和相关函数特性。三、实验设备1 .直流稳压电源2 .示波器3 .单片计算机实验电路装置四、实验原理1. m序列在通信系统中，为了研究随机噪声对系统的影响，往往需要人为生成随机噪声。而在20世纪60年代，人们发明了“伪随机噪声”，才真正满足了研究的需要。伪随机噪声具有类似 于随机噪声的某些统计特性，同时又能够重复产生，避免了随机噪声不可重现的缺点，因而 获得了广泛的

12、应用。而伪随机噪声是由周期性数字序列经滤波等处理后得到的，这种周期性 数字序列就是“伪随机序列”，有时也被称作伪随机信号或伪随机码。至今，最广为人知的二进制伪随机序列是“最长线性反馈移位寄存器序列”，简称m序列。m序列因其随机特性和预先可确定性及可重复实现的优点，在实际领域中得到广泛应用。2. m序列产生原理m序列是由带线性反馈的移存器产生的周期最长的序列。图4-1为一般的线性反馈移存器产生m序列的原理饭框图。图4-1线性反馈移位寄存器原理方框图由于带有线性反馈，在移位脉冲作用下，线性移位寄存器各级的状态将不断变化，通常 移位寄存器的最后一级做输出，当移位寄存器的级数及时钟一定时，输出序列就由

13、初始状态和反馈逻辑完全确定。输出序列为ak二ao.an.，是一个周期序列。经一次移位线性反馈，左端新得到的输入a n为nan 二 Ci nnd 二 C2anN 三 三 Cnao- Cian(模 2)y因此，一般说来，对于任意一个输入ak，有na k = ci ak _i.i 4上式称为“递推方程”，它给出移位输入与移位前各级状态的关系。按照递推方程，就可 以用软件产生 m序列。当初始状态为全零状态是，移位寄存器输出全零序列，因此，因避免 设置全0电路。在递推方程中，Ci的取值决定了序列的结构，所以Ci是一个很重要的参量。C的取值情况可以用特征多项式(或特征方程)f (x)描述：nf (x) =

14、 co C1x - . qxn 二、qxi=0对于n阶m序列，最大能产生的周期为 2n -1。其中能满足此要求的最简 f (x)也称“本 源多项式”。其中38阶序列的本源多项式如表 4-1所示。表4-1本源多项式表级数n8进制表示f (x)代数式313x3+x+1423x4+x+154552.x +x +16103x6+x+17211x7+x3+18435x +x +x +x +13. m序列的性质m序列具有几个有趣的性质和统计特性。1 对n级m序列周期为P =2n -12在一个周期中，“0”出现2n4 -1次，“1”出现2n4次3 在一个周期中，共有 2n4个游程，其中:2k长度为k的游程有

15、 个，1 . k n - 2 ；长度为n _ 1的有1个“ 0”游程;长度为n的有1个“ 1”游程。4 .归一化自相关函数为1 (j=0)-1/P (1 j P-1)且（j）的周期为P五、设计要求1 利用单片机编程实现，设计n级（n=38） m序列的发生器，每级可选择实现1种序列码型。2 .序列码元速率：（推荐）100Baud。3 .为便于示波器观测 m序列波形，给出序列周期同步信号作为示波器触发源，如图4-2示意。m序列m序列输出以n=3为例:序列周期同步脉冲图4-2 m序列发生和序列周期同步信号六、实验内容1. m序列生成本次试验由MCS5係列的AT89S52单片机完成 m序列机器周期同步

16、脉冲的生成。其主要 部分的程序设计如下：a）m序列的移位寄存器生成法m序列由移位寄存器法产生（详细原理见上面“实验原理”部分），反馈环路的设置采用了表4-1中的本源多项式的取值方法。每当定时器T0中断输出m序列的一个二进制位， 便把全局变量v置零。在主程序中，每当检测到v的值为零，则根据当前要求的阶数移位产生新的位，将 v置1,并等待定时器中断输出。为保证在38阶都可生成，并避免移位寄存器出现全零状态，移存器（实际上是一个数组）的初值设置为101,0,1,0,1,0n阶m序列则取其前n个作为m序列的初值。（详见“程序清单”主程序中无限循环部分）b）定时中断输出m序列我们选择m序列的波特率为 4

17、00Baud,而单片机晶振为12.0592MHz，所以定时器的计 数量为 pwm_time=11059200/12/400=2304。我们采用定时器 T0作为m序列输出定时器，工作于模式1 （即16位定时器）。中断服务程序如下：timer0（） interrupt 1 using 1/T0中断，发送 m序列EA = 0;TH0 = vth0;TL0 = vtl0;TRO = 1;pwm_out=pwm_status;syn c=s yn c_sta;v=O; EA = 1;2. 实验结果及m序列性质验证a）实验结果n = 3级m序列n = 4级m序列n = 5级m序列n = 6级m序列n =

18、7级m序列n = 8级m序列b）序列性质1.M序列的周期实验观察到的m序列为:1010011101011001000111101011001000111101011001000111101010110011011101101001001110 101010100110011101110100101100 101010100110010111010001110110 由表一可得生成的 M序列的周期p = 2n - 1，如n=3时，p=7; n=4时，p=152. M序列的游程4 1由上表所得的M序列所得，当n=4时，共有2 上8个游程，分别为1111, 0, 1, 0, 11,00, 1 , 0

19、00;其中长度为4的游程有1个（1111）,长度为3的游程有1个（000）,长度为2 的游程有2个（11, 00）,长度为1的游程有4个（0, 1 , 0 , 1）综上长度为1 （ k=1）的游程占了总游程数的一半（1/2 ）,长度为2 （ k=2）的游程占了1 1长度为3 （ n-1 ）的游程是连续0的游程（000）,长度为4 （ n）的游程为连续12* k 4的游程（1111 ）。同理分析n=3 , 5 , 6 , 7 , 8时得到结论：m序列的一个周期（p=2n-1）中，游程总数为2n-1。其中长度为1的游程个数占游程总数的1/2 ;长度为2的游程个数占游程总数的1/22=1/4 ;长度

20、为3的游程个数占游程总数的1/23=1/8 ;一般地，长度为 k的游程个数占游程总数的 1/2k=2-k ,其中1 k （n-2）。而且，在长度为k游程中，连1游程与连0游程各占一半，长为（n-1） 的游程是连0游程，长为n的游程是连1游程。七、程序清单#include #include /*/*/*/常量定义 /*/*/*/*定时器初值*/#define V_TH1 0xee /* 5ms , ! notice: crystal 11.0592MHz */#define V_TL1 0x00/*数码管段驱动寄存器地址*/#define ADDR_8SEG XBYTE0x2000/*数码管位驱

21、动和指示灯驱动寄存器地址*/ #define ADDR_SEL XBYTE0x4000/*按键*/sbit KEY1 =P1A0;sbit KEY2=P1A1;sbit KEY3=P1A2;sbit KEY4=P1A3;/*/*/*/变量定义/*/*/*/*数码管位驱动和指示灯驱动信号输岀缓存，定义了一个可位寻址的变量*/unsigned char bdata output_sel;sbit led_1 = output_selA5;sbit led_2 = output_seL6;sbit led_3 = output_seL7;sbit led_4 = output_seL4;/*数码管扫

22、描驱动指针，为测试外部存储器（ U3 6264）,特使用xdata类型*/ unsigned char xdata digi_scaner;/*测试用计数值十进制表示，为测试外部存储器（ U3 6264）,特使用xdata类型*/unsigned char xdata digi4=0,0,0,3;unsigned char temp4;/ 按键消抖/*pwm控制的变量*/unsigned int pwm_time;unsigned char vth0,vtl0;bit pwm_status;/ 信号bit sync_sta;/ 同步int power;/2的乘方的结果，等于2的n次方unsig

23、ned char j;/ 用于 for 循环int prd;/period，由于序列周期的计数sbit pwm_out = P3A5;/m序列输出sbit sync = P3A4;/序列周期同步信号unsigned char v=0;发送完一位置零/*/*/*/函数定义 /*/*/*/* 7 段数码显示译码参数：DATA需要显示的数字或符号；返回值：7段译码结果(D70 = PGFEDCBA )*/unsigned char NUMT0SEG7(unsigned char DATA) unsigned char AA;switch (DATA) case 0: AA=0xc0;break; /

24、*0 */case 1: AA=0xf9;break; /*T */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 6: AA=0x82;break; /*6 */case 7: AA=0xf8;break; /*7 */case 8: AA=0x80;break; /*8 */case 9: AA=0x90;break; /*9 */case 10: AA=0x88;break; /* A */case 1

25、1: AA=0x83;break; /* B */case 12: AA=0xc6;break; /*C */case 13: AA=0xa1;break; /* D */case 14: AA=0x86;break; /* E */case 15: AA=0x8e;break; /*F*/case -:AA=0xbf;break; /*破折号，此处原误为 0xdf，系04级王资凯同学指正*/case _:AA=0xf7;break; /*下划线*/case :AA=0xff;break; /*消隐*/default: AA=Oxff;return(AA); timer1() interrup

26、t 3 using 0 /* T1时钟中断服务程序在本例中，每 5ms被执行一次*/EA=O; /*关中断*/*重新对计数器赋初值，并启动定时计数*/TH1=V_TH1;TL仁 V_TL1;TR 1=1;EA=1;output_sel = OxfO;/* 初值，令数码管驱动位无效，指示灯全灭*/*按键4消抖和按键4确认的处理*/if (KEY4=0) led_4=0; temp4+; if(temp420 & KEY4=1) temp4=0;if(+digi38) digi3=3;power=1;prd=1;for(j=0;j=5) digi_scaner = 1;switch (digi_s

27、caner)case 1: /*驱动第一个数码管 */output_sel |= 0x01;ADDR_8SEG = NUMTOSEG7(pwm_status);/*输出到锁存器 U5 (在电路图中找)*/ break;case 2: /*驱动第二个数码管 */output_sel |= 0x02;ADDR_8SEG =0xff; /* 输出到锁存器 U5*/break;case 3: /*驱动第三个数码管 */output_sel |= 0x04;ADDR_8SEG = 0xff; /* 输出到锁存器 U5*/break;case 4: /*驱动第四个数码管 */output_sel |= 0

28、x08;ADDR_8SEG = NUMTOSEG7(digi3); /* 输出到锁存器 U5*/break;ADDR_SEL = output_sel; /*输岀到锁存器U6（在电路图中找）*/timerO() interrupt 1 using 1/T0 EA = 0;TH0 = vth0;TL0 = vtl0;TR0 = 1;pwm_out=pwm_status; sync=sync_sta;v=0;EA = 1;中断，发送m序列/*主程序*/main()int n;int i;int count=0;int temp;int yd8=1,0,1,0,1,0,1,0;/*单片机系统初始化*

29、/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;/数码管计数 pwm_time=2304;机器周期数，对应 400Hz vth0=(0xffff-pwm_time)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;/*主循环，本例中，在T1中断服务程序未被执行的空余时间里，处理机在以下程序中不断循环*/while (1)n