随机信号分析实验报告
更新时间:2023-10-27 10:14:01 阅读量: 综合文库 文档下载
Harbin Institute of Technology
实验报告
课程名称: 随机信号分析 院 系: 电子与信息工程学院 班 级: 姓 名: 学 号: 指导教师: 实验时间:
实验一、各种分布随机数的产生
(一)实验原理
1.均匀分布随机数的产生原理
产生伪随机数的一种实用方法是同余法,它利用同余运算递推产生伪随机数序列。最简单的方法是加同余法
yn?1?yn?c(modM)
xn?1?yn?1 M为了保证产生的伪随机数能在[0,1]内均匀分布,需要M为正整数,此外常数c和初值y0亦为正整数。加同余法虽然简单,但产生的伪随机数效果不好。另一种同余法为乘同余法,它需要两次乘法才能产生一个[0,1]上均匀分布的随机数
yn?1?ayn(modM)
xn?1?yn?1 Myn?1 M式中,a为正整数。用加法和乘法完成递推运算的称为混合同余法,即
yn?1?ayn?c(modM)
xn?1?用混合同余法产生的伪随机数具有较好的特性,一些程序库中都有成熟的程序供选择。
常用的计算语言如Basic、C和Matlab都有产生均匀分布随机数的函数可以调用,只是用各种编程语言对应的函数产生的均匀分布随机数的范围不同,有的函数可能还需要提供种子或初始化。
Matlab提供的函数rand()可以产生一个在[0,1]区间分布的随机数,
rand(2,4)则可以产生一个在[0,1]区间分布的随机数矩阵,矩阵为2行4列。Matlab提供的另一个产生随机数的函数是random('unif',a,b,N,M),unif 表示均匀分布,a和b是均匀分布区间的上下界,N和M分别是矩阵的行和列。 2.随机变量的仿真
根据随机变量函数变换的原理,如果能将两个分布之间的函数关系用显式表达,那么就可以利用一种分布的随机变量通过变换得到另一种分布的随机变量。
若X是分布函数为F(x)的随机变量,且分布函数F(x)为严格单调升函数,令Y=F(X),则Y必为在[0,1]上均匀分布的随机变量。反之,若Y是在[0,1]上
?1均匀分布的随机变量,那么X?FX(Y)
?1即是分布函数为FX(x)的随机变量。式中FX(?)为FX(?)的反函数。这样,欲
求某个分布的随机变量,先产生在[0,1]区间上的均匀分布随机数,再经上式变换,便可求得所需分布的随机数。 3.高斯分布随机数的仿真
广泛应用的有两种产生高斯随机数的方法,一种是变换法,一种是近似法。 如果X1,X2是两个互相独立的均匀分布随机数,那么下式给出的Y1,Y2
??Y1???2lnX1cos(2πX2)?m ?Y???2lnXsin(2πX)?m?12?2便是数学期望为m,方差为?2的高斯分布随机数,且互相独立,这就是变换法。 另外一种产生高斯随机数的方法是近似法。在学习中心极限定理时,曾提到n个在[0,1]区间上均匀分布的互相独立随机变量Xi (i=1,2…,n),当n足够大时,其和的分布接近高斯分布。当然,只要n不是无穷大,这个高斯分布是近似的。由于近似法避免了开方和三角函数运算,计算量大大降低。当精度要求不太高时,近似法还是具有很大应用价值的。 4.各种分布随机数的仿真
有了高斯随机变量的仿真方法,就可以构成与高斯变量有关的其他分布随机变量,如瑞利分布、指数分布和?2分布随机变量。
(二)实验目的
在很多系统仿真的过程中,需要产生不同分布的随机变量。利用计算机可以很方便地产生不同分布的随机变量,各种分布的随机变量的基础是均匀分布的随机变量。有了均匀分布的随机变量,就可以用函数变换等方法得到其他分布的随机变量。
(三)实验结果
附:源程序 subplot(2,2,1);
x=random('unif',2,5,1,1024);
plot(x);
title('均匀分布随机数')
subplot(2,2,2);
G1=random('Normal',0,1,1,20000); plot(G1);
title('高斯分布随机数')
subplot(2,2,3);
G2=random('Normal',0,1,1,20000); R=sqrt(G1.*G1+G2.*G2); plot(R);
title('瑞利分布随机数')
subplot(2,2,4);
G3=random('Normal',0,1,1,20000); G4=random('Normal',0,1,1,20000); X=G1.*G1+G2.*G2+G3.*G3+G4.*G4; plot(X);
title('x^2分布随机数')
实验二、随机变量检验(一)实验原理
1、均值的计算 在实际计算时,如果平稳随机序列满足各态历经性,则统计均值可用时间均值代替。这样,在计算统计均值时,并不需要大量样本函数的集合,只需对一个样本函数求时间平均即可。甚至有时也不需要计算N??时的极限,况且也不可能。通常的做法是取一个有限的、计算系统能够承受的N求时间均值和时间方差。根据强调计算速度或精度的不同,可选择不同的算法。
设随机数序列{x1,x2,?,xN},一种计算均值的方法是直接计算下m?1N?xnNn?1式中,xn为随机数序列中的第n个随机数。
另一种方法是利用递推算法,第n次迭代的均值也亦即前n个随机数的均值为mn?n?111mn?1?xn?mn?1?(xn?mn?1)迭代结束后,便得到随机数序列的均值 nnnm?mN
递推算法的优点是可以实时计算均值,这种方法常用在实时获取数据的场
合。
当数据量较大时,为防止计算误差的积累,也可采用m?m1?中,m1是取一小部分随机数计算的均值。
2、方差的计算
计算方差也分为直接法和递推法。仿照均值的做法
?2?1N1N222 (x?m)???n?xn?m2
Nn?1Nn?11n?1212(xn?mn?1) ?n?[?n?1?(xn?mn?1)2] nnn1N?(xn?m1)式Nn?1方差的递推算法需要同时递推均值和方差
mn?mn?1?2迭代结束后,得到随机数序列的方差为?2??N
其它矩函数也可用类似的方法得到。 3、统计随机数的概率密度直方图
假定被统计的序列x(n)的最大值和最小值分别为a和b。将(a,b)区间等分M(M应与被统计的序列x(n)的个数N相适应,否则统计效果不好。)份后的区间
(b?a)(b?a)2*(b?a)),(a?,a?),… , MMM(b?a)(i?1)2*(b?a)*i(b?a)(M?1),b)。用f(i),表示(a?,a?),… , (a?MMM(b?a)(i?1)2*(b?a)*i,a?)区间里的个数,统计序列序列x(n)的值落在(a?MM为(a,a?x(n)的值在各个区间的个数f(i),i?0,2,?,M?1,则f(i)就粗略地反映了随机
序列的概率密度的情况。用图形方式显示出来就是随机数的概率密度直方图。
(二)实验目的
正在阅读:
随机信号分析实验报告10-27
【优选精文】酒店销售部管理制度05-12
2018化学中考题分类汇编:微粒结构示意图与元素周期表综合03-27
禁塑从我做起作文500字06-30
基于VHDL的智能抢答器05-24
九年级语文单元检测题参考答案04-06
健全人格标准05-28
端面式高密度聚乙烯(HDPE)柔性抗震超级静音排水管材08-19
请别乱说话作文600字06-30
- Win7 安装MySql图示
- 计算器课程设计报告
- 部编版八年下语文第三单元第六单元古诗文理解默写练习及答案
- 13质量通病防治方案和施工措施
- 土力学试题~~~~
- 公务员打印资料
- 传热膜系数测定实验报告 - 图文
- 新时期煤矿协管安全工作的创新与实践
- 第五章 习题及参考答案
- 220kV架空线路强条执行记录表
- 音乐欣赏读后感
- 高炉
- 劳动教育需要新的时代内涵
- 10建筑地面工程施工质量验收规范GB50209-20021
- 银行会计练习题2答案
- 2013年七年级地理上册知识点复习提纲湘教版
- 人教版三年级语文上册第四单元测试题(A卷)(有答案)
- 营养师第九章练习题
- 湖北省武汉市2018届高三毕业生二月调研 理综化学
- 行业分析2018-2023年中国男性护肤品行业市场发展分析及投资前景
- 信号
- 随机
- 实验
- 报告
- 分析
- 青岛大学生职业生涯规划大赛书面作品评分标准 - 图文
- 小学五年级语文自我提高III(74)
- 2009-2010学年第二学期无机化学试卷
- 医院人力资源管理办法
- 如何更好地发挥教材中各栏目设置的功能
- 2016年兵团干部在线学习“互联网+”与公务员培训答案
- 起重机司机题库及答案
- 2015年电大公共关系学期末考试试题题库及答案 - 图文
- 银行排队系统实验
- 国家林业局关于加强林木采伐许可证核发管理工作的通知
- 2011西安世界园艺博览会考察报告
- 会计制度设计案例分析及答案
- 4.太阳能电池测试分选工序介绍
- ANSYS非稳态热分析及实例详解
- 财务管理复习题(附答案)
- 2017版步步高初高中化学(通用)衔接教材:第12讲 摩尔质量和气体摩尔体积
- 软件工程课程设计报告
- 2000湖南师范大学年真题
- 昆虫记测试题以及答案(经典阅读理解)
- 2019届高考语文作文经典写作素材名人名言分类及相应事例