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的流体质量为在时刻的流体质量为时间内微元内的流体净增量为时间内微元内的流体净增量为由于流动的连续性和质量守恒因此由于流动的连续性和质量守恒因此上面的方程称为一维的连续性方程

偏微分方程讲义

偏微分方程讲义

习题3.4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . §3.5 极坐标系下的分离变量法 . . . . . . . . . . . . . . . . . . . . . . . . . . 3.5.1 3.5.2 由射线和圆弧所界定区域中问题的解法 . . . . . . . . . . . . . . . 周期边界条件问题的解法 . . . . . . . . . . . . . . . . . . . . . .
iv 3.6.3 3.6.4 3.6.5 Legendre方程的级数解、 Legendre多项式 . . . . . . . . . . . . . . Bessel方程的级数解、 Bessel函数 . . . . . . . . . . . . . . . . . . 圆盘中热传导方程的解 . . . . . . . . . . . . . . . . . . . . . . . .
习题1.4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . §1.5 线性偏微分方程的叠加原理,定解问题的适定性 1.5.1 叠加原理 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ii
习题3.5 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . §3.6 高维曲线坐标系下的分离变量法、球函数和柱函数 . . . . . . . . . . . . 3.6.1 3.6.2 Bessel方程和Legendre方程的导出 . . . . . . . . . . . . . . . . . . 二阶线性齐次常微分方程的级数解法 . . . . . . . . . . . . . . . .

数学物理方程ch

数学物理方程ch

cut ku F
令 a2 k , f F , 得
c
c
ut a2u f
(6)
称(6)为三维热传导方程。如果物体内部没有热源
,即 f ≡0,则得齐次热传导方程
ut a2u
(7)
注1:在前面所讨论的热传导问题中,作为特例,如
果所考虑的物体是一根细杆或一块薄板,或者即使不
2u a2 f (x at) a2g(x at) t 2
u f (x at) g(x at) x
2u x2

f
(x at)
g(x at)

2u t 2

a2
2u x2
,
移项即证。
§2 三类典型方程的导出
一、弦振动方程
Q1
t2 t1

V
k
u n
dS dt
根据散度定理得,
Q1
[t2 kudv]dt t2
t1 V
t1
kudvdt
V
(3)
如果物体内有热源,设在单位时间内单位体积
所产生的热量为F(x,y,z,t), 则在 [t1,t2] 内热源放出的
热量为:
Q2
2u t 2
要比
g大得多,所以又可以把g略去。经过这样逐
步略去一些次要的量,抓住主要的量。最后得到 u(x,t)
应近似地满足的方程
2u t 2

a2
2u x 2
(1)
这里 a2 T 。
(1)式称为一维波动方程
如果弦还在横向(位移 u的方向)受到外力的作用。 设在时刻 t弦上 x点处的外力密度为 F(x,t)。仿照前面 的推导,有

第一章 偏微分方程和一阶线性偏微分方程解

第一章 偏微分方程和一阶线性偏微分方程解
3) (弦的一端固定在弹性支承上)(Robin条件)
在高维空间,相应的边界条件为
1)Dirichlet条件: ( 是边界)
2)Neumann条件:
3)Robin条件:
§1.4一阶线性偏微分方程解的特征线方法
对一阶齐次线性偏微分方程

从几何观点看,如果 满足该方程,则由函数 确定的平面上的向量场 ,与方程系数构成的向量场 正交。称由向量场 作为切向所确定的曲线
为方程的特征线。
例如,当 , 为常数,则过任意给定的点 的特征线为直线,方程为

之所以称其为特征线,是因为沿该直线函数 取常数值。以 为常数为例,
特征线上的任意一点可表示为 ,其中 是参数,由此

即 。
利用特征线的该性质,在给定适当的初始或边界条件后就可确定方程的解。
例1求解方程 。
解特征线 ,即 ,沿该直线, 是常数。所以,
解如图,在 内的流体,经过时间 ,一定处于 。所含污染物应相同,即

由此

从而,

【End】
可见偏微分方程是一个至少为两元的函数及其偏导数所满足的方程。
例2(扩散方程)假设水流静止,在 时间内,流经 处的污染物质(不计高阶无穷小)与该处浓度的方向导数(浓度变化)成正比,比例系数为 :

所以,在时间段 内,通过 的污染物为

