随机过程的模拟与特征估计

1、随机过程的模拟与特征估计,许可 ,手机民网Email 军网Email xukes4.gfkd.mtn,讲座提纲,Matlab简介 随机过程的模拟与特征估计 现场答疑,Matlab简介,什么是Matlab 如何获取和安装 参考书籍 我的一点经验,什么是Matlab,Matlab = Matrix + Laboratory 新墨西哥大学 Cleve Moler教授发明 目的：减轻学生编程负担 MathWorks公司的产品 最新版本是Matlab 2006a,Matlab常用的3个功能,*.m文件 GUI程序 Simulink仿真,如何获取和安装,在军网上下载 合训软件园，

2、Lab311，航天学院 安装之前请安装“java虚拟机” msjavx86 安装在线帮助文档,参考书籍,Matlab自带的帮助文档（best） 张志涌 精通MATLAB6.5版北京航空航天大学出版社 苏金明 MATLAB工程数学电子工业出版社,我的一点经验,没有必要追求高版本 多看help 多动手，多模仿 容易上手，比C语言简单,随机过程的模拟与特征估计,实验目的 实验原理 实验内容 程序详解 实验报告 补充,实验目的,了解随机过程特征估计的基本概念和方法 学会运用Matlab软件产生各种随机过程，对随机过程的特征进行估计 通过实验了解不同估计方法所估计出来的结果之间的差异。,实验原理,高斯白

3、噪声的产生 均值的估计 方差的估计 自相关函数的估计 功率谱的估计,高斯白噪声的产生,N=1000; x=randn(N,1); 高斯：概率分布 白噪声：功率谱,均值的估计,数学原理 mx=mean(x);,方差的估计,数学原理 vx=var (x);,自相关函数的估计,数学原理 rx=xcorr(x,x);,功率谱的估计,数学原理 px=pyulear(x,64);,实验内容,详见随机过程的模拟与特征估计 一共3+2个题目,程序详解,按如下模型产生一组随机序列， 其中w(n)为均值为1，方差为4的高斯白噪声序列，试估计随机过程的自相关函数和功率谱。,如何产生均值为1，方差为4的高斯白噪声序列

4、？ 如何产生x(n)？ 如何估计自相关函数？ 如何估计功率谱？ 如何画出美观，简洁的图像？,产生高斯白噪声序列,clear; clc; mu=1; %均值 sigma=2; %方差 N=1024; %产生的白噪声点数 w=mu+sigma*randn(N,1);,w的波形图,w的概率分布图像,产生x(n),x(1)=0; %初始化 for i=2:N x(i)=0.8*x(i-1)+w(i); end,x(n)的波形图,估计自相关函数,rx=xcorr(x,x);,自相关函数波形图,估计功率谱,px=pyulear(x,2);,功率谱的波形图,画出美观，简洁的图像,subplot(m,n,k)

5、; 在窗口中画m行，n列图像，当前图像是第k个 plot(n,x); 以为横坐标，为纵坐标连续画图 stem(n,x); 以为横坐标，为纵坐标离散画图,实验报告,实验内容和原理的简单阐述，分析； 得到的实验结果图形及简要分析，比较； 对“实验思考题”的详细分析和回答； 自己的实验心得；,补充,在*.m文件里面编程，不要在窗口编写 习惯使用help 习惯看看see also 学会看错误报告 学会debug,在*.m文件里面编程,习惯使用help, help stem STEM Discrete sequence or stem plot. STEM(Y) plots the data seque

6、nce Y as stems from the x axis terminated with circles for the data value. STEM(X,Y) plots the data sequence Y at the values specified in X. STEM(.,filled) produces a stem plot with filled markers. STEM(.,LINESPEC) uses the linetype specified for the stems and markers. See PLOT for possibilities. H = STEM(.) returns a vector of line handles. See also PLOT, BAR, STAIRS.,学会看错误报告,%test clear; clc; n=0:100; x=rand(100,1); plot(n,x);,? Error using = plot Vectors m