弹性力学方程有限差分法matlab,泊松方程的有限差分法的MATLAB实现

冯立伟　徐涛　屈福志

摘要：泊松方程是物理及工程应用领域中一类非常重要的方程，研究其数值求解方法具有重要意义。给出了使用有限差分法求解泊松方程的计算方法，并讨论了使用MATLAB编写计算程序，使用数值算例和静电场实例进行了数值实验，实验结果与理论一致，检验了算法的有效性。

关键词：泊松方程；五点差分格式；有限差分法

中图分类号：TP311 文献标识码：A 文章编号：1009-3044(2017)13-0233-03

1概述

物理过程，都可用椭圆型方程来描述。其中最典型的方程是泊松(Poisson)方程。传热学中带有稳定热源或内部无热源的稳定温度场的温度分布、流体动力学中不可压缩流体的稳定无旋流动、弹性力学中平衡问题及电磁学中静电场的电势等均满足泊松方程，泊松方程也是数值网格生成技术所遵循的基本方程。因此，研究其数值求解方法具有重要意义。

MATLAB是目前应用最广泛的科学和工程计算软件。MATLAB基于矩阵运算，具有强大数值运算能力，是方便实用、功能强大的数学软件；同时，MATLAB具有强大的图形绘制功能，用户只需提供绘图数据和指定绘图方式，用很少的程序指令就可得到将计算结果转化为直观、形象的圖像。使用MAT-LAB求解微分方程已有大量的研究。因此，近些年来，越来越多的人开始使用MATLAB来求解泊松方程。利用MAT-LAB强大的数值计算能力和图形绘制技术，可以实现使用差分法求解泊松方程并绘制出数值解的二维、三维图像，从而可以更好地理解泊松方程解的物理意义。

本文讨论使用差分法通过MATLAB编程求解二维矩形区域上的泊松方程，并使用两个算例进行检验和对结果进行分析。

边界条件为

将未知解函数在内部节点上的值按行排列，组成解向量为：

3差分格式的求解

为了便于使用MATLAB编写程序，将差分方程转化为矩阵形式：

4数值实验

算例1：

为了分析和比较差分格式在不同步长下的结果，使用2范数意义下的绝对误差和相对误差作为评价指标，表1给出了步长h=0.01取不同值的绝对误差和相对误差

从表1可看出随着网格步长h的减小数值解的绝对误差和相对误差在变小。

算例2：

二维静电场问题。有一个横截面为矩形的无限长金属槽，槽的宽度为16 m，高度为10m，槽的上盖板和两侧面绝缘，盖板处电势为IOOV，两侧及下地面处电势均为0，电势函数满足泊松方程，无场源f(x)=0，使用差分格式计算槽的横截面内的电势分布，计算结果见下图。

5结论

给出了泊松方程的五点差分格式，并使用MATLAB编制求解程序，使用一个数值算例和静电场算例进行了实验，验证了方法的有效性并分析了在不同步长下的误差情况。

弹性力学方程有限差分法matlab,泊松方程的有限差分法的MATLAB实现相关推荐

matlab求高阶微分的命令,matlab求解高阶微分方程辅导
同学你好,我们可以辅导matlab高阶微分方程.关于matlab解决高阶微分方程的知识点如下,如果同学感觉还死有不明白的地方,可以直接联系我们,安排专业老师进行一对一授课讲解. 微分方程(Differ ...
matlab求函数极值教程,MATLAB程序设计教程(7)—MATLAB解方程与函数极值
MATLAB程序设计教程(7)--MATLAB解方程与函数极值第7章MATLAB解方程与函数极值 7.1 线性方程组求解 7.2 非线性方程数值求解 7.3 常微分方程初值问题的数值解法 7. ...
matlab选修结课作业,matlab结课作业
文档来源为:从网络收集整理.word 版本可编辑.欢迎下载支持. 山西大同大学 matlab 课程结课作业 MATLAB 程序应用姓名: 课程序号: 2 班级: 学号: 2013 ...... 山西 ...
matlab 自定义对象,面向对象：MATLAB的自定义类 [MATLAB]
这几天刚刚开始学习MATLAB的面向对象编程.以前做的事情都是用MATLAB写一些简单的脚本或者函数,这方面MATLAB成熟的函数和直截了当的矩阵运算方法和语法都很容易上手,方便人专注于算法本身.前几 ...
matlab负荷预测,负荷预测matlab
本文主要针对BP 神经网络应用于电力系统短期负荷预测做了进一步的研究,并通过MATLAB设计BP神经网络,仿真结果表明BP神经网络在短期负荷预测中的应用是可行的,能较好...... MATLAB ...
matlab的实用与仿真,MATLAB建模与仿真实用教程
前言第1章MATLAB建模基础 1.1 MATLAB简介 1.1.1 MATLAB的安装 1.1.2 MATLAB通用命令 1.1.3 MATLAB应用窗口简介 1.2 MATLAB数学建模概述 1 ...
matlab绘图z＝sin（x_「matlab画三维图」Matlab 应用之绘制三维图形（基础篇） - seo实验室...
matlab画三维图在Matlab中,三维图形的绘制包括三维曲线,三维网线图和三维曲面图.闲话不多说,直接进入正题.首先介绍几个函数: 1.plot3(x,y,z,-) 其中,x,y,z为维数相同的 ...
matlab中的导函数驻点,Matlab用导数作定性分析
Matlab用导数作定性分析 5.1知识要点:函数作图 -用导数定性描述函数 [ clf,x=linspace(-8,8,30);f=(x-3).^2./(4*(x-1)); plot(x ...
MATLAB教程（1） MATLAB 基础知识(4)
第七部分:二.三维图二维图和三维图- MATLAB & Simulink- MathWorks 中国折线图 (1) 画图 x = 0:pi/1000:2*pi; y = sin(x); p ...
MATLAB程序详细解析,遗传算法——matlab代码解析
遗传算法--matlab代码解析本文为学习B站老哥数学建模课程之后的一点笔记,图片源自web,代码源自老哥程序包,侵权删. 详细的遗传算法原理不再赘述,百度即可找到. 算法定义遗传算法(GA)是模 ...

弹性力学方程有限差分法matlab,泊松方程的有限差分法的MATLAB实现

弹性力学方程有限差分法matlab,泊松方程的有限差分法的MATLAB实现相关推荐

最新文章

热门文章

弹性力学方程 有限差分法matlab,泊松方程的有限差分法的MATLAB实现

弹性力学方程 有限差分法matlab,泊松方程的有限差分法的MATLAB实现相关推荐

最新文章

热门文章

弹性力学方程有限差分法matlab,泊松方程的有限差分法的MATLAB实现

弹性力学方程有限差分法matlab,泊松方程的有限差分法的MATLAB实现相关推荐