随机过程-实验报告

  • 格式:doc
  • 大小:156.50 KB
  • 文档页数:10

随机过程试验报告

班级:信计09级01班

姓名:

学号:

实验一 1 实验题目

描绘出随机过程()cos()Xtxwt的图像

实验目的

利用MATLAB编程描绘出随机过程()cos()Xtxwt的图像

实验地点及时间 信息楼121机房 2012年5月31日

实验内容:描绘出随即过程X(t)=xcos(wt)的图像

程序如下:

x=0:0.1:2*pi;

t=0:0.1:2*pi;

y=x.*cos(4*t);

plot3(t,x,y);

axis square;

grid on;

结果:

实验总结

掌握应用随机过程的本质含义, 练习使用MATLAB描绘随机过程的三角函数图像,改变参数得到不同的随机过程图像。

实验成绩 评阅时间 评阅教师

实验二 2 实验题目 绘制随机相位正弦波()cos()Xtwt的均值,方差和自相关函数的图像

实验目的 通过绘制图像,深入理解随机相位正弦波的均值,方差和自相关函数

实验地点及时间 信息楼127机房 2012年6月1日

实验内容:绘制随机相位正弦波()cos()Xtwt的均值,方差和自相关函数的图像

实验习题

给出其程序与图像

>> w=1;a=1;

t=-pi:pi/100:pi;

b=-pi:pi/100:pi;

y=a.*cos(w*t+b);

>> plot3(t,b,y);

syms x y t;f=int(cos(t+pi/2),t,-pi,pi)

f =0

f=0;e=-pi:pi/100:pi;

plot(e,f) 3

s=-pi:pi/100:pi;

i=2;

plot(s,i);

k=-2*pi:pi/100:2*pi;a=1;

y=a.*cos(2*k);

plot(k,y);axis square;grid on; 4

实验总结

通过绘制()cos()Xtt 图像,由图像可知,随机相位正弦波的均值E(t)=0;方差Var(t)=2;由自相关函数图像的描述出不同时刻之间的相关程度。通过实验,我们可以知道,均值、方差、自相关函数是刻画随机过程的主要特性,因此,对于解决实际课题而言,常常能够起到重要作用。

实验成绩 评阅时间 评阅教师

实验三

实验题目 模拟Possion流

实验目的 用Matlab语言产生随机数,了解 Possion流

实验地点及时间 信息楼127机房 2012年6月4日 5 实验内容

用Matlab语言产生随机数,并编程实现possion流的模拟

程序:

U=rand(1,20);

a=2;

X=-a^(-1)*log(U);

S=zeros(1,22);

d=zeros(1,22);

S(1)=0;S(2)=X(1);

for n=3:21

S(n)=S(n-1)+X(n-1);

end

for i=0:21 %---

if 0<=i

d(i+1)=0;

else

for j=2:21

if (S(j)<=i)&(S(j+1)

d(i+1)=j;

end

end

end

end

plot(d)

结果 6

通过此次实验,模拟Possion流分布,运用MATLAB产生随机数,使我对泊松分布有了更深刻的理解。不过产生的随机数可以根据实际需要设置不同分布的随机数,便于解决实际问题。

实验成绩 评阅时间 评阅教师

实验四 7 实验题目 求Markov链的极限分布

实验目的 用Matlab语言求Markov遍历链的极限分布

实验地点及时间 信息楼127机房 2012年6月6日

实验内容

判定一个Markov链是否是遍历的,若是遍历的,求其极限分布。并能从实际问题中抽象出Markov链,并求出其极限分布,并理解其实际意义。

实验习题

1、已知齐次马氏链,0,1,2,Xnn的状态空间1,2,3E,状态转移矩阵为

1/21/31/61/31/31/31/31/21/6P

(1) 计算2步转移概率;(2) 已知初始分布为2/5,2/5,1/5P,求2X的分布律

(3) 求平稳分布,要求给出程序与结果。

程序:

S0=[2/5 2/5 1/5]; P=[1/2 1/3 1/6;1/3 1/3 1/3;1/3 1/2 1/6];

S2=S0*P.^2

P1=[P'-eye(3,3);1 1 1];

b=[0 0 0 1]';

T=P1\b

P2=P^2

结果:

S2 =

0.1667 0.1389 0.0611

T =

0.4000

0.3714

0.2286

P2 =

0.4167 0.3611 0.2222

0.3889 0.3889 0.2222

0.3889 0.3611 0.2500

8 (1) 2步转移概率

P2 =

0.4167 0.3611 0.2222

0.3889 0.3889 0.2222

0.3889 0.3611 0.2500

(2) 2X的分布律

S2 =

0.1667 0.1389 0.0611

(3) 平稳分布

T =

0.4000

0.3714

0.2286

2、为适应日益扩大的旅游事业的需要,某城市的A,B,C三个照相馆组成一个联营部,联合经营出租相机的业务,旅游者可由A,B,C三处任何一处租出相机,用完后还到A,B,C三处的任何一处即可.估计转移概率如表所示,今欲选择A,B,C之一附设租机维修点,问该点设在何处为好? (程序与结果)

还相机处

A B C

租相机处 A 0.2 0.8 0

B 0.8 0 0.2

C 0.1 0.3 0.6

程序:

p=[0.2 0.8 0;0.8 0 0.2; 0.1 0.3 0.6]; P2=p^2

a=[p'-eye(3);ones(1,3)];b=[0 0 0 1]';T=a\b

结果:

解:由题意可知,该问题的转移概率矩阵P为: 9 0.20.800.800.20.10.30.6P,20.680.160.160.180.70.120.320.260.42P

因为2P的所有元素都大于零,所以P为正规矩阵。当A,B,C三还相机处业务开展一定时期后,就会达到平衡条件,这样就可以得到一固定概率t,使得tPt成立,即

0.20.80,,10.800.2,,10.10.30.6xyxyxyxy 成立

上式展开,得

0.20.80.11xyxyx

0.80.31xxyy

0.20.611yxyxy

解上述联立方程式,得1741x,1641y

故 17168,,1,,414141xyxy

由上述计算可知,在稳定状态相机还到A处得概率为1741,在稳定状态相机还到B处得概率为1641,在稳定状态相机还到C处得概率为841,A处的概率最大,因此相机维修点设在A处是最佳得选择。

实验总结

熟悉理解markov链的概念及实际意义,并写出其转移概率矩阵,掌握markov链的定义及转移概率矩阵,掌握极限分布的存在的条件及求法,利用MATLAB求Markov遍历链的极限分布。

实验成绩 评阅时间 评阅教师