||
我的认识、错误和矩阵的行积定义等等(5)---用矩阵行积等运算求降水来源矩阵
张学文,20180701
本博客系列谈我在分析水分循环问题时的一些思路、感悟、错误以及感到需要定义一种矩阵的新计算规定等等。现在是第5部分,用列矢量与矩阵的行积语言表达地球的各个地区蒸发量、蒸发去向矩阵如何计算出来各个地区的降水来源矩阵。
上一个博客介绍了我们引入的列矢量与矩阵的行积计算的含义http://blog.sciencenet.cn/blog-2024-1121132.html 。本博客要利用包括这种行积运算的的语言去概括水分循环问题中从地球各大洲的蒸发量、蒸发去向矩阵求得各大洲的降水来源矩阵的计算过程。
l 我们现在以地球各纬度带的蒸发量情况为例说明它。下面的表5.1给出不同纬度带的蒸发量和蒸发的水分去了各个纬度带的百分比。
表5.1各个纬度带(xi)的蒸发量(列矢量,红字)与它们在各个纬度带(yi)变成降水的比例矩阵(黄底色)
(注:后面的表中的各个xi,表示地球上的互不相容的各个地区,各个yi 也是这个含义。我们仅在表5.1中具体列出它们代表的地区。现在的地区是指全球的各个纬度带)
水分降落的纬度带 | ||||||||||
y1 | y2 | y3 | y4 | y5 | y6 | y7 | ||||
蒸发量 | 90N-70N | 70N-40N | 40N-10N | 10N-10S | 10S-40S | 40S-70S | 70S-90S | |||
蒸发 | x1 | 1602 | 0.8 | 0.15 | 0.05 | 0 | 0 | 0 | 0 | 1 |
的 | x2 | 36819.6 | 0.03 | 0.8 | 0.17 | 0 | 0 | 0 | 0 | 1 |
纬度带 | x3 | 134921.6 | 0 | 0.18 | 0.62 | 0.2 | 0 | 0 | 0 | 1 |
x4 | 115586.1 | 0 | 0 | 0.15 | 0.67 | 0.18 | 0 | 0 | 1 | |
x5 | 159949 | 0 | 0 | 0 | 0.25 | 0.55 | 0.2 | 0 | 1 | |
x6 | 45637.9 | 0 | 0 | 0 | 0 | 0.09 | 0.9 | 0.01 | 1 | |
x7 | 667.8 | 0 | 0 | 0 | 0 | 0.04 | 0.11 | 0.85 | 1 |
l 列矢量与矩阵的行积运算:用蒸发量列矢量分别乘以各个行的矩阵元素以后得到的矩阵。表5.2
表5.2用行积运算获得的各个纬度带xi蒸发在各个纬度yi到形成的降水量矩阵
y1 | y2 | y3 | y4 | y5 | y6 | y7 | |
x1 | 1281.6 | 240.3 | 80.1 | 0 | 0 | 0 | 0 |
x2 | 1104.588 | 29455.68 | 6259.332 | 0 | 0 | 0 | 0 |
x3 | 0 | 24285.89 | 83651.39 | 26984.32 | 0 | 0 | 0 |
x4 | 0 | 0 | 17337.92 | 77442.69 | 20805.5 | 0 | 0 |
x5 | 0 | 0 | 0 | 39987.25 | 87971.95 | 31989.8 | 0 |
x6 | 0 | 0 | 0 | 0 | 4107.411 | 41074.11 | 456.379 |
x7 | 0 | 0 | 0 | 0 | 26.712 | 73.458 | 567.63 |
l 求各个纬度带yi的各种降水来源的合计值:即把上面的矩阵按列求合计值,而获得各个纬度带的降水总量表5.3(一个行矢量)
表5.3各个纬度带获得的降水量(行矢量)
y1 | y2 | y3 | y4 | y5 | y6 | y7 | |
列合计值 | 2386.188 | 53981.87 | 107328.7 | 144414.3 | 112911.6 | 73137.37 | 1024.009 |
l 用行矢量与矩阵的列积运算求得各个纬度带降水来源比例矩阵:即用行矢量表3的各个元素值的倒数(成为分母)与表2的对应列的各个值做乘法(列积运算)。其结果见表5.4
表5.4各个纬度带yi的降水中来自各个纬度带xi 的蒸发水分占有的比例(合计值=1)是列积运算的结果
y1 | y2 | y3 | y4 | y5 | y6 | y7 | |
x1 | 0.537091 | 0.004451 | 0.000746 | 0 | 0 | 0 | 0 |
x2 | 0.462909 | 0.545659 | 0.058319 | 0 | 0 | 0 | 0 |
x3 | 0 | 0.44989 | 0.779394 | 0.186854 | 0 | 0 | 0 |
x4 | 0 | 0 | 0.16154 | 0.536254 | 0.184264 | 0 | 0 |
x5 | 0 | 0 | 0 | 0.276893 | 0.779123 | 0.437393 | 0 |
x6 | 0 | 0 | 0 | 0 | 0.036377 | 0.561602 | 0.445679 |
x7 | 0 | 0 | 0 | 0 | 0.000237 | 0.001004 | 0.554321 |
1 | 1 | 1 | 1 | 1 | 1 | 1 |
l 表5.4已经是求得的各个纬度带的降水来源比例矩阵了。不过它是以列表达降水地区的,每列的合计值=1。这与我们的前面的习惯不同。我们把这个矩阵转置,就获得下面的矩阵,表5.5。这里左侧的各个yi是降水纬度带,而横的各个xi是水汽来源的纬度带。这个矩阵的每一行的合计值=1 。
表5.5各个纬度带yi的降水中来自各个纬度带xi的蒸发所占有的比例表(矩阵)
x1 | x2 | x3 | x4 | x5 | x6 | x7 | 合计值 | |
y1 | 0.537091 | 0.462909 | 0 | 0 | 0 | 0 | 0 | 1 |
y2 | 0.004451 | 0.545659 | 0.44989 | 0 | 0 | 0 | 0 | 1 |
y3 | 0.000746 | 0.058319 | 0.779394 | 0.16154 | 0 | 0 | 0 | 1 |
y4 | 0 | 0 | 0.186854 | 0.536254 | 0.276893 | 0 | 0 | 1 |
y5 | 0 | 0 | 0 | 0.184264 | 0.779123 | 0.036377 | 0.000237 | 1 |
y6 | 0 | 0 | 0 | 0 | 0.437393 | 0.561602 | 0.001004 | 1 |
y7 | 0 | 0 | 0 | 0 | 0 | 0.445679 | 0.554321 | 1 |
简单地说,当已经知道全球各个分区的蒸发量、蒸发去向矩阵(表5.1)以后通过;
1. 蒸发列矢量与蒸发转移矩阵的行积运算获得各地的降水来源数量矩阵,表5.2;
2. 把表2矩阵的各列元素分别做加法运算获得各地降水总量,表5.3,它是行矢量;
3. 用此行矢量的各个元素的倒数与表2矩阵做列积运算获得各地降水来源的比例矩阵;表5.4
4. 把表5.4这个矩阵做转置运算,获得各个地区yi(左侧,i=1,2,…,7)的降水中来自各个地区xi(i=1,2,…,7)的从左向右的降水来源比例矩阵(每行合计值=1)。
这样我们就用行积、列积等等矩阵的新运算的语言表达了各地蒸发去向求各地的降水来源的4个运算步骤。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-23 19:11
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社