热传导方程的初值问题

§2热传导方程的初值问题

一维热传导方程的初值问题(或Cauchy 问题)

??

???+∞<<∞-=>+∞<<∞-=??-??x x x u t x t x f x u a t

u ),()0,(0

,),,(2

2

2? （）

偏导数的多种记号xx x t u x

u

u x u u t u =??=??=??22,,. 问题也可记为

??

?+∞

<<∞-=>+∞<<∞-=-x x x u t x t x f u a u xx t ),()0,(0

,,),(2?.

Fourier 变换

我们将用Fourier 变换法求解热传导方程的柯西问题.为此我们将着重介绍Fourier 变换的基本知识.Fourier 变换在许多学科中是重要使用工具. 可积函数,设)(x f f =是定义在),(+∞-∞上的函数, 且对任意A B <，()f x 在[,]A B 上

可积，若积分

?

+∞

∞

-dx x f )(收敛,则称)(x f 在),(+∞-∞上绝对可积。

将),(+∞-∞上绝对可积函数形成的集合记为),(1

+∞-∞L 或),(+∞-∞L ， 即{

}

∞<=+∞-∞=+∞-∞?

+∞

∞

-dx x f f L L )(|

),(),(1

，称为可积函数空间.

连续函数空间: ),(+∞-∞上全体连续函数构成的集合,记为),(+∞-∞C ,

{}上连续在),(|),(+∞-∞=+∞-∞f f C ， {}上连续在),(,|),(1+∞-∞'=+∞-∞f f f C 。

定义 若),(+∞-∞∈L f ,那么积分

),(?)(21

λπ

λf dx e x f x i =?

+∞

∞

--

有意义,称为Fourier 变换, )(?

λf 称为)(x f 的Fourier 变式(或Fourier 变换的象). ?

+∞

∞

--=

=dx e x f f Ff x i λπ

λλ)(21)(?)(

定理 (Fourier 积分定理)若),(),(1

+∞-∞?+∞-∞∈C L f ,那么我们有

),()(?21lim

x f d e f N

N

x i N =?

+-∞

→λλπ

λ

公式称为反演公式.左端的积分表示取Cauchy 主值.

通常将由积分

)()(21

x g d e g x i ∨+∞

∞

-=?

λλπ

λ所定义的变换称为Fourier 逆变换.

因此亦可写成

()

f f =∨

?

即一个属于),(),(1

+∞-∞?+∞-∞C L 的函数作了一次Fourier 变换以后,再接着作一次Fourier 逆变换,就回到这个函数本身.

在应用科学中经常把)(?

λf 称为)(x f 的频谱.Fourier 变换的重要性亦远远超出求解偏微分方程的范围,它在其它应用科学中,如信息论,无线电技术等学科中都有着极为广阔的应

用.它是近代科学技术中得到广泛应用的重要数学工具.

定理的证明在经典书中都能查到(如姜礼尚,陈亚浙,<<数学物理方程讲义>>)

定理 设),(+∞-∞∈L f ，?

+∞

∞

--=dx e x f f

x i λπ

λ)(21

)(?，则)(?

λf 是有界连续函数,且 .0)(?lim =∞

→λλf

在运用Fourier 变换求解定解问题以前,我们先来介绍一些Fourier 变换的性质.

Fourier 变换的性质： 1.（线性性质） 若.2,1,),

