论文笔记|PINN解土壤水分运动方程

笔记内容

遇到的问题，模型的错误，参数的错误。

改进的思路

写在前面

在意识到已有代码和论文不是一个版本并且代码存在多处错误的客观条件下。经过我，一个编程小白弱鸡，的不懈努力，代码终于可以正常运行了。

不过新的问题又出现了，预测的各个参数值不在正常范围内……

好在模型能够正常运行，慢慢调整参数似乎是必经之路。

此次调整是在模型能够运行的基础上。

01调整模型错误

1.1设定参数

论文设定的参数如下

初始含水率theta = 0.03 对应的 psi =-0.215 m

上边界含水率theta = 0.3 对应的psi =-5.48 m

1.2模型参数

仔细对照代码文件，发现代码中设定的参数

K_s=13.524cm/h

时间/空间 步长设定

[[0,0],

[10,450]]

表示时间0-10 h，计算深度0- 450 cm。

但是，很奇怪，至少从结果看来，这里很有可能有问题。逻辑上0-10h，和0-600min意义相同，但是计算机只会带入数值进行计算。

以下是修正后的时间空间域

 在论文中，采用正交实验设计决定最优神经网络架构，结构如下

以下 在代码中将以上参数带入 

02总结

更改的参数有

K_s

时间/空间 域的数值

神经网络架构

此外，在此基础上我还调整了模型的学习率。原设定学习率为0.9

我更改为0.6后，含水率输出结果与论文趋势开始接近，更改学习率为0.09后，输出结果与论文趋势更加吻合。

这是所有参数未更改前的结果

这是更改参数，学习率为0.6

这是更改参数，学习率为0.09