随机过程-实验报告
离散时间随机过程建模实验报告
离散时间随机过程建模实验报告实验报告姓名：实验名称：离散时间随机过程建模学号：课程名称：统计信号处理基础班级：实验室名称：组号：实验日期：2012.10.10一、实验目的、要求本实验的目的是在了解了Matlab编程语言的编程和调试的基础上，利用Matlab本身自带的函数来验证随机信号建模，并掌握子函数的编写方法。
计算机根据理论模型生成随机数，学生需要根据观测的数据编程来计算随机过程的参数。
本实验主要是为了让学生在充分理解不同的随机过程建模的理论方法的基础上，用计算机来认识理论和仿真模型之间的差异。
要求包括以下几个部分：1．要求独立完成实验的内容所要求的各项功能，编制完整的Matlab程序，并在程序中注释说明各段程序的功能。
是一组脉冲序列，101()()()kk x n x k n n δ==-∑ 其中()k x k n 和的取值为 kn 25 40 55 65 85 95 110 130 140 155 ()x k 1 0.8 0.7 0.5 0.7 0.2 0.9 0.5 0.6 0.3 a. 根据上面的关系，画出观测数据()()()y n x n g n =*，并看看是否能通过()y n 的峰值来确定()x n 的幅度和位置。
（需要调用conv 函数）程序 n=[0:50];g=cos((n-25)/5).*exp(-(n-25).*(n-25)/100);g(51:200)=0; x=zeros(200,1);x(25)= 1; x(40)=0.8 ; x(55)=0.7 ; x(65)=0.5 ;x(85)= 0.7; x(95)=0.2; x(110)=0.9; x(130)=0.5;x(140)=0.6; x(155)=0.2;y=conv(x,g);figure(1)subplot(3,1,1),plot(g);title('滤波器冲击响应');xlabel('n');ylabel('响应幅值');subplot(3,1,2),plot(x);title('输入序列x');xlabel('n');ylabel('幅值');subplot(3,1,3),plot(y);title('滤波器输出');xlabel('n');ylabel('幅值');形b. 用教材中给出的spike.m函数来设计长度50N 的最小二乘逆滤波器()Nh n，并确定最佳的延迟。
随机实验报告1Poisson过程模拟
数学与计算科学学院实验报告实验项目名称随机数及Poisson过程的模拟所属课程名称随机过程实验类型综合实验日期班级学号姓名成绩一、实验概述： 【实验目的】通过模拟产生随机数，进一步编程实现对possion 过程样本轨道的模拟。
掌握生成随机变量的方法，深入了解poisson 过程的性质。
【实验原理】1、随机变量的生成（逆函数法）：利用均匀分布并结合分布函数的逆变换，生成分布函数为F （x ）的变换：若U 是[0,1]区间上的均匀分布，F （x ）为任一给定的分布函数，定义1()inf{:()}F x t F t x -=>，则随机变量1()Y F U -=的分布函数为F （x ）;2、Poisson 过程的模拟：（1）利用事件发生的间隔时间是独立同分布的随机变量序列，（2）给定事件发生次数的条件下，事件发生的时刻与该区间上对应的均匀分布的顺序统计量相同【实验环境】 硬件环境Windows 7 Microsoft Corporation Inter(R)Core(TM) i5-3210 软件环境 Matlab 7.0 二、实验内容： 【实验方案】1、利用求逆函数的方法生成指数分布随机变量；2、（a ）利用独立同分布的指数分布序列模拟强度为1的Poisson 过程； （b ）利用均匀分布的顺序统计量模拟强度为1的Poisson 过程 【实验过程】（实验步骤、记录、数据、分析） 1.利用求逆函数的方法生成指数分布随机变量；步骤一：我们知道一个指数分布的概率密度函数是：其中λ > 0是分布的一个参数，常被称为率参数（rate parameter ）。
即每单位时间发生该事件的次数。
指数分布的区间是[0,∞)。
如果一个随机变量X 呈指数分布，则可以写作：X ~ Exponential （λ）。
累积分布函数:累积分布函数可以写成：所以在 0≥x 时该分布函数的逆变换为：步骤二：生成均匀分布在[0,1]上的随机数Matlab 里生成[0，1]上的均匀随机数的语句是：rand(1,1); rand(n,m)。
随机过程实验报告
随机过程实验报告一.实验目的通过随机过程的模拟实验, 熟悉随机过程编码规律以及各种随机过程的实现方法, 通过理论与实际相结合的方式, 加深对随机过程的理解。
二. 实验原理及实现代码1.伪随机数的产生函数功能: 采用线性同余法, 根据输入的种子数产生一个伪随机数, 如果种子不变, 则将可以重复调用产生一个伪随机序列实现思路：利用CMyRand类中定义的全局变量：S, K, N, Y。
其中K和N为算法参数, S用于保存种子数, Y为产生的随机数, 第一次调用检查将seed赋值与S获得Y的初值, 之后调用选择rand（）函数赋值与Y。
代码如下:unsigned int CMyRand::MyRand(unsigned int seed){Y=seed;Y=K*seed%N;S=Y;return Y;}2.均匀分布随机数的产生在上面实验中, 已经产生了伪随机序列, 所以为了得到0~N 的均匀分布序列, 只需将其转化为min 到max 的均匀分布即可, 代码如下:double CMyRand::AverageRandom(double min,double max) {double dResult;dResult = (double(MyRand(S))/N)*(max-min)+min; dResult=(int(dResult*10000))/10000.0 ;return dResult; }3.正态分布随机数的产生由AverageRandom 函数获得0－1间隔均匀分布随机数U(0,1), i=1,2,…,n, 且相互独立, 由中心极限定理可知, 当n 较大时,()~(0,1)nU nE U Z N -=取n=12, 近似有, 也就是说, 只要产生12个伪随机数u1,u2,…u12, 将它们加起来, 再减去6, 就能近似得到标准正态变量的样本值。
代码如下:double CMyRand::NormalRandom(double miu, double sigma, double min, double max){double dResult;dResult = 0;for(int i=0;i<12;i++)dResult+=(double(MyRand(S))/N); //循环相加12次dResult-=6;dResult=(dResult*sigma+miu)*(max-min)+min;return dResult;}3.指数分布的随机数的产生用AverageRandom产生均匀分布随机数{ui}, 计算指数分布随机数: xi=-ln ui /λdouble CMyRand::ExpRandom(double lambda, double min, double max){double dResult = 0.0;dResult=-log(AverageRandom(min,max))/lambda;return dResult;}4.泊松分布的随机数产生unsigned int CMyRand::PoisonRandom(double lambda, double min, double max){unsigned int dResult = 0;double F=exp(-lambda);while(AverageRandom(0,1)>=F){F+=(lambda*F)/(dResult+1);dResult++;}return dResult;}5.计算任意分布的随机过程的均值根据大数定律, 调用任意函数加和求平均即为该分布的均值。
随机过程实验报告
一、实验目的1. 理解随机过程的基本概念和性质。
2. 掌握随机过程的基本运算和性质。
3. 通过实验验证随机过程的性质和规律。
二、实验原理随机过程是指一系列随机变量按照一定规则排列而成的序列。
在现实生活中，随机过程广泛存在于自然界和人类社会，如股票价格、气象变化、生物进化等。
随机过程的研究有助于我们更好地理解和预测这些现象。
随机过程可以分为两类：离散随机过程和连续随机过程。
本实验主要研究离散随机过程。
三、实验设备与材料1. 计算机2. 随机过程模拟软件（如Matlab）3. 纸笔四、实验内容1. 随机过程的基本概念（1）随机变量的概念随机变量是指具有不确定性的变量，它可以取多个值。
在随机过程中，随机变量是基本的研究对象。
（2）随机过程的概念随机过程是由一系列随机变量按照一定规则排列而成的序列。
2. 随机过程的基本性质（1）无后效性无后效性是指随机过程的前后状态相互独立。
（2）无记忆性无记忆性是指随机过程的状态只与当前时刻有关，与过去时刻无关。
（3）马尔可夫性马尔可夫性是指随机过程的状态只与当前时刻有关，与过去时刻无关。
3. 随机过程的运算（1）随机过程的和设{Xn}和{Yn}是两个随机过程，则它们的和{Zn}定义为Zn = Xn + Yn。
（2）随机过程的差设{Xn}和{Yn}是两个随机过程，则它们的差{Zn}定义为Zn = Xn - Yn。
（3）随机过程的乘积设{Xn}和{Yn}是两个随机过程，则它们的乘积{Zn}定义为Zn = Xn Yn。
4. 随机过程的模拟利用随机过程模拟软件（如Matlab）模拟随机过程，观察其性质和规律。
五、实验步骤1. 初始化随机数生成器2. 定义随机过程（1）根据随机过程的基本性质，定义随机过程{Xn}。
（2）根据随机过程的运算，定义随机过程{Yn}。
3. 模拟随机过程（1）使用随机过程模拟软件（如Matlab）模拟随机过程{Xn}和{Yn}。
（2）观察模拟结果，分析随机过程的性质和规律。
随机实验报告
随机信号实验报告课程：随机信号实验题目：随机过程的模拟与特征估计学院：学生名称：实验目的：1.学会利用MATLAB模拟产生各类随即序列。
2.熟悉和掌握随机信号数字特征估计的基本方法。
实验内容：1.模拟产生各种随即序列，并画出信号和波形。
（1）白噪声（高斯分布，正弦分布）。
（2）随相正弦波。
（3）白噪声中的多个正弦分布。
（4）二元随机信号。
（5）自然信号：语音，图形（选做）。
2.随机信号数字特征的估计（1）估计上诉随机信号的均值，方差，自相关函数，功率谱密度，概率密度。
（2）各估计量性能分析（选做）实验仪器：PC机一台MATLAB软件实验原理：随机变量常用到的数字特征是数字期望值、方差、自相关函数等。
相应地，随机过程常用到的数字特征是数字期望值、方差、相关函数等。
它们是由随机变量的数字特征推广而来，但是一般不再是确定的数值，而是确定的时间函数。
1.均值：m x(t)=E[X(t)]=；式中，p(x,t)是X（t）的一维概率密度。
m x(t)是随机过程X（t）的所有样本函数在时刻t的函数值的均值。
在matlab中用mea()函数求均值。
2.方差：（t）=D[X(t)]=E[]；（t）是t的确定函数，它描述了随机过程诸样本函数围绕数学期望m x(t)的分散程度。
若X（t）表示噪声电压，则方差（t）则表示瞬时交流功率的统计平均值。
在matlab中用var()函数求均值。
3.自相关函数：Rx(t1,t2)=E[X(t1)X(t2)]；自相关函数就是用来描述随机过程任意两个不同时刻状态之间相关性的重要数字特征。
在matlab中用xcorr（）来求自相关函数。
4.在matlab中可用函数rand、randn、normr、random即可生成满足各种需要的近似的独立随机序列。
实验步骤：（一）大体实验步骤（1）利用MATLAB编写程序。
（2）调试程序。
（3）得出各项输出结果，产生波形。
（4）分析各参数的物理意义，各个波形参数相比较。
随机过程上机实验报告讲解
2015-2016第一学期随机过程第二次上机实验报告实验目的：通过随机过程上机实验，熟悉Monte Carlo计算机随机模拟方法，熟悉Matlab的运行环境，了解随机模拟的原理，熟悉随机过程的编码规律即各种随机过程的实现方法，加深对随机过程的理解。
上机内容：(1 )模拟随机游走。
(2)模拟Brown运动的样本轨道。
(3)模拟Markov过程。
实验步骤：(1)给出随机游走的样本轨道模拟结果，并附带模拟程序。
①一维情形%—维简单随机游走% “从0开始，向前跳一步的概率为p，向后跳一步的概率为1-p”n=50;p=0.5;y=[0 cumsum(2.*(rand(1,n-1)v=p)-1)]; % n 步。
plot([0:n-1],y); %画出折线图如下。
w%一维随机步长的随机游动%选取任一零均值的分布为步长，比如，均匀分布。
n=50;x=rand(1,n)-1/2;y=[0 (cumsum(x)-l)];plot([0:n],y);②二维情形%在(u, v)坐标平面上画出点(u(k), v(k)), k=1:n,其中(u(k)) 和(v(k))是一维随机游动。
例%子程序是用四种不同颜色画了同一随机游动的四条轨道。
n=100000;colorstr=['b' 'r' 'g' 'y'];for k=1:4z=2.*(rand(2,n)<0.5)-1;x=[zeros(1,2); cumsum(z')];col=colorstr(k);plot(x(:,1),x(:,2),col);③%三维随机游走 ranwalk3dp=0.5;n=10000; colorstr=['b' 'r' 'g' 'y'];for k=1:4z=2.*(rand(3,n)v=p)-1; x=[zeros(1,3); cumsum(z')];col=colorstr(k);plot3(x(:,1),x(:,2),x(:,3),col);hold on end gridhold onendgrid4：04003?0-200-300-400-2OD20050、-100-200 -20D⑵给出一维，二维Brown运动和Poisson过程的模拟结果，并附带模拟程序，没有结果的也要把程序记录下来。
随机过程实验报告全
随机过程实验报告学院：专业：学号：姓名：一、实验目的通过随机过程的模拟实验，熟悉随机过程编码规律以及各种随机过程的实现方法，通过理论与实际相结合的方式，加深对随机过程的理解。
二、实验内容（1）熟悉Matlab工作环境，会计算Markov链的n步转移概率矩阵和Markov链的平稳分布。
（2）用Matlab产生服从各种常用分布的随机数，会调用matlab自带的一些常用分布的分布律或概率密度。
（3）模拟随机游走。
（4）模拟Brown运动的样本轨道的模拟。
（5）Markov过程的模拟。
三、实验原理及实验程序n步转移概率矩阵根据Matlab的矩阵运算原理编程，Pn = P ^n。
已知随机游动的转移概率矩阵为：P =0.5000 0.5000 00 0.5000 0.50000.5000 0 0.5000求三步转移概率矩阵p3及当初始分布为P{x0 = 1} = p{x0 = 2} = 0, P{x0 = 3} = 1 时经三步转移后处于状态3的概率。
代码及结果如下：P = [0.5 0.5 0; 0 0.5 0.5; 0.5 0 0.5] %一步转移概率矩阵P3 = P ^3 %三步转移概率矩阵P3_3 = P3(3,3) %三步转移后处于状态的概率1、两点分布x=0:1;y=binopdf(x,1,0.55);plot(x,y,'r*');title('两点分布');2、二项分布N=1000;p=0.3;k=0:N;pdf=binopdf(k,N,p);plot(k,pdf,'b*');title('二项分布');xlabel('k');ylabel('pdf');gridon;boxon3、泊松分布x=0:100;y=poisspdf(x,50);plot(x,y,'g.');title('泊松分布')4、几何分布x=0:100;y=geopdf(x,0.2);plot(x,y,'r*');title('几何分布');xlabel('x');ylabel('y');5、泊松过程仿真5.1 % simulate 10 timesclear;m=10; lamda=1; x=[];for i=1:ms=exprnd(lamda,'seed',1);x=[x,exprnd(lamda)];t1=cumsum(x);end[x',t1']5.2%输入：N=[];for t=0:0.1:(t1(m)+1)if t<t1(1)N=[N,0];elseif t<t1(2)N=[N,1];elseif t<t1(3)N=[N,2];elseif t<t1(4)N=[N,3];elseif t<t1(5)N=[N,4];elseif t<t1(6)N=[N,5];elseif t<t1(7)N=[N,6];elseif t<t1(8)N=[N,7];elseif t<t1(9)N=[N,8];elseif t<t1(10)N=[N,9];elseN=[N,10];endendplot(0:0.1:(t1(m)+1),N,'r-') 5.3% simulate 100 timesclear;m=100; lamda=1; x=[];for i=1:ms= rand('seed');x=[x,exprnd(lamda)];t1=cumsum(x);end[x',t1']N=[];for t=0:0.1:(t1(m)+1)if t<t1(1)N=[N,0];endfor i=1:(m-1)if t>=t1(i) & t<t1(i+1)N=[N,i];endendif t>t1(m)N=[N,m];endendplot(0:0.1:(t1(m)+1),N,'r-')6、泊松过程function I=possion(lambda,m,n)for j=1:mX=poissrnd(lambda,[1,n]); %参数为lambda的possion 过程N(1)=0;for i=2:nN(i)=N(i-1)+X(i-1);endt=1:n;plot(t,N)grid onhold onend7、布朗运动7.1一维布朗运动程序：function [t,w]=br1(t0,tf,h)t=t0:h:tf;t=t';x=randn(size(t));w(1)=0;for k=1:length(t)-1w(k+1)=w(k)+x(k);endw=sqrt(h)*w;w=w(:);end调用t0=1;tf=10;h=0.01;[t,w]=br1(t0,tf,h);figure;plot(t,w,'*');xlabel('t');ylabel('w');title('一维Brown运动模拟图'); 7.2二维布朗运动：function [x,y,m,n]=br2(x0,xf,y0,yf,h)x=x0:h:xf;y=y0:h:yf;a=randn(size(x));b=randn(size(y));m(1)=0;n(1)=0;for k=1:length(x)-1m(k+1)=m(k)+a(k);n(k+1)=n(k)+b(k);endm=sqrt(h)*m;n=sqrt(h)*n;end调用x0=0;xf=10;h=0.01;y0=0;yf=10;[x,y,m,n]=br2(x0,xf,y0,yf,h);figure;plot(m,n);xlabel('m');ylabel('n');title('二维Brown运动模拟图');7.3三维布朗运动：npoints =1000;dt = 1;bm = cumsum([zeros(1, 3); dt^0.5*randn(npoints-1, 3)]);figure(1);plot3(bm(:, 1), bm(:, 2), bm(:, 3), 'k');pcol = (bm-repmat(min(bm), npoints, 1))./ ...repmat(max(bm)-min(bm), npoints, 1);hold on;scatter3(bm(:, 1), bm(:, 2), bm(:, 3), ...10, pcol, 'filled');grid on;hold off;8、马尔科夫链离散服务系统中的缓冲动力学m=200;p=0.2;N=zeros(1,m); %初始化缓冲区A=geornd(1-p,1,m); %生成到达序列模型, for n=2:mN(n)=N(n-1)+A(n)-(N(n-1)+A(n)>=1);endstairs((0:m-1),N);9、随机数游走9.1 100步随机游走n = 100; %选取步数。
随机过程实验
实验名称：随机变量的仿真与实验实验内容：用MATLAB 分别产生服从（二项分布、泊松分布、正态分布、均匀分布、指数分布、瑞利分布）的随机变量，并分析他们的：1、分布函数或概率密度函数2、均值、方差1、服从二项分布的随机变量理论分析如果随机变量X 的分布律为k n k k n k q p C k X P p -===}{0<p<1, q=1-p, k=0,1,2，…n,则称X 服从参数为n ，p 的二项分布，记为X~B(n ，p)。
其期望和方差分别为E(X) = np ，D(X)=npq 。
随机变量X~B(20，0.4)，可以通过matla b 计算其期望和方差，绘制分布律和分布函数。
程序如下：n = 20;p = 0.4;[E,D] = binostat(n ,p); %计算期望和方差f = binopdf(1:21, n, p); %计算分布律F = binocdf(1:21, n, p); %计算分布函数subplot(2,2,1); stem(f); %绘制分布律title('二项分布理论分布律 n=20 p=0.4');xlabel('x');ylabel('p');subplot(2,2,3); stem(F); %绘制分布函数title('二项分布理论分布函数 n=20 p=0.4');xlabel('x');ylabel('f');计算得结果E(X) = 8，D(X) = 4.800，分布律和分布函数如图1。
图1 X~B(20，0.4)的分布律和分布函数样本分析利用matlab中binornd函数产生一个X~B(20，0.4)的样本，样本点总数为20000。
计算其均值和方差，计算分布律和分布函数，并与理论结果进行比较。
程序如下:n = 20;p = 0.4;R = binornd(n,p,1,20000);e = mean(R); %期望d = var(R); %方差f = zeros (1,21);F = zeros (1,21);for j = 1:21 %计算统计分布律for i=1:20000if j == R(i)f(1,j) = f(1,j) + 1;endendf(1,j) = f(1,j) / 20000;endsubplot(2,2,1);stem(f);title('二项分布样本分布律 n=20 p=0.4');xlabel('x');ylabel('p');for j = 1:21 %计算分布函数for i = 1:jF(1, j) = F(1, j) + f(1,i);endendsubplot(2,2,3);stem(F);title('二项分布样本分布函数 n=20 p=0.4');xlabel('x');ylabel('f');计算结果为e=8.0218，d=4.7760，与理论值(E(X)=8，D(X)=4.8)基本接近。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
结果
5
通过此次实验，模拟 Possion 流分布，运用 MATLAB 产生随机数，使我对泊松 分布有了更深刻的理解。 不过产生的随机数可以根据实际需要设置不同分布的随 机数，便于解决实际问题。
实验成绩
评阅时间
评阅教师
实验四
6
实验题目 实验目的 实验地点及时间
求 Markov 链的极限分布 用 Matlab 语言求 Markov 遍历链的极限分布 信息楼 127 机房 2012 年 6 月 6 日
1 / 2 P 1/3 1/ 3 1/3 1/3 1/ 2 1/ 6 1/3 1/ 6
(1) 计算 2 步转移概率;(2) 已知初始分布为 P 2 / 5, 2 / 5,1 / 5 ,求 X 2 的分布律 (3) 求平稳分布，要求给出程序与结果。 程序：
实验三
实验题目 实验目的 实验地点及时间 模拟 Possion 流 用 Matlab 语言产生随机数，了解 Possion 流 信息楼 127 机房 2012 年 6 月 4 日
4
实验内容 用 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<S(2) d(i+1)=0; else for j=2:21 if (S(j)<=i)&(S(j+1)<i) d(i+1)=j; end end end end plot(d)
P2 = 0.4167 0.3889 0.3889 0.3611 0.3889 0.3611 0.2222 0.2222 0.2500
7
(1) 2 步转移概率 P2 = 0.4167 0.3889 0.3889 0.3611 0.3889 0.3611 0.2222 0.2222 0.2500
(2) X 2 的分布律 S2 = 0.1667 (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
随机过程试验报告
班级：信计 09 级 01 班
姓名：
学号：
实验一
实验题目 描绘出随机过程 X ( t ) x cos( w t ) 的图像 实验目的 利用 MATLAB 编程描绘出随机过程 X ( t ) x cos( w t ) 的图像 实验地点及时间 信息楼 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
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;
3
实验总结 通过绘制 X ( t ) cos( t ) 图像， 由图像可知， 随机相位正弦波的均值 E （t） =0；方差 Var(t)=2；由自相关函数图像的描述出不同时刻之间的相关程度。通过 实验，我们可以知道，均值、方差、自相关函数是刻画随机过程的主要特性，因 此，对于解决实际课题而言，常常能够起到重要作用。 实验成绩 评阅时间 评阅教师
程序： 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 结果：
实验内容 判定一个 Markov 链是否是遍历的，若是遍历的，求其极限分布。并能从实际问 题中抽象出 Markov 链，并求出其极限分布，并理解其实际意义。 实验习题 1、已知齐次马氏链 X n , n 0,1, 2, 的状态空间 E 1, 2, 3 ,状态转移矩阵为
>> 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)
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
实验二
1
绘制随机相位正弦波 X ( t ) cos( w t ) 的均值，
实验题目
方差和自相关函数的图像
通过绘制图像，深入理解随机相位正弦波的均值，方差和自相 关函数 实验地点及时间 信息楼 127 机房 2012 年 6 月 1 日 实验目的 实验内容：绘制随机相位正弦波 X ( t ) cos( w t ) 的均值，方差和自相 关函数的图像 实验习题 给出其程序与图像