,(=∈+∞-∞∈j C L f j j α则

(),??22112211f f f f αααα+=+∧

2.（微商性质）

若),,(),()(),(+∞-∞?+∞-∞∈'L C x f x f 则.?

f i dx df λ=??

?

??∧

证明 由假设),,(),()(),(+∞-∞?+∞-∞∈'L C x f x f 故0)(lim =∞

→x f x ,

事实上由),()(+∞-∞∈'C x f ,则dt t f f x f x

?'+=0

)()0()(,

因为),()(+∞-∞∈'L x f ,故有

?

±∞

±±∞

→'+==0

)()0()(lim dt t f f a x f x

又因),()(+∞-∞∈L x f ,必有0=±a .

由0)(lim =∞

→x f x ,利用分部积分公式

?

∞

+∞

--∧

'=

??

?

??dx e x f dx df x i λπ)(21

??

????--=?

+∞

∞

--∞

+∞

--dx e i x f e x f x i x

i ))(()(21λλλπ

).(?)(2λλπ

λλf i dx e x f i x i ==?

+∞

∞

--

附注 这个性质说明微商运算经Fourier 变换转化为乘积运算,因此利用Fourier 变换可把常

系数微分方程简化为函数方程,或把偏微分方程简化为常微分方程,正是由于这个原因,Fourier 变换成为解微分方程的重要工具. 3.（乘多项式）

若),()(),(+∞-∞∈L x xf x f 则有[])(?

)(λλ

f d d i

x xf =∧

. 证明 由于),()(),(+∞-∞∈L x xf x f ,故)(?

λf 是λ的连续可微函数,且有 []∧

+∞

∞

---=-=?

)()())((21

)(?x xf i dx e ix x f f d d x i λπ

λλ

附注 作为性质2,3的推论,若),,(),()(),(),()

(+∞-∞?+∞-∞∈'L C x f

x f x f m Λ则 ())1(,)(?≥=???

? ??∧

m f i dx f

d m m m λλ 若),,()(),(),(+∞-∞∈L x f x x xf x f m

Λ则

[

]

)1(,)(?

)(≥=∧

m f d d i x f x m

m m

m

λλ

4．（平移性质）

若),,()(+∞-∞∈L x f 则

[])1()(?)(≥=--∧

m f e a x f a i λλ

证明

[])

(?)(21)(21

)()(λπ

π

λλλf e dy e y f y

a x dx e a x f a x f a i a y i x i -∞

+∞

-+-+∞

∞

--∧==--=

-?

?

5.（伸缩性质）

若),,()(+∞-∞∈L x f 则

[])0(,)(?1)(≠=

∧

k k

f k kx f λ

证明 无妨设,0

[]

)(?11)(1211)(21)(21)(k

f k dy k

e y

f k dy k e

y f y kx dx

e kx

f kx f k

y

i k

y

i x i λππ

π

λ

λ

λ=???

??-=

==

??

?

∞

+∞--∞

-∞+-+∞

∞

--∧

6.（对称性质）

若),,()(+∞-∞∈L x f 则 ,)(?

)(λλ-=∨

f f 证明?

+∞

∞

-∨

=dx e x f f x i λπ

λ)(21

)(?

+∞

∞

---=

dx

e x

f x i )()(21λπ

.)(?λ-=f

7.（卷积定理）

若),,()(),(+∞-∞∈L x g x f ?

+∞

∞

--=

*dt t g t x f x g f )()()(称为f 与g 的卷积,

则),()(+∞-∞∈*L x g f ,且有()).(?)(?2)(λλπλg

f g f =*∧

证明 由积分交换次序定理

??

?

+∞∞

-+∞

∞

-+∞

∞

--=*dx dt t g t x f dx x g f |)()(|)(??+∞

∞-+∞

∞-??

? ??-≤dt dx t g t x f )()(?

?+∞

∞

-+∞

∞-??

? ??-=dt dx t x f t g )()(??+∞∞-+∞∞-?=dt t g dx x f )()( 故),()(+∞-∞∈*L x g f ,又由积分交换次序定理

()()()().??2)(21)(212)()(21)()(21)(λλππππππ

λλλλλλg

f dy e y f dt e t

g dx e t x f dt e t g dt t g t x f dx e g f y

i t i t x i t

i x

i =??=-=-=*??????∞+∞-∞+∞

---∞+∞-∞+∞

----+∞∞-+∞∞

--∧

下面作为例子,我们根据Fourier 变换的定义与性质求一些具体函数的Fourier 变换.

例1 设 ????

?>≤=A

x A x x f ,0,

1)(1,(其中常数0>A ).

求)(?

1λf .

解 由定义

?

?

----=

=

A

A

x i A

A

x i dx e dx e x f f λλπ

π

λ21)(21

)(?11

A

A

x i e i --??

? ??-=

λλπ121λλπA sin 2=. 例2 设??

?<≥=-0

,

00

,

)(2x x e x f x , 求)(?

2λf . ?

+∞

--=

221)(?dx e

e f x

i x λπ

λ?

+∞

+-=

)1(21dx e x i λπ

∞

++-??

? ??+-

=

0)1(1121x i e i λλπλ

πi +=11

21

.

例3 设,)(3x

e

x f -=求)(?3λf

?

+∞

∞

---=

dx e e

f x i x

λπ

λ21)(?3????

??+=

??∞--+∞+-0)1(0)1(21dx e dx e x

i x i λλπ ??

?

??-++=

λλπi i 1111212

12

21λ

π+=

. 例4 设,)(2

4x e x f -=求)(?4λf

?

+∞

∞

---=

dx e

e

f x

i x λπ

λ2

21)(?4?

∞

+∞

---'

??

? ??-=

dx e i e

x i x λλπ

1212

??

????-??? ??-=

?∞+∞---∞+∞---dx e xe i e e i x i x x x i λλλλπ

222121

[]∧

-=

2

2x xe i

λ

)(?

24λλ

λf d d -

= ， 上面最后一个等式应用了性质3. 因为)(?4λf 作为λ的函数适合下面常微分方程初值问题：

???

??

??==-=?∞+∞--2121)0(?,)(?

2)(?2444dx e f f d f d x πλλλλ, 解之得

4

42

2

1)(?λλ-

=e

f .

例5 设,)(2

5Ax e x f -=（0>A ）,求)(?5λf .

由性质5

()()

A

e

A A f A x A f x f f 444552

21)(?1)()()(?λλλ-

∧∧====.

例6 ),(

)(462

2

B

x f e

e

x f B x B

x ===???

? ??--(0>B )

()4

4662

2

)/1(

?/11()(?λλ

λB e

B B

f B

x f f -

∨=

=

=.

()()?+∞

∞

-∨*=

*λλπ

λd e g f x g f x

i )(21)( ??+∞

∞-+∞

∞-??? ??-=

λλπ

λd e dy y g y f x i )()(21

dy d e y g y f x i ??+∞

∞-+∞

∞-??? ??-=

λλπ

λ)()(21

dy d e y f e y g x

y i iyx ??+∞

∞-+∞

∞--??

? ??-=

λλπ

λ)()()(21 )()(2x g x f ∨∨=π,

()()g f g

f

g f ?==

???

? ??*∨

∨

∨

∧∧??22121ππ

π,

于是()∧∧∧

*=

?g f g f π

21,

因为()g

f g f ??2?=

*∧π, 所以()

()[]

g f g f g f *=

*=

?∨

∧∨

π

π

2121??

.

最后我们简单地介绍一些有关多维Fourier 变换的基本知识

定义 设),(),,,()(21n

n R L x x x f x f ∈=Λ那么积分

())(?)(21

λπλf dx e x f n

R

x i n

=?

?-,

有意义,称为)(x f 的Fourier 变换,)(?

λf 称为)(x f 的Fourier 变式.

定理（反演公式）若)()()(1n

n R L R C x f ?∈,则有

(

))()(?21lim

x f d e f

N

x i n

N =?

≤?∞

→λλλλπ

. ()?

?∨=

n

R

x i n

d e g x g λλπλ)(21

)(称为)(λg 的Fourier 逆变换.

定理表明()()f f f f ＝∧

∨∨

=,?容易证明关于一维Fourier 变换的性质1—7对于多维Fourier

变换依然成立.根据上面Fourier 变换的定义,我们还有下面的结论： 8. 若),()()()(2211n n x f x f x f x f Λ=其中),,()(+∞-∞∈L x f i i 则有

)(?)(?1

i

i n

i f f λλ=∏= （） 利用这一性质,我们可求出函数2

2

1

)(i Ax n

i x

A e e

x f -=-∏==的Fourier 变式.

事实上(

)

A

Ax i i e

A

e

42

2

21λ-

∧

-=

,

()

(

)

A

n

A

n

i Ax n

i Ax n

i e

A

e A

e e

f i i

i 441

112

2

2

22121

)(?λ

λλ-

-

=∧

-=∧

-==

∏

=∏=??

? ??∏=.

Poisson 公式

在这一小节中我们应用Fourier 变换解初值问题

??

???+∞<<∞-=>+∞<<∞-=??-??x x x u t x t x f x u a t

u ),()0,(0

,),,(2

2

2? （）

在方程（）两边关于变量x 作Fourier 变换,

?

+∞

∞

--=dx e t x u t u

x i λπ

λ),(21),(? ,

利用性质1和性质2,得到

?????==+=),(??),,(???02

2λ?

λλt u t f u

a dt u d 其中 ?

+∞

∞

--=dx e

t x u t u

x

i λπ

λ),(21),(?,?

+∞

∞

--=dx e x x i λ?π

λ?

)(21

)(?

[]∧=),(),(?t x f t f λ.

解之得

?---+=t t a t a d e f e t u 0

)(2

222),(??),(?ττλ?λτλλ,

现在对上式两边求反演,由反演公式,得

(

)

(

)

?

∨

--∨

-+=t

t a t

a d e f e t x u 0

)(2222),(??),(ττλ?

τλλ （）

由()

,21

42

2

A

Ax e A

e

i

λ-

∧

-=

取t a A 241

=则t

a x t a e t

a e 2222241211λ-∧

-=???

? ??, 即t a x t a e

e t a 222241

21λ-∧

-=???

? ??, 令2

2

4121),(x t

a e

t

a t x g -

=

,[]t a e t x g 2

2),(λ-∧

=,

从而有(

)

()g g e t

a *21

???22?π

?

?λ==∨

∨

- ?

+∞

∞

--=

ξξξ?π

d x g )()(21

?

∞

+∞

---=

ξξ?πξd t a

t

a x 224)()(21 （）

同理我们有

(

)(

)

g f t g f e

f t a *21

),(?),(?),(?)(22π

τλτλτλτλ=

-=∨

∨

-- ?∞+∞

----

-=ξτξτπτξd e f t a t a x )

(4)(22),()(21

（）

于是得

???

∞

+∞

----

∞

+∞

---

-+=

ξτπτξτξξ?πτξξd e

t a f d d t a

t x u t a x t t

a x )

(4)(0

4)(2222)

(21

)

,()

(21),(

在一定条件下,可以证明上述表达式的函数是方程问题的解. 定理 若),()(+∞-∞∈C x ?,且)(x ?有界,则?

∞

+∞

---

=

ξξ?πξd e

t a

t x u t

a x 224)()(21),(

在),0(+∞?R 上连续,且在),0(+∞?R 上具有任意阶的连续偏导数，

),(t x u 是问题??

???+∞<<∞-=>+∞<<∞-=??-??x x x u t x x

u a t u ),()0,(0

,,022

2?的解,

即),(t x u 满足方程和)(),(lim 00

x t x u x x t ?=→→+

. ?

∞

+∞

---

=

ξξ?πξd e

t a

t x u t

a x 224)()(21),(

?

+∞

∞

--+-=ηη?π

ξηηd e t a x t

a x 2

)2(1

2/)(

特别说明：当)(x ?连续,)(x ?是某些无界函数时,),(t x u 的表达式亦是解（)(x ?无界时，也可以是解）.

例1 求解?????=??=??=x

u

x u a

t u t sin ,022

2

解 1、直接观察x e t x u t a sin ),(2

-=是解. 2、?

+∞

∞

--+=

ηη?π

ηd e t a x t x u 2

)2(1

),(

?

+∞

∞

--+=

ηηπ

ηd e t a x 2

)2sin(1

()

?

+∞

∞

---+=

ηηηπηηd e t a x e t a

x 2

22sin cos 2cos sin 1

?+∞∞

--=

ηηπ

ηd e

t a x 2

2cos sin 1

?

+∞

∞

---=ηπ

ηη

d e e x t a

i 2

2212

sin

4

422

12

sin t a e x -=4

422

12

sin t a e x -=x e t a sin 2

-=, ()

4

2

2

2

1λη

-

∧

-=

e e .

例2求初值问题?????=??=??=x u

x u a

t u t cos ,022

2的解x e t x u t a cos ),(2-=.

例3求初值问题?????+===1

,

2

02

x u u a u t xx t 的解. 解1 直接观察t a x t x u 2

221),(++= 2. []

?+∞∞

--++=

ηηπ

ηd e t a x t x u 2

1)2(1

),(2

[]

?+∞∞

--+++=

ηηηπ

ηd e t a t ax x 2

1441

222

t a x 2221++=

从这几个实例上,更直观明显的证明求解公式的正确,对模型方程的正确性,提供保证.

?????++===1cos ,

2

2

x x u u a u t xx t 定理 设)(x ?在),(+∞-∞上连续且有界，

),(t x f ，(,)x f x t 在],0[),(T ?+∞-∞上连续且有界，

令 ?∞

+∞

---

=

ξξ?πξd e

t

a

t x u t

a x 224)()(21),(??∞

+∞

----

-+

ξτ

τξτπ

τξd e t f d a t a x t )

(4)(0

221

)

,(21

，

其中常数0>a ，则有)(),(lim 00,0x t x u t x x ?=+

→→；(,)u x t 问题

??

??

?+∞<<∞-=>+∞<<∞-=??-??x x x u t x t x f x u a t u ),()0,(0

,),,(22

2?

的解。

证明 由于ηη?π

η?

+∞

∞

--+=

d e t a x t x u 2

)2(1

),(ητηττπ

η?

?+∞

∞

---++

d e t a x f d t 2

),2(1

，

利用控制收敛定理，得),(lim 0,0t x u t x x +

→→)(0)(1

002

x d e x ?η?π

η=+=

?

+∞

∞

--；

22()4(,)())x a t

u x t d t t ξ?ξξ

--

+∞

-∞??=??

?22()4()

(,))x t a t d f d t ξττξτξ--+∞

--∞

???

(,)f x t +；

22()2

2

422(,)())x a t

