R语言——数学计算,解决常微分方程式

一个有变量的等式称为方程式,其中带有微分符号的方程式称为微分方程式。在微分方程中,对一个变量微分称为常微分方程,对多个变量微分称为偏微分方程。

劳伦兹(Lorenz)方程式是由大气方程简化而来的,方程式如下:

R语言——数学计算,解决常微分方程式

其中,a,b,c为特定参数,在某些特定数值中可以观察到一些现象。下面通过R语言的"deSolve"软件包进行数值计算,并使用"scatterplot3d"软件包进行3D图形的绘制。

进入R之后,输入下面两行代码来安装软件包:

R语言——数学计算,解决常微分方程式

安装之后,加载这两个软件包并定义Lorenz函数,范例如下:

R语言——数学计算,解决常微分方程式

接下来就是进行两面映射图像的设置:

R语言——数学计算,解决常微分方程式

得到三面图像:

R语言——数学计算,解决常微分方程式

然后进行3D合成:

scatterplot3d(out[,-1],type = "l",color=\'red\')

效果如下:

R语言——数学计算,解决常微分方程式

通过上午的学习,大概了解了R和Python不一样的地方.继续努力!

期待自己把函数图像这一块儿学通,有利于了解接下来的偏微分方程学习!

By Z_Howe.01