如何解决随机选择代码中的不均匀数量错误
平台:Matematica
代码目标:在 0.1-1.0 的范围内更改传输概率 (TransmissionProb),并随着每次增量更改,找到将导致完全绿色动态图的相应恢复概率 (RecoveryProb)。
问题:将恢复概率增加到 1 以上会导致代码的 RandomChoice 部分崩溃。我不知道如何解决这个问题。 ):
CellN = 201
TimeN = 100
TransmissionProb = 0.2
RecoveryProb = 2
SpontaneousProb = 0
CurrentWorld =
Table[RandomInteger[{0,2}],{CellN},{CellN}]; Worldplot =
ArrayPlot[CurrentWorld,ColorFunction -> (Blend[{{0,Black},{1,Green},{2,Red}},#] &),ColorFunctionScaling -> False,Mesh -> False];
Dynamic[Worldplot]
NeighborStates = Table[0,{8}]
Tuples[{-1,1},2]
Select[{3,6,7,2},# > 3 &]
MoveBy = Select[Tuples[{-1,2],# != {0,0} &]
For[t = 1,t <= TimeN,t++,For[X = 1,X <= CellN,X++,For[Y = 1,Y <= CellN,Y++,For[n = 1,n <= 8,n++,{NeighborX,NeighborY} = {X,Y} + MoveBy[[n]];
If[NeighborX < 1,NeighborX = CellN];
If[NeighborX > CellN,NeighborX = 1];
If[NeighborY < 1,NeighborY = CellN];
If[NeighborY > CellN,NeighborY = 1];
NeighborStates[[n]] =
CurrentWorld[[NeighborX,NeighborY]];];(*close the For[n,...]*)
If[(CurrentWorld[[X,Y]] == 1) && (Count[NeighborStates,2] > 0),NextWorld[[X,Y]] =
RandomChoice[{1 - TransmissionProb,TransmissionProb} -> {1,2}]];
If[(CurrentWorld[[X,Y]] == 0) && (Count[NeighborStates,2] == 0),Y]] =
RandomChoice[{1 - RecoveryProb,RecoveryProb} -> {0,1}]];
If[(CurrentWorld[[X,Y]] =
RandomChoice[{1 - SpontaneousProb,SpontaneousProb} -> {1,2}]];
If[CurrentWorld[[X,Y]] == 2,Y]] = 0];] (*close the For[
Y,...]*)];(*close the For[X,...]*)CurrentWorld = NextWorld;
Worldplot =
ArrayPlot[CurrentWorld,Mesh -> False];
Pause[0];] (*close the For[t,...]*)
RandomChoice::wghtv: The weights given on the left-hand side of {-1,2}->{0,1} should be a list of positive numerical quantities having the same length as the list given on the right-hand side.
RandomChoice::wghtv: The weights given on the left-hand side of {-1,1} should be a list of positive numerical quantities having the same length as the list given on the right-hand side.
General::stop: Further output of RandomChoice::wghtv will be suppressed during this calculation.
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。