实验1 白噪声和M序列的产生.doc

实验1 白噪声和M序列的产生实验报告哈尔滨工业大学航天学院控制科学与工程系专业:           班级:            姓名:           学号:            : 白噪声和M序列的产生     、熟悉并掌握产生均匀分布随机序列方法以及进而产生高斯白噪声方法2、熟悉并掌握M序列生成原理及仿真生成方法 、混合同余法混合同余法是加同余法和乘同余法的混合形式,其迭代式如下:式中a为乘子,为种子,b为常数,M为模。混合同余法是一种递归算法,即先提供一个种子,逐次递归即得到一个不超过模M的整数数列。2、正态分布随机数产生方法由独立同分布中心极限定理有:设随机变量相互独立,服从同一分布,且具有数学期望和方差:则随机变量之和的标准化变量:近似服从分布。如果服从[0,1]均匀分布,则上式中,。即近似服从分布。3、M序列生成原理用移位寄存器产生M序列的简化框图如下图所示。该图表示一个由4个双稳态触发器顺序连接而成的4级移位寄存器,它带有一个反馈通道。当移位脉冲来到时,每级触发器的状态移到下一级触发器中,而反馈通道按模2加法规则反馈到第一级的输入端。 、生成均匀分布随机序列(1)利用混合同余法生成[0,1]区间上符合均匀分布的随机序列,并计算该序列的均值和方差,与理论值进行对比分析。要求序列长度为1200,推荐参数为a=65539,M=2147483647,0<x0<M。(2)将[0,1]区间分为不重叠的等长的10个子区间,绘制该随机序列落在每个子区间的频率曲线图,辅助验证该序列的均匀性。(3)对上述随机序列进行独立性检验。(该部分为选作内容)2、生成高斯白噪声利用上一步产生的均匀分布随机序列,令n=12,生成服从N(0,1)的白噪声,序列长度为100,并绘制曲线。3、生成M序列M序列的循环周期取为,时钟节拍,幅度,逻辑“0”为a,逻辑“1”为-a,特征多项式。生成M序列的结构图如下所示。要求编写Matlab程序生成该M序列,绘制该信号曲线,并分析验证M序列的性质。 -1实验1-2实验1-3 %实验1-1Xulie_Length=1200;A=65539;M=2147483647;b=1;R(Xulie_Length)=0;X(1)=199119;R(1)=X(1)/M;forn=1:(Xulie_Length-1)X(n+1)=mod((A*X(n)+b),M);R(n+1)=X(n+1)/M;endfigure(1);plot(R);figure(2);hist(R);mean_R=mean(R)var_R=var(R)---------------------------------------------------------------------------------------------------------------------%实验1-2clearall;N=12;SEQ_LENGTH=100;G(SEQ_LENGTH)=0;forn=1:SEQ_LENGTHfort=1:NG(n)=G(n)+R(N*(n-1)+t);endendG=G-6;figure(3);plot(G);figure(4);hist(G);---------------------------------------------------------------------------------------------------------------------%实验1-3clearall;number=2^6-1;a=1;%手动初始化MM(1)=1;M(2)=0;M(3)=0;M(4)=1;M(5)=1;M(6)=0;M_XuLie(number)=0;forn=1:numbertemp=xor(M(6),M(5));if(temp==0)M_XuLie(n)=a;elseM_XuLie(n)=-a;endM(6)=M(5);M(5)=M(4);M(4)=M(3);M(3)=M(2);M(2)=M(1);M(1)=temp;endstairs(M_XuLie);ylim([-]); 、生成均匀分布随机序列(1)生成的0-1均布随机序列如下所示:计算序列的均值和方差程序代码:mean_R=mean(R)var_R=var(R)均值和方差实际值:mean_R==