网格及元胞数据初始化：定义了元胞的边长a，元胞的初始尺寸M和N，时间步长dt，仿真终止时间ft等参数。其中，T表示温度，V表示风速，theta表示风向，Ks是可燃物配置格局的更正系数，n是可燃物植被类型，tm是单个元胞燃烧时间，x0和y0表示火源的位置。
初始化元胞状态A：使用一个M行N列的全零矩阵A0，表示初始状态全部未燃烧。然后，将火源的位置设置为1，表示该元胞处于完全燃烧状态。A是当前元胞状态，A1是更新后的元胞状态。
传播速度R：根据公式计算火蔓延速度R0，然后初始化一个M行N列的全一矩阵R，表示初始火蔓延速度。
地形坡度：使用高斯分布生成一个二维地形矩阵Z，表示不同位置处的地形高度。使用meshz函数将地形图绘制出来，并设置颜色映射。
标记火源：在地形图上标记火源位置，用红色星号表示。
4 I5 k0 h, H" m% H) F
