大连理工大学随机信号实验报告材料完整

1、大连理工大学实验预习报告学院（系）：信息与通信工程学院专业：电子信息工程班级：姓名：学号：组：_实验时间：2015.12.14实验室：C221实验台：指导教师：实验 I ：随机信号的产生、相关分析及其应用实验实验 1均匀分布随机数的产生，统计特性分析及计算机仿真一、实验目的和要求掌握均匀分布随机信号的基本产生方法二、实验原理和内容较简单的伪随机序列产生方法是采用数论中基于数环理论的线性同余法（乘同余法、混合同余法），其迭代公式的一般形式为f(x) = (r ·x+bM，其离散形式为sn) Mod(+1)=rs nbMsn为n时刻的随机数种子，r为扩展因子，·( ) + Mo

2、d。其中，( )b 为固定扰动项， M 为循环模， ModM 表示对 M 取模。为保证s( n) 的周期为 M，r的取值应满足 r= 4k+ 1 ，Mp， k 与 p 的选取应满足： r<M， r M31-1。s2， r， bM( -1)+1<2通常公式中参数常用取值为(0) =12357= 2045= 1，。=1048576三、实验步骤1. 编程实现产生 10000个在 (0, 1) 区间均匀分布随机数。2. 计算生成随机数的 14 阶矩，最大值，最小值，频度直方图。实验 2高斯分布随机数的产生，统计特性分析及计算机仿真一、实验目的和要求掌握高斯白噪声的基本产生方法二、实验原理和

3、内容1. 变换法2. 较简单的高斯白噪声产生方法是基于概率论中的中心极限定理。即无穷多个同分布随机变量之和构成随机变量服从高斯分布。方便起见，可以使用多个均匀分布随机变量之和近似高斯分布随机变量。三、实验步骤1. 编程实现产生10000个N(3,4)高斯随机数。2. 计算生成随机数的 14 阶矩，最大值，最小值，频度直方图。实验 3随机信号相关函数计算、相关分析及计算机仿真一、实验目的和要求掌握随机信号相关函数计算、相关分析及实现二、实验原理和内容根据自相关和互相关的定义， 自相关，互相关计算随机信号的自相关和互相关。三、实验步骤1. 产生高斯随机信号。2. 计算其自相关函数。3. 计算两个高

4、斯随机信号的互相关函数。大连理工大学实验报告学院（系）：信息与通信工程学院专业：电子信息工程班级：电子 1303姓名：实验时间：指导教师：李彤学号 ：201383081组：_2015.12.14实验室：C221实验台：实验 I ：随机信号的产生、相关分析及其应用实验实验 1均匀分布白噪声的生成一、实验目的和要求基于均匀分布伪随机数，掌握均匀分布白噪声典型生成方法。二、实验原理和内容较简单的伪随机序列产生方法是采用数论中基于数环理论的线性同余法（乘同余法、混合同余法），其迭代公式的一般形式为f(x) = (r ·x+bModM，其离散形式为sn)(+1)=rs nbMsn为n时刻的随机

5、数种子，r为扩展因子，·( ) + Mod。其中，( )b 为固定扰动项， M 为循环模， ModM 表示对 M 取模。为保证s( n) 的周期为 M，r的取值应满足r= 4k+ 1，M2p， k 与 p 的选取应满足：r<M， rM31-1。，r， b， M(-1)+1<2本实验中参数取值为s(1) =12357= 2025= 1。=1048576三、主要仪器设备微型计算机、 Matlab 开发环境四、实验步骤与操作方法1. 编程实现产生 10000 个在 (0, 1) 区间均匀分布随机数。2. 计算生成随机数的 14 阶矩，最大值，最小值，频度直方图五、实验数据记录和