u x t d x x ξφξξ

--

+∞

-∞

??=??

?22()2

4()

20

(,)

)x t

a t d f d x ξττξτξ--

+∞

--∞

???

，

显然成立2

22

(,)u u a f x t t x ??-=??，结论得证。

定理 假设函数),(t x f ，)(x ?关于x 都是解析的，则问题

??

???+∞<<∞-=>+∞<<∞-=??-??x x x u t x t x f x u a t

u ),()0,(0

,),,(2

2

2?

的解可以写成

∑?∑∞=∞

=-+=0

0)2(20)

2(2),(!)]([)(!)(),(n t n x n n n n d x f n t a x n t a t x u τττ?

其中)()

2(x n ?

和),()2(τx f n x 分别是)(x ?和),(τx f 关于x 的n 2阶导数。

证明：?

+∞

∞

--+=

ηη?π

ηd e

t a x t x u 2

)2(1

),(

??

+∞

∞

---++t d d e t a x f 02

),2(1

τητητπ

η

?∑

∞+∞

--∞

==

ηη?π

ηd e t a k x k k k 2

)()2(!

)

(1

??∑

∞+∞

--∞

=-+

t k k k x d d e t k x f 00

)(2

)(!

),(1

τηηττπ

η ∑?∑∞=∞

=-+=0

0)2(20)

2(2),(!)]([)(!)(n t n x n n n n d x f n t a x n t a τττ?。

例 求解定解问题2

22

,(,0),

(,0)cos ,().u u a Ax x t t x

u x x x θ???-=-∞<<+∞>?????=-∞<<+∞?

其中,A θ是常数。 解

方法一：2

(,)cos (2)u x t x e

d ηθη+∞

--∞

=

+?

??

+∞

∞

---++

t

d d

e t a x A 02

)2(1

τηητπ

η

2

cos cos2x a e d Axt ηθθη+∞

--∞

=

?+

22cos a t

e

x Axt θθ-=+；

方法二：2(2)0

()(,)(cos )!n

n n a t u x t x Axt n θ∞

==+∑ 220

()(1)cos !n

n n a t x Axt n θθ∞

==-+∑

22

cos a t e x Axt θθ-=+。

解的性质与物理解释（对齐次方程??

???==??-??),()0,(,022

2x x u x u a t u ?

1.（奇偶性与周期性）

若?是奇（偶,周期为l 的）函数,则解?

∞

+∞

---

=ξξ?πξd e

t a

t x u t

a x 224)()(21),(亦是x 的奇（偶,

周期为l 的）函数. 2.（无限传播速度）

如果杆的初始温度)(x ?只在小段),(00δδ+-x x 上不为零,不妨假设0)(>x ?,即

),(,0)(00δδ?+-∈>x x x x ,其它处0)(=x ?.那么当0>t ,杆上各点的温度0)(21)(21),(0022224)(4)(>=

=

?

?

+---

∞

+∞

---

δ

δ

ξξξξ?πξξ?πx x t

a x t

a x d e

t a

d e

t a

t x u

也就是说在顷刻之间,热量就传递到杆上的任意一点,当然在0x 附近的点所受到的影响较大（t

a x e

224)(ξ--来定）,而离0x 较远的点受到的影响较小.这与我们知道的物理现象一致.

初值问题解的渐近性态

讨论当+∞→t 时,热传导方程初值问题解的渐近性态.由前面的讨论可知,当)(x ?为有

界连续函数时,热传导方程的初值问题?????===)(,

2

x u u a u t xx t ?的解的唯一性,由下列Possion 积分给出

?

∞

+∞

---

=

ξξ?πξd e

t a

t x u t

a x 224)()(21),(.

为了讨论解的渐近性态,还需要对)(x ?加进一步的条件. 如果

?

+∞

∞

-dx x )(?收敛,则称),(1R L ∈?并记?+∞

∞

-=dx x R L )()

(1??

.

定理设?是有界连续函数,且),(1

R L ∈?则初值问题的唯一经典解（古典解）具有如下的渐近

性态：对一切0,>∈t R x 当+∞→t 时,一致地成立

0),(2

1→≤-

Ct

t x u ,（+∞→t ）.

其中C 为一个仅与a 与

)

(1R L ?

有关的正常数.

证明 由?

∞

+∞

---=

ξξ?πξd e

t a

t x u t

a x 224)()(21),(,

?

∞

+∞

---≤

ξξ?πξd e t a

t x u t

a x 224)()(21),(

?

+∞

∞

-≤

ξξ?πd t a

)(21

2

12

1)

(121--==Ct

t a R L ?π

证毕.

物理现象符合,高温→变低温,以至冷却到0. 热的不可逆性：

??

???=<+∞<<∞-=??-??),()0,(0

,,022

2x x u t x x u a

t

u ?

对一般的)(x ?,解不存在,说明热的逆现象是不确定的.

热传导方程的初值问题

§2热传导方程的初值问题 一维热传导方程的初值问题(或Cauchy 问题) ?? ???+∞<<∞-=>+∞<<∞-=??-??x x x u t x t x f x u a t u ),()0,(0 ,),,(2 2 2? （） 偏导数的多种记号xx x t u x u u x u u t u =??=??=??22,,. 问题也可记为 ?? ?+∞ <<∞-=>+∞<<∞-=-x x x u t x t x f u a u xx t ),()0,(0 ,,),(2?. Fourier 变换 我们将用Fourier 变换法求解热传导方程的柯西问题.为此我们将着重介绍Fourier 变换的基本知识.Fourier 变换在许多学科中是重要使用工具. 可积函数,设)(x f f =是定义在),(+∞-∞上的函数, 且对任意A B <，()f x 在[,]A B 上 可积，若积分 ? +∞ ∞ -dx x f )(收敛,则称)(x f 在),(+∞-∞上绝对可积。 将),(+∞-∞上绝对可积函数形成的集合记为),(1 +∞-∞L 或),(+∞-∞L ， 即{ } ∞<=+∞-∞=+∞-∞? +∞ ∞ -dx x f f L L )(| ),(),(1 ，称为可积函数空间. 连续函数空间: ),(+∞-∞上全体连续函数构成的集合,记为),(+∞-∞C , {}上连续在),(|),(+∞-∞=+∞-∞f f C ， {}上连续在),(,|),(1+∞-∞'=+∞-∞f f f C 。 定义 若),(+∞-∞∈L f ,那么积分 ),(?)(21 λπ λf dx e x f x i =? +∞ ∞ -- 有意义,称为Fourier 变换, )(? λf 称为)(x f 的Fourier 变式(或Fourier 变换的象). ? +∞ ∞ --= =dx e x f f Ff x i λπ λλ)(21)(?)( 定理 (Fourier 积分定理)若),(),(1 +∞-∞?+∞-∞∈C L f ,那么我们有

热传导+对流微分方程推导

热传导微分方程 导热又称热传导,是两个相互接触的物体或同一物体的各部分之间，由于温度不同而引起的热量传递现象.此时热量主要依靠分子、原子及自由电子等微观粒子的运动进行传递，没有明显的物质转移。热量可以通过固体、液体以及气体进行传导，但是严格来说,单纯的导热只发生在密实的固体物质中。 1 傅立叶定律 傅立叶定律是导热理论的基础。其向量表达式为： q gradT λ=-? （2-1) 式中：q ——热流密度,是一个向量,2/()Kcal m h gradT ——温度梯度，也是一个向量,℃／m . λ--导热系数，又称热导率，/()Kcal mh C ； 式中的负号表示q 的方向始终与gradT 相反。 2 导热系数（ｔh ｅrm ａl c ｏｎdu ｃt ｉv ｉty ）及其影响因素 导热系数λ（/()Kcal mh C ）是热传导过程中一个重要的比例常数,在数值上等于每小时每平方米面积上，当物体内温度梯度为1℃/m 时的导热量. 导热系数是指在稳定传热条件下,1m 厚的材料,两侧表面的温差为１度（Ｋ,°C),在1秒内，通过1平方米面积传递的热量，用λ表示,单位为瓦/米·度,w/m·k(W/m·K,此处的K 可用℃代替). 导热系数为温度梯度1℃/m ，单位时间通过每平方米等温面的热传导热流量.单位是:Ｗ/（m·K）。 在上述假设前提下，建立煤层瓦斯流动数学模型的控制方程. ３．热传导微分方程推导 在ｔ时刻ｗ界面的温度梯度为 x T ?? 在t 时刻e 界面的温度梯度为dx x T x T dx x x T x T 22??+??=???? +??

