- 1、本文档共21页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
实用标准文案
精彩文档
Harbin Institute of Technology
实验报告
课程名称: 随机信号分析
院 系: 电子与信息工程学院
班 级:
姓 名:
学 号:
指导教师:
实验时间:
实验一、各种分布随机数的产生
(一)实验原理
1.均匀分布随机数的产生原理
产生伪随机数的一种实用方法是同余法,它利用同余运算递推产生伪随机数序列。最简单的方法是加同余法
为了保证产生的伪随机数能在[0,1]内均匀分布,需要M为正整数,此外常数c和初值y0亦为正整数。加同余法虽然简单,但产生的伪随机数效果不好。另一种同余法为乘同余法,它需要两次乘法才能产生一个[0,1]上均匀分布的随机数
式中,a为正整数。用加法和乘法完成递推运算的称为混合同余法,即
用混合同余法产生的伪随机数具有较好的特性,一些程序库中都有成熟的程序供选择。
常用的计算语言如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]上均匀分布的随机变量,那么
即是分布函数为FX(x)的随机变量。式中为的反函数。这样,欲求某个分布的随机变量,先产生在[0,1]区间上的均匀分布随机数,再经上式变换,便可求得所需分布的随机数。
3.高斯分布随机数的仿真
广泛应用的有两种产生高斯随机数的方法,一种是变换法,一种是近似法。
如果X1,X2是两个互相独立的均匀分布随机数,那么下式给出的Y1,Y2
便是数学期望为m,方差为的高斯分布随机数,且互相独立,这就是变换法。
另外一种产生高斯随机数的方法是近似法。在学习中心极限定理时,曾提到n个在[0,1]区间上均匀分布的互相独立随机变量Xi (i=1,2…,n),当n足够大时,其和的分布接近高斯分布。当然,只要n不是无穷大,这个高斯分布是近似的。由于近似法避免了开方和三角函数运算,计算量大大降低。当精度要求不太高时,近似法还是具有很大应用价值的。
4.各种分布随机数的仿真
有了高斯随机变量的仿真方法,就可以构成与高斯变量有关的其他分布随机变量,如瑞利分布、指数分布和分布随机变量。
(二)实验目的
在很多系统仿真的过程中,需要产生不同分布的随机变量。利用计算机可以很方便地产生不同分布的随机变量,各种分布的随机变量的基础是均匀分布的随机变量。有了均匀分布的随机变量,就可以用函数变换等方法得到其他分布的随机变量。
(三)实验结果
附:源程序
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、均值的计算
在实际计算时,如果平稳随机序列满足各态历经性,则统计均值可用时间均值代替。这样,在计算统计均值时,并不需要大量样本函数的集合
您可能关注的文档
- 数学期望在经济决策中地运用.doc
- 数学实验5:线性代数方程组地数值解法.doc
- 数学中考模拟卷及问题详解.doc
- 数值代数实验报告材料(1)绝对经典.doc
- 数值分析报告(计算方法)期末试卷及参考问题详解.doc
- 数值分析报告课程设计.doc
- 数值分析报告模拟试卷1,2,3.doc
- 数值分析报告实验报告材料.doc
- 数值分析报告试卷及问题详解.doc
- 数值计算(二分法、简单迭代法、Newton迭代法、弦截法(割线法、双点弦法)).doc
- 2024年江西省寻乌县九上数学开学复习检测模拟试题【含答案】.doc
- 2024年江西省省宜春市袁州区数学九上开学学业水平测试模拟试题【含答案】.doc
- 《GB/T 44275.2-2024工业自动化系统与集成 开放技术字典及其在主数据中的应用 第2部分:术语》.pdf
- 中国国家标准 GB/T 44275.2-2024工业自动化系统与集成 开放技术字典及其在主数据中的应用 第2部分:术语.pdf
- GB/T 44285.1-2024卡及身份识别安全设备 通过移动设备进行身份管理的构件 第1部分:移动电子身份系统的通用系统架构.pdf
- 《GB/T 44285.1-2024卡及身份识别安全设备 通过移动设备进行身份管理的构件 第1部分:移动电子身份系统的通用系统架构》.pdf
- 中国国家标准 GB/T 44285.1-2024卡及身份识别安全设备 通过移动设备进行身份管理的构件 第1部分:移动电子身份系统的通用系统架构.pdf
- GB/T 44275.11-2024工业自动化系统与集成 开放技术字典及其在主数据中的应用 第11部分:术语制定指南.pdf
- 中国国家标准 GB/T 44275.11-2024工业自动化系统与集成 开放技术字典及其在主数据中的应用 第11部分:术语制定指南.pdf
- 《GB/T 44275.11-2024工业自动化系统与集成 开放技术字典及其在主数据中的应用 第11部分:术语制定指南》.pdf
文档评论(0)