随机信号分析实验报告基于MATLAB语言Word文档格式.docx

1、泊松分布poissrnd正态分布normrnd离散均匀分布unidrnd瑞利分布raylrnd均匀分布unifrndchi2rnd3.随机序列的数字特征估计对于遍历过程，可以通过随机序列的一条样本函数来获得该过程的统计特征。这里我们假定随机序列X（n）为遍历过程，样本函数为x（n），其中n=0,1,2,N-1。那么，X（n）的均值、方差和自相关函数的估计为 利用MATLAB的统计分析函数可以分析随机序列的数字特征。（1）均值函数meanm = mean（x）返回按1.3式估计X（n）的均值，其中x为样本序列x（n）。（2）方差函数varsigma2 = var（x）返回按（1.4）式估计X（n

2、）的方差，其中x为样本序列x（n），这一估计为无偏估计。（3）互相关函数xcorrc = xcorr（x,y）c = xcorr（x）c = xcorr（x,y,opition）c = xcorr（x,xcorr（x,y）计算X（n）与Y（n）的互相关，xcorr（x）计算X（n）的自相关。option 选项可以设定为：biased 有偏估计unbiased 无偏估计coeff m = 0 时的相关函数值归一化为1none 不做归一化处理实验内容及实验结果1.采用线性同余法产生均匀分布随机数1000个，计算该序列均值和方差与理论值之间的误差大小。改变样本个数重新计算。程序代码：y=1;k=7;

3、N=1010;xn=;for i=1:1000 y=mod（y*k,N）; x=y/N; xn=xn x;endm=mean（xn）n=var（xn）me=0.5-mne=1/12-n实验结果：m = 0.4813n = 0.0847me= 0.0187ne= -0.0013改变样本数量重新计算：（理论值 m=0.5 n=1/12）样本数量m误差n1000.41640.08360.09010.00685000.46680.03320.08590.00260.48130.01870.08470.0013100000.49730.00270.08480.0015500000.50090.00090

4、.08370.00041000000.49960.00032.参数为的指数分布的分布函数为利用反函数法产生参数为0.5的指数分布随机数1000个，测试其方差和相关函数。j=1:1999; y=mod（y*k,N）; x=y/N; xn=xn x;y=（-2）*log（1-xn）;n=var（y）c=xcorr（y,）;plot（j-1000,c）;方差 n=3.7596自相关函数：3.产生一组N（1,4）分布的高斯随机数（1000个样本），估计该序列的均值、方差、和相关函数。i=1:1000;x=normrnd（1,2,1,1000）;m=mean（x）n=var（x）c=xcorr（x,su

5、bplot（211）;plot（i,x）;title（随机序列）;subplot（212）;title（自相关函数）;均值 m=1.0082方差 n=3.8418实验小结本次实验对随机数的生成做了练习。具体来说，包括线性同余法，生成已知分布函数的随机数，rand函数等，还有就是有关均值、方差、相关的调用函数。实验二 随机过程的模拟与数字特征1.学习利用 MATLAB模拟产生随机过程的方法。2.熟悉和掌握特征估计的基本方法及其 MATLAB实现。1.正态分布白噪声序列的产生MATLAB提供了许多产生各种分布白噪声序列的函数，其中产生正态分布白噪声序列的函数为randn。产生 mn的标准正态分布随

6、机数矩阵。如果要产生服从 如果 2.相关函数估计MATLAB提供了函数 xcorr用于自相关函数的估计。c = xcorr（x,y） c = xcorr（x） ） xcorr（x,y）计算 X （n）与 Y（n）的互相关，xcorr（x）计算 X （n）的自相关。Option 选项可以设定为： 有偏估计。 无偏估计。 m =0时的相关函数值归一化为1。 不做归一化处理。3.功率谱估计对于平稳随机序列 X （n），如果它的相关函数满足那么它的功率谱定义为自相关函数的傅里叶变换：功率谱表示随机信号频域的统计特性，有着重要的物理意义。我们实际所能得到的随机信号的长度总是有限的，用有限长度的信号所得的

