土壤养分运移对流-弥散模型数值解的MATLAB实现

摘要

本文介绍对流-弥散模型数值解的一种简洁编程方法。应用MATLAB语言编写土壤氮磷钾养分运移基太模型求数值解程序，结果表明，与专业软件Uptake相比，供试植物的氮磷钾总吸收量的相对误差分别只有+2．94％、-1.35％和+2．83％。程序能绘制生动美观的二维图形；数值解采用有限元法,较好地克服了Uptake软件有限差分法收敛较幔，积分区问内各节点解值相对精度不均匀的缺陷；程序简短高效，易于使用和维护。该程序只要稍加修改.就可用于求解目前各种CDE改进模型和土壤水分运动方程，从而使该领域复杂数学问题的求解交得简单易行。