实验一-随机信号产生、相关分析与其应用实验.docVIP

电子1401 ***** ***** PAGE \* MERGEFORMAT1 2016-11-04 大连理工大学实验预习报告 学院（系）： 信息与通信工程学院 专业： 电 子 信 息 工 程 班级： 1401 姓 名： ***** 学号： ***** 组： ___ 实验时间： 2016/10/31 实验室： C227 实验台： 指导教师： 李 小 兵 实验I：随机信号的产生、相关分析及其应用实验 实验目的和要求 通过实验理解掌握随机信号样本生成的原理和方法、掌握随机过程相关函数的计算原理和方法。训练MATLAB程序代码编写能力，要求完成以下工作，并将实验结果与理论分析对照。 1.基于均匀分布伪随机数，掌握均匀分布白噪声典型生成方法。 2.基于均匀分布伪随机数，掌握高斯分布白噪声典型生成方法。 3.掌握随机信号相关函数计算、相关分析及实现方法。 实验原理和内容 1 实验原理 较简单的伪随机序列产生方法是采用数论中基于数环理论的线性同余法（乘同余法、混合同余法），其迭代公式的一般形式为f(x) = (r*x + b) Mod M，其离散形式为s(n + 1) = [r*s(n)+ b] Mod M。其中，s(n)为n 时刻的随机数种子，r 为扩展因子，b 为固定扰动项，M 为循环模，Mod M 表示对M 取模。为保证s(n)的周期为M，r 的取值应满足r = 4k + 1， p M 2 ? ，k与p 的选取应满足：r < M，r(M-1) + 1< 231-1。通常公式中参数常用取值为s(0) =12357，r = 2045，b = 1，M =1048576。 较简单的高斯白噪声产生方法是基于概率论中的中心极限定理。即无穷多个同分布随机变量之和构成随机变量服从高斯分布。方便起见，可用 N 个（通常 N=12）均匀分布随机变量之和 Xi近似高斯分布随机变量。若 Xi，i = 0, 1, …, 11 在[0, 1)上服从均匀分布，则 Y= i=011 随机信号相关函数计算、相关分析及实现方法原理： 离散随机序列自相关函数定义为 Rx(m)=E[x(n)x(n+m)]。对于各态历经随机过程，统计平均可用时间平均代替，即Rx(m)=limn→∞12N+1n=-NNxnx(n+m)工程实践中，无法获得无限长数据，只能用有限平均来近似，即 Rx(m)=1N-m 2实验内容 (1) 编程实现产生10000个在(0, 1)区间均匀分布随机数。计算生成随机数的1~4阶矩，最大值，最小值，频度直方图。 (2) 编程实现产生10000个N(3, 4) 高斯随机数。计算生成随机数的1~4阶矩，最大值，最小值，频度直方图。 (3) 编程实现产生10000个N(1, 2) 高斯随机数和10000个N(3, 4) 高斯随机数。计算其自相关函数， 计算两个高斯随机信号的互相关函数。 (4) 生成一个10000点高斯随机信号x(n)：N(1, 4)。将该信号延迟N点，N < 200，得到随机信号y(n)。计算两个信号的互相关函数，搜索峰值位置估算延迟。 3实验仪器设备 微型计算机、Matlab开发环境（本报告采用MATLAB2016a） 实验步骤 实现产生?10000?个在(0, 1)区间均匀分布随机数。? 计算生成随机数的1~4阶矩，最大值，最小值，频度直方图。 编程实现产生10000个N(3, 4)高斯随机数。 计算生成随机数的1~4阶矩，最大值，最小值，频度直方 编程实现产生10000个N(1, 2)高斯随机数和10000个N(3, 4)高斯随机数。 计算其自相关函数，计算两个高斯随机信号的互相关函数。 生成一个10000点高斯随机信号x(n)：N(1,4)。 将该信号延迟N点，N<200，得到随机信号y(n)。、 计算两个信号的互相关函数，搜索峰值位置估算延迟。 实验数据记录表格 表1-1 实验内容（1）数据处理记录表 均值 均方值 三阶原点矩 四阶原点矩 最小值 最大值 表1-2 实验内容（2）数据处理记录表 均值 均方值 三阶原点矩 四阶原点矩 最小值 最大值 表1-3 实验内容（4）数据处理记录表 峰值位置 延迟估算 大连理工大学实验报告 学院（系）： 信息与通信工程学院 专业： 电 子 信 息 工 程 班级： 1401 姓 名： ***** 学号： ***** 组： ___ 实验时间： 2016/10/31 实验室