有限差分法例题

有限差分法及其应用

1有限差分法简介

有限差分法（FDM）是计算机数值模拟最早采用的方法，至今仍被广泛运用。该方程将解域划分为差分网格，用有限个网络节点代替连续的求解域。有限差分法通过泰勒级数展开等方法，把控制方程中的导数用网格节点上的函数值得差商代替进行离散，从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似值解法，数学概念直观，表达简单，是发展较早且比较成熟的数值方法。

2有限差分法的数学基础

有限差分法的数学基础是用差分代替微分，用差商代替微商而用差商代替微商的意义是用函数在某区域内的平均变化率来代替函数的真是变化率。而根据泰勒级数展开可以看出，用差商代替微商必然会带来阶段误差，相应的用差分方程代替微分方程也会带来误差，因此，在应用有限差分法进行计算的时候，必须注意差分方程的形式，建立方法及由此产生的误差。

3有限差分解题基本步骤

有限差分法的主要解题步骤如下： 1） 建立微分方程

根据问题的性质选择计算区域，建立微分方程式，写出初始条件和边界条件。 2） 构建差分格式

首先对求解域进行离散化，确定计算节点，选择网格布局，差分形式和步长；然后以有限差分代替无线微分，以差商代替微商，以差分方程代替微

penglining 发表于 2007-5-16 08:26 有限元法、边界元法、有限差分法的区别和各自的优点

请问：有限元法、边界元法、有限差分法等方法有哪些区别和各自的优点？尤其是在声学方面。 谢谢！

fossiler 发表于 2007-5-19 14:00 网格的跑分上不同,差分要求模型规则,有限元可以是任意不规则模型,

hillyuan 发表于 2007-5-21 17:45 FEM: irregular grid-> easy to describe complex shape, hard in mesh generation

\\.a4hj

FDM: regular mesh -> easy in grid generation, hard to describe complex shape=> less accurate than FEM

BEM: irregular mesh in boundary -> mesh generation much easier than that of FEM. need much less computation resource than the above two. BUT ne

华北电力大学本科毕业设计（论文）

二维抛物方程的有限差分法

摘要

二维抛物方程是一类有广泛应用的偏微分方程，由于大部分抛物方程都难以求得解析解，故考虑采用数值方法求解。有限差分法是最简单又极为重要的解微分方程的数值方法。本文介绍了二维抛物方程的有限差分法。

首先，简单介绍了抛物方程的应用背景，解抛物方程的常见数值方法，有限差分法的产生背景和发展应用。讨论了抛物方程的有限差分法建立的基础，并介绍了有限差分方法的收敛性和稳定性。其次，介绍了几种常用的差分格式，有古典显式格式、古典隐式格式、Crank-Nicolson隐式格式、Douglas差分格式、加权六点隐式格式、交替方向隐式格式等，重点介绍了古典显式格式和交替方向隐式格式。进行了格式的推导，分析了格式的收敛性、稳定性。并以热传导方程为数值算例，运用差分方法求解。通过数值算例，得出古典显式格式计算起来较简单，但稳定性条件较苛刻；而交替方向隐式格式无条件稳定。

关键词：二维抛物方程；有限差分法；古典显式格式；交替方向隐式格式

I

华北电力大学本科毕业设计（论文）

FINITE DIFFERENCE METHOD FOR TWO-DIMENSIONAL PARABOLIC

EQUATION

第五次作业（前三题写在作业纸上）

一、用有限差分方法求解一维非定常热传导方程，初始条件和边界条件见说明.pdf文件，热扩散系数α=const，

?T?2T??2 ?t?x1. 用Tylaor展开法推导出FTCS格式的差分方程

2. 讨论该方程的相容性和稳定性，并说明稳定性要求对求解差分方程的影响。 3. 说明该方程的类型和定解条件，如何在程序中实现这些定解条件。

4. 编写M文件求解上述方程，并用适当的文字对程序做出说明。（部分由网络搜索得

到，添加，修改后得到。） function rechuandaopde

