随机噪声的产生和性能测试.docxVIP

第 第 PAGE \* MERGEFORMAT 1 页 共 NUMPAGES \* MERGEFORMAT 7 页 成绩 信息与通信工程学院实验报告 （软件仿真性实验） 课程名称：随机信号分析 实验题目：随机噪声的产生与性能测试 实验目的和任务 1．编程完成上述随机数和随机信号的产生以及相应特征量的计算，并加注释； 2．理论计算上述特征量，并与编程计算结果比较分析； 3．完成实验报告，说明随机数产生、特征量计算的理论依据。 实验设备 计算机, Matlab软件 实验内容及步骤 1．产生满足均匀分布、高斯分布、指数分布、瑞利分布、的随机数，长度为N=1024，并计算这些书的均值、方差、自相关函数、概率密度函数、概率分布函数、功率谱密度，画出时域、频域特性曲线； 2．编程分别确定当五个均匀分布过程和5个指数分布分别叠加时，结果是否是高斯分布； 3．采用幅度为2， 频率为25Hz 的正弦信号为原信号，在其中加入均值为2 ， 方差为0.04 的高斯噪声得到混合随机信号，编程求 的均值、相关函数、协方盖函数和方差，并与计算结果进行比较分析。 实验原理 平稳时间序列的模型估计与预测原理 1、样本自协方差函数： 样本自相关函数： 2、样本偏相关函数 实验数据分析及处理 实验代码： N=1024; fs=1000; n=0:N-1; signal=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,unbiased); 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(1) subplot(2,3,1); hist(signal); title(均匀随机数直方图); subplot(2,3,2); plot(signal); title(时域特性曲线); subplot(2,3,3); f=n*fs/N; %频率序列 plot(abs(signal_power)); title(功率幅度谱特性曲线); subplot(2,3,4); plot(angle(signal_power)); title(功率谱相频特性曲线 subplot(2,3,6); plot(n1,k1);); subplot(2,3,5); plot(1:2047,signal_corr); title(自相关函数图); title(概率密度); figure(2) plot(n2,k2); 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]=ksdensity(signal); figure(1) subplot(1,2,1); hist(signal); title(叠加均匀分布随机数直方图); subplot(1,2,2); plot(n1,k1); title(叠加均匀分布的概率密度); clear all; clc; t=0:320; x=sin(2*pi*t*25); x1=wgn(1,321,0); 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