从动态平衡分布求转移矩阵--《气象随机场-22》

张学文，2014/8/28-31

2014,9,3注：本稿的第4段在9月3日发现错误，这也影响了后面的结论。我已经用一个新版代替本稿，见

http://blog.sciencenet.cn/blog-2024-824687.html 。但是此稿依然保留在这里。

经过前面的分析与讨论，我们逐步认识到，面积很大的气象场（这通常理解是全球、半球或者面积很大的一个空间、区域）中的某气象变量在某时刻的分布函数经常具有稳定性。而前面还说明可以存在一种转移矩阵，把分布函数与转移矩阵做乘法而获得的新分布函数（一个时间步长）依然与原分布函数相同。而此时的分布函数称为极限分布函数。

我们后面着重分析这种具有时间不变性的极限分布函数与其转移矩阵应当是什么关系。或者说我们是否可以从分布函数去推求它要求的转移矩阵。即可以使原分布函数在经过一个（以致充分多个）时间步长的转移矩阵作用以后依然没有变化。

上一讲我们把这归结为公式：

xA=x

现在探索已经知道x求A。这里的x是一个矢量，它的各个分量是分布函数的各个函数值（自变量是规定的有限个离散状态的相空间），而A是待求的转移矩阵。

也许应当交代依据笔者的线性代数知识很差，以下的探索属于个人的经验摸索。承望高手指教。而如果这里有一定的创新性，也希望指出，而引用者请注明出处。

1.      矩阵中的未知数数量：我们这里的x是已经知道的有n个离散值的n维矢量。而待求的转移矩阵A自然n²个未知数（方阵中的元素数量）。我们要从x的n个已知数据中求得n²个未知数数据。

2.      矩阵中的“0”很多：注意我们在第16讲（状态转移矩阵的时间步长问题-《气象随机场-16》http://blog.sciencenet.cn/blog-2024-816969.html ），单独讨论过转移矩阵的时间步长问题，而特别指出我们可以把时间步长取得很短，以使在这么短的时间内气象场内各点（有时是指气象场内的元面积，有时是指空气微团）的气象状态仅可能是依然留在本相格（状态区间）内或者移到它的左右两个临近相格，而转移到离它更远的相格的百分率是0。所以对于有n²个元素转移矩阵A，它最多有3n-2个未知数。即此转移矩阵至少有[n²-(3n-2）]个“0”。具体到矩阵的每一行，则其第1行和最后一行，仅有两个未知数，其他的行有三个未知数。

3.      转移矩阵的每行的合计值=1：转移矩阵的每行的各个元素的合计值代表了从本状态转移到各个状态的百分比，它的合计值自然=1。这样我们就可以获得n个合计值=1的方程式，从而使未知数再减少n个。于是我们未知数就仅有（2n-2）个了。

4.      已知分布函数解决了（n-1）个未知数：由于本问题中的极限分布函数矢量本身就代表着n个已知数，它们显然参与（限制）着转移矩阵的各个元素值。所以n维已知矢量本身就获得(n-1)个方程，从而可以获得(n-1)个未知数（注：由于分布函数的各个分量的合计值应当=1的限制，所以它仅具有(n-1)个独立已知数）。这样我们的未知数就从（2n-2）个减少为只有（n-1）个了

5.      细致平衡原则：在已经达到极限分布的情况下（我们现在就是），在马尔科夫过程的理论中有一条所谓细致平衡原则。即这时各个相格中的存在量（也就是分布函数矢量的每个分量）与它转移到指定相格的转移率（矩阵的元素）的乘积，与对方转移回来的对应量是相等的。写成为公式就是

          x_ia_i,j=x_ja_j,i       （i.j=1,2，…，n）

     这里的x_i或者x_j表示分布函数的对这个矢量的对应分量，而a_i,j，a_j,i表示转移矩阵中的对应元素值。它体现在达到动态平衡时任何一个相格中从另外相格的流入量与自己流出到对方的量是相等的。对此在马尔科夫过程的理论中有论及。我们这里就不再多说了。

       在转移矩阵的非0元素中，这个细致平衡关系又使未知数减少很多（本段待完善）。

6.    给一个规定速度值：我们用限定一个时间步长相当短的规定已经使很多转移矩阵的元素值=0了。但是在没有具体规定时间步长的情况下，我们还是需要在转移矩阵的各个元素中，人为规定一个元素的值（一个时间步长的转移百分比），以具体体现转移的速度。在此我们一般规定转移矩阵中的a₂₂的值是我们事先给定的，a₂₂的含义是处于第2个相格中的空气状态在下一个时间步长依然是留在本相格的百分比。而且我们一般令a₂₂=0.98，即一个时间步长内处于第2个相格的空气维持原状态的百分率是98%。这样我们就又减少了一个未知数。

7.      转移矩阵的“解”不是唯一的：前面的分析指出如何在已知极限分布函数的情况下，并且在规定了a₂₂=0.98的情况下，可以（？）求得对应的转移矩阵。显然，如果修改的a₂₂值，我们还会获得另外的转移矩阵。这说明转移矩阵并不是唯一的。即一个极限分布函数在不同转移速度下，对应这不同的转移矩阵。

8.      在极限分布函数已知的情况下、在限定一次转移最多是可以到达最近的邻态的假设下、在a₂₂=0.98的规定下，在细致平衡原理的配合下，我们在理论上分析出可以从极限分布函数求得它要求的转移矩阵。极限分布函数也就是分布函数保持动态平衡（不变化）的分布函数。而它对应的转移矩阵则体现在动态平衡下的气象状态的变化规律。

以上是在代数学、马尔科夫过程知识、气象变量知识特点的基础上的从极限分布函数求其转移矩阵的一般分析。我们在下一讲中要以具体的例子具体计算其转移矩阵。