没有合适的资源?快使用搜索试试~ 我知道了~
数学建模_优化问题_偏微分方程的数值解
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 81 浏览量
2022-06-04
23:33:33
上传
评论 1
收藏 385KB PDF 举报
温馨提示
试读
35页
数学建模_优化问题_偏微分方程的数值解
资源推荐
资源详情
资源评论
-360-
第二十章 偏微分方程的数值解
自然科学与工程技术中种种运动发展过程与平衡现象各自遵守一定的规律。这些规
律的定量表述一般地呈现为关于含有未知函数及其导数的方程。我们将只含有未知多元
函数及其偏导数的方程,称之为偏微分方程。
方程中出现的未知函数偏导数的最高阶数称为偏微分方程的阶。如果方程中对于未
知函数和它的所有偏导数都是线性的,这样的方程称为线性偏微分方程,否则称它为非
线性偏微分方程。
初始条件和边界条件称为定解条件,未附加定解条件的偏微分方程称为泛定方程。
对于一个具体的问题,定解条件与泛定方程总是同时提出。定解条件与泛定方程作为一
个整体,称为定解问题。
§1 偏微分方程的定解问题
各种物理性质的定常(即不随时间变化)过程,都可用椭圆型方程来描述。其最典
型、最简单的形式是泊松(Poisson)方程
),(
2
2
2
2
yxf
y
u
x
u
u =
∂
∂
+
∂
∂
=Δ (1)
特别地,当
0),( ≡yxf 时,即为拉普拉斯(Laplace)方程,又称为调和方程
0
2
2
2
2
=
∂
∂
+
∂
∂
=Δ
y
u
x
u
u (2)
带有稳定热源或内部无热源的稳定温度场的温度分布,不可压缩流体的稳定无旋流动及
静电场的电势等均满足这类方程。
Poisson 方程的第一边值问题为
⎪
⎩
⎪
⎨
⎧
Ω∂=Γ=
Ω∈=
∂
∂
+
∂
∂
Γ∈
),(|),(
),(),(
),(
2
2
2
2
yxyxu
yxyxf
y
u
x
u
yx
ϕ
(3)
其中
Ω
为以 Γ 为边界的有界区域,
Γ
为分段光滑曲线,
Γ
Ω
U 称为定解区域,
),(),,( yxyxf
ϕ
分别为
ΓΩ,
上的已知连续函数。
第二类和第三类边界条件可统一表示成
),(
),(
yxu
n
u
yx
ϕα
=
⎟
⎠
⎞
⎜
⎝
⎛
+
∂
∂
Γ∈
(4)
其中
n 为边界 Γ 的外法线方向。当 0
=
α
时为第二类边界条件, 0
≠
α
时为第三类边界
条件。
在研究热传导过程,气体扩散现象及电磁场的传播等随时间变化的非定常物理问
题时,常常会遇到抛物型方程。其最简单的形式为一维热传导方程
)0(0
2
2
>=
∂
∂
−
∂
∂
a
x
u
a
t
u
(5)
方程(5)可以有两种不同类型的定解问题:
初值问题(也称为 Cauchy 问题)
-361-
⎪
⎩
⎪
⎨
⎧
+∞<<∞−=
+∞<<∞−>=
∂
∂
−
∂
∂
xxxu
xt
x
u
a
t
u
)()0,(
,00
2
2
ϕ
(6)
初边值问题
⎪
⎪
⎪
⎩
⎪
⎪
⎪
⎨
⎧
≤≤==
=
<<<<=
∂
∂
−
∂
∂
Tttgtlutgtu
xxu
lxTt
x
u
a
t
u
0),(),(),(),0(
)()0,(
0,00
21
2
2
ϕ
(7)
其中
)(),(),(
21
tgtgx
ϕ
为已知函数,且满足连接条件
)0()(),0()0(
21
glg
=
=
ϕ
ϕ
问题(7)中的边界条件
)(),(),(),0(
21
tgtlutgtu
=
=
称为第一类边界条件。第二类和
第三类边界条件为
Tttgut
x
u
Tttgut
x
u
lx
x
≤≤=
⎥
⎦
⎤
⎢
⎣
⎡
+
∂
∂
≤≤=
⎥
⎦
⎤
⎢
⎣
⎡
−
∂
∂
=
=
0),()(
0),()(
22
1
0
1
λ
λ
(8)
其中
0)(,0)(
21
≥≥ tt
λ
λ
。当
0)()(
21
≡
=
tt
λ
λ
时,为第二类边界条件,否则称为第三
类边界条件。
双曲型方程的最简单形式为一阶双曲型方程
0=
∂
∂
+
∂
∂
x
u
a
t
u
(9)
物理中常见的一维振动与波动问题可用二阶波动方程
2
2
2
2
2
x
u
a
t
u
∂
∂
=
∂
∂
(10)
描述,它是双曲型方程的典型形式。方程(10)的初值问题为
⎪
⎪
⎪
⎩
⎪
⎪
⎪
⎨
⎧
+∞<<∞−=
∂
∂
+∞<<∞−=
+∞<<∞−>
∂
∂
=
∂
∂
=
xx
t
u
xxxu
xt
x
u
a
t
u
t
)(
)()0,(
,0
0
2
2
2
2
2
φ
ϕ
(11)
边界条件一般也有三类,最简单的初边值问题为
-362-
⎪
⎪
⎪
⎩
⎪
⎪
⎪
⎨
⎧
≤≤==
≤≤=
∂
∂
=
<<>
∂
∂
=
∂
∂
=
Tttgtlutgtu
lxx
t
u
xxu
lxt
x
u
a
t
u
t
0)(),(),(),0(
0)(),()0,(
0,0
21
0
2
2
2
2
2
φϕ
如果偏微分方程定解问题的解存在,唯一且连续依赖于定解数据(即出现在方程
和定解条件中的已知函数),则此定解问题是适定的。可以证明,上面所举各种定解问
题都是适定的。
§2 偏微分方程的差分解法
差分方法又称为有限差分方法或网格法,是求偏微分方程定解问题的数值解中应用
最广泛的方法之一。它的基本思想是:先对求解区域作网格剖分,将自变量的连续变化
区域用有限离散点(网格点)集代替;将问题中出现的连续变量的函数用定义在网格点
上离散变量的函数代替;通过用网格点上函数的差商代替导数,将含连续变量的偏微分
方程定解问题化成只含有限个未知数的代数方程组(称为差分格式)。如果差分格式有
解,且当网格无限变小时其解收敛于原微分方程定解问题的解,则差分格式的解就作为
原问题的近似解(数值解)。因此,用差分方法求偏微分方程定解问题一般需要解决以
下问题:
(i)选取网格;
(ii)对微分方程及定解条件选择差分近似,列出差分格式;
(iii)求解差分格式;
(iv)讨论差分格式解对于微分方程解的收敛性及误差估计。
下面我们只对偏微分方程的差分解法作一简要的介绍。
2.1 椭圆型方程第一边值问题的差分解法
以 Poisson 方程(1)为基本模型讨论第一边值问题的差分方法。
考虑 Poisson 方程的第一边值问题(3)
⎪
⎩
⎪
⎨
⎧
Ω∂=Γ=
Ω∈=
∂
∂
+
∂
∂
Γ∈
),(|),(
),(),(
),(
2
2
2
2
yxyxu
yxyxf
y
u
x
u
yx
ϕ
取
τ
,h 分别为
x
方向和
y
方向的步长,以两族平行线
τ
jyykhxx
jk
==
=
=
,
),2,1,0,( L±±=jk 将定解区域剖分成矩形网格。节点的全体记为
},,,|),{( 为整数jijykhxyxR
jkjk
τ
=
=
= 。定解区域内部的节点称为内点,记内点
集
ΩIR 为
τ
h
Ω 。边界 Γ 与网格线的交点称为边界点,边界点全体记为
τ
h
Γ 。与节点
),(
jk
yx 沿
x
方向或
y
方向只差一个步长的点 ),(
1 jk
yx
±
和 ),(
1±jk
yx 称为节点
),(
jk
yx
的相邻节点。如果一个内点的四个相邻节点均属于
Γ
Ω
U ,称为正则内点,正
则内点的全体记为
)1(
Ω
,至少有一个相邻节点不属于
Γ
Ω
U 的内点称为非正则内点,
非正则内点的全体记为
)2(
Ω 。我们的问题是要求出问题(3)在全体内点上的数值解。
为简便记,记
),(),,(),(),,(),(
, jkjkjkjk
yxffyxujkuyxjk
=
=
= 。对正则内点
-363-
)1(
),( Ω∈jk ,由二阶中心差商公式
)(
),1(),(2),1(
2
2
),(
2
2
hO
h
jkujkujku
x
u
jk
+
−+−+
=
∂
∂
)(
)1,(),(2)1,(
2
2
),(
2
2
τ
τ
O
jkujkujku
y
u
jk
+
−+−+
=
∂
∂
Poisson 方程(1)在点
),( jk 处可表示为
)(
)1,(),(2)1,(),1(),(2),1(
22
,
22
τ
τ
++=
−
+
−
+
+
−
+
−
+
hOf
jkujkujku
h
jkujkujku
jk
(12)
在式(12)中略去
)(
22
τ
+hO ,即得与方程(1)相近似的差分方程
jk
jkjkjkjkjkjk
f
uuu
h
uuu
,
2
1,,1,
2
,1,,1
22
=
+
−
+
+−
−+−+
τ
(13)
式( 13)中方程的个数等于正则内点的个数,而未知数
jk
u
,
则除了包含正则内点处
解
u 的近似值,还包含一些非正则内点处 u 的近似值,因而方程个数少于未知数个数。
在非正则内点处 Poisson 方程的差分近似不能按式(13)给出,需要利用边界条件得到。
边界条件的处理可以有各种方案,下面介绍较简单的两种。
(1)直接转移。用最接近非正则内点的边界点上的
u 值作为该点上 u 值的近似,
这就是边界条件的直接转移。例如,点
),( jkP 为非正则内点,其最接近的边界点为 Q
点,则有:
)()(
,
QQuu
jk
ϕ
==
,
)2(
),( Ω∈jk
(ii)线性插值。这种方案是通过用同一条网格线上与点
P
相邻的边界点
R
与内点
T
作线性插值得到非正则内点 ),( jkP 处 u 值的近似。由点
R
与
T
的线性插值确定
)(Pu
的近似值为
)()(
,
Tu
dh
d
R
dh
h
u
jk
+
+
+
=
ϕ
其中
RPd = , PTh = ,其截断误差为 )(
2
hO 。
由式(13)所给出的差分格式称为五点菱形格式,实际计算时经常取
τ
=h ,此时
五点菱形格式可化为
jkjkjkjkjkjk
fuuuuu
h
,,1,1,,1,1
2
)4(
1
=−+++
−+−+
(14)
简记为
jkjk
fu
h
,,
2
1
=◊
(15)
其中
jkjkjkjkjkjk
uuuuuu
,1,1,,1,1,
4
−
+
++=◊
−+−+
。
求解差分方程组最常用的方法是同步迭代法,同步迭代法是最简单的迭代方式。除
边界节点外,区域内节点的初始值是任意取定的。
例 1 用五点菱形格式求解 Laplace 方程第一边值问题
-364-
⎪
⎩
⎪
⎨
⎧
Ω∂=Γ++=
Ω∈=
∂
∂
+
∂
∂
Γ∈
])1lg[(|),(
),( 0
22
),(
2
2
2
2
yxyxu
yx
y
u
x
u
yx
其中 }1,0|),{( ≤≤=Ω yxyx 。取
3
1
==
τ
h
。
图 1 网格划分图
解 节点编号为
),( jk
,
3,2,1,0=k
,
3,2,1,0
=
j
。网格中有四个内点,均为正则
内点。由五点菱形格式,得方程组:
⎪
⎪
⎪
⎪
⎩
⎪
⎪
⎪
⎪
⎨
⎧
=−+++
=−+++
=−+++
=−+++
0)(
1
0)4(
1
0)4(
1
0)4(
1
2,22,11,22,33,2
2
1,22,01,12,23,1
2
2,11,10,21,32,2
2
1,11,00,11,22,1
2
uuuuu
h
uuuuu
h
uuuuu
h
uuuuu
h
(16)
代入边界条件
2,31,33,23,12,01,00,20,1
,,,,,,, uuuuuuuu 的值,(16)式可以化成
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
+
+
+
+
−=
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
−
−
−
−
2,33,2
2,03,1
0,21,3
1,00,1
2,2
1,2
2,1
1,1
4110
1401
1041
0114
uu
uu
uu
uu
u
u
u
u
(17)
解非齐次线性方程组求得
6348.0
1,1
=u , 06.1
2,1
=u , 7985.0
1,2
=
u , 1698.1
2,2
=
u
计算的 Matlab 程序如下:
clc,clear
f1=@(x) 2*log(1+x);f2=@(x) log((1+x).^2+1);
f3=@(y) log(1+y.^2); f4=@(y) log(4+y.^2);
u=zeros(4);m=4;n=4;h=1/3;
u(1,1:m)=feval(f3,0:h:(m-1)*h)';
u(n,1:m)=feval(f4,0:h:(m-1)*h)';
u(1:n,1)=feval(f1,0:h:(n-1)*h);
u(1:n,m)=feval(f2,0:h:(n-1)*h);
b=-[u(2,1)+u(1,2);u(4,2)+u(3,1);u(2,4)+u(1,3);u(3,4)+u(4,3)];
a=[-4 1 1 0;1 -4 0 1;1 0 -4 1;0 1 1 -4];
剩余34页未读,继续阅读
资源评论
Sherry_shiry
- 粉丝: 2
- 资源: 1097
下载权益
C知道特权
VIP文章
课程特权
开通VIP
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功