%泊松過程的模擬、檢驗及參數(shù)估計
syms Un X S;
n=10;%生成n*n個隨機數(shù)
r=1;%參數(shù)
temp=0;
tem=0;
Un=rand(n,1);%共產(chǎn)生n*n個隨機數(shù)
for i=1:1:n
? ? X(i)=-log(Un(i))/r;
end
X=subs(X);
for i=1:1:n
? ? for j=1:1:i
? ? ? ? temp=temp+X(j); ? ?
? ? end
? ? S(i)=temp;
? ? temp=0;
end
S=subs(S);
%檢驗泊松過程使用第四條
for i=1:1:n
? ? tem=tem+S(i);
end
sigmaN=tem;
T=S(n);
alpha=0.05;%置信水平
p=sigmaN/T;
p1=(1/2)*(n-1.96*(n/3)^(1/2));
p2=(1/2)*(n+1.96*(n/3)^(1/2));
c1=subs(p-p1)
c2=subs(p-p2)
if (c1<=0&c2>=0)|(c1>=0&c2<=0)
? ? disp('這是一個泊松過程！')
? ? %參數(shù)估計使用極大似然估計
? ? r_=subs(n/T);
? ?if abs(subs(r_-r))<0.1
? ? ? disp('參數(shù)估計正確！')
? ? ? disp('參數(shù)估計值為：')
? ? ? r_
? ?end ??
? ?%繪制軌跡
? ?y=0;
? ?x=0:0.001:subs(S(1));
? ?plot(x,y)
? ?for k=1:1:n
? ? ? y=k;
? ? ? x=subs(S(k)):0.001:subs(S(k+1));
? ? ? hold on
? ? ? plot(x,y)
? ?end
else
? ? disp('這不是一個泊松過程！')
end
D116

總結(jié)

以上是生活随笔為你收集整理的使用matlab模拟、检验和估计泊松过程的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。