记 ,则非齐次的波动方程为

【end】
§1.2平面和空间上的偏微分方程
例1(三维空间中的扩散方程)假设污染流体充满三维空间的某区域, 是其密度。任取简单区域 ,相应的边界 。假设,在 时间内,流出 的流与密度关于 处的法向导数成正比,即 ,因此在 流出曲面 的流量为

Ch1第一章

Ch1第一章

..
. ..
. ..
.
. . . . . . . . . . . . .. .. .. .. .. .. .. .. .. .. .. .. ..
. ..
. ..
. ..
. ..
.
引论
数值分析
12/72
例:1/3 + 1/6 = 1/2(符号) 0.333333 + 0.166666 = 0.499999(数值) 无穷大和无穷小 函数的表达——节点的值 f ( x1 ) = y1 , · · · , f ( xn ) = yn 。 插值理论、微分方程的数值解 理论和算法都是围绕这个特点进行的。
. ..
. ..
.
引论
数值分析
13/72
2. 可靠的理论分析:收敛性、稳定性、 误差分析 3. 好的计算复杂性:时间、空间 4. 数值试验(有效性)
..
. ..
. ..
.
. . . . . . . . . . . . .. .. .. .. .. .. .. .. .. .. .. .. ..
. ..
. ..
. ..
. ..
. ..
.
引论
数值分析
13/72
2. 可靠的理论分析:收敛性、稳定性、 误差分析 3. 好的计算复杂性:时间、空间
..
. ..
. ..
.
. . . . . . . . . . . . .. .. .. .. .. .. .. .. .. .. .. .. ..
. ..
. ..
.. . .. . .. . . . . . . . . . . . . . .. .. .. .. .. .. .. .. .. .. .. .. .. . .. . .. . .. . .. .

偏微分内容及作业

偏微分内容及作业

实验和实地勘察结果表明,如果巷道内有障碍物或巷道结构的改变都可更快地导致 DDT 的产生。一旦爆轰产生其强大的破坏能力远超人们的想象。
因此,瓦斯爆炸之所以破坏强度如此大,就是因为产生了爆轰现象,而爆轰现象的产生又 是因为两个自由界面运动并发生重叠的结果,即火焰面追上了冲击波。所以,将瓦斯爆炸控制 在点火燃烧初期或爆燃转变为爆轰之前,具有重大的现实意义。
xk
=
k
·
l n
,k
=
1, · · ·
, n.
通过分析第
k
个质点的力,约翰

伯努利已经证明,如果
yk
是第
k
个质点的位移,则
d2y ( na )2 dt2 = l (yk+1 − 2yk + yk−1), k = 1, · · · , n − 1
其中 a2
=
lT M
,T
是弦中的张力(弦振动时它被当作常数),M
历史源头问题之三:作为实际问题,在工业上为了处理金属;作为科学问题,是企图确定地球 内部的温度
1807 年,傅里叶向巴黎科学院提交了一篇关于热传导的基本论文,这篇论文经拉格朗日, 拉普拉斯和勒让德评审后被拒绝了. 但科学院确实想鼓励傅里叶发展他的思想,所以把热传导 问题确定为将于 1812 年授予高额奖金的课题. 傅里叶在 1811 年提交了修改过的论文,受到上 述诸人和另外一些人评审,得到了奖金,但因受到缺乏严密性的批评而未发表在当时的科学院 的《报告》里. 傅里叶对他所受到的待遇感到愤恨,他继续对热的课题进行研究,在 1822 年发 表了数学的经典文献之一《热的解析理论》(Theoorie analytique de la chaleur),编入了他实际上 未作改动的 1811 年论文的第一部分. 此书是傅里叶思想的主要出处. 两年以后,他成为科学院 的秘书,于是他将 1811 年的论文原封不动地发表在《报告》里.

数学物理方程,偏微分方程答案第一章1-25 课后答案

