随机过程实验报告

随机过程试验报告班级:姓名:学号:实验一实验题目__MATLAB基本操作实验题目实验目的 了解MATLAB基础知识，掌握MATLAB基本操作实验地点及时间 信息楼实验内容1.MATLAB简介及基础知识；2.MATLAB的及矩阵的运算；3、MATLAB的基本流程结构；4、MATLAB的M文件、函数文件及常用函数实验习题1建立一个矩阵A，求矩阵A的转置B并求A与B的和与乘积；A的逆；解方程组AX=b，其中b为一个向量。(程序与结果)建立矩阵A：>>A=[234;567;891]A的转置B:>>B=A'B=TOC\o"1-5"\h\z5 86 97 1求解：b=[1;2;3];>>x=inv(A)*bx二-0.00000.3333-0.00002通过具体实例构文件说明建MATLAB基本流程结构：for,while,if等while--end循环语句的格式While条件％初值开始，终值结束%循环体中的执行语句%循环结束%循环体中的执行语句%循环结束end程序:>>x=1;>>while(x)<100;x=1*x"3+14;end>>x结果：x二33893列举一些常用函数，并举例说明其应用。如length,ones等MATLAB部分常用向量函数

函数说明函数说明max最大值eye对角线为1的矩阵min最小值ones全1矩阵length长度rand均匀分布的随机矩阵mean均值randn止态分布随机矩阵std标准差diag对角矩阵>>magic(4)ans=16 2 3 135 11 10 89 7 6 124 14 15 14通过实例说明常用画图命令。绘制数据点确定的曲线>>y=rands(10)';>>x=linspace(l,10,10);>>plot(x,y)0.8-\-〔 1 r 1 1 iiii'123456789 10在同一张图中绘制正余弦曲线t=O:O.1:1Oy1=sin(t);y2=cos(t);plot(t,y1,'r',t,y2,'b--');x=[1.7*pi;1.6*pi];y=[-0.3;0.8];s=['sin(t)';'cos(t)'];text(x,ys);title('正弦和余弦曲线');legend('正弦','余弦')xlabel(时间t'),ylabel('正弦、余弦')gridaxissquare

正苍和余荘fijfl'0 2 J S B IDRtfQi实验总结经过这次认真的实验，我了解到了MATLAB基础知识，掌握了MATLAB基本操作，对MATLAB有了进一步的了解及认识，友好的工作平台和环境编程是MATLAB的优点。实验成绩评阅时间评阅教师实验二实验题目描绘出随机过程X(t)=xcos(wt)的图像实验目的利用MATLAB编程描绘出随机过程X(t)=Xcos(wt)的图像实验地点及时间实验内容实验习题给出其程序与结果>>t=l;w=2;>>y=t*cos(w.*x);>>plot(x,y':')

实验总结明确了遇到有些问题可以通过自己定义参数来解决问题。学会了利用MATLAB编程描绘出随机过程数学表达式实验成绩评阅时间评阅教师实验三实验题目绘制随机相位余弦波X(t)=acos(wt+9)的均值，方差和自相关函数的图像实验目的通过绘制图像，深入理解随机相位正弦波的均值，方差和自相关函数实验地点及时间实验内容：绘制随机相位正弦波X(t)=acos(wt+9)的均值，方差和自相关函数的图像实验习题给出其程序与图像余弦波的图像：>t=0:0.01:4*pi;>>w=4;a=5;>>r=pi/4;>>y=a*cos(w*t+r);>>plot(t,y,'*')

自相关的图像：>>symst1t2k>>yl=5*sin((pi/10)*tl+k);>>y2=5*sin((pi/10)*t2+k);>>y=y1*y2;>>r=1/(2*pi)*int(y,k,0,2*pi);>>ezmeshc(r)TicosfTi^(-1.0/1.0e1)+7it20.O/1.Oe1))0.433540284805665e17Z3.602879701896397e16)

方差的图像:>>symstk;>>y=2*cos((pi/5)*t+k);>>h=3/(2*pi)*int(y,k,0,2*pi);>>m=(y-h)A3;>>d=l/(2*pi)*int(m,k,0,2*pi);>>ezplot(d)o0.80.60.40.20-02-04-0..6-08-20246j-6o0.80.60.40.20-02-04-0..6-08-20246j-6实验总结通过绘制图像，深入理解随机相位正弦波的均值，方差和自相关函数评阅时间评阅教师实验成绩实验四实验题目求Markov链的极限分布实验目的用Matlab语言求Markov遍历链的极限分布实验地点及时间实验内容判定一个Markov链是否是遍历的，若是遍历的，求其极限分布。并能从实际问题中抽象出Markov链，并求出其极限分布，并理解其实际意义。实验习题1、已知齐次马氏链｛x(n),n=0,1,2, ｝的状态空间E=｛1,2,3)，状态转移矩阵为(1/21/31/6)P=1/31/31/3、1/31/21/6丿计算2步转移概率;程序:p=[1/21/31/6;1/31/31/3;1/31/21/6];>>P=pA2结果：P=0.41670.36110.22220.38890.38890.22220.38890.36110.2500(2)>>P=[2/52/51/5]P=0.40000.40000.2000>>F=P*EF=0.40000.37220.2278(3)>>P=[1/21/3l/6;l/31/31/3;1/31/21/6]>>M=[100;010;001];>>N=p-MN=-0.50000.33330.16670.3333-0.66670.33330.33330.5000-0.8333>>O=[L;111]-0.50000.33330.16670.3333-0.66670.33330.33330.5000-0.83331.00001.00001.0000>>b=[0001];>>rl二rank(O)r1=3>>r2=rank([O,b'])r2=3>>formatrat;>>x=O\b'x二1/31/31/32、1、为适应日益扩大的旅游事业的需要，某城市的A,B,C三个照相馆组成一个联营部，联合经营出租相机的业务，旅游者可由A,B,C三处任何一处租出相机，用完后还到A,B,C三处的任何一处即可.估计转移概率如表所示，今欲选择A,B,C之一附设租机维修点，问该点设在何处为好？(程序与结果)还相机处ABC租相机处A0.20.80B0.800.2C0.10.30.6解：由于旅客还相机的情况只与该租机点有关，而与相机以前所在的店无关，所以可用Xn表示相机第n次被租时所在的店址;“Xn=1、2、3”分别表示相机第n次被租时在甲，已丙馆。贝0{Xn，n=1,2……}是一个马氏链，其转移矩阵p由上表给出，考虑维修点的设置地点问题，实际上要计算这个马氏链的极限概率分布。解方程组：P1=0.2p1+0.8p2+0.1p3P2=O.8p1+O.3p3P3=0.2p2+0.6p3P1+p2+p3=1程序如下：>>p=[0.20.80;0.800.2;0.10.30.6];>>a=[p'-eye(3);ones(1,3)];>>b=[zeros(3,1);1];>>p_limit=a\