大连理工大学随机信号试验报告完整

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

2、对M取模。为保证s(n)的周期为Mr的取值应满足r=4k+1,M2p,k与p的选取应满足：rMr(M1)+1工(m)x(n.m)&5)=根据自相关和互相关的定义，自相关，互相关+in)-计算随机信号的自相关和互相关。三、实验步骤1 .产生高斯随机信号。2 .计算其自相关函数3 .计算两个高斯随机信号的互相关函数大连理工大学实验报告学院(系)：信息与通信工程学院专业:电子信息工程班级:电子1303姓名：李彤学号:201383081组：实验时间：2015.12.14实验室：C221实验台：指导教师：实验 I：随机信号的产生、相关分析及其应用实验实验1均匀分布白噪声的生成实验目的和要求基于均

3、匀分布伪随机数，掌握均匀分布白噪声典型生成方法。实验原理和内容较简单的伪随机序列产生方法是采用数论中基于数环理论的线性同余法(乘同余法、混合同余法)，其迭代公式的一般形式为f(x)=(r-x+b)ModM其离散形式为s(n+1)=r-s(n)+bModM其中，s(n)为n时刻的随机数种子，r为扩展因子，b为固定扰动项，M为循环模，ModM表示对M取模。为保证s(n)的周期为Mr的取值应满足r=4k+1,M2p,k与p的选取应满足：rM,r(M1)+1231-1。本实验中参数取值为s(1)=12357,r=2025,b=1,M=104857&三、主要仪器设备微型计算机、Matlab开发环

4、境四、实验步骤与操作方法1 .编程实现产生10000个在(0,1)区间均匀分布随机数。2 .计算生成随机数的14阶矩，最大值，最小值，频度直方图五、实验数据记录和处理程序如下：M=1048576;b=1;r=2025;s=zeros(1,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,40)t

