资源描述：
《随机过程实验四_线性系统参数估计与随机过程预测》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、.......HUNANUNIVERSITY课程实验报告题目：线性系统参数估计及随机过程预测学习参考.......实验四线性系统参数估计及随机过程预测一、实验目的通过本仿真实验了解基于随机过程的线性系统参数的估计方法以及基于线性系统模型的随机过程预测方法；培养计算机编程能力。二、实验要求采用MATLAB或VB语言进行编程(1)运用正态分布随机数产生函数产生均值为零、根方差s=1的白色噪声样本序列[或可参考实验1的正态分布产生方法]{u(n)
2、n=1,2,…,2000}；画出噪声u(n)的波形图(2)设离散时间线性系统的差分方程为画出x(n)的波形图。(3)假设已知线
3、性系统为二阶全极点系统，参数未知，满足以x(n)(n=3,4,…,1500)为已知数据，估计系统参数观察a,b与0.9、-0.2的相近性及估计误差(4)利用系统参数的估计值以及已获取的数据，采用单步递推预测方法对随机过程x(n)在区间nÎ[1501,2500]的值进行预测在x(n)的波形图上用不同的颜色画出y(n)的波形图，，观察和比较在[1501,2000]区间上二者的相近性及差异性。三、程序代码clc;%(1)画出噪声u(n)的波形图a=1;N=2000;学习参考.......u=normrnd(0,a,1,N);%-------------normrnd为正态
4、分布随机数产生函数figure(1);%----------------------画出第一个图stem([1:N],u,'.r');%-------------用红色画出1到N范围内的随机数分布图gridon;%------------------------画出网格title('u(n)波形');%(2)画出x(n)的波形图x=zeros(1,1998);x(1)=u(1);x(2)=0.9*x(1)+u(2);forn=3:Nx(n)=0.9*x(n-1)-0.2*x(n-2)+u(n);endfigure(2);stem([1:N],x,'.r');grid
5、on;title('x(n)波形');%(3)以x(n)(n=3,4,…,1500)为已知数据，估计系统参数m1=0;m2=0;m3=0;m4=0;n1=0;n2=0;%--赋初值forn=3:1500m1=m1+x(n-1)^2;m2=m2+x(n-1)*x(n-2);m3=m3+x(n-1)*x(n-2);m4=m4+x(n-2)^2;%-------------第一个矩阵的四个表达式n1=n1+x(n)*x(n-1);n2=n2+x(n)*x(n-2);%----------第二个矩阵的两个表达式endA=[m1,m2;m3,m4];A1=inv(A);%--
6、--------------------求A的逆矩阵B=[n1;n2];C=A1*B;%-------------------------将两个矩阵相乘a=C(1)b=C(2)%--------------------------将所求得值赋值给a,b%(4)在x(n)的波形图上用不同的颜色画出y(n)的波形图y=zeros(1,2000);n=[1501:2000];y(n)=a*x(n-1)+b*x(n-2);figure(3);stem(x,'.r');gridon;holdon;stem(y,'.b');%--------------------在x(n)的
7、波形图上用蓝色画出y(n)的波形图gridon;三、实验结果（1）学习参考.......（2）学习参考.......（3）（4）五、实验心得通过本次试验，我了解基于随机过程的线性系统参数的估计方法以及基于线性系统模型的随机过程预测方法，进一步熟悉了MATLAB的编程语言。1.若不给自己设限，则人生中就没有限制你发挥的藩篱。2.若不是心宽似海，哪有人生风平浪静。在纷杂的尘世里，为自己留下一片纯静的心灵空间，不管是潮起潮落，也不管是阴晴圆缺，你都可以免去浮躁，义无反顾，勇往直前，轻松自如地走好人生路上的每一步3.花一些时间，总会看清一些事。用一些事情，总会看清一些人。有
8、时候觉得自己像个神经病。既纠结了自己，又打扰了别人。努力过后，才知道许多事情，坚持坚持，就过来了。4.岁月是无情的，假如你丢给它的是一片空白，它还给你的也是一片空白。岁月是有情的，假如你奉献给她的是一些色彩，它奉献给你的也是一些色彩。你必须努力，当有一天蓦然回首时，你的回忆里才会多一些色彩斑斓，少一些苍白无力。只有你自己才能把岁月描画成一幅难以忘怀的人生画卷。学习参考