单位时间内六面体在x 方向流入的热流量为：dydz x T ??-λ ; 单位时间内六面体在x 方向流出的热流量为:dydz dx x T x T ?? ? ? ????+??-22λ； 单位时间内六面体在x 方向流入的净热量为：dxdydz x T 22??λ 图3-1 微分单元体各面上进出流量示意图 同理，单位时间内六面体在y 方向流入的净热量为:dxdydz y T 22??λ 单位时间内六面体在y 方向流入的净热量为：dxdydz z T 22??λ 单位时间内流入六面体的总热量为： dxdydz z T y T x T ??? ?????+??+??222222λ （3－1)

导热方程求解matlab

使用差分方法求解下面的热传导方程 2 (,)4(,) 0100.2t xx T x t T x t x t =<<<< 初值条件：2(,0)44T x x x =- 边值条件：(0,)0(1,)0 T t T t == 使用差分公式 1,,1,2 2 2 (,)2(,)(,) 2(,)()i j i j i j i j i j i j xx i j T x h t T x t T x h t T T T T x t O h h h -+--++-+= +≈ ,1,(,)(,) (,)()i j i j i j i j t i j T x t k T x t T T T x t O k k k ++--= +≈ 上面两式带入原热传导方程 ,1,1,,1,2 2i j i j i j i j i j T T T T T k h +-+--+= 令2 24k r h =，化简上式的 ,1,1,1,(12)()i j i j i j i j T r T r T T +-+=-++ i x j t j