6、处理程序如下：M=1048576;b=1;r=2025;n(2)=n(2)+(s(i)-avr)2;%求 二阶s=zeros(1,10000);矩s(1)=12357;n(3)=n(3)+(s(i)-avr)3;%求三阶s=zeros(1,10000);n(4)=n(4)+(s(i)-avr)4;%求四阶for i=2:10000ends(i)=mod(s(i-1)*r+b,M);m=m/10000;endn=n/10000;s=s/M;% 均匀分布随机生成10000 数据disp('均值 =',num2str(m(1);figure,plot(s)%全部数据画线disp(&#

7、39;二阶原点矩title('全部数据连线 ')=',num2str(m(2);figure,plot(s,'.')%全部数据画点disp('三阶原点矩title('全部数据画点 ')=',num2str(m(3);%画直方图disp('四阶原点矩hist(s,40)=',num2str(m(4);title('40个区间 ')disp('均值 =',num2str(n(1);sum=0;disp('二阶中心矩for i=1:10000=',num2str(n

8、(2);sum=sum+s(i);% 求所有数的总和disp('三阶中心矩end;=',num2str(n(3);avr=sum/10000;% 求所有数的平均数disp('四阶中心矩m=zeros(1,4);=',num2str(n(4);for i=1:10000disp('方差 =',num2str(var(s);m(1)=m(1)+s(i); %求均值c=0;m(2)=m(2)+s(i)2;%求二阶矩d=1;m(3)=m(3)+s(i)3;%求三阶for i=1:10000%求最大值m(4)=m(4)+s(i)4;%求四阶if c<

9、s(i)endc=s(i);n=zeros(1,4);end;for i=1:10000end;n(1)=n(1)+(s(i)-avr); %求均值for i=1:10000%求最小值if d>s(i)end;d=s(i);cend;d六、实验结果与分析运行程序， command窗口中显示的结果如下：均值 =0.49699二阶原点矩 =0.32916三阶原点矩 =0.24551四阶原点矩 =0.19546均值 =-8.3666e-017二阶中心矩 =0.08217三阶中心矩 =0.00024392四阶中心矩 =0.012193方差 =0.082179c =1.0000d =0实验得到图表

10、如下：图 1.1图 1.2七、讨论、建议、质疑本实验中编写了生成随机序列的程序，通过设定不同的参数值可以得到不同的随机序列，通过计算多阶原点矩、中心矩和绘制直方图可以更加清楚的看到产生的随机序列的特点，对随机序列的理解更加深刻。实验 2 高斯分布白噪声的生成一、实验目的和要求基于均匀分布伪随机数，掌握高斯分布白噪声典型生成方法。二、实验原理和内容1. 变换法2. 较简单的高斯白噪声产生方法是基于概率论中的中心极限定理。即无穷多个同分布随机变量之和构成随机变量服从高斯分布。方便起见，可以使用多个均匀分布随机变量之和近似高斯分布随机变量。三、主要仪器设备微型计算机、 Matlab 开发环境四、实验

11、步骤与操作方法1. 编程实现产生10000个N(3,4)高斯随机数。2. 计算生成随机数的 14 阶矩，最大值，最小值，频度直方图。五、实验数据记录和处理实验程序如下：m=3;n(3) = n(3) + (s(i)-avr)3;%求三n=4;阶n1=sqrt(n);n(4) = n(4) + (s(i)-avr)4;%求四pi=3.1416;阶s=zeros(1,10000);endfor i=1:10000m=m/10000;a=sqrt(-2*log(rand);n=n/10000;b=2*pi*rand;disp('均值 = ',num2str(m(1) );s(i)=n

12、1*a*cos(b)+m;%生成10000 个disp('二阶原点矩=N(3, 4) 高斯随机数',num2str(m(2) );enddisp('三阶原点矩=figure',num2str(m(3) );plot(s)disp('四阶原点矩=sum=0;',num2str(m(4) );for i=1:10000disp('一阶中心矩 = ',num2str(n(1) );sum=sum+s(i);% 求所有数总数disp('二阶中心矩=end;',num2str(n(2) );avr=sum/10000;% 求

13、平均数disp('三阶中心矩=m = zeros(1,4);',num2str(n(3) );for i = 1 : 10000disp('四阶中心矩=m(1) = m(1) + s(i);%求均值',num2str(n(4) );m(2) = m(2) + s(i)2;%求二阶矩hist(s,100)m(3) = m(3) + s(i)3;%求三阶title('100个区间 ')% 显示频率m(4) = m(4) + s(i)4;%求四阶c=0;endd=1;n = zeros(1,4);for i=1:10000%求最大for i = 1 :

14、 10000if c<s(i)n(1) = n(1) + (s(i)-avr);%求一c=s(i);阶矩end;n(2) = n(2) + (s(i)-avr)2;%求end;二阶矩for i=1:10000%求最小if d>s(i)end;d=s(i);cend;d六、实验结果与分析运行程序， command窗口中显示的结果如下：均值 = 2.993二阶原点矩 = 12.9888三阶原点矩 = 62.522四阶原点矩 = 341.0162一阶中心矩 =-5.8249e-015二阶中心矩 = 4.0306三阶中心矩 = -0.48106四阶中心矩 = 49.8857c =10.21

15、22d =-4.5628实验得到图像如下：图 1.3七、讨论、建议、质疑本实验和上一个实验类似， 但需要用到高斯函数的编写方法，利用循环结构可以求得若干个随机数，通过直方图可以对高斯分布有更直观实际的了解。实验 3随机信号相关函数估计一、实验目的和要求掌握随机信号相关函数计算、相关分析及实现二、实验原理和内容根据自相关和互相关的定义， 自相关，互相关计算随机信号的自相关和互相关。三、主要仪器设备微型计算机、 Matlab 开发环境四、实验步骤与操作方法4. 产生高斯随机信号。5. 计算其自相关函数。6. 计算两个高斯随机信号的互相关函数。五、 实验数据记录和处理实验程序如下：m=3;N(3,

16、4)高斯随机数 sn=4;endn1=sqrt(n);q=zeros(1,10000);pi=3.1416;for i=1:10000Fs=1000;c=sqrt(-2*log(rand);s=zeros(1,10000);d=2*pi*rand;for i=1:10000q(i)=n*c*cos(d)+m;%生成10000 个a=sqrt(-2*log(rand);N(3, 4)高斯随机数 qb=2*pi*rand;ends(i)=n1*a*cos(b)+m;%生成 10000 个figure(1);subplot(211),plot(s),title('s(n)');%co

17、r2 lag2=xcorr(q);画出 s 的图像 sfigure(3);plot(lag2/Fs,cor2),title('q的自相关subplot(212),plot(q),title('q(n)');%函数 ');% 画出 q 的自相关函数画出 q 的图像cor3 lag3=xcorr(s,q);cor1 lag1=xcorr(s);figure(4);figure(2);plot(lag3/Fs,cor3),title('s与 q 的互相关函数 ');% 画出 s 与 q 的互相关函数plot(lag1/Fs,cor1),title(&#

18、39;s的自相关函数 ');% 画出 s 的自相关函数六、实验结果与分析运行程序，得到各图像如下：图 1.4 s和 q 的数据图像图 1.5 s的自相关函数图 1.6 q的自相关函数图 1.7 s与 q 的互相关函数七、讨论、建议、质疑本实验与实验2 相似，也要产生高斯随机数，关键在于自相关程序以及互相关程序的编写。通过对自相关以及互相关图像的观察，可以让我们对于自相关以及互相关有更全面的了解。大连理工大学实验预习报告学院（系）：信息与通信工程学院专业：电子信息工程班级：电子 1303姓名：实验时间：指导教师：李彤学号 ：201383081组：_2015.12.15实验室：C221实验

19、台：实验 II ：系统对随机信号响应的统计特性分析、功率谱分析及应用实验实验 4 随机信号的功率谱分析1一、实验目的和要求掌握直接法估计随机信号功率谱。二、实验原理和内容根据随机信号的 Fourier变换结果，计算信号的功率谱。功率谱与频谱关系：三、实验步骤1. 生成高斯白噪声的随机信号。2. 通过 Fourier 变换计算高斯白噪声的功率谱，并绘图。实验 5 随机信号的功率谱分析方法2一、实验目的和要求掌握间接法估计随机信号功率谱。二、实验原理和内容根据维纳 - 辛钦定理，计算信号的功率谱。即三、实验步骤1. 生成高斯白噪声的随机信号。2. 计算高斯白噪声的自相关函数。3. 通过计算高斯白噪

20、声自相关函数的 Fourier 变换，得到噪声功率谱并绘图。实验 6 系统对随机信号响应的统计特性分析及仿真1一、实验目的和要求掌握时域上系统对随机信号响应的统计特性分析及仿真实现。二、实验原理和内容根据系统卷积性质，计算系统输出信号的统计特性。有如下性质：三、实验步骤1. 生成均匀分布的随机信号。2. 计算均匀分布随机信号通过平均低通滤波器。3. 计算输出信号均值、方差等统计特性。实验 7 系统对随机信号响应的统计特性分析及仿真2一、实验目的和要求掌握频域上系统对随机信号响应的统计特性分析及仿真实现。二、实验原理和内容根据卷积定理，计算系统输出信号的统计特性。即：三、实验步骤1. 生成均匀分

21、布的随机信号。2. 频域上计算均匀分布随机信号通过平均低通滤波器。3. 计算输出信号均值、方差等统计特性。大连理工大学实验报告学院（系）：信息与通信工程学院专业：电子信息工程班级：电子 1303姓名：实验时间：指导教师：李彤学号 ：201383081组：_2015.12.15实验室：C221实验台：实验 II ：系统对随机信号响应的统计特性分析、功率谱分析及应用实验实验 4 随机信号的功率谱分析1一、实验目的和要求掌握直接法估计随机信号功率谱。二、实验原理和内容根据维纳 - 辛钦定理，计算信号的功率谱。即三、主要仪器设备微型计算机、 Matlab 开发环境四、实验步骤与操作方法1. 生成高斯白

22、噪声的随机信号。2. 通过 Fourier 变换计算高斯白噪声的功率谱，并绘图。五、实验数据记录和处理程序如下：n=65536;+a(2)+randn(1,n);fs=2000;pi=3.14159;b1=fft(xn);t=(0:n-1)/fs;%b2=cos(2*pi/n)a=random('unif',0,1,1,2)*2*pi;b2=abs(b1).2/n;xn=cos(2*pi*30*t+a(1)+3*cos(2*pi*100*tb3=log(b2(1:n/2);f=(0:n/2-1)*fs/n;title('功率谱密度 ');figure,plot(

23、f,10*b3);六、实验结果与分析实验得到图表如下：图 2.1七、讨论、建议、质疑通过本实验，我掌握了直接估计随机信号功率谱的方法。实验 5 随机信号的功率谱分析2一、实验目的和要求掌握间接法估计随机信号功率谱。二、实验原理和内容根据随机信号的 Fourier变换结果，计算信号的功率谱。功率谱与频谱关系：三、主要仪器设备微型计算机、 Matlab 开发环境四、实验步骤与操作方法1. 生成高斯白噪声的随机信号。2. 计算高斯白噪声的自相关函数。3. 通过计算高斯白噪声自相关函数的 Fourier 变换，得到噪声功率谱并绘图。五、实验数据记录和处理程序如下：n=65536;b1=xcorr(xn

24、,'biased');fs=2000;b2=fft(b1);t=(0:n-1)/fs;b3=abs(b2);a=random('unif',0,1,1,2)*2*pi;f=(0:n-1)*fs/n/2;xn=cos(2*pi*30*t+a(1)+3*cos(2*pi*1figure,plot(f,10*log(b3(1:n);00*t+a(2)+randn(1,n);title('功率谱密度 ');六、实验结果与分析实验得到图表如下：图 2.2七、讨论、建议、质疑本实验与上一个实验结果相同，但是用的方法却不同。通过本实验，我掌握了间接法估计随机信

25、号功率谱。实验 6 系统对随机信号响应的统计特性分析及仿真1一、实验目的和要求掌握时域上系统对随机信号响应的统计特性分析及仿真实现。二、实验原理和内容根据系统卷积性质，计算系统输出信号的统计特性。有如下性质：三、主要仪器设备微型计算机、 Matlab 开发环境四、实验步骤与操作方法1. 生成均匀分布的随机信号。2. 计算均匀分布随机信号通过平均低通滤波器。3. 计算输出信号均值、方差等统计特性。五、实验数据记录和处理程序如下：m=1048576;b=1;r=2045;x=zeros(1,10000); s(1)=12357;for i=2:10000s(i)=mod(s(i-1)*r+b,m)

26、;endmy=ms*sum(s); for i=1:10000h(i)=0.9i;endn1=1024;rs=zeros(1,1024);for m=0:(n1-1)for k=0:(n1-1)s=s/m;hist(s)r=zeros(1,1001);t=zeros(1,1001);ms=sum(s)/length(s);rs(m+1)=rx(m+1)+s(k+1)*s(k+m+1);endendp=0;rs=rs/1000;endfigureplot(rs(1:1000)ry(i+1)=ry(i+1)+rx(a+1)*h(i+1)*h(k+ry=zeros(1,100);1);for i=0

27、:99endfor j=0:99endfor k=0:99enda=i+j-k;ry=ry/200;if(a<0)figurea=-a;plot(ry(1:100)均值=2.4725e+03方差 =-6.1134e+06六、实验结果与分析实验得到图表如下：图 2.3图 2.4图 2.5七、讨论、建议、质疑通过本实验我对时域上系统对随机信号响应的统计特性分析的了解更加深刻，并掌握了对随机信号响应的仿真。实验 7 系统对随机信号响应的统计特性分析及仿真2一、实验目的和要求掌握频域上系统对随机信号响应的统计特性分析及仿真实现。二、实验原理和内容根据卷积定理，计算系统输出信号的统计特性。即：三、

28、主要仪器设备微型计算机、 Matlab 开发环境四、实验步骤与操作方法1. 生成均匀分布的随机信号。2. 频域上计算均匀分布随机信号通过平均低通滤波器。3. 计算输出信号均值、方差等统计特性。五、实验数据记录和处理程序如下：n=500;w=(1:n)/n;xt=random('norm',0,1,1,n);t=(-n:n-1)/n*(n/20000);ht=fir1(500,0.3 0.4);subplot(4,1,1);plot(w,abs(sxx(1hw=fft(ht,2*n);:n);rxx=xcorr(xt,'biased');subplot(4,1,2

29、);plot(w,abs(hw2(1sxx=abs(fft(xt,2*n).2)/(2*n);:n);hw2=abs(hw).2;subplot(4,1,3);plot(w,abs(syy(1syy=sxx.*hw2;:n);ryy=fftshift(ifft(syy);subplot(4,1,4);plot(t,ryy);六、实验结果与分析实验得到图表如下：图 2.6七、讨论、建议、质疑通过本实验我对频域上系统对随机信号响应的统计特性分析的了解更加深刻，并掌握了对随机信号响应的仿真。大连理工大学实验预习报告学院（系）：信息与通信工程学院专业：电子信息工程班级：电子 1303姓名：实验时间：指

30、导教师：李彤学号 ：201383081组：_2015.12.17实验室：C219实验台：实验 III ：通信信号统计特性分析、仿真实现及应用实验 8 Hilbert变换算法实现1一、实验目的和要求掌握Hilbert变换的时域实现。二、实验原理和内容根据信号的卷积性质，进行Hilbert变换。Hilbert变换等效系统单位响应为三、实验步骤生成正弦波信号，卷积得到Hilbert变换信号。生成余弦波信号，卷积得到Hilbert变换信号。实验 9 Hilbert变换算法实现2一、实验目的和要求掌握Hilbert变换的频域实现。二、实验原理和内容根据信号的卷积定理，在频域进行Hilbert变换。 Hi

31、lbert变换等效系统函数为三、实验步骤生成正弦波信号，在频域相乘得到Hilbert变换信号。生成余弦波信号，在频域计算得到Hilbert变换信号。实验10窄带平稳随机信号统计特性分析一、实验目的和要求掌握时域窄带平稳随机信号统计特性分析方法。二、实验原理和内容根据窄带平稳随机信号的特点，进行信号分析。窄带过程可以表示为 Ac(n)cos( 0n)-As(n)sin( 0n) 。可以通过提取载波，相干解调得到x(n)=Ac(n) 和 As(n) ，进一步求得包络和相位三、实验步骤1. 按照实验 2中方法得到高斯白噪声 xc(n) 和 xs(n) ，设定这两个信号采样频率为10MHz。2. 用M

32、ATLAB的fdatool 函数生成采样频率为 10MHz，通带频率为 1KHz、阻带频率为10KHz、通带波动为 1dB、阻带衰减为 60dB 的低通滤波器 h(n) 。3. 用低通滤波器 h1(n) 分别对 xc(n) 和 xs(n) 滤波得到 Ac(n) 和As(n) 。4. 用载波频率为 0.2MHz的余弦、正弦波 ( 等效数字角频率 0 = 0.04 ) 分别对 Ac(n) 和As(n) 调制，并组合为 x(n) 。5. 用载波频率为 0.2MHz 的余弦、正弦波对 x(n) 进行解调，并分别通过 h(n) 得到 Ac(n)和 As(n) 的估计值 Ac?(n) As?(n)。6.根

33、据，计算包络、相位。7. 计算包络分布和相位概率分布。大连理工大学实验报告学院（系）：信息与通信工程学院专业：电子信息工程班级：电子 1303姓名：实验时间：指导教师：李彤学号 ：201383081组：_2015.12.17实验室：C219实验台：实验 III ：通信信号统计特性分析、仿真实现及应用实验 8 Hilbert变换算法实现1六、实验目的和要求掌握Hilbert变换的时域实现。七、实验原理和内容根据信号的卷积性质，进行Hilbert变换。Hilbert变换等效系统单位响应为八、主要仪器设备微型计算机、 Matlab开发环境九、实验步骤与操作方法生成正弦波信号，卷积得到Hilbert变

34、换信号。生成余弦波信号，卷积得到Hilbert变换信号。十、实验数据记录和处理程序如下：h=zeros(1,15);x=cos(pi/10*n);g=zeros(1,15);for i=1:15;n=1:400;h(i)=1/(i*pi)-(-1)*i*1/(i*pi);endg(j)=1/(j*pi)-(-1)*j*1/(j*pi);y=conv(x,h);endsubplot(2,1,1)b=conv(a,g);plot(x);subplot(2,1,1)title('余弦信号 ');figuresubplot(2,1,2)plot(a);plot(y);title(

35、9;正弦信号 ');title('余弦 Hilbert变换 ');subplot(2,1,2)a=sin(pi/10*n);plot(b);for j=1:15;title('正弦 Hilbert变换 ');十一、实验结果与分析生成的余弦信号及其Hilbert变换信号如下图：图 3.1余弦信号及其Hilbert变换生成的正弦信号及其Hilbert变换信号如下图：图 3.2 正弦信号及其Hilbert变换七、讨论、建议、质疑本次实验加深了我对希尔伯特变换的理解，并且编程实现了正弦和余弦信号希尔伯特变换的时域实现。实验 9 Hilbert变换算法实现2一、实

36、验目的和要求掌握Hilbert变换的频域实现。二、实验原理和内容根据信号的卷积定理，在频域进行Hilbert变换。 Hilbert变换等效系统函数为三、主要仪器设备微型计算机、 Matlab开发环境四、实验步骤与操作方法生成正弦波信号，在频域相乘得到Hilbert变换信号。生成余弦波信号，在频域计算得到Hilbert变换信号。五、实验数据记录和处理程序如下：Y=zeros(1,400);H=zeros(1,400);n=1:400;a=sin(pi/10*n);x=cos(pi/10*n);subplot(2,1,1)subplot(2,1,1)figure,plot(a);plot(x);t

37、itle('正弦信号 ');title('正弦 Hilbert变换 ');G=fft(a,400);X=fft(x,400);for k=2:200for k=2:200H(k)=-1i*G(k);Y(k)=-1i*X(k);endendfor k=202:400for k=202:400H(k)=1i*G(k);Y(k)=1i*X(k);endendH(1)=G(1);Y(1)=X(1);H(201)=G(201);Y(201)=X(201);b=ifft(H);y=ifft(Y);subplot(2,1,2)subplot(2,1,2)plot(b);plot(y);title('正弦 Hilbert变换 ');title('余弦 Hilbert变换 ');六、实验结果与分析生成的余弦信号及其Hilbert变换信号如下图：图 3.3余弦信号及其Hilbert变换生成的正弦信号及其Hilbert变换信号如下图：图 3.4正弦信号及其Hilbert变换七、讨论、建议、质疑本实验与上一个实验类似，也是实现正弦和余弦信号的