7、功率谱只是真实功率谱的估计，称为谱估计或谱分析。功率谱估计的方法有很多种，以下是两种通用谱估计方法。（1）自相关法先求自相关函数的估计，然后对自相关函数做傅里叶变换。 其中N表示用于估计样本序列的样本个数。（2）周期图法先对样本序列 x（n）做傅里叶变换其中，则功率谱估计为MATLAB函数 periodogram实现了周期图法的功率谱估计。periodogram Pxx,w = periodogram（x） Pxx,w = periodogram（x,window） Pxx,w = periodogram（x,window,nfft） Pxx,f = periodogram（x,window,

8、nfft,fs） periodogram（.）实现周期图法的功率谱估计。其中：Pxx为输出的功率谱估计值；f为频率向量；w为归一化的频率向量；window代表窗函数，这种用法种对数据进行了加窗，对数据加窗是为了减少功率谱估计中因为数据截断产生的截断误差，表 2.1列出了产生常用窗函数的 MATLAB函数窗函数MATLAB函数矩形窗boxcarBlackman窗blackman三角窗triangChebyshev窗chebwinHanning窗hannBartlett窗bartlettHamminghammingKaiser窗kaiser nfft设定 FFT算法的长度；fs表示采样频率；如果不

9、指定输出参数（最后一种用法），则直接会出功率谱估计的波形。1.按如下模型产生一组随机序列是均值为1，方差为4的正态分布白噪声序列。估计过程的自相关函数和功率谱。w=normrnd（1,4,1,1024）;x（1）=w（1）;i=2;while i1025x（i）=0.8*x（i-1）+w（i）;i=i+1;R=xcorr（x）;S,W=periodogram（x）;subplot（3,1,1）; plot（x）;title（x（n）axis tight;subplot（3,1,2）; plot（R）;R（m）subplot（3,1,3）; plot（S）;S（W）2.设信号为，为正态分布白噪声

10、序列，试在N=256和N=1024点时，分别产生随机序列，画出的波形并估计的相关函数和功率谱。N=256：N=256;n=1:1:N;w=randn（1,N）;f1=0.05;f2=0.12;x=sin（2*pi*f1*n）+2*cos（2*pi*f2*n）+w（n）;N=1024：N=1024;x（n） N=256R（m） N=256S（W） N=256本次实验对随机序列的产生进行了复习，对自相关函数与功率谱密度的产生进行了练习。并且验证了自相关函数与功率谱密度呈傅立叶变换关系。实验三 随机过程通过线性系统的分析1.理解和分析白噪声通过线性系统后输出的特性。2.学习和掌握随机过程通过线性系统

11、后的特性，验证随机过程的正态化问题。1.白噪声通过线性系统设连续线性系统的传递函数为或，输入白噪声的功率谱密度为，那么系统输出的功率谱密度为输出自相关函数为输出相关系数为输出相关时间为输出平均功率为上述式子表明，若输入端是具有均匀谱的白噪声，则输出端随机信号的功率谱主要由系统的幅频特性决定，不再是常数。2.等效噪声带宽在实际中， 常常用一个理想系统等效代替实际系统的H（），因此引入了等效噪声带宽的概念，他被定义为理想系统的带宽。等效的原则是，理想系统与实际系统在同一白噪声的激励下，两个系统的输出平均功率相等，理想系统的增益等于实际系统的最大增益。实际系统的等效噪声带宽为3.线性系统输出端随机过

12、程的概率分布（1）正态随机过程通过线性系统若线性系统输入为正态过程，则该系统输出仍为正态过程。（2）随机过程的正态化随机过程的正态化指的是，非正态随机过程通过线性系统后变换为正态过程。任意分布的白噪声通过线性系统后输出是服从正态分布的；宽带噪声通过窄带系统，输出近似服从正态分布。设白噪声通过图 3.1所示的 RC电路，分析输出的统计特性。图3.1 RC电路（1）试推导系统输出的功率谱密度、相关函数、相关时间和系统的等效噪声带宽。经计算，结果如下：传输函数输出功率谱密度输出自相关函数为：相关时间为：等效噪声带宽为：若输入为高斯白噪声，则输出均值为0，方差，服从高斯分布，概率密度函数为：（2）采用

