大连理工大学随机信号实验二解说.doc

大连理工大学实验报告 学院（系）：信息与通信工程学院 专业： 电子信息工程 班级： 姓 名： 学号： 组： ___ 实验时间： 2015.12.15 实验室： C221 实验台： 指导教师： 实验II：系统对随机信号响应的统计特性分析、功率谱分析及应用实验 实验4 随机信号的功率谱分析1 实验目的和要求 掌握直接法估计随机信号功率谱。 实验原理和内容 根据维纳-辛钦定理，计算信号的功率谱。即 主要仪器设备 微型计算机、Matlab开发环境 实验步骤与操作方法 生成高斯白噪声的随机信号。 通过Fourier变换计算高斯白噪声的功率谱，并绘图。 实验数据记录和处理 程序如下 n=65536; fs=2000;pi=3.14159; 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); b1=fft(xn); %b2=cos(2*pi/n) b2=abs(b1).^2/n; b3=log(b2(1:n/2)); f=(0:n/2-1)*fs/n; figure,plot(f,10*b3); title(功率谱密度); 实验结果与分析 实验得到图表如下 图.1 七、讨论、建议、质疑 通过本实验，我掌握了直接估计随机信号功率谱的方法。 实验5 随机信号的功率谱分析2 实验目的和要求 掌握间接法估计随机信号功率谱。 实验原理和内容 根据随机信号的Fourier变换结果，计算信号的功率谱。功率谱与频谱关系： 主要仪器设备 微型计算机、Matlab开发环境 实验步骤与操作方法 生成高斯白噪声的随机信号。 计算高斯白噪声的自相关函数。 通过计算高斯白噪声自相关函数的Fourier变换，得到噪声功率谱并绘图。 实验数据记录和处理 程序如下 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); 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(功率谱密度);实验结果与分析 实验得到图表如下 图.2 七、讨论、建议、质疑 本实验与上一个实验结果相同，但是用的方法却不同。通过本实验，我掌握了间接法估计随机信号功率谱。 实验6系统对随机信号响应的统计特性分析及仿真1 实验目的和要求 掌握时域上系统对随机信号响应的统计特性分析及仿真实现。 实验原理和内容 根据系统卷积性质，计算系统输出信号的统计特性。?有如下性质： 主要仪器设备 微型计算机、Matlab开发环境 实验步骤与操作方法 生成均匀分布的随机信号。 计算均匀分布随机信号通过平均低通滤波器。 计算输出信号均值、方差等统计特性。 实验数据记录和处理 程序如下m=1048576; b=1; r=2045; x=zeros(1,10000); s(1)=12357; for i=2:10000 s(i)=mod(s(i-1)*r+b,m); end s=s/m; hist(s) r=zeros(1,1001); t=zeros(1,1001); ms=sum(s)/length(s); my=ms*sum(s); for i=1:10000 h(i)=0.9^i; end n1=1024; rs=zeros(1,1024); for m=0:(n1-1) for k=0:(n1-1) rs(m+1)=rx(m+1)+s(k+1)*s(k+m+1); end end p=0; rs=rs/1000; figure plot(rs(1:1000)) ry=zeros(1,100); for i=0:99 for j=0:99 for k=0:99 a=i+j-k; if(a0)