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-2i%print running processs1=0;for j=0:i-1s1=s1+(-1)^j/factorial(j);%i=1 or nend s2=0;for j=0:n-is2=s2+(-1)^j/factorial(j);%i=1 or nend Psum=Psum+s1*s2; end Psum+2*exp(-1)

matlab運算結果：

?

而計算器結果(n+2)/(e^2)=270.941237

書上結果應該是正確的。

總結

以上是生活随笔為你收集整理的matlab验证Ross随机过程(第二版)P19页的结果的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。