数学物理方程,偏微分方程答案第一章1-25 课后答案
o
x
u ( x, t )
1 [(h x at ) ( x at ) (h x at ) ( x at )] 2(h x)
u v u u 2v [(h x) 2 (u ) (h x) (h x) 2 (h x)(u 2 ) x x x x x x
( x) s( x)u tt ( ESu x ) x
ww w.
利用微分中值定理,消去 x ,再令 x 0 得
u ∣ x 0 k[u (0, t ) v(t )] x u ( u ) ∣ x 0 f (t ). x
E
x 2 u x 2 2u 3. 试证:圆锥形枢轴的纵振动方程为 E [(1 ) ] (1 ) x h x h t 2
t有
即对任何 x, G(x) C 0 又 G(x)=
(1) 如果初始条件在 x 轴的区间[x 1 ,x 2 ]上发生变化,那末对应的解在区间[ x1 ,
1 1 x C ( x) ( )d x 2 2a 0 2a
x 2 ]的影响区域以外不发生变化;
w.
2u x 2 t x
2

u (t 2 x 2 y 2 ) 2 x x
da

同理 所以 运动方程为:
2

且 T ( x) 的方向总是沿着弦在 x 点处的切线方向。仍以 u ( x, t ) 表示弦上各点在时刻 t 沿垂直于 x 轴 方向的位移,取弦段 ( x, x x), 则弦段两端张力在 u 轴方向的投影分别为
+
x at 1 (h ) ( )d . 2a(h x) x at
即为初值问题的解散。 2.问初始条件 ( x) 与 ( x) 满足怎样的条件时,齐次波动方程初值问题的解仅由右传播波

ch1偏微分方程第一章

ch1偏微分方程第一章

在从事热流动的研究中,1822年发表了《热的解析 理论》,在文章中他提出了三维空间的热传导方程,也 就是一种偏微分方程。他的研究对偏微分方程的发展 的影响是很大的。
1.4 偏微分方程的发展
现在偏微分方程相关理论及其方法已经应用到各个 自然科学,工程技术领域和社会科学领域中。
由于其特殊的地位,偏微分方程现在是数学领域中 最活跃,最核心的领域之一。在菲尔兹奖获得者中与 偏微分方程研究相关的,就有十位左右的数学家。
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 ]内产生的冲量
PDE的解
区域内有m阶连续偏导数的函数.
广义解 (弱解)
PDE的分类 线性PDE 非线性PDE
半线性PDE 拟线性PDE 完全非线性PDE
线性PDE: PDE中关于未知函数及其各阶偏导数都 是线性的。例如:
∑ ∑ n aij ( x1 ,
i , j=1
∂2u , xn ) ∂xi∂x j
+
n
bj ( x1,
AB
x1 x2
下面利用 Fourier 热力学定律和能量守恒定律来建 立热传导方程。
由Fourier 热力学定律,单位时间单位面积内通过A 端的热量为
Qx1

偏微分方程 PDE-Ch1

偏微分方程  PDE-Ch1

(*)
在适当情况下, 方程中描述空间坐标的自变量数目可以减少. 例如当物体是各向同性的均匀细杆时, 如果它的侧面不产生热交 换(即绝热), 且在同一截面上温度的分布是相同的, 则温度函数u 仅与坐标x及时间t有关, 这时得到的就是一维热传导方程
ut a2uxx 0
21
机动 目录 上页 下页 返回 结束
9
机动 目录 上页 下页 返回 结束
《偏微分方程》第一章 绪论 第10页
1.1.4. 线性偏微分方程 如果方程中关于未知函数及其各阶偏导数都是线性的, 则称 它为线性偏微分方程。 例子:
一阶线性偏微分方程 二阶线性偏微分方程
在线性偏微分方程中, 不含有u及它的偏导数的项称为自由项; 当自由项为零时, 称方程为线性齐次方程。 当自由项不为零时, 称方程为线性非齐次方程。
《偏微分方程》第一章 绪论 第19页
类似地可导出二维波动方程和三维波动方程, 它们 的形式分别为
utt a2 (uxx uyy ) f ( x, y, t ) utt a2 (uxx uyy uzz ) f ( x, y, z, t )
二维波动方程可视为薄膜的振动所满足的运动规律, 即在平 面上放置一个框架, 对于固定在该框架上作微小横振动的薄膜上 各点的运动规律. 三维波动方程表示的是声波、电磁波的传播所 满足的规律. 类似地,我们可考虑函数 u u( x1 , x2 ,
Du (ux1 , ux2 ,
, uxn )
则偏微分方程的一般形式为
5
机动 目录 上页 下页 返回 结束
《偏微分方程》第一章 绪论
第 6页
其中是F自变量x,未知函数u及u的有限多个偏导数的已知函数. 例如关系式

【数理方程】91偏微分方程的建立

【数理方程】91偏微分方程的建立
2
(3)利用前面的推导方法,并略去弦本身的 重量,可得弦的, t ) 2 2 t x
2 2
这里a
2
T

,f(x,t)
1

F(x,t)
表示t时刻单位质量的弦在x点处所受的外力。 上面的方程也称为一维 波动方程。
u 2 u a f ( x, t ) 2 2 t x
偏微分方程的特解
即为满足偏微分方程的函数,同时又满足附加 条件的解。
偏微分方程的定解条件
在偏微分方程的问题中,附加条件可分为初 始条件和边值条件,统称为定解条件。
导出数学物理方程的一般方法:
(1) 确定所研究的物理量; (2) 建立适当的坐标系; (3) 划出研究单元,根据物理定律和实 验 资料写出该单元与邻近单元的相互作用, 分析这种相互作用在一个短时间内对所研 究物理量的影响,表达为数学式; (4) 简化整理,得到方程。
u 2 u a 2 2 t x T 2 这里a ,上面的方程称为一维 波动方程。
2 2
2. 受外力作用时弦振动方程的推导
若弦是柔软的线,弦上任何一点处的张力不随 时间而变,总是沿着该点处弦的切线方向。弦 上另外还受一个与振动方向平行的外力影响。 (1)假定在时刻t弦上x点处的外力密度为 F(x,t). (2)考虑 MM 弧段在t时刻的受力情况:
其中ds为小弧段上的质量, gds为小弧段上的重力。
u 可得: T sin T sin gds ds 2 t
2
由 0, 0 可得:
u( x , t ) sin tan 2 x 1 tan
tan
u( x dx, t ) sin tan 2 x 1 tan

常微分方程-精讲自学版 ch1

常微分方程-精讲自学版 ch1

第一章 基本概念当我们描述实际对象的某些特性随时间(或空间)而演变的过程、分析它的变化规律、预测它的未来性态,研究它的控制手段时,通常要建立对象的动态模型.建模时首先要根据建模目的和对问题的具体分析作出简化假设,然后按照对象内在的或可以类比的其他对象的规律列出微分方程,求出方程的解并将结果翻译回实际对象,就可以进行描述、分析、预测或控制了,微分方程的基本问题在于求解和研究解的各种属性.而微积分的产生和发展,与人们求解微分方程的需要有密切的关系.所谓微分方程,就是联系着自变量,未知函数,及其导数在内的方程.物理学,化学,生物学,工程技术和某些社会科学中的大量问题一旦加以精确的数学描述,往往会出现微分方程.在本书的各章中,将举出引导到微分方程的各种例子.一个实际问题只要转化为微分方程,那么问题的解决就有赖于对微分方程的研究.就是在数学本身的理论探讨中,微分方程也是常用的工具.本讲义主要介绍常微分方程的一些最基本的理论和方法.我们在第一章首先给出微分方程及其解的定义,并予以相应的几何解释.实际上,这也是为以后各章进一步的学习所作的必要准备.1.1 微分方程及其解的定义微分方程是一门十分活跃的数学分支,利用数学手段研究自然现象和社会现象,或解决工程技术问题,一般需要对问题建立数学模型,再对它进行分析求解或近似计算,然后按实际的要求对所得的结果做出分析和探讨.数学模型最常见的表达方式,是包含自变量和未知函数的函数方程.在很多情形这类方程还包含未知函数的导数,它们就是微分方程.例如,人口定量分析、生物种群的发展变化、在人文学科领域及交通、环境以及用牛顿第二运动定律列出的质点运动方程都是微分方程,而质点运动方程中的未知函数代表质点的坐标,它们对自变量(时间)的一阶导数和二阶导数分别表示质点的运动速度和加速度.现在,我们给出如下的定义.定义 1.1 凡是联系自变量x ,与这个自变量的未知函数)(x y y =,和它的导数)(x y y '='以及直到n 阶导数)()()(x y y n n =在内的方程0),,,,()(='n y y y x F Λ (1.1)叫作常微分方程,其中导数实际出现的最高阶数n 叫作常微分方程(1.1)的阶.注 1.1 这里F 是一个关于变元(),,,,n x y y y 'L 的给定的已知函数。

  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 =
相关文档
最新文档