欢迎来到天天文库
浏览记录
ID:44190044
大小:104.50 KB
页数:12页
时间:2019-10-19
《大连理工大学随机信号实验一实验报告》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、大连理工大学实验预习报告学院(系):信息与通信工程学院专业:电子信息工程班级:电子1202班姓名:孙相国学号:201281617组:实验时间:2014/11/5实验室:创新园C219实验台:指导教师签字:成绩:实验一随机信号的产生、相关分析及其应用实验实验1均匀分布随机数的产生,统计特性分析及计算机仿真一、实验目的和要求掌握均匀分布随机信号的基本产生方法二、实验原理和内容较简单的伪随机序列产生方法是采用数论中基于数环理论的线性同余法(乘同余法、混合同余法),其迭代公式的一般形式为f(x)=(r·x+b)ModM,其离散形式为s(n+1)=[r·s(n)+b]ModM。
2、其中,s(n)为n时刻的随机数种子,r为扩展因子,b为固定扰动项,M为循环模,ModM表示对M取模。为保证s(n)的周期为M,r的取值应满足r=4k+1,pM2=,k与p的选取应满足:r3、,10000);s(1)=12357;s=zeros(1,10000);fori=2:10000s(i)=mod(s(i-1)*r+b,M);ends=s/M;%均匀分布随机生成10000数据%%画线figure,plot(s)%全部数据画线title('全部数据连线')%%画点figure,plot(s,'.')%全部数据画点title('全部数据画点')%%画直方图hist(s,100)title('100个区间')sum=0;fori=1:10000sum=sum+s(i);%求所有数的总和end;avr=sum/10000;%求所有数的平均数m=zeros(14、,4);fori=1:10000m(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;%四阶endn=zeros(1,4);fori=1:10000n(1)=n(1)+(s(i)-avr);%均值n(2)=n(2)+(s(i)-avr)^2;%二阶矩n(3)=n(3)+(s(i)-avr)^3;%三阶n(4)=n(4)+(s(i)-avr)^4;%四阶endm=m/10000;n=n/10000;disp(['均值=',num2str(m(1))]);disp(['二阶原5、点矩=',num2str(m(2))]);disp(['三阶原点矩=',num2str(m(3))]);disp(['四阶原点矩=',num2str(m(4))]);disp(['均值=',num2str(n(1))]);disp(['二阶中心矩=',num2str(n(2))]);disp(['三阶中心矩=',num2str(n(3))]);disp(['四阶中心矩=',num2str(n(4))]);c=0;d=1;fori=1:10000%求最大值ifcs(i)d=s(i);end6、;end;cd五、实验数据记录以及结果与分析实验结果如下:均值=0.50527二阶原点矩=0.33873三阶原点矩=0.25499四阶原点矩=0.20455均值=-3.749e-016二阶中心矩=0.083429三阶中心矩=-0.00046705四阶中心矩=0.012513最大最小值如下:c=0.9998d=0图表如下:六、讨论、建议、质疑通过本实验让我了解了生成随机数的程序,并了解了编写多阶原点矩和中心矩的程序,通过实际编程锻炼了自己的使用Matlab编写随机信号实验的能力。实验2高斯分布随机数的产生,统计特性分析及计算机仿真一、实验目的和要求掌握高斯白噪声的基本产7、生方法二、实验原理和内容(1)变换法(2)较简单的高斯白噪声产生方法是基于概率论中的中心极限定理。即无穷多个同分布随机变量之和构成随机变量服从高斯分布。方便起见,可以使用多个均匀分布随机变量之和近似高斯分布随机变量。三、实验步骤(1)编程实现产生10000个N(3,4)高斯随机数。(2)计算生成随机数的1~4阶矩,最大值,最小值,频度直方图。实验程序如下:m=3;n=4;pi=3.1415;s=zeros(1,10000);fori=1:10000a=sqrt(-2*log10(rand));b=2*pi*rand;s(1,i)=n*a*cos(b)
3、,10000);s(1)=12357;s=zeros(1,10000);fori=2:10000s(i)=mod(s(i-1)*r+b,M);ends=s/M;%均匀分布随机生成10000数据%%画线figure,plot(s)%全部数据画线title('全部数据连线')%%画点figure,plot(s,'.')%全部数据画点title('全部数据画点')%%画直方图hist(s,100)title('100个区间')sum=0;fori=1:10000sum=sum+s(i);%求所有数的总和end;avr=sum/10000;%求所有数的平均数m=zeros(1
4、,4);fori=1:10000m(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;%四阶endn=zeros(1,4);fori=1:10000n(1)=n(1)+(s(i)-avr);%均值n(2)=n(2)+(s(i)-avr)^2;%二阶矩n(3)=n(3)+(s(i)-avr)^3;%三阶n(4)=n(4)+(s(i)-avr)^4;%四阶endm=m/10000;n=n/10000;disp(['均值=',num2str(m(1))]);disp(['二阶原
5、点矩=',num2str(m(2))]);disp(['三阶原点矩=',num2str(m(3))]);disp(['四阶原点矩=',num2str(m(4))]);disp(['均值=',num2str(n(1))]);disp(['二阶中心矩=',num2str(n(2))]);disp(['三阶中心矩=',num2str(n(3))]);disp(['四阶中心矩=',num2str(n(4))]);c=0;d=1;fori=1:10000%求最大值ifcs(i)d=s(i);end
6、;end;cd五、实验数据记录以及结果与分析实验结果如下:均值=0.50527二阶原点矩=0.33873三阶原点矩=0.25499四阶原点矩=0.20455均值=-3.749e-016二阶中心矩=0.083429三阶中心矩=-0.00046705四阶中心矩=0.012513最大最小值如下:c=0.9998d=0图表如下:六、讨论、建议、质疑通过本实验让我了解了生成随机数的程序,并了解了编写多阶原点矩和中心矩的程序,通过实际编程锻炼了自己的使用Matlab编写随机信号实验的能力。实验2高斯分布随机数的产生,统计特性分析及计算机仿真一、实验目的和要求掌握高斯白噪声的基本产
7、生方法二、实验原理和内容(1)变换法(2)较简单的高斯白噪声产生方法是基于概率论中的中心极限定理。即无穷多个同分布随机变量之和构成随机变量服从高斯分布。方便起见,可以使用多个均匀分布随机变量之和近似高斯分布随机变量。三、实验步骤(1)编程实现产生10000个N(3,4)高斯随机数。(2)计算生成随机数的1~4阶矩,最大值,最小值,频度直方图。实验程序如下:m=3;n=4;pi=3.1415;s=zeros(1,10000);fori=1:10000a=sqrt(-2*log10(rand));b=2*pi*rand;s(1,i)=n*a*cos(b)
此文档下载收益归作者所有