ch1偏微分方程第一章

合集下载

偏微分方程

偏微分方程

第一章. 波动方程§1 方程的导出。

定解条件1.细杆(或弹簧)受某种外界原因而产生纵向振动,以u(x,t)表示静止时在x 点处的点在时刻t 离开原来位置的偏移,假设振动过程发生的张力服从虎克定律,试证明),(t x u 满足方程()⎪⎭⎫⎝⎛∂∂∂∂=⎪⎭⎫ ⎝⎛∂∂∂∂x u E x t u x t ρ 其中ρ为杆的密度,E 为杨氏模量。

证:在杆上任取一段,其中两端于静止时的坐标分别为 x 与+x x ∆。

现在计算这段杆在时刻t 的相对伸长。

在时刻t 这段杆两端的坐标分别为:),();,(t x x u x x t x u x ∆++∆++其相对伸长等于 ),()],([)],([t x x u xxt x u x t x x u x x x ∆+=∆∆-+-∆++∆+θ令0→∆x ,取极限得在点x 的相对伸长为x u ),(t x 。

由虎克定律,张力),(t x T 等于),()(),(t x u x E t x T x =其中)(x E 是在点x 的杨氏模量。

设杆的横截面面积为),(x S 则作用在杆段),(x x x ∆+两端的力分别为x u x S x E )()(x u x x S x x E t x )()();,(∆+∆+).,(t x x ∆+于是得运动方程tt u x x s x ⋅∆⋅)()(ρx ESu t x =),(x x x x x ESu x x |)(|)(-∆+∆+利用微分中值定理,消去x ∆,再令0→∆x 得tt u x s x )()(ρx∂∂=x ESu () 若=)(x s 常量,则得22)(tu x ∂∂ρ=))((x u x E x ∂∂∂∂即得所证。

2.在杆纵向振动时,假设(1)端点固定,(2)端点自由,(3)端点固定在弹性支承上,试分别导出这三种情况下所对应的边界条件。

解:(1)杆的两端被固定在l x x ==,0两点则相应的边界条件为.0),(,0),0(==t l u t u(2)若l x =为自由端,则杆在l x =的张力xux E t l T ∂∂=)(),(|l x =等于零,因此相应的边界条件为x u∂∂|lx ==0 同理,若0=x 为自由端,则相应的边界条件为x u∂∂∣00==x (3)若l x =端固定在弹性支承上,而弹性支承固定于某点,且该点离开原来位置的偏移由函数)(t v 给出,则在l x =端支承的伸长为)(),(t v t l u -。

《偏微分方程》课件

《偏微分方程》课件

非线性偏微 分方程:方 程中含有偏 导数,且偏 导数项的系 数不是常数
椭圆型偏微 分方程:方 程中只含有 二阶偏导数, 且二阶偏导 数项的系数 是常数
抛物型偏微 分方程:方 程中只含有 二阶偏导数, 且二阶偏导 数项的系数 不是常数
双曲型偏微 分方程:方 程中只含有 二阶偏导数, 且二阶偏导 数项的系数 是常数,但 方程的解不 是实数
边界条件:确定求解区域和边界条件,如Dirichlet边界条件、 Neumann边界条件等
初值条件:确定求解区域的初值条件,如Cauchy问题、初边值问题等
稳定性和收敛性:分析求解方法的稳定性和收敛性,确保解的准确性和 可靠性
应用实例:通过具体实例,展示求解方法的应用和效果
课件结构
课件目录
偏微分方程的应用
物理领域:描述 流体力学、热力 学、电磁学等现 象
工程领域:解决 结构力学、材料 力学、电子工程 等问题
生物领域:模拟 生物系统的生长、 扩散、反应等过 程
经济领域:用于 金融、经济模型、 风险管理等方面
偏微分方程的求解方法
分析法:通过分析方程的性质,寻找解的性质和形式
数值法:通过数值计算,求解偏微分方程的数值解
偏微分方程的求解方法:展示偏微分方程的求解方法,如分离变量法、积分因子法等
公式素材
偏微分方程的 定义和性质
偏微分方程的 应用实例
偏微分方程的 求解方法
偏微分方程的 扩展和研究进

动画素材
动画类型:2D动画、3D动画、Flash动画等 动画内容:偏微分方程的求解过程、应用实例等 动画风格:简洁明了、生动有趣、易于理解 动画时长:根据课件内容需要,控制在5-10分钟以内
偏微分方程PPT课件

偏微分方程讲义

偏微分方程讲义

G(x) = φ(x/2) − F (0), F (0) + G(0) = φ(0) = ψ (0). ( ) ( ) x − at x + at ⇒ u(x, t) = ψ +φ − φ(0). 2 2
例 2.4 对非齐次波动方程的初值问题 (2.5)、(2.6), 证明: 当 f (x, t) 不变时, (1) 如果初始条件在 x 轴的区间 [x1 , x2 ] 上发生变化, 那么对应的解在区间 [x1 , x2 ] 的影响区域外不发生变化; (2) 在 x 轴区间 [x1 , x2 ] 上所给的初始条件唯一确定区间 [x1 , x2 ] 的决定区域中解的 数值. 解: 弄清影响区域、决定区域的定义. 例 2.5 求解 utt − a2 uxx = 0, x > 0, t > 0, u|t=0 = φ(x), ut |t=0 = 0, ux − kut |x=0 = 0,
其中 k 为正常数. 齐海涛 htqi2008@ 4
山东大学威海分校
2
达朗贝尔公式、波的传播
解: 波动方程的通解为 u = F (x − at) + G(x + at), 由初始条件得 F (x) + G(x) = φ(x), −aF ′ (x) + aG′ (x) = 0
1 C 1 C F (x) − G(x) = C, F (x) = φ(x) + , G(x) = φ(x) − , 2 2 2 2
1 其中 C = F (0) − G(0). 由于 x + at ≥ 0, G(x + at) = 2 φ(x + at) − C 2 . 当 x − at ≥ 0 1 C 1 时, F (x − at) = 2 φ(x − at) + 2 . 此时 u(x, t) = 2 [φ(x + at) + φ(x − at)]. 当 x − at < 0 时, 由边界条件知

一阶偏微分方程教程

一阶偏微分方程教程

N (t) 0 p(a,t)da
18
若不考虑死亡,则在时刻 t+t,年龄在[a, a+a] 中的人口数量 p(a, t+t)a,应等于在时刻 t,年龄 在区间[a−t, a+a−t]中的人口数量p(a−t, t)a, 即
p(a,t t) p(a t,t)
因此 p(a, t)应满足
dx
dy du
1 u x y 1 2
首次积分为 u 2 y, 2 u x y y
于是原方程的隐式通解为
u 2y, 2 u x y y 0
其中 为任意二元连续可微函数。
16
例5. 求解hy问题
u
u x

xz u y

xy u z
0

u yy0 f (x, z)
11
解:特征方程组为 dx dy dz yz xz xy
首次积分为 x2 y2, x2 z2
于是原方程的通解为 u x2 y2, x2 z2 ,其中
为任意二元连续可微函数。
研究的数据包括50根圆柱组织样本中每一根所含 药物的测量值(见表1、表2及图1)。每一圆柱的长度 为0.76mm,直径为0.66mm。这些平行圆柱的中心 位于1mm×0.76mm×1mm的网格点上。因此,圆

a 0, t 0

p(a,
0)

p0 (a),
a0
(4)


p(0,
t
)


(a,t, N (t)) p(a,t)da,
0
t 0


N (t) 0 p(a,t)da, t 0

偏微分方程讲义

偏微分方程讲义

2 达朗贝尔公式、波的传播
例 2.1 证明方程 [ ] ∂ ( x )2 ∂u 1 ( x )2 ∂ 2 u 1− = 2 1− ∂x h ∂x a h ∂t2
(h > 0 常数)的通解可以写成 u= F (x − at) + G(x + at) , h−x
其中 F, G 为任意的具有二阶连续导数的单变量函数, 并由此求解它的初值问题: t=0: u = φ(x), ∂u = ψ (x). ∂t
其中 k 为正常数. 齐海涛 htqi2008@ 4
山东大学威海分校
2
达朗贝尔公式、波的传播
解: 波动方程的通解为 u = F (x − at) + G(x + at), 由初始条件得 F (x) + G(x) = φ(x), −aF ′ (x) + aG′ (x) = 0
1 C 1 C F (x) − G(x) = C, F (x) = φ(x) + , G(x) = φ(x) − , 2 2 2 2
G(x) = φ(x/2) − F (0), F (0) + G(0) = φ(0) = ψ (0). ( ) ( ) x − at x + at ⇒ u(x, t) = ψ +φ − φ(0). 2 2
例 2.4 对非齐次波动方程的初值问题 (2.5)、(2.6), 证明: 当 f (x, t) 不变时, (1) 如果初始条件在 x 轴的区间 [x1 , x2 ] 上发生变化, 那么对应的解在区间 [x1 , x2 ] 的影响区域外不发生变化; (2) 在 x 轴区间 [x1 , x2 ] 上所给的初始条件唯一确定区间 [x1 , x2 ] 的决定区域中解的 数值. 解: 弄清影响区域、决定区域的定义. 例 2.5 求解 utt − a2 uxx = 0, x > 0, t > 0, u|t=0 = φ(x), ut |t=0 = 0, ux − kut |x=0 = 0,

一阶偏微分方程讲义

一阶偏微分方程讲义

偏微分方程(Partial Differential Equations)許多物理規律、物理過程和物理狀態都可以用微分方程描述。

當物理過程和狀態只由一個因素決定時,往往提出常微分方程。

例如質點的運動,通過解常微分方程就能得到質點的運動規律。

例:()()()()mu t cu t ku t P t++=(質點運動方程式)當物理問題由多個因素決定時,就會涉及到偏微分方程,偏微分方程為應用數學中重要的課題之一,物理問題之數學模式與偏微分方程式有關,許多數學理論與方法的發展往往肇因於求解偏微分方程。

例:222u uat x∂∂=∂∂(熱傳導方程式)22222u uat x∂∂=∂∂(波動方程式)2222u ux y∂∂+=∂∂(拉普拉斯方程式)u uut x∂∂+=∂∂(衝擊波方程式)33u u uut x xσ∂∂∂++=∂∂∂(KdV方程式)1. 偏微分方程的定義與解設()12,,,n u x x x =為自變數12,,,n x x x 之函數,任何包含其偏導數之關係式21211,,,,,,,0n n u u f x x x x x x ⎛⎫∂∂= ⎪∂∂∂⎝⎭(1)稱為偏微分方程式(簡稱P.D.E.)。

基本名詞:(1) 階數(order):P.D.E.中所含最高階偏導數之階數。

(2) 線性(linear):P.D.E.中,其未知函數以及其偏導數均滿足 (i)次數均為一次。

(ii)無互相的乘項。

(iii)無非線性函數。

則稱為線性P.D.E.。

(3) 擬線性(quasi-linear):P.D.E.中,其最高階的偏導數之次數為1次,且彼此無互乘項,則稱為擬線性P.D.E.。

(4) 非線性(non-linear):若P.D.E.不為線性或擬線性,則稱為非線性P.D.E.。

解之分類:(1) 通解(general solution):滿足P.D.E.且包含任意函數之解。

(2) 全解(complete solution):滿足P.D.E.且包含任意常數之解。

偏微分方程课件-Ch1

偏微分方程课件-Ch1
《偏微分方程》第一章 绪论
第一章 绪论
1.1
《偏微分方程》第一章 绪论
在偏微分方程中最高阶偏微商的阶数叫做该偏微分方 程的阶。
在偏微分方程组中最高阶偏微商的阶数叫做该偏微分 方程组的阶。
如果一个函数在其自变量 (x1, x2 , , xn )的某变化范 围内连续并且具有方程(方程组)的一切连续偏微商 将它代入方程后使其成为恒等式,则称该函数是方程
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
有三种方程的划分(1.3.16)可知,在椭圆型区域内不存 在实的特征方向;在双曲型区域内存在两族实的特征方 向;而在抛物型的点上仅有一个实的特征方向。因此, 方程双曲型区域被两族实特征曲线网覆盖,在抛物型的 点集被一簇实特征线网覆盖
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
考虑方程
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
《偏微分方程》第一章 绪论
考虑(1.3.8)的线性主部
《偏微分方程》第一章 绪论
支覆盖了整个下半平面4,并且每个分支都与 x 轴相切。
(3)在椭圆型区域 y 0 中特征方程的其中的一个复解为
x 2i y c
去实部和虚部作变换
x 2 y
《偏微分方程》第一章 绪论
经过计算便得到方程在上半平面的标准型
u
u
1

u

0,
y

0.
《偏微分方程》第一章 绪论

1 ch1 数理方程第一章1

1 ch1 数理方程第一章1

∂u ( x2 , t ) Qx2 = − k ∇u • n( x2 ) = − k ∂x
24
• 在 dt 时段内通过微元的两端流入的热量
∂u ( x2 , t ) ∂u ( x1 , t ) dQ1 = −(Qx1 + Qx2 ) dt = k ( ) dt − ∂x ∂x x2 2 ∂ u ( x, t ) = k∫ dxdt 2 ∂x x1
i =1

7
数学物理方程的导出
• 波动方程
– 均匀弦的微小横振动方程 – 推广
• 扩散方程
– 一维热传导方程 – 推广
• 稳定场方程
8
• 弦振动方程
• 弦的特点:匀、细、软、紧的一根弹性细线。 • 振动特性:微小的、横向振动:在一个平面内弦上各点
的运动方向垂直于最初的平衡位置. “微小的”是指弦上各 点的位移与弦的长度相比很小, 弦的纵向伸长可以忽略不 计
数理方程的基本概念
一. 偏微分方程的基本概念
偏微分方程:凡含有多元未知函数及未知函数关于自变量 的偏导数的等式。 自变量 1 2 n
x = (x , x ,
,x )
u ( x) = u ( x1 , x2 ,
, xn )
未知函数
1
偏微分方程: Partial Differential Equation, 简写 为: PDE
在流体柱上任意取一微元在流体柱上任意取一微元处两个截面处两个截面任取一个时段任取一个时段流体在流体在这段时间间隔内从x这段时间间隔内从x处截面流入的质量为处截面流入的质量为处截面流出的质量为处截面流出的质量为所以流体在所以流体在时间间隔内微元中流体净增量为时间间隔内微元中流体净增量为由于在时刻t的流体质量为在时刻的流体质量为由于在时刻t的流体质量为在时刻的流体质量为时间内微元内的流体净增量为时间内微元内的流体净增量为由于流动的连续性和质量守恒因此由于流动的连续性和质量守恒因此上面的方程称为一维的连续性方程
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
AB
x1 x2
下面利用 Fourier 热力学定律和能量守恒定律来建 立热传导方程。
由Fourier 热力学定律,单位时间单位面积内通过A 端的热量为
Qx1
=
−k ∂u ∂n
x= x1
=
−k∇u1, t ) ∂x
单位时间单位面积内通过 B 端的热量为
Qx2
=
−k
∂u ∂n
各点的运动方向垂直于最初的平衡位置. “微小的”是 指弦上各点的位移与弦的长度相比很小, 弦的纵向伸 长可以忽略不计. y 考虑一根拉紧的长为l 的弦,线密度ρ , 以弦的平衡 位置所在直线为 x 轴,并以弦的左端点为坐标原 点,则右端点的坐标为 l。求它在平衡位置附近做微 小的横向振动的规律。
动量守恒律
在从事热流动的研究中,1822年发表了《热的解析 理论》,在文章中他提出了三维空间的热传导方程,也 就是一种偏微分方程。他的研究对偏微分方程的发展 的影响是很大的。
1.4 偏微分方程的发展
现在偏微分方程相关理论及其方法已经应用到各个 自然科学,工程技术领域和社会科学领域中。
由于其特殊的地位,偏微分方程现在是数学领域中 最活跃,最核心的领域之一。在菲尔兹奖获得者中与 偏微分方程研究相关的,就有十位左右的数学家。
∫ ∫t2
b
F ( x,t)dxdt
t1 a
由动量守恒律有
∫b a
ρ
(
ut
(
x,
t2
)

ut
(
x
,
t1
)
)
dx
∫ ( ) ∫ ∫ =
T t2
t1
ux (b, t) − ux (a,t)
dt +
t2 t1
b
F ( x,t)dxdt
a
由Newton-Leibnitz 公式,有
ρutt ( x,t) = Tuxx ( x,t) + F ( x,t),
(4). (ut )2 + (ux )2 = u2
拟线性PDE: 在非线性方程中, 如果关于未知函数的 所有最高阶偏导数是线性的. 例如:
(1 + u2y )uxx − 2uxuyuxy + (1 + ux2 )uyy = 0
utt − uxx + u3 = 0
ut + uux + uxxx = 0 拟线性PDE的一般形式:
y T 在[ t1, t2 ]内产生的冲量:
∫ ( ) T t2 t1
ux (b,t) − ux (a,t)
dt
y [ a, b ]的动量变化为:
∫b a
ρ
(
ut
(
x
,
t2
)

ut
(
x,
t1
))
dx
y 在点 x 处 t 时刻外力密度为F(x, t), 则F(x, t)在微弦段
[ a, b ]上[ t1, t2 ]内产生的冲量
utt − uxx + u3 = 0 ut + uux + uxxx = 0
半线性PDE的一般形式:
∑ aα ( x)Dα u + b( x, u, Du, Dm−1u) = 0,
|α |= m
其中 Dα
=
∂α ∂x1α1 ∂xnαn
, Dk u = (Dα u :| α
|=
k ).
1.5 叠加原理
即弦振动方程 (又称为一维波动方程) utt ( x, t ) − a2uxx ( x, t ) = f ( x, t ),
其中 f ( x,t ) = F ( x,t ) / ρ 表示单位质量所受的力。
二维、三维波动方程的形式分别为
( ) utt = a2 uxx + uyy + f ( x, y, t), ( ) utt = a2 uxx + uyy + uzz + f ( x, y, z, t ).
偏微分方程
主讲:赵志红
本课程要求
y 平时成绩占30分,期末考试占70分; y 每周周一交作业,由课代表课前收齐放到讲台上 y 答疑安排:双周周一晚7:00-9:00,学楼103 y E-mail: zzh_math@
第一章 方程的导出及定解问题 的提法
y 基本概念 y 几个经典方程 y 定解问题
以二阶线性偏微分方程的为例来说明叠加原理.
偏微分方程可用偏微分算符来表示.
一般含n个自变量的二阶线性偏微分方程可写为以 下形式:
∑ ∑ n aij
i, j=1
∂2u ∂xi∂x j
+
n
bi
i =1
∂u ∂xi
+ cu =
f
∑ ∑ 引入以下算符
L=
n
aij
i , j=1
∂2 ∂xi∂x j
+
n
bi
i =1
和欧拉同时代的瑞士数学家丹尼尔·贝努利也研究了数 学物理方面的问题,提出了解弹性系振动问题的一般方 法,对偏微分方程的发展起了比较大的影响。拉格朗日 也讨论了一阶偏微分方程,丰富了这门学科的内容.
偏微分方程得到迅速发展是在十九世纪,那时候,数 学物理问题的研究繁荣起来了,许多数学家都对数学物 理问题的解决做出了贡献。这里应该提一提法国数学家 傅立叶,他年轻的时候就是一个出色的数学学者。
x= x2
=
−k∇uin( x2 )
=
−k
∂u( x2 , t ) ∂x
在 dt 时段内通过A, B的两端流入的热量
dQ1
=
(Qx1

Qx2
)dt
=
k( ∂u( x2 ,t ) ∂x
小结
叠加原理使得以后在使用分离变量法时能够将分离 变量法得到的线性无关的解叠加在一起, 然后去构造 原问题的解.
二、几个经典方程
y 波动方程 y 均匀弦的微小横振动方程 y 推广
y 热传导方程 y 一维热传导方程 y 推广
y 稳定场方程
2.1 弦振动方程 y 弦的特点:匀、细、软、紧的一根弹性细线。 y 振动特性:微小的、横向振动:在一个平面内弦上
1.3 偏微分方程的起源
微积分方程这门学科产生于十八世纪,欧拉在他的 著作中最早提出了弦振动的二阶方程,随后不久, 法国数学家达朗贝尔也在他的著作《论动力学》中 提出了特殊的偏微分方程。这些著作当时没有引起多 大注意。
1746年,达朗贝尔在他的论文《张紧的弦振动时形成 的曲线的研究》中,提议证明无穷多种和正弦曲线不同 的曲线是振动的模式。这样就由对弦振动的研究开创了 偏微分方程这门学科。
3. 两个非齐次方程的解的线性组合为一个新的非齐次方 程的解,新方程的自由项为原方程自由项的同样组合。
若 L u1 = f1 , L u2 = f2, 则:L (au1+ bu2)= af1 + bf2
线性方程的叠加原理
设ui (i = 1, 2, 3, ) 满足方程 Lui = fi (i = 1, 2, 3, ), 的解,ci (i = 1, 2, 3, ) 为常数,而级数
总结
建立数学物理方程是一个辩证分析的过程。由于客观 事物的复杂性,要求对所研究的对象能够抓住事物发展 的主要因素,摈弃次要因素,使问题得到适度的简化。
在上面的推导过程中,我们作了一些假设。 ¾弦是完全柔软的,张力才会沿着弦的切线方向;
¾弦的横振动是很小的, 所以才可用sinθ 代替 tanθ .
¾弦的纵向伸长可以忽略不计,不然各点张力的不同,
千禧年大奖难题,又称世界七大数学难题,其中之 一就是纳维叶-斯托克斯(Navier-Stokes)方程的存在 性与光滑性。
偏微分方程: Partial Differential Equation (PDE)
PDE的阶 所含有的未知函数最高阶导数的阶数. m = m1 + m2 + + mn
古典解 是指满足方程,并且在所考虑的
j=1
主部
,
∂u xn ) ∂x j
+
c( x1,
, xn )u = f ( x1,
, xn )
其中aij ,bj ,c, f 是给定的函数.
线性PDE的主部: 具有最高阶数偏导数组成的部分.
常系数线性PDE: 齐次线性PDE:
系数aij ,bj ,c均为常数. 否则称为变系数的PDE. f ≡ 0. 否则称为非齐次的.
PDE的解
区域内有m阶连续偏导数的函数.
广义解 (弱解)
PDE的分类 线性PDE 非线性PDE
半线性PDE 拟线性PDE 完全非线性PDE
线性PDE: PDE中关于未知函数及其各阶偏导数都 是线性的。例如:
∑ ∑ n aij ( x1 ,
i , j=1
∂2u , xn ) ∂xi∂x j
+
n
bj ( x1,

∑ u = ciui i =1
一致收敛且能够逐项微分两次,则u满足方程
Lu = f ,

∑ 此处要求级数 f = ci fi 一致收敛。 i =1
特别是,如果 ui (i = 1,2,3, ) 是二阶线性齐次方程 Lu = 0

∑ 的解, 则只要 u = ciui 一致收敛,且可以逐项微分两次, i =1 则u一定也是此方程的解.
t2时的动量- t1时的动量=[t1, t2]内外力产生的冲量
取弦的平衡位置为ox 轴,运动平面为 xou
在时刻 t ,弦线在 x 点的位移为 u(x, t)
o
a, b两点受力图示
利用微元法建立方程:
在任一时刻 t,任取一小段弦 [a,b], 它弧长为
∫ ∫ s =
相关文档
最新文档