采用动态的流动意向平均值

src/Constants.java

@@ -18,6 +18,11 @@ public class Constants {
      * <p>
      * -0.99 人群流动最慢速率，甚至完全控制疫情传播
      * 0.99为人群流动最快速率, 可导致全城感染
+     * 流动意向平均值随感染(包括死亡和潜伏和隔离)人数占比的函数，可以使用双曲正切函数拟合得到，函数解析式为
+     * 使用f(x)将感染比例值域从[0, 1]转化为[-1 , 1], f(x) = -2 * (感染人数+死亡人数+隔离人数+潜伏人数)/城市总人数 -1，
+     * 简写为f(x) = -2x + 1
+     * 意向流动平均值g(x) = tanh(2 * f(x)) 简写为 y = tanh(3x)
+     * g(x)函数的定义域为[0, 1], 值域为(-1,1)
      */
     public static float u = 0.99f;
     public static int CITY_PERSON_SIZE = 5000;//城市总人口数量

src/MyPanel.java

@@ -97,8 +97,10 @@ public void paint(Graphics g) {
         g.drawString("急需病床：" + (needBeds > 0 ? needBeds : 0), captionStartOffsetX, captionStartOffsetY + 6 * captionSize);
         g.setColor(new Color(0xccbbcc));
         g.drawString("病死人数：" + PersonPool.getInstance().getPeopleSize(Person.State.DEATH), captionStartOffsetX, captionStartOffsetY + 7 * captionSize);
-        g.setColor(new Color(0xffffff));
+        g.setColor(new Color(0xdd00cc));
         g.drawString("世界时间（天）：" + (int) (worldTime / 10.0), captionStartOffsetX, captionStartOffsetY + 8 * captionSize);
+        g.setColor(new Color(0xffffff));
+        g.drawString("意向流动率：" + String.format("%.2f", Constants.u), captionStartOffsetX, captionStartOffsetY + 9 * captionSize);
 
     }
 

src/Person.java

@@ -70,6 +70,13 @@ public Person(City city, int x, int y) {
      * @return
      */
     public boolean wantMove() {
+        // 计算感染人数占比
+        double nfectedRate = (double) (PersonPool.getInstance().getPeopleSize(State.SHADOW)
+                + PersonPool.getInstance().getPeopleSize(State.CONFIRMED)
+                + PersonPool.getInstance().getPeopleSize(State.DEATH)
+                + PersonPool.getInstance().getPeopleSize(Person.State.FREEZE)) / Constants.CITY_PERSON_SIZE;
+        // 使用由双曲正切函数缩两倍变化而来的函数，拟合流动意向平均值随感染率的变换
+        Constants.u = (float) Math.tanh(2 * ((-2 * nfectedRate) + 1));
         return MathUtil.stdGaussian(sig, Constants.u) > 0;
     }