MATLAB在有限差分法中应用.PDFVIP

维普资讯 第21卷 第2期 桂 林 工 学 院 学 报 Ⅷ 21 2 2001年 4月 JOURNALOFGUILIN INSTITUTEOFTECHNOLOGY ADr． 2001 文章编号：1006—544X(2001)02—0104—06 MATLAB在有限差分法中的应用 熊 彬，阮百尧 (桂林工学院资振与环境工程系．广西桂林 541004) 摘 要：将MATLM3应用于有限差分法的正演计算中．充分发挥了其强大而方便的功能。 通过对二维稳定电流场模型的试算表明，MATLAB在解决实际的工程和数学问题中，与其 它计算机程序设计语言 FORTRAN相比，具有使用更为简便、语句功能更强，用户界面 良好的特点，适台于在工程计算．尤其是物探数据处理领域推广应用。 关键词：MATLM3；有限差分；正演计算；二维电场 中图分类号：P631322 文献标识码 ：A 1 MATLAB概述… MATLAB是一个为科学和工程计算而专门设计的高级交互式的软件包。它集数值分析、矩阵运 算、信号处理和图形显示于一体，构成 了一个方便的、界面友好的用户环境。在这个环境下，对所要 求解的问题，用户只需简单的列出数学表达式，其结果便 以数值或图形方式显示出来。例如，为解用 矩阵／向量形式表达的形如 A．z：b的线性方程组，首先生成系数矩阵 A和右侧的b．然后简单地敲八 z=A／b即可。 MAⅡAB中有大量的命令和事先定义的可用函数集，也可通称为MATLABM文件，这就使得用 它来求解问题通常比传统编程快得多；另外一点，也是它最重要的特点，易于扩展。它允许用户 自行 建立完成指定功能的M文件．从而构成适合于其它领域的工具箱。 MATLAB既是一种编程环境，又是一种程序设计语言。它与其它高级程序设计语言C、Fortran等 一 样，也有其内定的规则，但其规则更接近于数学表示，使用起来更为方便，避免了诸如C、Fortrna 语言的许多限制，比方说，变量、矩阵无须事先定义；其次，它的语句功能之强大，是其它语言所无 法比拟的，例如，一条 fft语句就可完成对指定数据的快速傅氏变换．这就相当于几十条甚至于几百条 c或Fortrna语言语句的功能；再者，也是广大程序爱好者所梦寐以求的，MATLAB提供 了良好的用 户界面，许多函数本身会 自动绘制出图形，而且会 自动选取坐标刻度。 2 稳定电流场的边值问题 在数学上解各种物探方法的正演问题．通常需要求解相应的物理场所满足的数学物理方程。其中， 最常遇到的是下列椭圆型偏微分方程 ： - (pVU) qU= f 收{斋日期：2000 06 17；修订日期：2000 07—14 基金项目：教育部高等学校骨干教师资助计划资助项 目 (桂 [2000]413)；广西 自然科学基金资助项 目 (桂科 白 0007024)；广西培养新世纪学术和技术带头人专项资金资助项目 (99216) 作者简介：熊 彬 (1974 )．男．湖北剂州人，在读硕士研究生．地球探测与信息技术专业，主要从事电阻率法的正演 模拟研究。 维普资讯 第 21卷 第2期 熊 彬等：M．TLAB在有限差分法中的应用 式中：u一物理场的标量位或矢量位； ，q一表征介质性的参数；严 场源项。它们一般都是空间坐 标的函数。微分算子加在标量之前表示梯度 而加在矢量之前，则表示取散度。 椭圆型方程最简单的形式 就是Laplace方程 a2u／Ox +a2¨／a =0 和 Poisson方程 a2u／a +a2u／a，=f(，) (2) 其中：为 z．)已知函数。本文主要以(1)和(2)为模型研究椭圆型方程边值 问题