%以下所用数据，除了t的范围我根据题目要求取到了20000，其余均从pdf中得来 a=0.00001;%a的取值 xspan=[0 1];%x的取值范围 tspan=[0 20000];%t的取值范围

ngrid=[100 10];%分割的份数，前面的是t轴的，后面的是x轴的 f=@(x)0;%初值

g1=@(t)100;%边界条件一 g2=@(t)100;%边界条件二

[T,x,t]=pdesolution(a,f,g1,g2,xspan,tspan,ngrid);%计算所调用的函数 [x,t]=meshgrid(x,t);

mesh

第五次作业（前三题写在作业纸上）

一、用有限差分方法求解一维非定常热传导方程，初始条件和边界条件见说明.pdf文件，热扩散系数α=const，

?T?2T??2 ?t?x1. 用Tylaor展开法推导出FTCS格式的差分方程

2. 讨论该方程的相容性和稳定性，并说明稳定性要求对求解差分方程的影响。 3. 说明该方程的类型和定解条件，如何在程序中实现这些定解条件。

4. 编写M文件求解上述方程，并用适当的文字对程序做出说明。（部分由网络搜索得

到，添加，修改后得到。） function rechuandaopde

%以下所用数据，除了t的范围我根据题目要求取到了20000，其余均从pdf中得来 a=0.00001;%a的取值 xspan=[0 1];%x的取值范围 tspan=[0 20000];%t的取值范围

ngrid=[100 10];%分割的份数，前面的是t轴的，后面的是x轴的 f=@(x)0;%初值

g1=@(t)100;%边界条件一 g2=@(t)100;%边界条件二

[T,x,t]=pdesolution(a,f,g1,g2,xspan,tspan,ngrid);%计算所调用的函数 [x,t]=meshgrid(x,t);

mesh

商法例题

第一章 国际商法导论

香港同进公司为了其在中国内地合作经营所需设备的购买与安装，与香港美善公司在香港签订了供应和安装设备的合同，规定由美善公司供应并负责安装同进公司所需设备。合同签订之后，同进公司用港币预付了合同的部分价款，美善公司则在设备安装所在地的中国内地某市向工商管理部门办理了安装登记证，同时提供进口设备进行安装。后来，因美善公司安装的部分设备与会谈规定的品名不符，并且还有部分设备未进行安装，同进公司便拒绝支付所欠价款。于是，美善公司在设备安装所在地某人民法院对同进公司提起诉讼，要求其支付所欠合同价款并赔偿利息损失。在开庭审理时，被告以双方当事人是香港公司并在香港签订合同为由，要求依香港法律确认合同无效，原告主张应依中国法律确认合同有效，以令被告按合同规定支付欠款。

问：根据合同法律适用的规定，我国法院应该适用我国的法律，还是适用香港的法律？

答：根据我国《民法通则》和《合同法》关于合同法律适用的规定，如果合同当事人没有选择适用的法律，则以最密切联系的国家或地区的法律为合同准据法。 具体到本案，应适用我国的法律（20页）

第二章 合伙企业法

[案例1] 刘某与A、B两人欲设立一普通合伙企

专题：层次分析法

一般情况下，物流系统的评价属于多目标、多判据的系统综合评价。如果仅仅依靠评价者的定性分析和逻辑判断，缺乏定量分析依据来评价系统方案的优劣，显然是十分困难的。尤其是物流系统的社会经济评价很难作出精确的定量分析。

层次分析法（Analytical Hierarchy Process）由美国著名运筹学家萨蒂（T．L．Saaty）于1982年提出，它综合了人们主观判断，是一种简明、实用的定性分析与定量分析相结合的系统分析与评价的方法。目前，该方法在国内已得到广泛的推广应用，广泛应用于能源问题分析、科技成果评比、地区经济发展方案比较，尤其是投入产出分析、资源分配、方案选择及评比等方面。它既是一种系统分析的好方法，也是一种新的、简洁的、实用的决策方法。

◆ 层次分析法的基本原理

人们在日常生活中经常要从一堆同样大小的物品中挑选出最重的物品。这时，一般是利用两两比较的方法来达到目的。假设有n个物品，其真实重量用w1，w2，…wn表示。要想知道w1，w2，…wn的值，最简单的就是用秤称出它们的重量，但如果没有秤，可以将几个物品两两比较，得到它们的重量比矩阵A。

