数理方程基于matlab的数值解法.docVIP

数理方程数值解法与其在matlab软件上的实现 张体强 1026222 廖荣发 1026226 [摘要] 数学物理方程的数值解在实际生活中越来越使用，首先基于偏微分数值解的思想上，通过matlab软件的功能，研究其数学物理方程的数值解，并通过对精确解和数值解进行对比，追究其数值解的可行性，在此，给出相关例子和程序代码，利于以后的再次研究和直接使用。 [关键字] 偏微分方程 数值解 matlab 数学物理方程的可视化 一：研究意义 在我们解数学物理方程，理论上求数学物理方程的定解有着多种解法，但是有许多定解问题却不能严格求解，只能用数值方法求出满足实际需要的近似解。而且实际问题往往很复杂，这时即便要解出精确解就很困难，有时甚至不可能，另一方面，在建立数学模型时，我们已作了很多近似，所以求出的精确解也知识推导出的数学问题的精确解，并非真正实际问题的精确解。因此，我们有必要研究近似解法，只要使所求得的近似解与精确解之间的误差在规定的范围内，则仍能满足实际的需要，有限差分法和有限元法是两种最常用的求解数学物理方程的数值解法,而MATLAB在这一方面具有超强的数学功能，可以用来求其解。 二：数值解法思想和步骤 2.1:网格剖分 为了用差分方法求解上述问题，将求解区域作剖分。将空间区间作等分，将时间区间作等分，并记 。分别称和为空间和时间步长。用两簇平行直线将分割成矩形网格。 2.2：差分格式的建立 ………………………………（1）是平面任一有界域，据Green公式(参考数学物理方程第五章)： 其中。于是可将（1）式写成积分守恒形式： …………………………（2）Lax格式设网格如下图所示 图1 取为以，，和为顶点的矩形。 为其边界，则 …………（3）Lax-Friedrich格式 截断误差为 所以Lax-Friedrich格式的截断误差的阶式 令：则可得差分格式为 2.3差分格式的求解 差分格式 写成如下矩阵形式： 则需要通过对k时间层进行矩阵作用求出k+1时间层。 对上面的矩阵形式我通过C++（或matlab）编出如附录的程序求出数值解、真实解和误差。 例1：如下方程 利用 matlab的数值解法结果并填入表中关系对比如下： 表 j k （x,t） 数值解 真实解 误差 450 100 (0.5,0.1) -0.308981 -0.309017 0.000036 450 200 (0.5,0.2) -0.587647 -0.587785 0.000138 450 300 (0.5,0.3) -0.808731 -0.809017 0.000286 450 400 (0.5,0.4) -0.950609 -0.951056 0.000448 450 500 (0.5,0.5) -0.999409 -1.000000 0.000591 450 600 (0.5,0.6) -0.950496 -0.951057 0.000560 450 700 (0.5,0.7) -0.808539 -0.809017 0.000478 450 800 (0.5,0.8) -0.587437 -0.587785 0.000348 450 900 (0.5,0.9) -0.308833 -0.309017 0.000184 450 1000 (0.5,1.0) 0.000002 -0.000000 0.000002 从上面可以看出，数值解精度相当高的 三：matlab的在数学物理方程上简单的应用 Matlab是一个强大的计算工具，超强的计算能力和绘图能力，下面几个例题来说明matlab数学物理上的应用 例1：将函数1/（1－a）2在z=0 处展成幂级数。 解：syms a； Taylor (1/(1-a)^2,0) ans =1+2*a+3*a^2+4*a^3+5*a^4+6*a^5 例2：写出函数f（x）＝ 1/(x2＋p2 )（a0）的Fourier 变换式。 解：syms x w; syms a potitive f=1/(x^2+p^2); F=Fourier (f,x,w) F=pi*(signum(0,Re(p),0)*cosh(p*w)-2*Heaviside(w)*sinh(p*w)+sinh(p*w))/p 例2：已知函数f（x）＝x3 e－x ，试求取该函数的Laplace 变换,并对结果进行Laplace 反变换。 解：syms x w; f=x^3*exp(-x); F=laplace(f,x,w) F=6/(w+1)^4 对得出的结果进行Laplace 反变换,从而有 ilaplace(F) ans=x^3*exp(-x