5、itle(40个区间)sum=0;fori=1:10000sum=sum+s(i);%求所有数的总和end;avr=sum/10000;%求所有数的平均数m=zeros(1,4);fori=1:10000m(1)=m(1)+s(i);%求均值m(2)=m(2)+s(i)A2;%求二阶矩m(3)=m(3)+s(i)A3;%求三阶m(4)=m(4)+s(i)A4;%求四阶endn=zeros(1,4);fori=1:10000n(1)=n(1)+(s(i)-avr);%求均值n(2)=n(2)+(s(i)-avr)A2;%求二阶矩n(3)=n(3)+(s(i)-avr)A3;%求三阶n(4)=n(

6、4)+(s(i)-avr)A4;%求四阶endm=m/10000;n=n/10000;disp(均值=,num2str(m(1);disp(二阶原点矩=,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);disp(方差=,num2str(var(s);c=0;d=1;fori=1:10000%求最大值ifcs(i)d=s(i);

7、end;end;cd六、实验结果与分析运行程序，command窗口中显示的结果如下:均值=0.49699二阶原点矩=0.32916三阶原点矩=0.24551四阶原点矩=0.19546均值=-8.3666e-017二阶中心矩=0.08217三阶中心矩=0.00024392四阶中心矩=0.012193方差=0.082179c=1.0000d=0实验得到图表如下:全部数据连线全部数据连线图 1.13004件区间件区间七、讨论、建议、质疑本实验中编写了生成随机序列的程序，通过设定不同的参数值可以得到不同的随机序列，通过计算多阶原点矩、中心矩和绘制直方图可以更加清楚的看到产生的随机序列的特点，对随机序列

8、的理解更加深刻。实验2高斯分布白噪声的生成一、实验目的和要求基于均匀分布伪随机数，掌握高斯分布白噪声典型生成方法。二、实验原理和内容1 .变换法2 .较简单的高斯白噪声产生方法是基于概率论中的中心极限定理。 即无穷多个同分布随机变量之和构成随机变量服从高斯分布。方便起见，可以使用多个均匀分布随机变量之和近似高斯分布随机变量。三、主要仪器设备微型计算机、Matlab开发环境四、实验步骤与操作方法1.编程实现产生10000个N(3,4)高斯随机数2.计算生成随机数的14阶矩，最大值,五、实验数据记录和处理实验程序如下：m=3;n=4;n1=sqrt(n);pi=3.1416;s=zeros(1,1

9、0000);fori=1:10000a=sqrt(-2*log(rand);b=2*pi*rand;s(i)=n1*a*cos(b)+m;%生成10000个N(3,4)高斯随机数endfigureplot(s)sum=0;fori=1:10000sum=sum+s(i);%求所有数总数end;avr=sum/10000;%求平均数m=zeros(1,4);fori=1:10000m=m+s(i);%求均值m(2)=m(2)+s(i)A2;%求二阶矩m(3)=m(3)+s(3;%求三阶m(4)=m(4)+s(i)A4;%求四阶endn=zeros(1,4);fori=1:10000n(1)=n(

10、1)+(s(i)-avr);%求一阶矩n(2)=n(2)+(s(i)-avr)A2;%求二阶矩六、实验结果与分析最小值，频度直方图n(3)=n(3)+(s(i)-avr)A3;%求三阶n(4)=n(4)+(s(i)-avr)A4;%求四阶endm=m/10000;n=n/10000;disp(均值=,num2str(m(1);disp(二阶原点矩=,num2str(m(2);disp(三阶原点矩=,num2str(m(3);disp(四阶原点矩=,num2str(m(4);disp(一阶中心矩=,num2str(n(1);disp(二阶中心矩=,num2str(n(2);disp(三阶中心矩=

11、,num2str(n(3);disp(四阶中心矩=,num2str(n(4);hist(s,100)title(,100个区间)显示频率c=0;d=1;fori=1:10000%求最大ifcs(i)d=s(i);end;end;cd运行程序，command窗口中显示的结果如下:均值=2.99310.2122-4.5628实验得到图像如下:七、讨论、建议、质疑本实验和上一个实验类似，但需要用到高斯函数的编写方法，利用循环结构可以求得若干个随机数，通过直方图可以对高斯分布有更直观实际的了解。二 阶 原 点 矩三 阶 原 点 矩四 阶 原 点 矩一 阶 中 心 矩二 阶 中 心 矩三 阶 中 心 矩

12、四阶中心矩=12.9888=62.522=341.0162=-5.8249e-015=4.0306=-0.48106=49.8857实验3随机信号相关函数估计、实验目的和要求掌握随机信号相关函数计算、相关分析及实现.、实验原理和内容Ri(n)=Vx(m)X(U-Fni)氏 S)二根据自相关和互相关的定义，自相关彳，互相关工刀(帆)卜。十加)-计算随机信号的自相关和互相关。三、主要仪器设备微型计算机、Matlab开发环境四、实验步骤与操作方法4 .产生高斯随机信号。5 .计算其自相关函数。6 .计算两个高斯随机信号的互相关函数。五、实验数据记录和处理实验程序如下：m=3;n=4;n1=sqrt(

13、n);pi=3.1416;Fs=1000;s=zeros(1,10000);fori=1:10000a=sqrt(-2*log(rand);b=2*pi*rand;s(i)=n1*a*cos(b)+m;%生成10000个N(3,4)高斯随机数sendq=zeros(1,10000);fori=1:10000c=sqrt(-2*log(rand);d=2*pi*rand;q(i)=n*c*cos(d)+m;%生成10000个N(3,4)高斯随机数qendfigure(1);subplot(211),plot(s),title(s(n);%画出s的图像ssubplot(212),plot(q),t

14、itle(q(n);%画出q的图像cor1lag1=xcorr(s);figure(2);plot(lag1/Fs,cor1),title(s的自相关函数);%画出s的自相关函数六、实验结果与分析cor2lag2=xcorr(q);figure(3);plot(lag2/Fs,cor2),title(q的自相关函数,);%画出q的自相关函数cor3lag3=xcorr(s,q);figure(4);plot(lag3/Fs,cor3),title(s与q的互相关函数);%画出s与q的互相关函数运行程序,图 1.6q 的自相关函数图 1.7s 与 q 的互相关函数七、讨论、建议、质疑本实验与实验

15、2相似，也要产生高斯随机数，关键在于自相关程序以及互相关程序的编写。通过对自相关以及互相关图像的观察，可以让我们对于自相关以及互相关有更全面的了解。大连理工大学实验预习报告学院（系）：信息与通信工程学院专业:电子信息工程班级:电子1303姓名：李彤学号:201383081组：实验时间：2015.12.15实验室：C221实验台：指导教师：实验 II:系统对随机信号响应的统计特性分析、功率谱分析及应用实验实验4随机信号的功率谱分析1一、实验目的和要求掌握直接法估计随机信号功率谱。二、实验原理和内容根据随机信号的Fourier变换结果，计算信号的功率谱。功率谱与频谱关系：二|尸（或三、实验步骤1

16、.生成高斯白噪声的随机信号。2 .通过Fourier变换计算高斯白噪声的功率谱，并绘图。实验5随机信号的功率谱分析方法2一、实验目的和要求掌握间接法估计随机信号功率谱。二、实验原理和内容根据维纳-辛钦定理，计算信号的功率谱。即尸）二网七1三、实验步骤1 .生成高斯白噪声的随机信号。2 .计算高斯白噪声的自相关函数。3,通过计算高斯白噪声自相关函数的Fourier变换，得到噪声功率谱并绘图。实验6系统对随机信号响应的统计特性分析及仿真1实验目的和要求掌握时域上系统对随机信号响应的统计特性分析及仿真实现实验原理和内容根据系统卷积性质，计算系统输出信号的统计特性。有如下性质:*X啊啊= =用用于于火

17、火(m)=Wj-k)h(j)li(k)三、实验步骤1,生成均匀分布的随机信号。2,计算均匀分布随机信号通过平均低通滤波器。3,计算输出信号均值、方差等统计特性。实验7系统对随机信号响应的统计特性分析及仿真2一、实验目的和要求掌握频域上系统对随机信号响应的统计特性分析及仿真实现。二、实验原理和内容根据卷积定理，计算系统输出信号的统计特性。即：网h(t)*x(t)=H(w)X(w)三、实验步骤1,生成均匀分布的随机信号。2.频域上计算均匀分布随机信号通过平均低通滤波器。3,计算输出信号均值、方差等统计特性。大连理工大学实验报告学院（系）：信息与通信工程学院专业:电子信息工程班级:电子1303姓名：

18、李彤学号:201383081组：实验时间：2015.12.15实验室：C221实验台：指导教师：实验 II： 系统对随机信号响应的统计特性分析、 功率谱分析及应用实验实验4随机信号的功率谱分析1实验目的和要求掌握直接法估计随机信号功率谱。二、实验原理和内容根据维纳-辛钦定理，计算信号的功率谱。即尸（如）二1川%（3三、主要仪器设备微型计算机、Matlab开发环境四、实验步骤与操作方法1 .生成高斯白噪声的随机信号。2 .通过Fourier变换计算高斯白噪声的功率谱，并绘图五、实验数据记录和处理程序如下：n=65536;fs=2000;pi=3.14159;t=(0:n-1)/fs;a=rand

19、om(unif,0,1,1,2)*2*pi;)+randn(1,n);b1=fft(xn);%b2=cos(2*pi/n)b2=abs(b1).A2/n;b3=log(b2(1:n/2);xn=cos(2*pi*30*t+a(1)+3*cos(2*pi*100*t+a(2)f=(0:n/2-1)*fs/n;title(功率谱密度,);figure,plot(f,10*b3);六、实验结果与分析七、讨论、建议、质疑通过本实验，我掌握了直接估计随机信号功率谱的方法。实验5随机信号的功率谱分析2一、实验目的和要求掌握间接法估计随机信号功率谱。二、实验原理和内容根据随机信号的Fourier变换结果，计

20、算信号的功率谱。功率谱与频谱关系:49)=|尸3三、主要仪器设备微型计算机、Matlab开发环境四、实验步骤与操作方法1 .生成高斯白噪声的随机信号。2 .计算高斯白噪声的自相关函数。3 .通过计算高斯白噪声自相关函数的五、实验数据记录和处理程序如下:n=65536;fs=2000;t=(0:n-1)/fs;a=random(unif,0,1,1,2)*2*pi;xn=cos(2*pi*30*t+a(1)+3*cos(2*pi*100*t+a(2)+randn(1,n);六、实验结果与分析实验得到图表如下：七、讨论、建议、质疑本实验与上一个实验结果相同，但是用的方法却不同。通过本实验，我掌握了

21、间接Fourier变换，得到噪声功率谱并绘图b1=xcorr(xn,biased);b2=fft(b1);b3=abs(b2);f=(0:n-1)*fs/n/2;figure,plot(f,10*log(b3(1:n);title(功率谱密度,);法估计随机信号功率谱。实验6系统对随机信号响应的统计特性分析及仿真实验目的和要求掌握时域上系统对随机信号响应的统计特性分析及仿真实现实验原理和内容根据系统卷积性质，计算系统输出信号的统计特性。有如下性质:三、主要仪器设备微型计算机、Matlab开发环境四、实验步骤与操作方法1 .生成均匀分布的随机信号。2 .计算均匀分布随机信号通过平均低通滤波器。3

22、 .计算输出信号均值、方差等统计特性。五、实验数据记录和处理程序如下：m=1048576;b=1;r=2045;x=zeros(1,10000);s(1)=12357;fori=2:10000s(i)=mod(s(i-1)*r+b,m);ends=s/m;hist(s)r=zeros(1,1001);t=zeros(1,1001);my=ms*sum(s);fori=1:10000h(i)=0.9Ai;endn1=1024;rs=zeros(1,1024);form=0:(n1-1)fork=0:(n1-1)rs(m+1)=rx(m+1)+s(k+1)*s(k+m+1);endendms=su

23、m(s)/length(s);p=0;endry(i+1)=ry(i+1)+rx(a+1)*h(i+1)*h(k+1);endendendry=ry/200;figureplot(ry(1:100)均值=2.4725e+03方差=-6.1134e+06六、实验结果与分析实验得到图表如下：图 2.3rs=rs/1000;figureplot(rs(1:1000)ry=zeros(1,100);fori=0:99forj=0:99fork=0:99a=i+j-k;if(a0)a=-a;0.360.340.320.30.2S图 2.4图 2.5七、讨论、建议、质疑通过本实验我对时域上系统对随机信号响

24、应的统计特性分析的了解更加深刻，并掌握了对随机信号响应的仿真。实验7系统对随机信号响应的统计特性分析及仿真2一、实验目的和要求掌握频域上系统对随机信号响应的统计特性分析及仿真实现。二、实验原理和内容根据卷积定理，计算系统输出信号的统计特性。即：Flh(t)*x(t)=H(w)X(w)三、主要仪器设备微型计算机、Matlab开发环境四、实验步骤与操作方法1 .生成均匀分布的随机信号。2 .频域上计算均匀分布随机信号通过平均低通滤波器。3 .计算输出信号均值、方差等统计特性。五、实验数据记录和处理程序如下:ryy=fftshift(ifft(syy);xt=random(norm,0,1,1,n)

25、;n=500;w=(1:n)/n;hw=fft(ht,2*n);subplot(4,1,1);plot(w,abs(sxx(1:n);rxx=xcorr(xt,biased);subplot(4,1,2);plot(w,abs(hw2(1:n);sxx=abs(fft(xt,2*n).A2)/(2*n);subplot(4,1,3);plot(w,abs(syy(1:n);hw2=abs(hw).A2;subplot(4,1,4);plot(t,ryy);syy=sxx.*hw2;六、实验结果与分析实验得到图表如下:ht=fir1(500,0.30.4);t=(-n:n-1)/n*(n/200

26、00);图 2.6七、讨论、建议、质疑通过本实验我对频域上系统对随机信号响应的统计特性分析的了解更加深刻，并掌握了对随机信号响应的仿真。大连理工大学实验预习报告学院（系）：信息与通信工程学院专业:电子信息工程班级:电子1303姓名：李彤学号:201383081组：实验时间：2015.12.17实验室：C219实验台：指导教师：实验 III:通信信号统计特性分析、仿真实现及应用实验8Hilbert变换算法实现1一、实验目的和要求掌握Hilbert变换的时域实现。二、实验原理和内容根据信号的卷积性质，进行Hilbert变换。Hilbert变换等效系统单位响应为0n=2k%,开=2上+】IH7T三、

27、实验步骤生成正弦波信号，卷积得到Hilbert变换信号。生成余弦波信号，卷积得到Hilbert变换信号。实验9Hilbert变换算法实现2一、实验目的和要求掌握Hilbert变换的频域实现。二、实验原理和内容根据信号的卷积定理，在频域进行Hilbert变换。Hilbert变换等效系统函数为三、实验步骤生成正弦波信号，在频域相乘得到Hilbert变换信号。生成余弦波信号，在频域计算得到Hilbert变换信号。实验10窄带平稳随机信号统计特性分析一、实验目的和要求掌握时域窄带平稳随机信号统计特性分析方法。二、实验原理和内容根据窄带平稳随机信号的特点，进行信号分析。窄带过程可以表示为x(n)=Ac(

28、n)cos(0n)-As(n)sin(0n)。可以通过提取载波，相干解调得到Ac(n)和As(n),进一=yJ；(n)+乂式)Otn)=arc4步求得包络和相位,4三、实验步骤1 .按照实验2中方法得到高斯白噪声xc(n)和xs(n),设定这两个信号采样频率为10MHz02 .用MATLAB的fdatool函数生成采样频率为10MHz,通带频率为1KHz、阻带频率为10KHz、通带波动为1dB、阻带衰减为60dB的低通滤波器h(n)。3 .用低通滤波器h1(n)分别对xc(n)和xs(n)滤波得到Ac(n)和As(n)。4 .用载波频率为0.2MHz的余弦、正弦波(等效数字角频率CD0=0.0

29、4冗)分别对Ac(n)和As(n)调制，并组合为x(n)05 .用载波频率为0.2MHz的余弦、正弦波对x(n)进行解调，并分别通过h(n)得到Ac(n)和As(n)的估计值Ac?(n)As?(n)。5)t叫打卜arcUtiJ6.根据4S)，计算包络、相位。7 .计算包络分布和相位概率分布。大连理工大学实验报告学院(系)：信息与通信工程学院专业:电子信息工程班级:电子1303姓名：李彤学号:201383081组：实验时间：2015.12.17实验室：C219实验台：指导教师：实验 III：通信信号统计特性分析、仿真实现及应用实验8Hilbert变换算法实现1六、实验目的和要求掌握Hilbert

30、变换的时域实现七、实验原理和内容根据信号的卷积性质，进行Hilbert变换。Hilbert变换等效系统单位响应为八、主要仪器设备微型计算机、Matlab开发环境九、实验步骤与操作方法生成正弦波信号，卷积得到Hilbert变换信号。生成余弦波信号，卷积得到Hilbert变换信号。十、实验数据记录和处理程序如下：h=zeros(1,15);fori=1:15;g=zeros(1,15);h(i)=1/(i*pi)-(-1)*i*1/(i*pi);n=1:400;endx=cos(pi/10*n);y=conv(x,h);subplot(2,1,1)plot(x);title(余弦信号);subpl

31、ot(2,1,2)plot(y);title(余弦Hilbert变换)；a=sin(pi/10*n);forj=1:15;g(j)=1/(j*pi)-(-1)*j*1/(j*pi);十一、实验结果与分析生成的余弦信号及其Hilbert变换信号如下图:余弦信号余弦信号图 3.1 余弦信号及其Hilbert变换endb=conv(a,g);subplot(2,1,1)figureplot(a);title(正弦信号);subplot(2,1,2)plot(b);title(正弦Hilbert变换)2?06010015020D250300350400450生成的正弦信号及其Hilbert变换信号如下

32、图:正弦信号正弦信号o50Too150200-25030035040U正弦正弦Hilbert变换变换410501001502002503003S0400450图 3.2 正弦信号及其 Hilbert 变换七、讨论、建议、质疑本次实验加深了我对希尔伯特变换的理解， 并且编程实现了正弦和余弦信号希尔伯特变换的时域实现。实验9Hilbert变换算法实现2一、实验目的和要求掌握Hilbert变换的频域实现。二、实验原理和内容根据信号的卷积定理，在频域进行Hilbert变换。Hilbert变换等效系统函数为碌。卜厂”二j一刀0三、主要仪器设备微型计算机、Matlab开发环境四、实验步骤与操作方法生成正弦

33、波信号，在频域相乘得到生成余弦波信号，在频域计算得到五、实验数据记录和处理程序如下：Y=zeros(1,400);n=1:400;x=cos(pi/10*n);subplot(2,1,1)plot(x);title(正弦Hilbert变换)；X=fft(x,400);fork=2:200丫(k)=-1i*X(k);endfork=202:400Y(k)=1i*X(k);endY(1)=X(1);Y(201)=X(201);y=ifft(Y);subplot(2,1,2)plot(y);title(余弦Hilbert变换)；六、实验结果与分析生成的余弦信号及其Hilbert变换信号如下图:Hil

34、bert变换信号Hilbert变换信号H=zeros(1,400);a=sin(pi/10*n);subplot(2,1,1)figure,plot(a);title(正弦信号);G=fft(a,400);fork=2:200H(k)=-1i*G(k);endfork=202:400H(k)=1i*G(k);endH(1)=G(1);H(201)=G(201);b=ifft(H);subplot(2,1,2)plot(b);title(正弦Hilbert变换)余寇信号余寇信号图 3.3 余弦信号及其Hilbert变换生成的正弦信号及其Hilbert变换信号如下图:正弦信号正弦信号正弦正弦Hilhert变变换换七、讨论、建议、质疑本实验与上一个实验类似，也是实现正弦和余弦信号的希尔伯特变换，但是本实验是在频域实现实验10窄带平稳随机信号统计特性分析实验目的和