||
我的认识、错误和矩阵的行积定义等等(2)
张学文,20180624
这里谈我在分析水分循环问题时的一些思路、感悟、错误以及感到需要定义一种矩阵的新计算规定等等。这可能要几个博客去表达。现在是第2部分,即蒸发去向矩阵的逆矩阵不是降水来源矩阵。
l 在本系列的第一个博客http://blog.sciencenet.cn/blog-2024-1120503.html 中,我们给出了一个过去认可的在气象学、物理上比较可以接受的各大洲的水分蒸发去向矩阵。也给出了我另外计算出来的各大洲降水来源矩阵。还提到蒸发去向矩阵的逆矩阵不是降水来源矩阵。
l 现在我们通过计算给出的蒸发去向矩阵的逆矩阵看看它是什么样子。手工做这个求逆矩阵计算显然很费力,而excel中有现成的命令完成这个计算。于是我就用 EXCEL的命令计算了它。下面是原始的蒸发去向矩阵(前面博客中用到的)以及计算的它的逆矩阵。
表2.1 地球各大洲的蒸发水分去向矩阵
(即水分次左侧各大洲水分去了上面第一行的各个大洲的水分的百分比值)
1欧洲 | 2亚洲 | 3非洲 | 4澳洲 | 5北美洲 | 6南美洲 | 7南极洲 | 8海洋 | ||
1欧洲 | 0.07 | 0.32 | 0.12 | 0 | 0 | 0 | 0 | 0.49 | 1 |
2亚洲 | 0.05 | 0.33 | 0.06 | 0.03 | 0.03 | 0 | 0 | 0.5 | 1 |
3非洲 | 0.05 | 0.14 | 0.21 | 0 | 0 | 0 | 0 | 0.6 | 1 |
4澳洲 | 0 | 0.025 | 0 | 0.015 | 0 | 0 | 0 | 0.96 | 1 |
5北美洲 | 0 | 0.005 | 0 | 0 | 0.17 | 0.16 | 0 | 0.665 | 1 |
6南美洲 | 0 | 0 | 0 | 0 | 0.07 | 0.3 | 0.015 | 0.615 | 1 |
7南极洲 | 0 | 0 | 0 | 0 | 0 | 0.02 | 0.02 | 0.96 | 1 |
8海洋 | 0.011 | 0.049 | 0.037 | 0.016 | 0.028 | 0.05 | 0.005 | 0.804 | 1 |
表2.2 用excel 里MINVERSE 命令求得的逆矩阵如下。
这个具有很多负值的矩阵显然不可能是各大洲的降水来源矩阵
1欧洲 | 2亚洲 | 3非洲 | 4澳洲 | 5北美洲 | 6南美洲 | 7南极洲 | 8海洋 | |
1欧洲 | 76.17381 | -73.8039 | -61.0563 | -86.1708 | -11.478 | -28.1616 | -33.6707 | 219.1674 |
2亚洲 | -9.16754 | 12.93087 | 8.141197 | 14.07767 | 1.887026 | 4.8527 | 5.721273 | -37.4432 |
3非洲 | -12.8878 | 9.841445 | 13.7909 | 7.41733 | 1.035405 | 3.430266 | 3.778915 | -25.4065 |
4澳洲 | -4.04913 | -1.62292 | -15.3776 | -39.9991 | -12.9884 | -8.97808 | -19.027 | 103.0421 |
5北美洲 | -1.32208 | 1.230057 | -0.46799 | -7.7394 | 6.726919 | -4.55799 | 2.181159 | 4.94933 |
6南美洲 | 0.435987 | -0.41684 | 0.125359 | 2.37987 | -1.59556 | 4.633403 | -3.11266 | -1.44955 |
7南极洲 | -14.9323 | 15.36323 | -1.48203 | -64.7821 | -5.78694 | -5.30109 | 45.99403 | 31.92714 |
8海洋 | 0.302006 | -0.31138 | 0.028264 | 1.300046 | 0.153802 | 0.01391 | 0.148305 | -0.63495 |
l 显然,用蒸发去向矩阵通过求逆矩阵而得到的降水来源(水分来自各个大洲的百分比)矩阵中出现了很多负值。这显然是错误的。
l 所以我在2017年11月14日的博客http://blog.sciencenet.cn/blog-2024-1085247.html 中的这个认识是错误的。借此机会也就此说明该博客的这个结论不对。
l 关于如何从各大洲蒸发去向矩阵和蒸发量计算合理的降水来源的事,我们在后面的博客中再具体讨论。
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-12-23 23:53
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社