随机信号分析实验报告

《随机信号分析实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、WORD文档下载可编辑随机信号分析实验报告班级：13050141姓名：学号：日期：2016年5月24日专业资料整理分享WORD文档下载可编辑实验一随机噪声的产生与性能测试一、实验内容1．产生满足均匀分布、高斯分布、指数分布、瑞利分布的随机数，长度为N=1024，并计算这些数的均值、方差、自相关函数、概率密度函数、概率分布函数、功率谱密度，画出时域、频域特性曲线；2．编程分别确定当五个均匀分布过程和5个指数分布分别叠加时，结果是否是高斯分布；3．采用幅度为2，频率为25Hz的正弦信号为原信号，在其中加入均值为2，方差为0.04的高斯

2、噪声得到混合随机信号，编程求的均值、相关函数、协方差函数和方差，并与计算结果进行比较分析。二、实验步骤1.程序N=1024;fs=1000;n=0:N-1;signal=chi2rnd(2,1,N);%rand(1,N)均匀分布,randn(1,N)高斯分布,exprnd(2,1,N)指数分布,raylrnd(2,1,N)瑞利分布,chi2rnd(2,1,N)卡方分布signal_mean=mean(signal);signal_var=var(signal);signal_corr=xcorr(signal,signal,'un

3、biased');signal_density=unifpdf(signal,0,1);signal_power=fft(signal_corr);%[s,w]=periodogram(signal);[k1,n1]=ksdensity(signal);[k2,n2]=ksdensity(signal,'function','cdf');figure;hist(signal);title('频数直方图');figure;plot(signal);title('均匀分布随机信号曲线');f=n*fs/N;%频率序列figure;pl

4、ot(abs(signal_power));title('功率幅频');figure;plot(angle(signal_power));title('功率相频');figure;plot(1:2047,signal_corr);title('自相关函数');figure;专业资料整理分享WORD文档下载可编辑plot(n1,k1);title('概率密度');figure;plot(n2,k2);title('分布函数');结果（1）均匀分布专业资料整理分享WORD文档下载可编辑（2）高斯分布专业资料整理分享WORD文档下载可编

5、辑（3）指数分布专业资料整理分享WORD文档下载可编辑（4）瑞利分布专业资料整理分享WORD文档下载可编辑专业资料整理分享WORD文档下载可编辑（5）卡方分布专业资料整理分享WORD文档下载可编辑2.程序N=1024;signal_1=rand(1,N);signal_2=rand(1,N);signal_3=rand(1,N);signal_4=rand(1,N);signal_5=rand(1,N);signal=signal_1+signal_2+signal_3+signal_4+signal_5;[k1,n1]=ksde

6、nsity(signal);figure(1)subplot(1,2,1);hist(signal);title('叠加均匀分布随机数直方图');subplot(1,2,2);plot(n1,k1);title('叠加均匀分布的概率密度');专业资料整理分享WORD文档下载可编辑结果指数分布叠加均匀分布叠加结果：五个均匀分布过程和五个指数分布分别叠加时，结果是高斯分布。专业资料整理分享WORD文档下载可编辑3.程序clearall;clc;t=0:320;x=2*sin(2*pi*t*25);z=imnoise(x,'gaussi

7、an',2,0.04);%z=x+x1;y=trapz(t,z);%y=int(z,x,0,t);subplot(3,2,1),plot(z);title('随机信号序列')meany=mean(z);figure(1)subplot(3,2,3),plot(t,meany,'.');title('随机信号均值')vary=var(y);%方差subplot(3,2,4),plot(t,vary,'.');title('随机信号方差')cory=xcorr(z,'unbiased');%自相关函数subplot(3,2,2),pl

8、ot(cory);title('随机信号自相关函数')covv=cov(y);subplot(3,2,5),plot(t,covv,'.');title('随机信号协方差')t=[0:0.0005:0.045];X1=sin(2*pi*25*t);