matlab验证Ross随机过程(第二版)P19页的结果

Ross书上的式子为:


下面用Matlab验证最终结果ΣPi,n

取n=2000进行验证,运行时间大概是两三分钟。

代码主要思路是:

先处理2<i<n-1的情况,

代码中分别用一个for循环计算Σ,然后在[3,n-2]范围内求和,最后再加上i=1,n时的值1/e+1/e=2/e

s=0;
n=5;
for j=0:n-1
s=s+(-1)^j/factorial(j);%i=1 or n
end 
P4=s%The result of P4 is right.

% ----------------the above is to verify when i=2,n-1-----------



Psum=0
n=2000
for i=3:n-2
    i%print running process
    s1=0;
    for j=0:i-1
    s1=s1+(-1)^j/factorial(j);%i=1 or n
    end 

    s2=0;
    for j=0:n-i
    s2=s2+(-1)^j/factorial(j);%i=1 or n
    end 

Psum=Psum+s1*s2;
end
Psum+2*exp(-1)

matlab运算结果:


而计算器结果(n+2)/(e^2)=270.941237

书上结果应该是正确的。

阅读更多
版权声明:本文为博主原创文章,可以随便转载 https://blog.csdn.net/appleyuchi/article/details/80457895
个人分类: 概率论与数理统计
上一篇virtualpc设置共享文件夹
下一篇泊松分布的来源—公式推导—应用
想对作者说点什么? 我来说一句

没有更多推荐了,返回首页

关闭
关闭
关闭