编程MA TLAB 程序，运行结果如下 1 x t T function mypdesolution c=1; xspan=[0 1]; tspan=[0 0.2]; ngrid=[100 10]; f=@(x)4*x-4*x.^2; g1=@(t)0; g2=@(t)0; [T,x,t]=rechuandao(c,f,g1,g2,xspan,tspan,ngrid); [x,t]=meshgrid(x,t); mesh(x,t,T); xlabel('x') ylabel('t') zlabel('T') function [U,x,t]=rechuandao(c,f,g1,g2,xspan,tspan,ngrid) % 热传导方程：

热传导方程及其定解问题的导出

第一章 热传导方程 本章介绍最典型的抛物型方程—热传导方程,在研究热传导,扩散等物理现象时都会遇 到这类方程. §1 热传导方程及其定解问题的导出 1.1热传导方程的导出 物理模型 在三维空间中,考虑一均匀,各向同性的物体Ω,假定它内部有热源,并且与周围介质有热交换,需要来研究物体内部温度的分布和变化. 以函数),,,(t z y x u 表示物体Ω在位置),,(z y x 及时刻t 的温度.物体内部由于各部分温度不同,产生热量的传递,它们遵循能量守恒定律. 能量守恒定律 物体内部的热量的增加等于通过物体的边界流入的热量与由物体内部的热源所生成的热量的总和 . 在物体Ω内任意截取一块D .现在时段],[21t t 上对D 使用能量守恒定律. 设),,,(t z y x u u =是温度(度),c 是比热(焦耳∕度·千克),ρ是密度(千克/米3), q 是热流密度(焦耳/秒·米2),0f 是热源强度(焦耳/千克·秒). 注意到在dt 时段内通过D 的边界D ?上小块dS 进入区域D 的热量为dSdt n q ?-(n 是 D ?的外法向),从而由能量守恒律,我们有 ,)||(21 21 120??????????+?-=-?==t t D t t D D t t t t dxdydz f dt ds n q dt dxdydz u u c ρρ （1.1） 大家知道,热量流动的原因是因为在物体内部存在温差.依据传热学中的傅立叶实验定律,在一定条件下,热流向量与温度梯度成正比 ,u k q ?-= （梯度? ?? ? ????????==?z u y u x u gradu u ,,） （1.2） 这里负号表明热量是由高温向低温流动,k 是物体的导热系数.

热传导方程的求解

应用物理软件训练 前言 MATLAB 是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。 MATLAB是矩阵实验室（Matrix Laboratory）的简称，和Mathematica、Maple 并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其

他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。本部分主要介绍如何根据所学热传导方程的理论知识进行MATLAB数值实现可视化。本部分主要介绍如何根据所学热传导方程的理论知识进行MATLAB数值实现可视化。本部分主要介绍如何根据所学热传导方程的理论知识进行MATLAB数值实现可视化。 本部分主要介绍如何根据所学热传导方程的理论知识进行MATLAB数值实现可视化。

题目：热传导方程的求解 目录 一、参数说明 (1) 二、基本原理 (1) 三、MATLAB程序流程图 (3) 四、源程序 (3) 五、程序调试情况 (6) 六、仿真中遇到的问题 (9) 七、结束语 (9) 八、参考文献 (10)

一、参数说明 U=zeros(21,101) 返回一个21*101的零矩阵 x=linspace(0,1,100);将变量设成列向量 meshz(u)绘制矩阵打的三维图 axis([0 21 0 1]);横坐标从0到21，纵坐标从0到1 eps是MATLAB默认的最小浮点数精度 [X,Y]=pol2cart(R,TH);效果和上一句相同 waterfall(RR,TT,wn)瀑布图 二、基本原理 1、一维热传导问题 （1）无限长细杆的热传导定解问题 利用傅里叶变换求得问题的解是： 取得初始温度分布如下 这是在区间0到1之间的高度为1的一个矩形脉冲，于是得 （2）有限长细杆的热传导定解问题

一维热传导MATLAB模拟

昆明学院2015届毕业设计（论文） 设计（论文）题目 一维热传导问题的数值解法及其MATLAB模拟子课题题目无 姓名伍有超 学号201117030225 所属系物理科学与技术系 专业年级2011级物理学2班 指导教师王荣丽 2015 年 5 月

摘要 本文介绍了利用分离变量法和有限差分法来求解一维传导问题的基本解，并对其物理意义进行了讨论。从基本解可以看出，在温度平衡过程中，杠上各点均受初始状态的影响，而且基本解也满足归一化条件，表示在热传导过程中杆的总热量保持不变。通过对一维杆热传导的分析，利用分离变量法和有限差分法对一维热传导进行求解，并用MATLAB 数学软件来对两种方法下的热传导过程进行模拟，通过对模拟所得三维图像进行取值分析，得出由分离变量法和有限差分法绘制的三维图基本相同，且均符合热传导过程中温度随时间、空间的变化规律，所以两种方法均可用来解决一维热传导过程中的温度变化问题。 关键词：一维热传导；分离变量法；有限差分法；数值计算；MATLAB 模拟

Abstract In this paper, the method of variable separation and finite difference method are introduced to solve the problem of one-dimensional heat conduction problems, and the physical significance of numerical methods for heat conduction problems are discussed. From the basic solution, we can see the temperature on the bar are affected by the initial state during the process of temperature balance, and basic solution also satisfy the normalization condition which implied the invariance of the total heat in the bar during the heat conduction process. Through the analysis of the one-dimensional heat conduction, by taking use of variable separation method and finite difference method, we simulated the one-dimensional heat conduction problem by MATLAB. The three-dimensional images of the simulation results obtained by the method of separation of variables and finite difference method are similar to each other, and the temperature curve is in accordance with the law of temperature variation during heat conduction. Thus, we can go to the conclusion that both methods can be used to deal with the one-dimensional heat conduction problems. Keywords: One-dimensional heat conduction; method of variable separation; finite difference method; numerical method; MATLAB simulation

一维非稳态导热问题的数值解

计算传热学程序报告 题目：一维非稳态导热问题的数值解 ： 学号： 学院：能源与动力工程学院 专业：工程热物理 日期：2014年5月25日

一维非稳态导热问题数值解 求解下列热传导问题： ? ?? ????=====≤≤=??- ??1,10),(,1),0(0)0,()0(01T 22ααL t L T t T x T L x t T x 1.方程离散化 对方程进行控制体积分得到： dxdt t T dxdt x T t t t e w t t t e w ? ?? ??+?+??=??α 1 2 2 ? ? -=??-???+?+e w t t t w e t t t dx T T dt x T x T )(1])()( [α 非稳态项：选取T 随x 阶梯式变化，有 x T T dx T T t p t t p e w t t t ?-=-?+?+? )()( 扩散项：选取一阶导数随时间做显示变化，有 t x T x T dt x T x T t w t e w e t t t ???-??=??-??? ?+])()[(])()[( 进一步取T 随x 呈分段线性变化，有 e P E e x T T x T )()( δ-=?? ， w W P w x T T x T )()(δ-=?? 整理可以得到总的离散方程为： 2 21x T T T t T T t W t P t E t P t t E ?+-=?-?+α 2.计算空间和时间步长 取空间步长为： h=L/N 网格Fourier 数为： 2 2 0x t x t F ??= ??= α（小于0.5时稳定）

热传导方程傅里解

热传导方程傅里解

————————————————————————————————作者：————————————————————————————————日期：

热传导在三维的等方向均匀介质里的传播可用以下方程表达： 其中： ?u =u(t, x, y, z) 表温度，它是时间变量t 与空间变量(x,y,z) 的函数。 ?/是空间中一点的温度对时间的变化率。 ?, 与温度对三个空间座标轴的二次导数。 ?k决定于材料的热传导率、密度与热容。 热方程是傅里叶冷却律的一个推论（详见条目热传导）。 如果考虑的介质不是整个空间，则为了得到方程的唯一解，必须指定u 的边界条件。如果介质是整个空间，为了得到唯一性，必须假定解的增长速度有个指数型的上界，此假定吻合实验结果。 热方程的解具有将初始温度平滑化的特质，这代表热从高温处向低温处传播。一般而言，许多不同的初始状态会趋向同一个稳态（热平衡）。因此我们很难从现存的热分布反解初始状态，即使对极短的时间间隔也一样。 热方程也是抛物线偏微分方程最简单的例子。 利用拉普拉斯算子，热方程可推广为下述形式

其中的是对空间变量的拉普拉斯算子。 热方程支配热传导及其它扩散过程，诸如粒子扩散或神经细胞的动作电位。热方程也可以作为某些金融现象的模型，诸如布莱克-斯科尔斯模型与 Ornstein-Uhlenbeck 过程。热方程及其非线性的推广型式也被应用于影像分析。量子力学中的薛定谔方程虽然有类似热方程的数学式（但时间参数为纯虚数），本质却不是扩散问题，解的定性行为也完全不同。 就技术上来说，热方程违背狭义相对论，因为它的解表达了一个扰动可以在瞬间传播至空间各处。扰动在前方光锥外的影响通常可忽略不计，但是若要为热传导推出一个合理的速度，则须转而考虑一个双曲线型偏微分方程。 以傅里叶级数解热方程[编辑] 以下解法首先由约瑟夫·傅里叶在他于1822年出版的著作Théorie analytique de la chaleur（中译：解析热学）给出。先考虑只有一个空间变量的热方程，这可以当作棍子的热传导之模型。方程如下： 其中u = u(t, x) 是t和x的双变量函数。 ?x是空间变量，所以x∈[0,L]，其中L表示棍子长度。

一维热传导方程

一维热传导方程Last revision on 21 December 2020

一维热传导方程 一． 问题介绍 考虑一维热传导方程： （1） ,0),(22T t x f x u a t u ≤<+??=?? 其中a 是正常数，)(x f 是给定的连续函数。按照定解条件的不同给法，可将方程（1）的定解问题分为两类： 第一类、初值问题（也称Cauthy 问题）：求具有所需次数偏微商的函数),(t x u ，满足方程（1）（∞<<∞-x ）和初始条件： （2） ),()0,(x x u ?= ∞<<∞-x 第二类、初边值问题（也称混合问题）：求具有所需次数偏微商的函数),(t x u ，满足方程（1）（l x <<0）和初始条件： （3） ),()0,(x x u ?= l x <<0 及边值条件 (4) .0),(),0(==t l u t u T t ≤≤0 假定)(x ?在相应区域光滑，并且在l x ,0=满足相容条件，使上述问题有唯一充分光滑的解。 二． 区域剖分 考虑边值问题（1），（4）的差分逼近。去空间步长N l h /=和时间步长M T /=τ，其中N,M 都是正整数。用两族平行直线： 将矩形域}0;0{T t l x G ≤≤≤≤=分割成矩形网格，网格节点为),(k j t x 。以h G 表示网格内点集合，即位于开矩形G 的网点集合；h G 表示所有位于闭矩形G 的网点集合；Γ=G --G 是网格界点集合。

三． 离散格式 第k+1层值通过第k 层值明显表示出来，无需求解线性代数方程组，这样的格式称为显格式。 第k+1层值不能通过第k 层值明显表示出来，而由线性代数方程组确定，这样的格式称为隐格式。 1. 向前差分格式 （5） ,22111j k j k j k j k j k j f h u u u a u u ++-=--++τ )(j j x f f =， )(0 j j j x u ??==， 00==k N k u u ， 其中j = 1,2,…，N-1，k = 1,2,…,M-1。以2/h a r τ=表示网比。则方程（5）可以改写为： 易知向前差分格式是显格式。 2. 向后差分格式 （6） ,11111)21(j k j k j k j k j f u ru u u ru τ+=-++-+-+++ )(0 j j j x u ??==， 00==k N k u u ， 其中j = 1,2,…，N-1，k = 1,2,…,M-1，易知向前差分格式是显格式。 3. 六点对称格式（Grank-Nicolson 格式） 将向前差分格式和向后差分格式作算术平均，即得到六点对称格式： （7） 111112)1(2+-+++-++-k j k j k j u r u r u r =j k j k j k j f u r u r u r τ++-+-+112 )1(2 利用0j u 和边值便可逐层求到k j u 。六点对称格式是隐格式，由第k 层计算第k+1层时需解线性代数方程组（因系数矩阵严格对角占优，方程组可唯一求解）。

第四章导热题的数值解法

第四章导热问题的数值解法 1 、重点内容：①掌握导热问题数值解法的基本思路； ②利用热平衡法和泰勒级数展开法建立节点的离散方程。 2 、掌握内容：数值解法的实质。 3 、了解内容：了解非稳态导热问题的两种差分格式及其稳定性。 §4—1导热问题数值求解的基本思想及内节点方程的建立由前述 3 可知，求解导热问题实际上就是对导热微分方程在定解条件下的积分求解，从而获得分析解。但是，对于工程中几何形状及定解条件比较复杂的导热问题，从数学上目前无法得出其分析解。随着计算机技术的迅速发展，对物理问题进行离散求解的数值方法发展得十分迅速，并得到广泛应用，并形成为传热学的一个分支——计算传热学（数值传热学），这些数值解法主要有以下几种： （１）有限差分法（ 2 ）有限元方法（ 3 ）边界元方法 数值解法能解决的问题原则上是一切导热问题，特别是分析解方法无法解决的问题。如：几何形状、边界条件复杂、物性不均、多维导热问题。 一．分析解法与数值解法的异同点： ?相同点：根本目的是相同的，即确定① t=f(x ， y ， z) ；② 。 ?不同点：数值解法求解的是区域或时间空间坐标系中离散点的温度分布代替连续的温度场；分析解法求解的是连续的温度场的分布特征，而不是分散点的数值。 数值求解的基本思路及稳态导热内节点离散方程的建立 二．解法的基本概念 ?实质 对物理问题进行数值解法的基本思路可以概括为：把原来在时间、空间坐标系中连续的物理量的场，如导热物体的温度场等，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，来获得离散点上被求物理量的值。该方法称为数值解法。 这些离散点上被求物理量值的集合称为该物理量的数值解。 2 、基本思路：数值解法的求解过程可用框图 4-1 表示。 由此可见： 1 ）物理模型简化成数学模型是基础； 2 ）建立节点离散方程是关键； 3 ）一般情况微分方程中，某一变量在某一坐标方向所需边界条件的个数等于该变量在该坐标方向最高阶导数的阶数。 ?数值求解的步骤 如图 4-2 （ a ），二维矩形域内无内热源、稳态、常物性的导热问题采用数值解法的步骤如下：（1）建立控制方程及定解条件 控制方程：是指描写物理问题的微分方程 针对图示的导热问题，它的控制方程（即导热微分方程）为：（ a ） 边界条件： x=0 时， x=H 时， 当 y=0 时，

3热传导方程的初边值问题(精选、)

例4 周期初始温度分布 求解热传导方程t xx u u =,(,0)x t -∞<<+∞>给定初始温度分布 (,0)1cos 2,()u x x x =+-∞<<+∞。 解 4(,)1cos2t u x t e x -=+. 初始高斯温度分布 例 5求解定解问题22 22 0,(,0) (,0),()kx u u a x t t x u x e x -???-=-∞<<+∞>?????=-∞<<+∞? ， 其中常数0k >. 解 22()4(,)()x s a t u x t s e ds ?-- +∞ -∞ = ? 22 2()4x s ks a t e e ds -- +∞ --∞ = ? 222 2(41)24ka t s xs x a t e ds +-+- +∞ -∞ = ? 222 22224(41)()41414x ka t ka t s x ka t ka t a t e ds +- +++- +∞ -∞ = ? 22 2 222(41)()41 441 k ka t x x s ka t a t ka t e e ds +---+∞ ++-∞ = ? 2241 k x ka t e - += 2241 k x ka t - += . §3初边值问题 设长度为l ,侧表面绝热的均匀细杆,初始温度与细杆两端的温度已知,则杆上的温度分布 ),(t x u 满足以下初边值问题 ?? ? ??≤<==≤≤=<<<<=-T t t g t l u t g t u l x x x u T t l x t x f u a u xx t 0),(),(),(),0(,0), ()0,(0,0),,(212? 对于这样的问题,可以用分离变量法来求解. 将边值齐次化

偏微分方程数值解法

《偏微分方程数值解法》 课程设计 题目：六点对称差分格式解热传导方程的初边值 问题 姓名：王晓霜 学院：理学院 专业：信息与计算科学 班级： 0911012 学号： 091101218 指导老师：翟方曼 2012年12月14日

一、题目 用六点对称差分格式计算如下热传导方程的初边值问题 222122,01,01(,0),01 (0,),(1,),01x t t u u x t t x u x e x u t e u t e t +???=<<<≤?????=≤≤??==≤≤??? 已知其精确解为 2(,)x t u x t e += 二、理论 1．考虑的问题 考虑一维模型热传导方程 (1.1) )(22x f x u a t u +??=??，T t ≤<0 其中a 为常数。)(x f 是给定的连续函数。（1.1）的定解问题分两类： 第一，初值问题（Cauchy 问题）：求足够光滑的函数()t x u ,，满足方程（1.1）和初始条件： （1.2） ()()x x u ?=0,, ∞<<∞-x 第二，初边值问题（也称混合问题）：求足够光滑的函数()t x u ,，满足方程（1.1）和初始条件： ()13.1 ()()x x u ?=0,, l x l <<- 及边值条件 ()23.1 ()()0,,0==t l u t u ， T t ≤≤0 假定()x f 和()x ?在相应的区域光滑，并且于()0,0，()0,l 两点满足相容条件,则上述问题有唯一的充分光滑的解。 现在考虑边值问题（1.1），（1.3）的差分逼近 取 N l h = 为空间步长，M T =τ为时间步长，其中N ，M 是自然数， jh x x j ==, ()N j ,,1,0Λ=； τk y y k ==, ()M k ,,1,0Λ=

一维热传导方程(Richardson格式)

中南林业科技大学 偏微分方程数值解法学生姓名：周晓虹 学号：20083710 学院：理学院 专业年级：08信计1班 设计题目：一维热传导方程的Richardson格式 2011年06月

一． 问题介绍 考虑一维热传导方程： （1） ,0),(22 T t x f x u a t u ≤<+??=?? 其中a 是正常数，)(x f 是给定的连续函数。按照定解条件的不同给法，可将方程（1）的定解问题分为两类： 第一类、初值问题（也称Cauthy 问题）：求具有所需次数偏微商的函数),(t x u ，满足方 程（1）（∞<<∞-x ）和初始条件： （2） ),()0,(x x u ?= ∞<<∞-x 第二类、初边值问题（也称混合问题）：求具有所需次数偏微商的函数),(t x u ，满足方 程（1）（l x <<0）和初始条件： （3） ),()0,(x x u ?= l x <<0 及边值条件 (4) .0),(),0(==t l u t u T t ≤≤0 假定)(x ?在相应区域光滑，并且在l x ,0=满足相容条件，使上述问题有唯一充分光滑 的解。 二． 区域剖分 考虑边值问题（1），（4）的差分逼近。去空间步长N l h /=和时间步长M T /=τ，其中N,M 都是正整数。用两族平行直线： ),,1,0(N j jh x x j === ),,1,0(M k k t t k ===τ 将矩形域}0;0{T t l x G ≤≤≤≤=分割成矩形网格，网格节点为),(k j t x 。以h G 表示网格内点集合，即位于开矩形G 的网点集合；h G 表示所有位于闭矩形G 的网点集合； h Γ=h G --h G 是网格界点集合。

热传导+对流微分方程推导(精.选)

热传导微分方程 导热又称热传导，是两个相互接触的物体或同一物体的各部分之间，由于温度不同而引起的热量传递现象。此时热量主要依靠分子、原子及自由电子等微观粒子的运动进行传递，没有明显的物质转移。热量可以通过固体、液体以及气体进行传导，但是严格来说，单纯的导热只发生在密实的固体物质中。 1 傅立叶定律 傅立叶定律是导热理论的基础。其向量表达式为: q gradT λ=-? (2-1) 式中:q ——热流密度，是一个向量，2/()Kcal m h gradT ——温度梯度，也是一个向量，℃/m 。 λ——导热系数，又称热导率，/()Kcal mh C o ； 式中的负号表示q 的方向始终与gradT 相反。 2 导热系数（thermal conductivity ）及其影响因素 导热系数λ（ /()Kcal mh C o ）是热传导过程中一个重要的比例常数，在数值上等于每小时每平方米面积上，当物体内温度梯度为1℃/m 时的导热量。 导热系数是指在稳定传热条件下，1m 厚的材料，两侧表面的温差为1度（K ，°C）,在1秒内，通过1平方米面积传递的热量，用λ表示，单位为瓦/米·度，w/m·k （W/m·K,此处的K 可用℃代替）。 导热系数为温度梯度1℃/m ，单位时间通过每平方米等温面的热传导热流量。单位是：W/(m·K)。 在上述假设前提下，建立煤层瓦斯流动数学模型的控制方程。 3．热传导微分方程推导 在t 时刻w 界面的温度梯度为 x T ?? 在t 时刻e 界面的温度梯度为dx x T x T dx x x T x T 22??+??=???? +??

单位时间内六面体在x 方向流入的热流量为：dydz x T ??-λ ； 单位时间内六面体在x 方向流出的热流量为：dydz dx x T x T ?? ? ? ????+??-22λ； 单位时间内六面体在x 方向流入的净热量为：dxdydz x T 22??λ 图3-1 微分单元体各面上进出流量示意图 同理，单位时间内六面体在y 方向流入的净热量为：dxdydz y T 22??λ 单位时间内六面体在y 方向流入的净热量为：dxdydz z T 22??λ 单位时间内流入六面体的总热量为： dxdydz z T y T x T ??? ?????+??+??222222λ (3-1)

偏微分方程数值解法的MATLAB源码

[原创]偏微分方程数值解法的MATLAB源码【更新完毕】 说明：由于偏微分的程序都比较长，比其他的算法稍复杂一些，所以另开一贴，专门上传偏微分的程序谢谢大家的支持！其他的数值算法见：..//Announce/Announce.asp?BoardID=209&id=82450041、古典显式格式求解抛物型偏微分方程（一维热传导方程） function [U x t]=PDEParabolicClassicalExplicit(uX,uT,phi,psi1,psi2,M,N,C) %古典显式格式求解抛物型偏微分方程 %[U x t]=PDEParabolicClassicalExplicit(uX,uT,phi,psi1,psi2,M,N,C) % %方程：u_t=C*u_xx 0 <= x <= uX,0 <= t <= uT %初值条件：u(x,0)=phi(x) %边值条件：u(0,t)=psi1(t), u(uX,t)=psi2(t) % %输出参数：U -解矩阵，第一行表示初值，第一列和最后一列表示边值，第二行表示第2层…… % x -空间变量 % t -时间变量 %输入参数：uX -空间变量x的取值上限 % uT -时间变量t的取值上限 % phi -初值条件，定义为内联函数 % psi1 -边值条件，定义为内联函数 % psi2 -边值条件，定义为内联函数 % M -沿x轴的等分区间数 % N -沿t轴的等分区间数 % C -系数，默认情况下C=1 % %应用举例： %uX=1;uT=0.2;M=15;N=100;C=1; %phi=inline('sin(pi*x)');psi1=inline('0');psi2=inline('0'); %[U x t]=PDEParabolicClassicalExplicit(uX,uT,phi,psi1,psi2,M,N,C); %设置参数C的默认值 if nargin==7 C=1; end %计算步长 dx=uX/M;%x的步长 dt=uT/N;%t的步长 x=(0:M)*dx; t=(0:N)*dt;

热传导方程差分格式

热传导方程的差分格式第2页 一维抛物方程的初边值问题 分别用向前差分格式、向后差分格式、六点对称格式，求解下列问题： .:u ；：2u a 2，0 ：：： x ：： 1, .:t ；x u(x,0) =sin 二x, 0 ：： x :: 1 u(O,t) =u(1,t) =0, t 0 2 在t =0.05,0.1和0.2时刻的数值解，并与解析解u(x,t)=ef t s in (二x)进行比较 1差分格式形式 2 设空间步长h =1/ N ,时间步长? ? 0, T =M ?，网比r = ? / h . (1) 向前差分格式 该问题是第二类初边值问题(混合问题)，我们要求出所需次数的偏微商的函数 Eu c2u u(x,t)，满足方程—a—^, 0 ：：：x ：：：1,和初始条件u(x,0)= sin x , 0 x ：：： 1抚ex 及边值条件u(0,t) =u(1,t) =0, t 0。 已知sin二x在相应区域光滑，并且在x =0,丨与边值相容，使问题有唯一充分光滑的解。 取空间步长h =1/ N，和时间步长-T /M，其中N,M都是正整数。用两族平行 直线x= j X =( j h0Hj 1 ，和tlNt k =X(k=0,1,|||,M) 将矩形域 G二｛0 — x — 1,t — 0｝分割成矩形网络，网络格节点为(X j,t k)。以G h表示网格内点集合， 即位于矩形G的网点集合；G h表示闭矩形G的网格集合；丨h=G h-G h是网格界点的集合。 向前差分格式，即 k 1 k k k k u, -u, u, 1 -2比■ u, 4 - j二a 亠2亠t ( 1) ￡h