ImageVerifierCode 换一换
格式:DOCX , 页数:33 ,大小:557.73KB ,
资源ID:18299815      下载积分:3 金币
快捷下载
登录下载
邮箱/手机:
温馨提示:
快捷下载时,用户名和密码都是您填写的邮箱或者手机号,方便查询和重复下载(系统自动生成)。 如填写123,账号就是123,密码也是123。
特别说明:
请自助下载,系统不会自动发送文件的哦; 如果您已付费,想二次下载,请登录后访问:我的下载记录
支付方式: 支付宝    微信支付   
验证码:   换一换

加入VIP,免费下载
 

温馨提示:由于个人手机设置不同,如果发现不能下载,请复制以下地址【https://www.bdocx.com/down/18299815.html】到电脑端继续下载(重复下载不扣费)。

已注册用户请登录:
账号:
密码:
验证码:   换一换
  忘记密码?
三方登录: 微信登录   QQ登录  

下载须知

1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。
2: 试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓。
3: 文件的所有权益归上传用户所有。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 本站仅提供交流平台,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

版权提示 | 免责声明

本文(实验一 随机信号的产生相关分析及其应用实验Word文档下载推荐.docx)为本站会员(b****3)主动上传,冰豆网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知冰豆网(发送邮件至service@bdocx.com或直接QQ联系客服),我们立即给予删除!

实验一 随机信号的产生相关分析及其应用实验Word文档下载推荐.docx

1、 231-1。通常公式中参数常用取值为s(0) =12357,r = 2045,b = 1,M =1048576。较简单的高斯白噪声产生方法是基于概率论中的中心极限定理。即无穷多个同分布随机变量之和构成随机变量服从高斯分布。方便起见,可用 N 个(通常 N=12)均匀分 . 布随机变量之和 Xi近似高斯分布随机变量。若 Xi,i = 0, 1, , 11 在0, 1)上服从均匀分布,则 Y= 错误!未找到引用源。近似服从均值为 0,方差为 1 的高斯分布。随机信号相关函数计算、相关分析及实现方法原理:离散随机序列自相关函数定义为 Rx(m)=Ex(n)x(n+m)。对于各态历经随机过程,统计平

2、均可用时间平均代替,即Rx(m)=错误!工程实践中,无法获得无限长数据,只能用有限平均来近似,即 Rx(m)=错误!,m=0, 1, , K-1;KN。为保证估计质量,通常要求 K此时 Rx(m)也可以简化为错误!。同理,也类似地计算互相关函数。 2实验内容 (1) 编程实现产生10000个在(0, 1)区间均匀分布随机数。计算生成随机数的14阶矩,最大值,最小值,频度直方图。(2) 编程实现产生10000个N(3, 4) 高斯随机数。(3) 编程实现产生10000个N(1, 2) 高斯随机数和10000个N(3, 4) 高斯随机数。计算其自相关函数, 计算两个高斯随机信号的互相关函数。(4)

3、 生成一个10000点高斯随机信号x(n):N(1, 4)。将该信号延迟N点,N 200,得到随机信号y(n)。计算两个信号的互相关函数,搜索峰值位置估算延迟。3实验仪器设备 微型计算机、Matlab开发环境(本报告采用MATLAB2016a) 三、 实验步骤 (1) 实现产生 10000 个在(0, 1)区间均匀分布随机数。计算生成随机数的(2) 14阶矩,最大值,最小值,频度直方图。编程实现产生10000个N(3, 4)高斯随机数。 (3) 计算生成随机数的14阶矩,最大值,最小值,频度直方 (4) 编程实现产生10000个(5) N(1, 2)高斯随机数和10000个N(3, 4)高斯随

4、机数。计算其自相关函数,计算两个高斯随机信号的互相关函数。 (6) 生成一个10000点高斯随机信号x(n):N(1,4)。(7) 将该信号延迟N点,N200,得到随机信号y(n) (8)。、 (9) 四、 实验数据记录表格 )数据处理记录表实验内容(1表1-1 . 三阶原点矩均值 均方值 四阶原点矩 最大值最小值 )数据处理记录表2表1-2 实验内容( 三阶原点矩 均方值 均值 最大值最小值 四阶原点矩 )数据处理记录表实验内容(4表1-3 延迟估算峰值位置 大连理工大学实验报告 电 子 信 息 工 程班级: 1401 姓 名:* 学号*组: 2016/10/31 C227 1 实验原理 较

5、简单的伪随机序列产生方法是采用数论中基于数环理论的线性同余法(乘同余法、混合同余法),其迭代公式的一般形式为f(x) = (r*x + b) Mod M,其离散形式为s(n + 1) = r*s(n)+ b Mod M。 (1) 编程实现产生10000个在(0, 1)区间均匀分布随机数。计算生成随机数的14阶矩,最 . 大值,最小值,频度直方图。三、 主要仪器设备 (注:以下为代码内容已经有详细的说明已包括步骤和操作办法) 四、实验步骤与操作方法 实验内容(1):编程实现产生10000个在(0, 1)区间均匀分布随机数。计算生成随机数的14阶矩,最大值,最小值,频度直方图 代码如下:%清除以前

6、数据缓存以及命令行窗口的显示 clear clc % 均匀分布随机数的产生 % 生成 10000 数据 s = randuniform(1048576,1,2045,12357,10000);% 画线 figure,plot(s) %全部 title(全部数据连线) figure,plot(s(1:100) %前100个数据 前100个数据连线% 画点 figure,plot(s,.) %全部 全部数据画点100),*) %前100个数据 前100个数据画点% 画直方图 hist(s) 10个区间hist(s,100) 100个区间% 估计随机数的分布 f,xi=ksdensity(s);fi