如果用物品重量向量W=[w1，w2，…wn]右乘矩阵A，则有：

T

由上式可知，n是A的特征值

实验目的：

熟悉有关层次分析法模型的建立与计算，熟悉Matlab的相关命令。

实验准备：

1. 在开始本实验之前，请回顾教科书的相关内容；

2. 需要一台准备安装Windows XP Professional操作系统和装有Matlab的计算机。

实验内容及要求

试用层次分析法解决一个实际问题。问题可参考教材P296第4大题。

实验过程：

某物流企业需要采购一台设备，在采购设备时需要从功能、价格与可维护性三个角度进行评价，考虑应用层次分析法对3个不同品牌的设备进行综合分析评价和排序，从中选出能实现物流规划总目标的最优设备，其层次结构如下图所示。以A表示系统的总目标，判断层中B1表示功能，B2表示价格，B3表示可维护性。C1，C2，C3表示备选的3种品牌的设备。

购买设备A 目标层： 判断层： 功能B1 价格B2 维护性B3 方案层： 产品C1 产品C2 设备采购层次结构图

产品C3

解题步骤：

1、标度及描述

人们定性区分事物的能力习惯用5个属性来表示，即同样重要、稍微重要、较强重要、强烈重要、绝对重要，当需要较高精度时，可以取两个相邻属性之间的值，这样就得到9个数值，即9个标度。

为了便于将比较判断定量化，引入1～9比率标度方法，规定用

时域有限差分法对平面TE波的

MATLAB仿真

摘 要

时域有限差分法是由有限差分法发展出来的数值计算方法。自1966年Yee

在其论文中首次提出时域有限差分以来，时域有限差分法在电磁研究领域得到了广泛的应用。主要有分析辐射条线、微波器件和导行波结构的研究、散射和雷达截面计算、分析周期结构、电子封装和电磁兼容的分析、核电磁脉冲的传播和散射以及在地面的反射及对电缆传输线的干扰、微光学元器件中光的传播和衍射特性等等。

由于电磁场是以场的形态存在的物质，具有独特的研究方法，采取重叠的研究方法是其重要的特点，即只有理论分析、测量、计算机模拟的结果相互佐证，才可以认为是获得了正确可信的结论。时域有限差分法就是实现直接对电磁工程问题进行计算机模拟的基本方法。在近年的研究电磁问题中，许多学者对时域脉冲源的传播和响应进行了大量的研究，主要是描述物体在瞬态电磁源作用下的理论。另外，对于物体的电特性，理论上具有几乎所有的频率成分，但实际上，只有有限的频带内的频率成分在区主要作用。

文中主要谈到了关于高斯制下完全匹配层的差分公式的问题，通过MATLAB程序对TE波进行了仿真，模拟了高斯制下完全匹配层中磁场分量瞬态分布。得到了相应的磁场幅值效果图。

实验一：有限差分法研究静电场边值问题

实验报告人： 年级和班级： 学号：

1. 实验用软件工具： Matlab 2. 实验原理：电磁场课本P36-38

1) 差分方程

2) 差分方程组的解

简单迭代法

高斯-赛德尔迭代法 逐次超松弛法

3. 实验步骤： 1）简单迭代法 程序：

hx=41;hy=21; v1=zeros(hy,hx); v1(hy,:)=zeros(1,hx); v1(1,:)=ones(1,hx)*100; v1(:,1)=zeros(hy,1); v1(:,hx)=zeros(hy,1); v1

v2=v1;maxt=1;t=0; k=0;

while(maxt>1e-5) k=k+1; maxt=0; for i=2:hy-1 for j=2:hx-1

v2(i,j)=(v1(i,j+1)+v1(i+1,j)+v1(i-1,j)+v1(i,j-1))/4; t=abs(v2(i,j)-v1(i,j)); if(t>maxt) maxt=t;end end end v1=v2; end v2 k clf

subplot(1,2,1),m