利用MATLAB实现基于流域水热平衡和水量平衡方差径流变化归因分析

Jude · 发表于 2018-10-5 15:34:28

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

   本帖参考了杨大文老师的水热耦合方程，如图片所示。该方程其实有点像budyko假设下的傅抱璞公式。其中R表示流域多年径流深；P表示流域多年降雨深；n为下垫面指数，反映人类活动影响；E0是潜在蒸发量。
   基于该方程，首先对流域径流深突变点进行识别，利用弹性系数法，分解气候因素（P、E0）和人类活动因素n分别造成径流深改变的贡献率，来达到流域径流深归因分析的目的。
   根据图中公式，以及P值、R值、E0值，可以借由MATLAB的solve语句反解出n值，然后找到突变点，按照文献中弹性系数的定义，构造偏导公式，利用MATLAB中的syms变量和diff函数完成求解，相应的代码见附件。

Jude · 发表于 2018-10-5 15:37:51

附件中的程序，利用文献中数据亲测算过，基本无误差，有相关研究方向的大神们，可以下载来使用。另外，只要将代码中的水热平衡方程换成傅抱璞公式，也可以用傅抱璞公式完成径流变化的归因分析。

yaotianci · 发表于 2018-10-9 00:16:43

感谢分享，谢谢

ztf_qy · 发表于 2018-10-11 09:14:12

感谢分享，学习一下！

lxj7177 · 发表于 2018-10-12 10:13:43

感谢分享，多多学习~

早西v · 发表于 2018-10-20 11:48:33

你好我试了一下n的公式为什么返回是空值呢？

Jude · 发表于 2018-10-20 20:16:59

早西v 发表于 2018-10-20 11:48
你好我试了一下n的公式为什么返回是空值呢？

我觉得是数据问题。你的降雨是84，产生径流深是150，这显然是不可能的。如果你把Q=150换成Q=15，就可以计算出结果。

早西v · 发表于 2018-10-21 15:59:12

Jude 发表于 2018-10-20 20:16
我觉得是数据问题。你的降雨是84，产生径流深是150，这显然是不可能的。如果你把Q=150换成Q=15，就可以计 ...

嗯嗯对谢谢啦

sciren · 发表于 2019-4-27 17:51:59

就是输入P值、R值、E0值然后用代码算的n值？然后还得到P E0的弹性系数吗

sciren · 发表于 2019-4-27 17:56:39

早西v 发表于 2018-10-20 11:48
你好我试了一下n的公式为什么返回是空值呢？

代码怎么用啊，，，您的1961.txt是什么数据啊为啥你的Q A E都注释掉了呢

		自动登录	找回密码
密码			立即注册

[源程序] 利用MATLAB实现基于流域水热平衡和水量平衡方差径流变化归因分析

登录后查看更多精彩内容~

评分