MATLAB的答案

如何模拟泊松分布过程?

97(30天)
萨菲亚”class=
萨菲亚 2022年9月5日
评论道: Farshid R2022年10月1日
大家好!
我有许多的车辆经过一段公路,到达率遵循泊松分布与λ= 3(汽车/分钟)
T = 60分钟,dt = 1分钟。
我怎么能模拟这个过程呢?

在回答这个问题。

接受的答案

明星黾”class=
明星黾 2022年9月5日
看看 poissrnd 函数将做你想做的事情
λ= 3;
v = poissrnd(λ1 60)
v = 1×60
4 3 1 2 2 4 5 4 3 5 6 5 6 1 4 1 5 1 2 0 6 2 1 1 3 4 4 3 6
。
14日的评论
显示13年长的评论
明星黾”class=
明星黾 2022年9月29日
我不确定我理解。
泊松过程决定了汽车到达时,并不会有多少。的要求 3 车辆/分钟意味着汽车的数量将是相同的(或几乎相同的)任何特定的测量时间间隔,提供时间间隔是一样的长度,长度是,例如,几分钟的时间。

登录置评。

更多的答案(5)

威廉•罗斯”class=
威廉•罗斯 2022年9月5日
@Safia
是的。你可以使用二项式(还是伯努利?)近似:每分钟分解成N =段,λ/ N < < 1。在每个segmeent, 1辆车通过的几率是p =λ/ N。两辆车通过的机会如果N是足够大可以忽略不计。这就是为什么这是一个近似。所以你“抛硬币”通过借鉴rand()在每一段。汽车计数加1,每次rand()小于p。
然后重新启动你的车与零,再做下一分钟。
以上只是一个大纲;我没有时间检查或providee代码细节,但我认为这将给一个很好的近似。
试试看,祝你好运。
图像分析”class=
图像分析 2022年9月5日
你想模拟一个吗 随机过程 吗?像有1000车道的高速公路的汽车旅行。这些是矩阵的行。和列《纽约时报》,从1到60(每一列是一分钟)。你填充第一列。然后在下一次迭代(一分钟)移动列1 /第2列在第一列(因为汽车驾驶),现在做一个新列1。然后在下一次迭代你改变所有列一列和几辆车装满第一列。
4评论
显示3年长的评论

登录置评。

布鲁诺陈德良”class=
布鲁诺陈德良 2022年9月5日
编辑:布鲁诺陈德良 2022年9月5日
如果你没有staistics工具箱
%生成n遵循泊松随机整数参数λ(= =)
λ= 3;
n = 1 e6;% 60在你的情况中
%的最大可能值
k = 0;
托尔=每股收益(exp(λ));
p = 1;
而p >托尔
k = k + 1;
p = p *λ/ k;
结束
kmax = k;
%累计金额
K = 0: kmax;
c = [0 cumsum(λ。^ K /阶乘(K))) * exp(λ);
c = c / c(结束);
r =离散化(兰德(1,n), c) 1;
意思是(右)% ~λ
ans = 2.9968
性病(r) ^ 2% ~λ
ans = 2.9934
柱状图(右)
5个评论
显示4年长的评论
布鲁诺陈德良”class=
布鲁诺陈德良 2022年9月6日
泰勒指数计算的部分和不完整的γ函数:
%生成n遵循泊松随机整数参数λ(= =)
λ= 3;
n = 1 e6;% 60在你的情况中
%的最大可能值
k = 0;
托尔=每股收益(exp(λ));
p = 1;
而p >托尔
k = k + 1;
p = p *λ/ k;
结束
kmax = k;
%累计金额
c = gammainc(λ,0:kmax + 1,“上”);
c = c / c(结束);%确定最后一本是1。
r =离散化(兰德(1,n), c) 1;
意思是(右)% ~λ
ans = 3.0034
性病(r) ^ 2% ~λ
ans = 2.9986
柱状图(右)

登录置评。

布鲁诺陈德良”class=
布鲁诺陈德良 2022年9月29日
@Safia 与固定的金额是你想要的整数,我提出这个问题,它不再是泊松法律严格来说吗
λ= 3;
T = 60;%周期长度,在几分钟内
n = 1 e3;%的组合数的T
% r n x T,随机求和(r, 2)
%的“频率”r s / T ~λ
s =圆(λ* T);
r =兰迪([0],n, t - 1);
r = diff ([0 (n, 1),排序(r, 2), s + 0 (n, 1)), 1, 2);
% r并不遵循泊松,但不是遥远的显示
直方图(r (:))
意思是(r (:))
ans = 3
性病(r (:))
ans = 3.0077
布鲁诺陈德良”class=
布鲁诺陈德良 2022年9月6日
编辑:布鲁诺陈德良 2022年9月6日
另一种方法基于exponetial分布,这两个事件之间的时间
%生成n遵循泊松随机整数参数λ(= =)
λ= 3;
n = 1 e6;
r = [];
tlast = 0;
而tlast < n + 1
%,我们通常只需要一个迭代
m =装天花板((n + 1-tlast) *λ* 1.1);% 10%的保证金
s = cumsum(日志(1-rand (m, 1)) /λ);
r = [r;tlast + s];% #好吧
tlast = r(结束);
结束
r = accumarray ((r),即:1);
r (n + 1:结束)= [];%将超过尾巴
%检查
意思是(右)
ans = 3.0006
性病(r) ^ 2
ans = 3.0012
柱状图(右)
1评论
Farshid R”class=
Farshid R 2022年10月1日
你好 @Bruno陈德良 ,
我有一个优化问题(共识优化问题),不幸的是,我在MATLAB站点上发布的问题,
@Sam翟告诉我,我可以发送一个消息给你,我可以和你讨论我的问题。
我的链接问题:
//www.tianjin-qmedu.com/matlabcentral/answers/1812615-optimization-with-fmincon-command-in-万博1manbetxsimulink?s_tid=mlc_ans_men_view&mentions=true comment_2390935
的问题,我给模型和关系完全和我讨论了@Sam翟但没有达到期望的结果和@万博1manbetxSam翟说给你发送一条消息。
请指导我。我期待着你的肯定的答复。
最好的问候,

登录置评。

类别

找到更多的在泊松分布在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!