随机信号马尔可夫仿真.docVIP

随即信号分析仿真 专业：电子信息工程 实验一、离散时间马尔科夫链仿真 1.源程序：见附录shiyanyi.m 文件 2.实验步骤： (a)设置初始状态X0，令时刻n=1，当前状态i=1； (b)生成u，利用条件分布{Pi1，Pi2,,,,PiI}，计算Yi并令Xi=Yi； (c)如果nN，令i=Xn，生成u并利用条件分布{Pi1，Pi2,,,,PiI} 计算Yi，令n=n+1并令Xn=Yi；如果n=N则停止。（调用while循环两次） (d)回到第三步继续循环。 最后算法运行得到的Xn与n之间的关系即为马尔可夫链的一个样本函数。 3.程序运行截图： 4.实验一分析总结：仿真图得到与实际理论图一样，满足离散时间马尔科夫链。总结：这个实验用到了两次循环，内循环是找出符合要求的j值，外重循环是判断n是否小于N，然后重复。 实验二、泊松过程仿真 1.源程序：见附录 shiyaner.m 文件 2.实验步骤： (a)令当前时刻t=0，泊松事件计数值N=0。 (b) 生成U，令 。 (c) 令t=t+E，N=N+1并且设tN=t，如果tT，则停止。 （此步处用一个while循环实现，判断条件是tT） (d)回到第2步继续循环。 最后，算法运行得到的N与tN之间的关系即为泊松过程的一个样本函数。 3.运行截图： 4.实验二分析总结：程序每产生一个指数分布随机数E，行坐标时间t就增加E，纵坐标N也相应的增加1.得到如上图的泊松过程的一个仿真图，实验现象与实际相符合。程序中也是将时间t和N各存在一个数组中，然后调用画图阶梯函数stairs()。 实验三、关于布朗运动的仿真 1.源程序：见附录 shiyansan.m 文件 2.实验步骤： (a)设初始状态W=0，t0=0，n=1。 (b) 如果nN，生成 如果n=N，停止。 （用一个while循环实现，判断条件是nN） (c) 回到第二步继续循环。 3.运行截图： 4.实验三分析总结：得到的实验现象与理论分析的一致。具体，程序在每次循环开始就调用normrnd产生一个正态分布随机数，根据所给公式产生一个增量W作为纵坐标的数值，横坐标为以一定间隔dt增加的时间t。其中程序中也是将时间t和W各存在一个数组中，然后调用画图函数plot()。 实验心得：通过本次随机信号分析课后实验，不仅自己在课后主动的翻阅复习了离散时间马尔科夫链、泊松过程、布朗运动的相关内容，通过实验现象加深对它们在宏观上的理解，而且通过应用matlab软件编程，也提高了自己对matlab编程的能力，对自己的很多方面的都有很大的提高。 附件： %实验二 关于泊松过程仿真 %************************************** %样本函数一 源程序 T=10; %设置步长T为10 t=0; %初始化变量时间t为0 N=0; %计数器N初始化为0 n=0; %计数器初始化为0 x=2; %设置参数 while tT %循环判断条件 u=rand(); %产生一个随机数u E=-log(u)/x; %获得一个指数分布随机数 n=n+1; %计数器加1 t=t+E; %t加e，作为今后的判断条件 t1(n)=t; %将时间t存在一个数组中 N=N+1;%计数器N加1 N1(n)=N; %将N存在一个数组中 end %************************************** %样本函数二 源程序 T=10; %设置步长T为10 t=0; %初始化变量时间t为0 N=0; %计数器N初始化为0 n=0; %计数器初始化为0 x=2; %设置参数 while tT %循环判断条件 u=rand(); %产生一个随机数u E=-log(u)/x; %获得一个指数分布随机数 n=n+1; %计数器加1 t=t+E; %t加e，作为今后的判断条件 t2(n)=t; %将时间t存在一个数组中 N=N+1;%计数器N加1 N2(n)=N; %将N存在一个数组中 end %************************************** %样本函数三 源程序 T=10; %设置步长T为10 t=0; %初始化变量时间t为0 N=0; %计数器N初始化为0 n=0; %计数器初始化为0 x=2; %设置参数 while tT %循环判断条件 u=rand(); %产生一个随机数u E=-log(u)/x; %获得一个指数分布随机数 n=n+1; %计数器加1 t=t+E; %t加e，作为今后的判断条件 t3(n)=t; %将时间t存在一个数组中