7、gure,plot(xi,f);利用MATLAB函数ksdensity估计得概率密度% 随机数的均值、均方值、一维三阶原点矩、一维四阶原点矩 m = zeros(1,4); %生成0 0 0 0 for i = 1 : 10000 m(1) = m(1) + s(i); % 均值 m(2) = m(2) + s(i)2; % 二阶矩 m(3) = m(3) + s(i)3; %三阶 m(4) = m(4) + s(i)4; %四阶 end m5=max(s);%最大值 m6=min(s);%最小值 m=m/10000;disp(均值 = ,num2str(m(1) );均方值 = ,num2s

8、tr(m(2) );三阶原点矩 = ,num2str(m(3) );四阶原点矩 = ,num2str(m(4) );最大值 = ,num2str(m5) );最小值 = ,num2str(m6) );% 对比:Matlab函数 disp(=以下对比Matlab库函数生成的数据=按任意键继续=pause;sm = rand(1,10000);% 画线 figure,plot(sm) %全部 figure,plot(sm(1:% 画点 figure,plot(sm,% 画直方图 hist(sm) hist(sm,100) % 估计随机数的分布 f,xi=ksdensity(sm);% 求随机数的均

9、值、均方值、一维三阶原点矩、一维四阶原点矩 m(1) = m(1) + sm(i); m(2) = m(2) + sm(i)2; m(3) = m(3) + sm(i)3; m(4) = m(4) + sm(i)4;m5=max(sm); %最大值 m6=min(sm); %最小值 实验内容(2): 编程实现产生10000个N(3, 4) 高斯随机数。%高斯分布随机数的产生 m=0;a=1;num=10000;s = rnd(m,a,num); %rnd()自己编写的产生高斯随机数的函数 % 求随机数的均值、均方值、一维三阶原点矩、一维四阶原点矩、最大值和最小值 % 对比:=按任意键继续=%

10、pause;% 生成 10000 数据 sm = random(normal,0,1,1,num);实验内容(3):(3) 编程实现产生10000个N(1, 2) 高斯随机数和10000个N(3, 4) 高斯随机数。%定义高斯随机的方差和均值及产生的数量 m1=1;a1=2;m2=3;a2=4;%利用自带函数产生高斯随机数 d1=random(norm,1,2,1,num);d2=random(,3,4,1,num);%计算自相关和互相关函数 Rx1=xcorr(d1);Rx2=xcorr(d1);Rx1x2=xcorr(d1,d2);m=-num+1:num-1;figure,plot(m,

11、Rx1);figure,plot(m,Rx2);figure,plot(m,Rx1x2);实验内容(4):a1=4;fs=500;N=100;long=1000;%由于数据过多,延时图比较不明显,显示1000个数据 x1=random(,1,4,1,num);x2=random(x=random(,1,4,1,N); %产生一个1维N个随机高斯随机数,用来替换延时后的空位 %延时N for i=1:10000 if(imax11) max11=cor11(1,i); delay11=lag11(1,i);max11 delay11 if (iN) x2(i)=x1(i-N); x2(i)=x(

12、i);%将延时后的前100用产生的高斯随机数代替 用高斯随机数替换后的延时信号cor12, lag12=xcorr(x1,x2);figure,plot(lag12/fs,cor12) max12=0; if(cor12(1,i)max12) max12=cor12(1,i); delay12=lag12(1,i);max12 delay12 附:本次作业中用到的两个随机信号源代码 1. 产生(0, 1)区间均匀分布随机数的函数:function s = randuniform(M,b,r,first,num) %调用示例: s = randuniform(1048576,1,2045,123

13、57,10000) % 其中,s(n)为n时刻的随机数种子,r为扩展因子,b为固定扰动项,M为循环模,Mod M表示对M取模。% 通常公式中参数常用取值为s(0) =12357,r = 2045,b = 1,M =1048576。s = zeros(1,num);s(1) = first;for i = 2 : num s(i) = mod(s(i-1)*r+b,M);s = s/M;2. 高斯随机序列产生函数:function s1,s2 =rnd( m,a,num ) %其中a是均值,m是方差 ,num变量则是数据量 s = rnd(0,1,1000), x1=rand(1,num);x2=rand(1,num);y1=sqrt(-2*log(x1).*cos(2*pi*x2);y2=sqrt(-2*log(x1).*sin(2*pi*x2);s1=a*y1+m;s2=a*y2+m; 五、实验数据记录和处理 2)结果:实验内容(图(自编(自编函数)100001-1.1:个均匀随机数的连线图 :图1-1.410000个均匀随机变量频度直方图 函数)图 图1-1.5:估得的概率密度图(自编函数) 1001-1.2:前个均匀随机数连线图(自编函数)图(库函数)10000个均匀随机数数据点图10000(自编函数)1-1.6 图:1-1.3 .

copyright@ 2008-2022 冰豆网网站版权所有

经营许可证编号:鄂ICP备2022015515号-1