13、Matlab模拟正态分布白噪声通过上述RC电路，观察输入和输出的噪声波形以及输出噪声的概率密度。R=10;C=0.1;b=1/（R*C）;128;h=b*exp（-n*b）;x=randn（1,4096）;y=conv（x,h）;fy yi=ksdensity（y）;subplot（3,1,1）plot（x）;subplot（3,1,2）plot（y）;y（n）subplot（3,1,3）plot（fy）;fy运行结果：（3）模拟产生均匀分布的高斯白噪声通过上述RC电路，观察输入和输出的噪声波形以及输出噪声的概率密度。256;x=rand（1,1024）;x=x-0.5;y=conv（h,x）

14、; plot（y）; title（ plot（fy）; axis tight;（4）改变RC电路的参数（电路的RC值），重做（2）和（3），与之前的结果进行比较。a）改变RC的值为10000 0.1高斯输入：均匀分布输入：b）改变RC的值为 10 0.001非正态随机过程通过线性系统后变换为正态过程。低通系统通带过窄时，输出与输入差别很大，因为只有低频可以通过，高频量被抑制了。实验四 窄带随机过程的产生及其性能测试2.掌握窄带随机过程的特性，包括均值（数学期望）、方差、相关函数及功率谱密度等。1.窄带随机过程的莱斯表达式任何一个实平稳窄带随机过程X （t）都可以表示为上式称为莱斯表达式，根据上

15、式可以模拟产生窄带随机过程，具体过程如图4.1 所示。图4.1 窄带随机过程的产生2.窄带随机过程包络与相位的概率密度包络的概率密度为，服从瑞利分布。相位的概率密度为，呈均匀分布。3.窄带随机过程包络平方的概率密度包络平方的概率密度为，为指数概率密度函数。1.按图4.1所示结构框图，基于随机过程的莱斯表达式，用MATLAB产生一满足条件的窄带随机过程。程序代码1024;h=exp（-n）;c1=randn（1,1024）;a=conv（c1,h）;c2=randn（1,1024）;b=conv（c2,h）;fc=10000;i=1;x=zeros（1,1024）;x（i）=a（i）*cos（2

16、*pi*fc*i）-b（i）*sin（2*pi*fc*i）;x实验结果2.画出该随机过程的若干次实现，观察其形状。3.编写MATLAB程序计算该随机过程的均值函数、自相关函数、功率谱、包络、包络平方及相位的一维概率密度，画出相应的图形并给出解释。subplot（3,2,1）;%-均值&自相关函数-subplot（3,2,2）;plot（R）;R（x）%-功率谱密度-S,w=periodogram（x）;subplot（3,2,3）;plot（S）;S（w）%-包络-A=zeros（1,1024）;while（i1025） 利用循环，根据得到包络 A（i）=sqrt（a（i）2+b（i）2）;f

17、A j=ksdensity（A）;subplot（3,2,4）;plot（fA）;f（A）%-包络平方-1025） %利用循环得到包络的平方A2（i）=A（i）2;fA2 j=ksdensity（A2）;subplot（3,2,5）;plot（fA2）;f（A2）%-相位- i=1; while（i1025） %根据得到包络相位fai（i）=atan（b（i）/a（i）;fp j=ksdensity（fai）;subplot（3,2,6）;plot（fp）;f（p）m =0.0038 均值约等于0包络的概率密度近似为瑞利分布包络平方的概率密度近似为指数概率密度包络的相位近似为均匀分布，此时样本

18、数为1024如果加大样本数至32768，则相位的概率密度函数就更接近于均匀分布了通过本次实验我学会了基于随机过程的莱斯表达式产生窄带随机过程的方法。掌握了窄带随机过程的特性，包括均值、方差、相关函数及功率谱密度等。通过实验验证了窄带随机过程包络的概率密度分布为瑞丽分布，相位的概率密度分布为均匀分布，包络平方的概率密度分布为指数分布。实验总结总结这四次随机信号分析实验，我觉得收获最大的地方便是学会利用MATLAB进行随机信号分析。经过此次实验之后，我掌握了随机信号的初始化，信号的各种分析方法，同时也通过MATLAB的画图功能直观的了解了随机信号处理之后的结果。通过这次实验，我收获了很多，同时也非常感谢老师的教导，才有了这样的成果。