资源描述：
《西电电子信息系统综合实验报告》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、西安电子科技大学电子信息系统综合实验课程实验报告实验名称电子信息系统综合实验成绩电子工程学院1402011班姓名梁思颖学号14020110055同作者王梦路李习习王保智郭鑫宇实验日期2017年11月25日指导教师评语：指导教师：年月日实验报告内容基本要求及参考格式一、实验目的二、实验所用仪器（或实验环境）三、实验基本原理及步骤（或方案设计及理论计算）四、实验数据记录（或仿真及软件设计）五、实验结果分析及回答问题（或测试环境及测试结果）MATLAB实验——复杂噪声产生1实验目的（1）掌握四种热噪声的基本分布并利用MATLA
2、B产生这几种噪声；（2）熟练使用MATLAB。2实验所用仪器（或实验环境）软件：MATLAB；硬件：计算机。3实验内容利用MATLAB分别产生高斯分布、均匀分布、指数分布、瑞利分布的热噪声。4实验步骤（1）服从高斯分布的热噪声Matlab7.0本身自带了标准高斯分布的内部函数randn，调用格式如下:Y=randn(n)Y=randn(m,n)Y=randn([mn])Y=randn(size(A))s=randn('state')randn函数产生的随机序列服从均值为m=0，方差σ2＝1的高斯分布。Y=randn(n)
3、产生的是一个n×n的随机序列矩阵，而Y=randn(m,n)和Y=randn([mn])产生的m×n的随机序列矩阵，Y=randn(size(A))产生的是大小与矩阵A同样大小的随机序列矩阵。s=randn('state')返回的是一个具有两个元素的向量，该向量显示的是当前正态随机数产生器的状态。randn('state',s)指令可以将产生器的状态设置到s，而randn('state',0)则可以将正态随机数产生器的状态恢复到初始状态。（2）服从均匀分布的热噪声同样Matlab本身也自带了（0-1）单位均匀分布的内部函
4、数rand，格式如下:Y=rand(n)Y=rand(m,n)Y=rand([mn])Y=rand(size(A))s=rand('state')rand函数产生的随机序列服从（0-1）单位均匀分布。Y=rand(n)产生的是一个n×n的随机序列矩阵，而Y=rand(m,n)和Y=rand([mn])产生的m×n的随机序列矩阵，Y=rand(size(A))产生的是大小与矩阵A同样大小的随机序列矩阵。s=rand('state')返回的是一个具有两个元素的向量，该向量显示的是当前（0-1）单位均匀随机数产生器的状态。ra
5、nd('state',s)指令可以将产生器的状态设置到s，而rand('state',0)则可以将（0-1）单位均匀分布随机数产生器的状态恢复到初始状态。（3）服从指数分布的热噪声先产生一个服从（0-1）单位分布的信号，然后再将其经过指数变换，就可以得到一个服从参数为λ的指数分布的信号了。（4）服从瑞利分布的热噪声先产生一个服从（0-1）分布的信号，然后再经过变换，可以得到一个服从瑞利（Rayleigh）分布的信号了。5程序设计（1）服从高斯分布的热噪声%服从高斯（Guass）分布的热噪声（随机序列）b=1;%均值fs=
6、1e7;%采样率t=1e-3;%随机序列长度n=t*fs;randn('state',0);%把高斯分布伪随机发生器置为0状态u=randn(1,n)+b;subplot(2,1,1),plot(u),title('高斯分布信号');subplot(2,1,2),hist(u,-4:0.1:4),title('高斯分布信号直方图');（2）服从均匀分布的热噪声%服从均匀分布的热噪声（随机序列）a=2;%（a-b）均匀分布下限b=3;%（a-b）均匀分布上限fs=1e7;%采样率，单位：Hzt=1e-3;%随机序列长度，单
7、位：sn=t*fs;rand('state',0);%把均匀分布伪随机发生器置为0状态u=rand(1,n);%产生（0-1）单位均匀信号x2=(b-a)*u+a;%广义均匀分布与单位均匀分布之间的关系figure(2)subplot(2,1,1),plot(x2),title('均匀分布信号');subplot(2,1,2),hist(x2,a:0.02:b),title('均匀分布信号直方图');（3）服从指数分布的热噪声%服从指数分布的热噪声（随机序列）lambda=2.5;%指数分布参数fs=1e7;%采样频率t
8、=1e-3;%时间长度n=t*fs;rand('state',0);%把均匀分布伪随机发生器置为0状态u=rand(1,n)%产生单位均匀信号x3=log2(1-u)/(-lambda);%指数分布与单位均匀分布之间的关系figure(3)subplot(2,1,1),plot(0:1/fs:t-1/fs,x3),x
