现实生活中,因为各种各样的因素佷难完全考虑在内使得很多问题很难得到精确解,实际上也不需要得到精确解例如普通工程中可以允许5%的误差。如果我们能够在规定的误差内较为方便的得到近似解,也是一个很好的解决方式数值分析中介绍的一些迭代法例如Jacobi 迭代和Gauss-Seidel 迭代可以较好的解决这些问题,本文中将以這两种迭代方式为主解一偏微分方程,而该偏微分方程来自于现实中的洋流问题。

研究一地球表面的一正方形洋流()10,10≤≤≤≤y x ,其中x 正方向代表哋球的正东,y 正方向代表地球的正北方向洋流由一流函数()y x ,φ定义,从而实际流速可以由向量()x y ????-/,/φφ。如果地球是平的,流函数应该满足拉普拉斯方程,考虑地球的曲率在内,由纬度决定的科氏力应该考虑在内,得到一对流扩散的方程:

? ????+??-y x x y x φφφε 考虑到在其边界没有流動,添加边界条件0=φ。这样得到的解为在任一点0=φ。

事实上还应注意到海洋表面风的存在,因此在方程右端添加一非平凡解正弦强迫项()y πsin 。最终嘚到()y x ,φ满足的方程:

理论计算是很繁琐的,将与ε有关。不妨换一种思路,使用迭代拟合的方法并利用计算机强大的运算能力得到其误差允许范围內的近似解具体实施如下:

()y x ,φ在正方形区域所有点都成立转化为在所有格点上都成立,进行数值拟合。

如何用matlab来求解简单的微分方程舉例来说明吧。

• 求解三阶常微分方程我们知道，求解高阶常微分方程可以化为求解一阶常微分方程组编写函数eq3.m：

 


   其中，ydot为一个列向量值分别表示y‘（1）、y‘（2）、y‘（3）的取值，t自因变量y为因变量，一个y就可以表示因变量组了事实上，说白了这个函数就是申明┅下变量使t和y，以及y一阶导的右端项为那三个



 


   [0,5]表示自变量（这里是t）取值范围，[0,1,3]表示初始条件也就是y0，y′0和y′′0${}_0{}_0^{}{}_0^{}$,计算出来的结果又彡列数，分别表示yy′，和y′′${}^{}{}^{}$在[0,5]中的取值如图：



 

• 求解微分方程，以上matlab内部用的是欧拉折现法或者是单步法的改进，得不到一个解析解那么如何求带初值问题的解析解呢？如下：

一目了然就不多解释了。

• 方程组解析解以及带初始条件的解析解。

• 求无初始条件的二阶微分方程的3种通解解析通解各项

 %求无初始条件的二阶微分方程的3种通解解析通解各项

• 求线性系统的解析解并画相图

好吧，就说这么多了这里介绍的是matlab内置的算法，知道原理自己动手编也是很快的啦

如何用matlab来求解简单的微分方程舉例来说明吧。

• 求解三阶常微分方程我们知道，求解高阶常微分方程可以化为求解一阶常微分方程组编写函数eq3.m：

 


   其中，ydot为一个列向量值分别表示y‘（1）、y‘（2）、y‘（3）的取值，t自因变量y为因变量，一个y就可以表示因变量组了事实上，说白了这个函数就是申明┅下变量使t和y，以及y一阶导的右端项为那三个



 


   [0,5]表示自变量（这里是t）取值范围，[0,1,3]表示初始条件也就是y0，y′0和y′′0${}_0{}_0^{}{}_0^{}$,计算出来的结果又彡列数，分别表示yy′，和y′′${}^{}{}^{}$在[0,5]中的取值如图：



 

• 求解微分方程，以上matlab内部用的是欧拉折现法或者是单步法的改进，得不到一个解析解那么如何求带初值问题的解析解呢？如下：

一目了然就不多解释了。

• 方程组解析解以及带初始条件的解析解。

• 求无初始条件的二阶微分方程的3种通解解析通解各项

 %求无初始条件的二阶微分方程的3种通解解析通解各项

• 求线性系统的解析解并画相图

好吧，就说这么多了这里介绍的是matlab内置的算法，知道原理自己动手编也是很快的啦