数理方程参考答案4第四章 积分变换法
复变函数与积分变换第四章习题解答
2!
3!
2!
3!
3!
5!
2
4
而收敛半径 R=扛'fJ •
而收敛半径 R=+oo;
(7)
z
而收敛半径 R=l 。
cos 土 ==1- 上 (z+z2 + z3 + .. 一上 (z+z2 +z3 + ...r +... =1-2. z2 - z 3 +...' I zI< 1 I 1-z 2 4! 2
'
In n
1
n
1
4)因 cos in= cbn,
( 1)每一个幕级数在它的收敛圆周上处处收敛 ;
4. 下列说法是否正确?为什么?
而lim—-=1=0,
II�")
chn
2"
故
cosm 2 — " 发散。 11=2 2
00
(2)每一个幕级数的和函数在收敛圆内可能有奇点;
解 (1)不对。如Iz"在收敛圆lzl < 1内收敛, 但在收敛圆周日=l上井不收敛; (3)不对。如八 z) 三在全平面上连续, 但它在任何点的邻域内均不能展开成 Taylor 级 5幕级数LC11 (z-2)" 能否在z=0收敛而在z=3发散?
=
=早-(于)2 f ()
11=]
一I
干是收敛半径 R=2 。 (2)因
(-1t z-1 "' "
2
+ ... + ( -1 y,-1
(早厂
lz-11<2
l
及
飞(z�2 一言) = z�2 一士 2 = = 1-'� 厂; J- J [ =』 z�2 4 +(:-2i ± + � 2 �
工程数学_积分变换
1 T2 j n t jn f (t ) lim T fT ( )e d e T T n 2
2p 当n取一切整数时, n =n n 所对应的点便 T 均匀分布在整个数轴上,
如图
2p 2p 2p T T T 2p T
m 1
2
T 2
T an T cos nt d t an 2 2 T 2 2 即 an T fT (t )cos nt d t T 2
T 2
2
同理, 为求bn, 计算[fT(t), sin nt], 即
T 2 T 2
a0 fT (t )sin nt d t T sin nt d t 2 2
t
最常用的一种周期函数是三角函数 fT(t)=Asin(t+j) 其中 A 称为振幅,=2p/T 称为角频率,j 称为初相角
t
而Asin(t+j)又可以看作是两个周期函数sint和cost 的线性组合 Asin(t+j)=asint+bcost
人们发现, 所有的工程中使用的周期函数都可以用一 系列的三角函数的线性组合来逼近.
T 2
T 2
一. Fourier级数
1.Dirichlet条件
若函数在区间[-T/2,T/2]上满足: 1, 连续或只有有限个第一类间断点; 2, 只有有限个极值点
则称函数满足Dirichlet条件. 注: 这两个条件实际上就是要保证函数是可积函数.
第一类间断点和第二类间断点的区别:
第二类间断点
a0 fT (t ) (an cos nt bn sin nt ) (1.1) 2 n 1 为求出a0 , 计算[ fT ,1], 即
数学物理方程课件 积分变换法
设F[ f1(x)] F1(), F[ f2 (x)] F2 (),
则F[ f1(x) f2 (x)] F1() F2 ()
(5)
其中,为常数,逆变换也成立,即
F-1[ F1() F2 ()] f1(x) f2 (x)
(6)
试证明Fourier正弦变换和Fourier余弦变换的公式分别为
Fs1[Fs ()]
f (x)
2
0 fs (x) sin xdx
Fc1[Fc ()]
f
(x)
2
0 fc (x) cos xdx
§4.1.1 Fourier变换法
证明:F () F[ f (x)] f (x)eixdx
i
2
0
Fs
(
)
ei
x
d
(欧拉公式)
即Fourier正弦变换的公式为
f (x) 2
0 Fs () cos xd
§4.1.1 Fourier变换法
例9:证明
x 0 1 x2
sin xdx
2
e
(
0)。
证明:本题直接积分不易计算,考虑到fs
1 l
l l
f (x) cos n
l
xdx, n 0,1, 2,...
bn
1 l
l l
f (x) sin n
l
xdx, n 1, 2,...
§4.1.1 Fourier变换法
二、Fourier变换
设f (x)在(-, )上满足
i)逐段光滑(可导);
数学物理方程——8 积分变换法
下午9时10分
数学物理方法
第五章
积分变换法
拉普拉斯逆变换
1 σ + i∞ f (t ) = F ( p )e pt dp, 2π i ∫σ −i∞
p = σ + iω
又称 f (t )为 原函数 ⇔ F ( p )
为像函数
13
下午9时10分
数学物理方法
第五章
积分变换法
例2
(1) 求 L[1]
1 L[1] = ∫ 1 ⋅ e − pt dt = − e − pt 0 p
∞ ∞ 0
=
1 . p
(2) 求 L[t ]
1 ∞ 1 − pt ∞ 1 ∞ − pt 1 ∞ − pt 1 − pt L[t] = ∫ t ⋅ e dt = − ∫ t ⋅ d(e ) = − [t ⋅ e ] 0 + ∫ e dt = ∫ e dt = 2 . 0 p 0 p p 0 p0 p
− pt ∞
数学物理方法
第五章
积分变换法
1. Fourier变换 1.1 Fourier变换的定义
+∞ +∞
1 f ( x) = 2π
∫ ∫
−∞
(
−∞
f (τ )e −iωτ dτ )e iω x dω ,
(*)
傅里叶积分定理:设f 在 (−∞,+∞) 内满足下面两个条件:
+∞
(1)积分
−∞
∫
f ( x) dx 存在;
⎧ d 2U (ω , t ) t>0 = − a 2ω 2U (ω , t ), ⎪ ⎪ dt 2 ⎨ ⎪U (ω ,0) = Φ (ω ), dU (ω ,0) = Ψ (ω ), ⎪ dt ⎩ U (ω , t ) = A cos aωt + B sin aωt Ψ (ω ) B= U (ω , 0) = A = Φ (ω ) aω Ψ (ω ) U (ω , t ) = Φ (ω ) cos aωt + sin aωt aω
积分变换课后习题答案
积分变换课后习题答案积分变换是数学分析中的一个重要概念,它涉及到对函数的积分进行变换以简化问题或求解特定的数学问题。
以下是一些积分变换课后习题的答案示例:1. 习题一:求函数 \( f(x) = x^2 \) 在区间 \( [0, 1] \) 上的定积分的傅里叶变换。
答案:首先计算 \( f(x) \) 在给定区间上的定积分:\[\int_{0}^{1} x^2 dx = \left[ \frac{x^3}{3} \right]_{0}^{1} = \frac{1}{3}\]然后,根据傅里叶变换的定义,计算其傅里叶变换:\[F(k) = \int_{-\infty}^{\infty} f(x) e^{-2\pi i k x} dx\]由于 \( f(x) \) 只在 \( [0, 1] \) 上非零,我们可以将积分区间限制在这个区间内:\[F(k) = \int_{0}^{1} x^2 e^{-2\pi i k x} dx\]通过换元和积分计算,我们可以得到 \( F(k) \) 的表达式。
2. 习题二:证明拉普拉斯变换 \( \mathcal{L} \{t^n\} =\frac{n!}{s^{n+1}} \),其中 \( n \) 是非负整数,\( s > 0 \)。
答案:根据拉普拉斯变换的定义:\[\mathcal{L} \{f(t)\} = F(s) = \int_{0}^{\infty} e^{-st}f(t) dt\]对于 \( f(t) = t^n \),我们有:\[\mathcal{L} \{t^n\} = \int_{0}^{\infty} e^{-st} t^n dt\]通过分部积分法,我们可以逐步计算这个积分,最终得到:\[\mathcal{L} \{t^n\} = \frac{n!}{s^{n+1}}\]3. 习题三:求函数 \( f(x) = e^{-x} \) 在 \( x > 0 \) 时的傅里叶变换。
江苏省专转本高等数学第四章定积分核心知识点例题讲解(含答案)
第四章 定积分本章主要知识点● 定积分计算● 特殊类函数的定积分计算 ● 变限积分● 定积分有关的证明题 ● 广义积分敛散性 ● 定积分应用(1)面积 (2)旋转体体积一、定积分计算定积分计算主要依据牛顿—莱伯尼兹公式:设⎰+=C x F dx x f )()(,则()()()()bb aaf x dx F b F a F x =-=⎰。
其主要计算方法与不定积分的计算方法是类似的,也有三个主要方法,但需要指出的是对于第Ⅱ类直接交换法,注意积分限的变化:()111()()()()()(())x t bb aa t x f x dx f t t dt ϕϕϕϕϕϕ---=='=⎰⎰。
例4.1.111)edx x ⎰解:原式=e11)ln d x ⎰=32125((ln )ln )|33ex x +=例4.2.3dx ⎰ 解:原式t x t x =+-==11222 1121t tdt t -+⎰=32 121t t dt t -+⎰=322125()|33t t -= 例4.3.⎰22sin πxdx x解:原式=⎰-202cos 21πx xd =⎰+-2022cos 21|2cos 21ππxdx x x=20|2sin 414ππx +=4π 二、特殊类函数的定积分计算1.含绝对值函数利用函数的可拆分性质,插入使绝对值为0的点,去掉绝对值,直接积分即可。
例4.4.⎰--21|1|dx x解:原式=121 1(1)(1)x dx x dx --+-⎰⎰=212|)2(2x x -+=)121(02--+=25 例4.5.⎰--++22|)1||1(|dx x x解:原式=112211(|1||1|)(|1||1|)(|1||1|)x x dx x x dx x x dx ---++-+++-+++-⎰⎰⎰=112211(11)(11)(11)x x dx x x dx x x dx ------++++-+++-⎰⎰⎰=112211222xdx dx xdx ----++⎰⎰⎰=212122|4|x x ++---=)14(4)41(-++--=102.分段函数积分例4.6.⎩⎨⎧≤+>=0,10,)(2x x x x x f ,求⎰-11)(dx x f解:原式=⎰⎰-+0110)()(dx x f dx x f =⎰⎰-++01102)1(dx x dx x =103012|31|)2(x x x ++- =31)121(+--=65例4.7.⎩⎨⎧≤>+=1,1,12)(x x x x x f ,求⎰-+12)1(dx x f解:原式11221(1)()u x f x dx f u du =+--=+==⎰⎰1211()()f u du f u du -+⎰⎰1222111(21)0()udu u du u u -=++=++⎰⎰624=-=3.奇函数积分如果 ()f x 为定义在[],a a -的奇函数,则()0aaf x dx -≡⎰,这是一个很重要考点。
数学物理方法3-4积分变换法
§3.4.1
第三章 偏微分方程的定解问题 第四节 积分变换法
直线上的初值问题
例3.4.1求解热传导 问题
dU(, t) 2 2 a U(, t), t 0 解:利用傅立 dt 叶变换的性质 U(, 0) (), t a22 a22t C () U(, t) e C F(, ) e d
思考 利用积分变换方法求解问题的好处是什么?
第三章 偏微分方程的定解问题 第四节 积分变换法
傅立叶变换的定义
U ( , t ) u ( x, t )e
j x
1 dx , u ( x , t ) 2
U ( , t )e j x d
傅立叶变换的性质 微分性 位移性 f ( n ) (x) ( j ) n F ( )
e
d d
1 2a
t
( )e
2 x
4 a 2t
d
第三章 偏微分方程的定解问题 第四节 积分变换法
§3.4.2
半无界直线上的问题
半无界区域上的热传导(扩散)问题 2 u 2 u 0 x , t 0 t a x 2 0, 例3.4.4 求解 t 0 u (0, t ) u0 , u ( x, 0) 0, 0 x 做代换 u ( x, t ) v( x, t ) u0 转化为直线上热传导方程 2 v v 2 对称延拓法(奇延拓) a , 0 x , t 0 2 x t x0 u0 , v(0, t ) 0, t0 ( x) u0 , x0 v( x, 0) u0 , 0 x 考虑与无界区域上 波传播问题的差别
积分变换习题解答1-5
1.求微分方程()()(),()x t x t t t δ'+=-∞<<+∞的解. 分析:求解微分、积分方程的步骤:1)对微分、积分方程取Fourier 变换得象函数的代数方程; 2)解代数方程得象函数;3)取Fourier 逆变换得象原函数(方程的解).解:设()(),x t X ω⎡⎤=⎣⎦F 对方程两边取Fourier 变换,得 ()()j 1.X X ωωω+= 即()1.1X j ωω=+其逆变换为()0,.e ,0tt x t t -⎧<⎪=⎨≥⎪⎩4.求解下列积分方程: 1)()()()222210;y a b t bt aτττ+∞-∞=<<+-+⎰d2)()22tt y τττ+∞----∞=⎰ed .解:1)利用卷积定理可以求解此类积分方程.显然,方程的左端是未知函数()y t 与221t a+的卷积,即()221y t t a+.设()(),y t Y ω⎡⎤=⎣⎦F 对方程两边取Fourier 变换,有()222211y t t a t b ⎡⎤⎡⎤=⎢⎥*+⎢⎥⎣⎦⎣+⎦F F即()222211y t t a t b ⎡⎤⎡⎤⎡⎤⋅=⎣⎦⎢+⎥⎢⎥⎣⎦⎣⎦+FFF易知:22cos 2ttβωωβωβ+∞-=+⎰πd e,有()222211ttY t t t at bωωω+∞+∞---∞-∞⋅=++⎰⎰j j ed ed()2222cos cos 22t t Y t t t a t bωωω+∞+∞⋅=++⎰⎰d d所以()()22b b a a a bY ba ωωωω----==πe eπe由上可知222201cos π2d e a t t t a t a a ωω+∞-⎡⎤=⎢⎦=⎥++⎣⎰F ,()()-1b a a y t e b ω--⎥=⎡⎤⎢⎣⎦F()-1-b a a b a b b a ω--=⋅-⎡⎤⎢⎥⎣⎦F πe π()()22--a b a b t b a =⎡⎤+⎣⎦π.2)设()(),y t Y ω⎡⎤=⎣⎦F 对方程两边取Fourier 变换,同理可得()22etty t --⎤⎡⎤=⎥⎣⎦⎥⎣⎦FF利用钟形脉冲函数的Fourier变换224e et A ωββ--⎡⎤=⎣⎦F 及由Fourier 变换的定义可求得:222e tβββω-⎡⎤=⎣⎦+F ,从而()22e tt y t --⎤⎡⎤⎡⎤⋅=⎥⎣⎦⎣⎦⎥⎣⎦FFF即()()2222222121Y ωωωωω--==++πe πe()22222ωωω--=-πeπj e从而()()222-1-122y t ωωω--⎡⎤⎡⎤=-⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦πe πj eFF , 其中,记()22ef t ω-⎡⎤=⎣⎦F ,则()221tf t -=,上式中第二项可利用微分性质()()()()2222f t f t ωωω-''⎡⎤⎡⎤==⎣⎦⎣⎦F F j j e,则()()2222-12221tf t t ωω--⎡⎤⎛⎫''== ⎪⎢⎥ ⎪⎢⎥⎣⎦⎝⎭Fd j ed 222t -=因此()22222tty t --=-ππ22221t t -⎫=-⎪⎭e. 5.求下列微分方程的解()x t :()()()()d ax t b x f t ch t τττ+∞-∞'+-=⎰其中()(),f t h t 为已知函数,,,a b c 均为已知常数.解:设()()()()()(),,.f t F h t H x t X ωωω⎡⎤⎡⎤⎡⎤===⎣⎦⎣⎦⎣⎦F F F 对方程两边取Fourier 变换,可得()()()()j a X bX F cH ωωωωω+=即()()(),j cH X a bF ωωωω=+从而()()()()-1.12tcH X a bF x t ωωωωωω+∞-∞⎡⎤==⎣⎦+⎰Fj πed j。
积分变换课后答案.docx
1-11.试证:若 f t 满足Fourier积分定理中的条件,则有f t a cos td b sin td001f cos d , b 1sin d .其中 a fππ分析:由 Fourier 积分的复数形式和三角形式都可以证明此题,请读者试用三角形式证明 .证明:利用 Fourier积分的复数形式,有f t1f e j t e j t d2π11f cos j sin d e j t d2π1j b cos t j sin t da2由于 aa, b b, 所以f1a cos td1b sin td t22a cos tdb sin t d002.求下列函数的 Fourier积分:1)f1t 2 ,t 212)f0,t0 tt 2;t;0,1 e t sin 2t, t00,t13)f1,1t0 t0t11,0,1t分析:由 Fourier积分的复数形式和三角形式都可以解此题,请读者试用三角形式解 .解: 1)函数f1t 2 , t 21tt 2为连续的偶函数,其 Fourier 变换为0,1F () F [ f (t )] f (t)e j t d t2 f (t )cos tdt 21t 2 )cos tdt (1—sin t2t cos t2sin t t 2 sin t 1cos )4(sin(偶函2233数)f(t)的 Fourier积分为f (t )1 F ()e j t d1 F ()cos td2ππ 04(sin cos)tdπ 03cos2) 所给函数为连续函数,其Fourier变换为F ω F f (t ) f (t )e j t dt e t sin 2te j t dt0e t e2tj e 2tj e j t dt1[e( 1 2jj ) te (1 2j j )t ]d t2j2j1e( 1 2j j)t e (1 2j j)t2j 1 2j j 1 2j j0j11 2 5 2 1 (2)j 1 (2)j25 622 j24(实部为偶函数,虚数为奇函数)f (t)的 Fourier变换为f t1 F ()e j t d2π12522jcos t jsin t d2π25624152 cos t2sin t152 sin t 2 cos tπ25624dπ25 624d252 cos t2sin tπ 025624d这里用到奇偶函数的积分性质 .3)所给函数有间断点 -1 ,0,1且 f(- t)= - f(t)是奇函数,其 Fourier变换为F F f ( t ) f ( t)e j t dt2jf (t )sin tdt—2j1tdt2j(cos1)(奇函数)1 sinf(t)的 Fourier 积分为f ( t) = 1Fe j t djFsin td2π 0π 02 1 costdsinπ 0其中 t-1 , , (在间断点 f t 00 f t 0代替).0 1 t 0 处,右边 f(t)应以23.求下列函数的 Fourier变换,并推证下列积分结果:1) f te t(0), 证明:cos tπ t22 d2 e;t 2 2cosπe t cos t;2) f ( t)e cost ,证明:0 4td42sin t, tπsin πsint ππ 3) f ( t)sin t , t0,t,证明:12d2π0,tπ证明: 1)函数 f te t 为连续的偶函数,其 Fourier变换为FFf tete j t dt2e t cos tdte tcostsin tt22222 2t 0再由 Fourier 变换得ft 1 Fej t d122 costdt2ππ 02即costπ e t22d22)函数f t e t cos 为连续的偶函数,其Fourier 变换为tF ( )f t e j t dte t cos te j t dt—e t e j te j te j t dt21 0j )tdt 0e (1 jj )tdte ( 1 j ) t dte (1 j j ) t dte ( 1 jj21e (1 jj )te (1 j j ) te ( 1 j j ) te (1 j j)t2 1 jj 1 jj1 j j0 1 j j1 11112 2421 jj 1 j j1 j j1 j j4再由 Fourier 变换公式得1e j t d11 2f ( t)FFcos td42ππ 0π 02 2cos tdπe t即0 4cos t422c os td43)给出的函数为奇函数,其 Fourier 变换为Ff t e j tdt πtdt πt jsintdtsin te jsin t cosπππtdt j π1 t cos1 t dt2j sin t sincossin 1 t π sin1 t π jsinsin2jsinj111121F-1F 1 F ej td1 2jsin πcos t jsint d2π2π212 sin πsin tsin t, t ππ21d0,tπ故sin πsin tπ π2 sin t , t12d0, t π4. 求函数 f te t0,t 0 的 Fourier 正弦积分表达式和 Fourier 余弦积分表达式 .解:根据 Fourier 正弦积分公式,并用分部积分法,有2 fsin d sin tdf t2e t sin d sin tdπ 02 π 2 πesincos t22sin td0 22 sin td .根据 Fourier 余弦积分公式,用分部积分法,有2 fcosdcos tdf tπ 0 02 e t cosdcostdπ 0 02esincostcos tdπ 0222 2 2 cos td .π 01-21.求矩形脉冲函数 f (t )A, 0t变换 .0,的 Fourier其他解:ej tA 1 e j tf (t )e j tdtAe j tdtF ( ) Ff ( t)Ajj2. 设 F 是函数 f t 的 Fourier 变换,证明 F与 f t 有相同的奇偶性 .证明: F与 f t 是一个 Fourier变换对,即Ff t e jtdt , f t1 Fe j t d2π如果 F为奇函数,即 FF,则ft1 Fe jtd1 Fe j t d(令u )1 F u e jut du2π(换积分变量 u 为)1F e j t d f t2π所以 f t 亦为奇函数 .如果 f t为奇函数,即 f t f t ,则F f t e j t dt f t e j t d t(令 t u ) f u e j u du(换积分变量 u 为t)f t e j t dt F所以 F亦为奇函数 .同理可证 f t 与 F同为偶函数 .4.求函数 f t e t t0的 Fourier正弦变换,并推证sin2dπe012解:由 Fourier正弦变换公式,有F s () F s f t f t sin t dt e t sin tdt00e t sin t cos t12012由 Fourier正弦逆变换公式,有f t F s1F s ()2F s ()sin td2sin2t dπ 0π 01由此,当 t0 时,可得sind ππ2f e1225.设 F f t F ( ) ,试证明:—1) f t为实值函数的充要条件是 F () F () ;2) f t为虚值函数的充要条件是 F () F () .证明:在一般情况下,记 f t f r t j f i t其中 f r t和 f i t 均为 t 的实值函数,且分别为f t的实部与虚部 .因此F f t e j t dt f r t j f i t cos t jsin t dtf r t cos t f i t sin t dt j f r t sin t f i t cos t dtRe F j Im F其中 Re F f r t cos t f i t sin t dt ,aIm F f r t sin t f i t cos t dt b1)若 f t 为 t 的实值函数,即 f t f r t , f i t0. 此时, a 式和 b 式分别为Re F f r t cos tdtIm F f r t sin tdt所以F Re F jIm FRe F jIm F F反之,若已知 F F,则有Re F jIm F Re F jIm F此即表明 F的实部是关于的偶函数;F的虚部是关于的奇函数.因此,必定有F f r t cos tdt j f r t sin tdt亦即表明 f t f r t 为 t 的实值函数 . 从而结论 1)获证 .2)若 f t 为 t 的虚值函数,即 f t j f i t , f r t0 . 此时, a 式和 b 式分别为Re F f i t sin tdtIm F f i t cos tdt所以F Re F jIm FRe F jIm FRe F jIm FF反之,若已知 F F,则有Re F jIm F Re F jIm F此即表明 F的实部是关于的奇函数; F的虚部是关于的偶函数 . 因此,必定有F f i t sin tdt j f i t cos tdt ,亦即表明 f t jf i t为 t 的虚值函数 . 从而结论 2)获证 .6. 已知某函数的 Fouriersin,求该函数 f t .变换 F ( )sin解: F ( )为连续的偶函数,由公式有f tπF e j t d1sin cos td2π 01sin 1t1sin 1td 2π 0d2π 0但由于当 a0 时sin a sin ad( a)sin tdtπ0d0t2 0当 a0 时sin a sin(a)d π0d021,12t当 a0时,sin a0,所以得f1,104,10t7.已知某函数的Fourier变换为Fπ δδ,求该00函数 f t .解:由函数δ t t0g t dt g t 0,易知f t1F e j t d2π1πδ0e j t d1πδ0 e j t d2π2π1 e j t 1 e j t cos0t20208.求符号函数(又称正负号函数)sgn t1, t0变换 . 1, t的 Fourier解:容易看出 sgn t u t u t,而 F[u(t )] F ( )1jπδ( ).9.求函数 f t 1aδ taδ ta的 Fourier δ t a δ t222变换 .解:FFf t1 δ t aδ t aδ ta δ ta e j t d2221 e j tae jte jtae j tta2tt at22cosaacos.210 . 求函数 f t cost sin t 的 Fourier 变换 .解: 已知F sin0 tj π δδ由 ftcost sin t1sin 2t 有 Ff tπj δ2δ2 2211. 求函数 f t sin 3 t 的 Fourier 变换 .解: 已知 Fej 0t2πδ, 由e jt e jt3f tsin 3 t2jje 3j t 3e jt 3e -j t e 3j t8即得F f tπjδ3 3δ1 3δ1 δ3412. 求函数 f tsin 5tπ的 Fourier变换 .3解: 由于f tsin 5tπ 1sin5 t3cos5t322故 Ff tπj δ5 δ53πδ5 δ5 .2214. 证明:若 F ej tF ,其中t为一实数,则F cost1 FF2F sin t1 FF2j其中 F 为 F的共轭函数 .证明:因为 Fe jte j t dt F e j t e jtdte jte j t dt1 F Fe jte j te j t d tcos t e j t dt F cos t22同理可证另一等式 .17.求作如图的锯齿形波的频谱图 . (图形见教科书) .1t T解 : 02π, f tht ,0TT0,其他1 C 0TT 1T 1 hf t dtTTht dt0 021C nF nTTj n 0tdt1 Thtjn 0tdth Tj n 0 tdtf t eTeT 2teTh 1 ejn 0t T1 T j n 0 tdtj h T 2j n 0j n e2n π0 0Fh2πδj h 2πδ n 0πh δ2n2n πnn 0n 0j h δn 0 .n1- 31.若 F 1( ) F [ f 1( t )], F 2 ( )F [ f 2 (t )],, 是常数,证明(线性性质):Ff 1 (t ) f 2 (t )F 1 ( )F 2 ( )F -1F 1 ( )F 2 ( )f 1 (t )f 2 (t )分析:根据 Fourier 变换的定义很容易证明 .证明:根据 Fourier 变换与逆变换的公式分别有Ff 1 (t )f 2 (t )f 1 ( t) f 2( t ) e j t dtf 1( t )e j tdtf 2 (t )e j tdtF 1 () F 2 ( )F-1F 1 ( )F 2 ( )1 F 1 ( )F 2 ( ) e j t d2π1 F 1 ( ) e jtd1 F2 (e jt d2π2π)f 1( t ) f 2 ( t)6.若 F ( ) F [ f (t)] ,证明(翻转性质): F () F [ f ( t )]分析:根据 Fourier 变换的定义,再进行变量代换即可证明 .证明: F[ f ( t )]ft ejtdt(令 tu )f u e judu(换 u 为 t )f t e jtdtF ()9.设函数 f t1, t 1 sin t π, 10, t,利用对称性质,证明:Ft0,.11证明: F [ f (t )]f t etdt1tdtje j11 costdt1sin t dt由对称性质: F[ f (t )] F,则 F [ F (t )]2πf, 有F[ F ( t)] Fsin t 2πftF sin t πf π, 1 t0,112.利用能量积分f t21 F2,求下列积分的值:d t2πd1)1cos xd x ; 2 )sin 4xd x ;x 2x 23)12 d x ;4)x 22 d x . 1x 21x 21cos x d x2sin 2x解: 1)x22dxx 2(令xsin t2t ) d t2t12Fsin t2πtd112ππd2π1422sin x d sin x1cos x2)x d xx 2x 2sin x2sin x cos x2d x d xx xπ 1sin t2d t 2tπ-π=π22112112 3)2 d x dt F d ,其1 x2 1 t 22π 1 t 2中F1212e j t d t2cos2td t 2πeπe1t1t01t2从而1d x 1πe2122dπ1e2π2πdππe02022 1x24)x 22 d x x 21 12 d x12 d x12 dx1x21x21x1x2arctan x ππ π π π222221- 41.证明下列各式:2) f1 t f 2 t f 3 t f 1 t f 2 t f 3 t;6)df1 t f 2 t df1 t f 2 t f 1 tdf 2 t ;dt dt dt10) f t u t td f分析:根据卷积的定义证明.证明: 2) f 1t f 2t f 3t f1 f 2t f 3 t df 1 f 3 u f 2 t u du df1 f 3 u f 2 t u du df1 f 2 t u d f 3 u duf 1 t u f 2 t u f 3 u duf1 t f 2t f 3t6)df 1 t f 2 t d f1 f 2 t ddt dtd f1 f 2 t d f1 tdf 2 t,d t d td1 tf 2 tdf 1 t f 2dfdt dtdf 1 t f 2d df1 t f 2 t.dt dt 10)f t u t f u t du t1,t tf d.0,t2.若 f1 t e t u t, f 2t sin tu t ,求 f 1t f 2t .注意:不能随意调换f1 t和 f 2 t的位置 .e t, t, f 2sin t , t0解:由 f 1t u t t sin tut e t,0,t00, t0所以 f1 t f 2 t f 2 t f 1 t f2f1 td要确定 f 2 f 1 t0 的区间,采用解不等式组的方法 . 因为0, f20; t0, f 1t0 . 即必须满足, 即t, 因此0tf1 t f 2 t f 2 t f 1 tf 2 f 1 t dttsin e de tte d sin(分部积分法)e te sin cos t210e te sin cos12121sin cos e t214 . 若 F1F f1t, F2F f 2t,证明 :F f 1 t f 2 t1F1* F 2 2π证明 : 1F1F21F1 u F2u du2π2π1F2u f1t e j ut dt d u2π1F2u f 1 t e j ut dt du2π1F2u e j ut f1 t du dt2π1 f 1 tF 2u e j ut du dt2π1 f 1 tF 2s e jste j t ds dt2πf 1 te j tf 2 t dt F f 1 t f 2 t5. 求下列函数的 Fourier 变换:1) f t sin0tu t ;2) f t e t sin0tu t ;5) f te j 0t u t t 0;解: 1 )已知 Fu tπδ1, 又jf tsin0 tu t1 e j 0 tut e j0 tu t .2j由位移性质有Ff t1 πδ1πδ12jjjπ δδ2j0 22.2)由 Fourier 变换的定义,有Fe t sin0tu tetsin0tu t e jtdtsin 0 tejtd te j t j sin0 t0 cos0tj22 00 j 2 25)利用位移性质及 u t 的 Fourier 变换,有Fu t t 0e jt 0Fu te jt 01 πδj—再由象函数的位移性质,有F e j0t u t t0e j0 t 01πδ0j07.已知某信号的相关函数R1e 2a,求它的能量谱密度 S,其4中 a0 .解由定义知S R e j d1e 2 a e j d4102ae j d 1e2a j4e4e d1 e2 a j01e 2a j4 2a j42a j0111a42a j2a j4a22 9.求函数 f t e t u t,0的能量谱密度 .解: 因为f t e t u t e t , t0, 0,t0t e t,tf t e tu0,t当0 时, f t f t0 的区间为 0,,所以R f t f t dt e t tee e 2 t dt e1e 2 t020dt1e2当0 时, f t f t0 的区间为,,所以R f t f t dte t etdte e 2 t dt e 1 e2te1 e2 1 e22因此, R1 e ,现在可以求得 f t 的能量谱密度,即2SRe j d1 ee jd21 0jdej2ed11j 01j2ejej1 1 12jj1221- 51.求微分方程 x tx t t ,( t ) 的解 .分析:求解微分、积分方程的步骤:1)对微分、积分方程取 Fourier 变换得象函数的代数方程; 2)解代数方程得象函数;3)取 Fourier 逆变换得象原函数(方程的解) .解:设 Fx tX, 对方程两边取 Fourier 变换,得j XX1.即X1j .1其逆变换为 x 0, t 0 tt , t.e 04.求解下列积分方程:1)yd12 0 a b ;2t 2ta 2 bt 22)e ty d2πe 2 .解: 1)利用卷积定理可以求解此类积分方程 . 显然,方程的左端是未知函数 y t与1的卷积,即 y t1 . 设 F y tY, 对方程两边取2a 22a 2ttFourier 变换,有F yt1F1t 2 a 22b 2t即Fy tF1F1t2a2t2b2易知:cos t dπ e t,有222Y1 e j t d t t2 1b 2 e j t d t t 2a 2即Y2 cos t d t2cos t d t 0 t 2 a 2t 2 b 2πe ba所以 Y2beb aπab2a e由上可知1cos tπ aFt 2a 220 t 2a 2dta e ,y tF-1 a e b aba b - a F-1π e bab abπa b - a2 .π t 2b - a b2)设 F y tY , 对方程两边取 Fourier 变换,同理可得t2 F y t e t F2πe 2利用钟形脉冲函数的Fourier 变换F Ae t 2πAe24及由 Fourier变换的定义可求得: F e t 2,从而22F y t F e tt 2 F2πe 2即2πe Y2122 2π 12e2 2πe 从而22π j22e222y t πF -1 e 2πF-1j 22,e21t 2其中,记 F f t e 2,则f t e 2,上式中第二项可利用微分性质2π2F f t j 2f t j22 ,则F e22d2t 22F -1j f t1e t 1 ee2d t 22π2π2因此t2 2y tπ 1e2πt 2t212eπ2πt222πt 2t 21 e 2.25.求下列微分方程的解 x t :ax t b x f t d ch t其中 f t ,h t 为已知函数,a,b, c均为已知常数 .解:设F f t F,F h t H,F x t X. 对方程两边取—Fourier 变换,可得aj X bX F cH即cH X,ajbF从而x t F -1X 1 cH e j t d . π j 2 a bF2- 11.求下列函数的 Laplace 变换,并给出其收敛域,再用查表的方法来验证结果 .1) f tsint.2分析:用 Laplace 变换的定义解题 .t t st1jjs t解: Ls tsinsindt22 dt22 e2j 0e e1 112Re(s) 0 .2j sjs j4s 21222) f t e 2 t .解:2t2tst2 s t 1Lee edtedt s 2 Re(s).3) f t t 2 .解:22st1 2st1 t2eststLttedtst d(e)s 02te dt2td(est)2 te stestdt2 Re( s) 0 .223sss4) f t sin t cost .解: Lsin t cost0 sin t cos te st dt1sin 2te stdt2 0—1 21 .2 s 24 s 247) f t cos 2 t .解 : Lcos 2 tcos 2 te st dt1 cos2t e st dt21 1 cos2te st dt 1 1e 2j s t2s 2 0 2s4 01 111s 2 22s4 s 2js 2js s 2 .42.求下列函数的 Laplace 变换:3, 0 t21) ft1, 2 t 4.0, t4st2st4st解: Lf tf t edt 3 0 e dt2 edt3st 2 1 st 4 12 s4ss e 0s e 2s 3 4ee .3, tπ2) f2 . tcost , tπ2e st dtπstdtπcoste st dt解: Lf tf t 3 2 e23πj t -j t 3πsst 2πee e stdt1ee2s2s2e 2js tdtj s t j s t1 eeπ2 j sj s 2πsj s πj s ππsπs31e1 e 2e2 31 e1e Re(s) 022s 2 2s2 j sj ss13)f t2 tδ e5 t2t2 s te1 5e st t51Re( s) 2 .s2s24) f t cost δ t sin t u t解:Lf t0 δ t cos tu t sin te stdtcostest1j tj tstdtte ee2j111 1 112j j s j s 1 s 22- 21.求下列函数的 Laplace 变换式: 1) f t t 2 3t2 .解 : 由 Lmm ! 及 L1有 L2ts m11st3t22) ft1 te t .1t1t解 : Lts 2 ,Ltes+ 12 ,Lte3) f tt2e t .1 解:Lt-1 2Lt 2e t 2te t e te t2 32 21 s2 4s3 5 .s-1s-1s- 1s-1δ t coste st dtsin te st dtj s tj s t11 ee2j j sj ss 2Re(s).1 s 22 3 2 . s 3 s 2 s112 .ss- 15) f t t cosat .解 : 由微分性质有 :t cosatd d s s 2 a 2 LdsLcosata 22 2ds s 2s2 a6) f t 5sin2 t 3cos2 t解 : 已知 L sin t2 , Lcos tss2 s 22 , 则Lsin 2t3cos 2t5 22 10 3s434s 24s 2s 2 8) f te 4 t cos4t .解 : 由 L cos4t=s及位移性质有4tLecos4t s 216s+ 42( s +4).163.若 Lf tF s ,证明(象函数的微分性质) :F nsnLt n f t , Re s c1特别地, Ltf tF s ,或 f t1 L 1F s,并利用此结论计t算下列各式:1) f tte 3t sin2 t ,求 F s .3 t2解: L e sin 2t222s + 32,s + 34Lte3t d 2 2 2 s34 s3sin 2td s22222s+ 34s+ 34s+ 342) f t tt 3tsin 2tdt ,求 F s .et3 t1 3t1 2解: Lesin 2tdtsLe sin 2ts s 3 2 4 ,t3t22 3s 2 12s 13eL tsin 2tdt23 22s s 3 4s 2 s 43) Fs ln s1,求 f t .s 1解: F slns1, 令 LF sf t ,s 1F ' s2 s 11 s 11tts 21 Le e L tf tLtf t故 LF s f2sinh tt.t4.若 L f tF s ,证明(象函数的积分性质) :f tF s ds ,或 f ttL1F s dsLs ts并利用此结论计算下列各式:1) ftsin kt ,求 F s .t解 :L sin ktk,s 22k2k 2sLsin kts s 2k ds1 2dsarctansπarctan stk 2s 1skk s2kk2)f te 3 t sin 2t ,求 F s .t解 : Le 3 t sin 2ts2 2 ,3 4Le 3t sin 2t2dsπ arctans 3tss324 222- 31.设 f 1 t , f 2 t 均满足 Laplace 变换存在定理的条件 (若它们的增长指数均 为 c ), 且 L f 1 tF 1 s , Lf 2 tF 2 s ,则 乘积 f 1 tf 2 t 的Laplace 变换一定存在,且L f 1 t f 21t2πj—jj F1 q F2 s q dq其中c,Re s c.证明:已知 f 1 t , f 2t 均满足 Laplace变换存在定理的条件且其增长指数均为 c ,由Laplace变换存在定理知f1t f 2t也满足 Laplace变换存在定理的条件且f1 t f2 t f 1t f 2 t M e ct M e ct M 2 e2ct ,0t表明f 1 t f2 t 的增长指数为 2c . 因此 f 1t f 2t的 Laplace 变换F sf1 t f 2 t e st dt在半平面 Re s2c 上一定存在,且右端积分在Re s c c 上绝对且一致收敛,并且在Re s2c 的半平面内, F s 为解析函数 .根据 L f1t F1s ,则 f1 t 的 Laplace 反演积分公式为1jf1t2πj jF1 q e qt dq从而L f 1 t f 2 t f1 t f2 t e st dt1jst dtF1 q e q t dq f 2 t e 2πj j1(交换积分次序)2πjjf 2 t e s q t dt dq F1 qj01 2πjjj F1 q F2s q dq2.求下列函数的 Laplace 逆变换(象原函数);并用另一种方法加以验证.11)F s s2a2.2) F sss .a s b—3) F ss c2 .s as b10) F s s.s21 s24解: 1 ) L111sin at .s 2 a 2a2)s1ab ,s a s ba b s as bL1a sb1 b e at a e bt b . s s a3) F ss cc a1 1 b c 1 ,22s a s bb a2s a s bb as b故L 1s c b 2c a 2 e at b c t a c 2 e bts a s b ab aa b10)由 F ss 2s= 1s - s ,有1 s 243 s 2 1 s 2 4f tL1Fs1 cost cos2t.33.求下列函数的 Laplace 逆变换:1) F s12 .s246) F ss 21.lns 213) F s1e 2s. s 2解 : 1 )用留数计算法,由于s 1 2j, s 22j 均为 F s 的二级极点,所以f t L1F sL1 12Res Fst22s es 2js 2js s kk 1—de stlimde stlims 2j 2s 2 s2jdss2 jds2jst2 s 2jst2 s2jlimtee stlimtee st2424s2js 2js 2j s2 j s 2js2 jt e 2j t 8j e 2j t t e 2j t 8j e 2j t16 256 16256t e 2j t e 2j t1 e 2j t e 2j tsin 2t t cos2t82162j1686)令 F sln s 2 1 , F sln 2,s 2 s s 21F s1 s 12 L e t e t2L tf t,s 1 1 sL1 ln s2 1f t2cosht .s 21t13) L11 e2 sL11L1 e2 ss 2s 2s 22stt 2 u t22t 1, t 2t ,0 .t 22- 41. 求下列卷积:3) t m t n ( m, n 为正整数) .解: tmtntnkC n k t n k k dtdtnmm1k 0tnknkt1 m kdC n k t n k1 C n k t n k m k dk 0knk t m k 1 tn knkk t m n 1k1 m kC n1C n m k 1k 01km ! n ! t m n 1 m n 1 ! .注: 本小题可先用卷积定理求出t mt n 的 Laplace 变换 , 再由 Laplace 逆变—换求出卷积6) sin kt sin kt k0 .解:sin kt sin kt t1 sin k sin k t d02tcoskt cos 2k kt d 01t cos kt1t cosk 2 t d 2 t k 24k01t coskt sin2t k24k7) t sinh t解: t sinh t sinh t t sinh t dt0t1sin ktt coskt.22k1td 1td2e t2e t 001 2t1t t d(e )1t tt d(e )2t e e sinh t t 0202009)u t a f t a 0.解: u t a f t t a f t d 0,t au t.0 f t d , t aa10)δ t af t a0 .解:当 t a , δ t a f t0 .当 t a ,δ t a f t t a f t dδδ a f t d f t f t a .a2. 设 L f t F s,利用卷积定理,证明 : L tF s f t d t0s证明: L f t u t L f t L u t F s 1 ,s—L f tu tLt f tdLt dLt u f tf t dt3.利用卷积定理,证明 : L1s2tsin at .s2a 22a证明 : F sss1, 由2 2a2s 2a 2s2 a 2sL1s 2 s cosat , L111sin ata 2s 2a 2a有f tL 1 Fs 1 sin at1tcosa tdcosat asin aa1tsin atsin2aatd2a 0t sin at1 ttd a 2t2a4a2 sin a 2t sin at12 cosa 2 t2a4a tt sin at2a2-51. 求下列常系数微分方程的解:1) y y e 2 t , y 0 0 ;8) y3 y 3 yy 1, y 0 y 0y 0 0 ;12) y 4 2 yy 0, y 0 y 0y 0 0, y 0 1 ;16) y y 10sin 2t , y0 π 1 。
积分变换法
特别的,
f (x) (x)dx f (0)
(2) 对称性: (x) 为偶函数,则有
特别的,
(x x0 ) (x0 x) (x) (x)
自然也有
f (x) (x0 x)dx f (x0 )
7
例1 求函数 (x a) 的傅里叶变换,其中 a 是与
自变量 x 无关的数。
解 由定义知
F[ f (x)ei0x ] fˆ( 0 ) 傅里叶变换
L[ f (t)eat ] F (s a) 拉普拉斯变换
(6) 延迟定理
对变换的自变量而言
若 fˆ () F[ f (x)], F(s) L[ f (t)], 则有
F[ f (x x0 )] fˆ()eix0 傅里叶变换
L[ f (t t0 )u(t t0 )] F (s)est0 拉普拉斯变换
fˆ () F ( f ) f (x)eix dx
f (x) F 1 ( fˆ ) 1 fˆ ()eix d.
2
F (s) L( f ) f (t)est dt. 0
拉普拉斯逆变换记为
f (t) L1 (F (s)),
可用留数定理求得:设F(s) 除在半平面 Re s c内
20
例1 求解下列问题的解 ut a 2uxx f (x,t) ( x , t 0), (37)
u |t0 (x).
(38)
对方程(39)两端关于 t取拉氏变换,并结合条件
(40)得
sU (, s) () a22U (, s) G (, s),
U
(, s)
s
1
2a 2
()
s
1
2a 2
s 2U (s) k 2U (s) f (s)
高等教育出版社积分变换第四版答案
a cos t d
0
b sin t d
0, t0 2) f t ; t e sin 2 t , t 0
2.求下列函数的 Fourier 积分:
1 t 2 , t 2 1 1) f t ; 2 0, t 1
1 a jb cos t j sin t d 2
由于 a a , b b , 所以
f t 1 1 a cos t d b sin td 2 2
f(t)的 Fourier 积分为
f (t ) =
1 j F e j t d F sin td 2π 0 π 0 2 1 cos sin td π0
f t0 0 f t0 0 2
其中 t -1,0,1(在间断点 t0 处,右边 f(t)应以
1-1
1. 试证:若
f t 满足 Fourier 积分定理中的条件,则有
f t
0
a cos t d
0
b sin t d
其中 a
1 1 f cos d , b f sin d . π π
2 2 0 4
4
π sin π sin t sin t , t π sin t , t π 3) f ( t ) ,证明: d 2 2 0 1 t π 0, 0, t π
t 证明:1)函数 f t e 为连续的偶函数,其 Fourier 变换为
《数理方程》积分变换法解析
x2
x2
1 p2
dU dx
2x p
x2 p3
.
而 u |x1 cos y
变为
U
x,
p
|x1
1
p p2
,
解常微分方程得
U x, p
1 3 p3
x3
1 p
x2
p 1 p2
1 3 p3
1 p
.
取拉普拉斯逆变换,得
L(t n )
n! pn1 , n 0,1,
u
|x
0
f
t.
思考:需要对哪一个自变量进行哪一种积分变换?
对 t 进行拉普拉斯变换,设
u x,t U x, p, f t F p
于是方程变为
a2
d 2U x,
dx 2
p
pU
x,
p,
U x, p |x0 F p
这是二阶常微分方程的边值问题,它的通解为
根据傅里叶变换的微分性质,
方程转化为
dU ,
t
2U , t
dt
U , t |t0 F
于是 U ,t F e2t .
为了求出原方程的解,下面对 U ,t 关于 进行
傅立叶逆变换.
U ,t F e2t .
再由边值条件 U x, p |x0 F p 可知,C = F(p).
U
x,
p
F
pe
p a
x
.
为求出 u(x,t), 需要对 U(x,p) 进行拉普拉斯 逆变换。
数学物理方程-第四章积分变换法
第四章 积分变换法积分变换法是求解偏微分方程的一种基本方法. 不仅如此,在自然科学和工程技术的许多领域也有着广泛应用. 本章介绍Fourier 变换在求解偏微分方程定解问题中的应用. 主要以一维热传导方程,一维波动方程及平面上的Laplace 方程为主. 对于高维情形,由于计算过程要复杂一些,故只做简单介绍,也不做过多要求.§4⋅1 热传导方程Cauchy 问题4.1.1 一维热传导方程Cauchy 问题 考虑如下问题2(,), , 0 (1. 1)(,0)(), (1. 2)t xx u a u f x t x t u x x x ϕ⎧-=-∞<<∞>⎨=-∞<<∞⎩ 下面利用Fourier 变换求解该定解问题.设0>β为常数,函数2x e β-的Fourier 变换为()224()x F e e ωββπωβ--=(1.3)为书写方便起见,引入记号ˆ()(())(),f F f x ωω=,如果f 为二元函数),(t x f , ),))(,((),(ˆt t x f F t fωω=表示对),(t x f 中的空间变量x 作Fourier 变换的像函数,此时t 作为参数对待.对(1.1)—(1.2)关于空间变量x 作Fourier 变换得22ˆ(,)ˆˆ(,)(,), 0ˆˆ(,0)() .du t a u t f t t dtuωωωωωϕω⎧+=>⎪⎨⎪=⎩ 上面是一阶线性常微分方程的初值问题,解之可得2222()0ˆˆˆ(,)()(,)t a t a t ut e f e d ωωτωϕωωττ---=+⎰ (1.4) 利用(1.3)得),)( 21(22224t eta F e ta x t a ωπω--=),)()(21()(4)(2222τωτπττω--=----t e t a F e t a x t a记2241Γ(,)()2x a tx t eu t a πt-=(1.5)其中)(t u 为单位阶跃函数. 则有22ˆ((,))()(,)a t e F x t t ωωω-=Γ=Γ22()ˆ((,))()(,)a t e F x t t ωττωωτ--=Γ-=Γ-利用上面结果将(1.4)改写为ˆˆˆˆˆ(,)()(,)(,)(,)tu t t f t d ωϕωωωτωττ=Γ+Γ-⎰ (1.6) 对(1.6)两边取Fourier 逆变换,并利用Fourier 变换卷积公式 ))(()))((ˆ)(ˆ(21211x f f x f fF *=-ωω 便得0(,)()Γ(,(,)*Γ(,)tu x t x x t)f x x t d ϕτττ=*+-⎰()(,)(,)(,)t x t d d f x t d ϕξξξτξτξτξ∞∞-∞-∞=Γ-+Γ--⎰⎰⎰2222()()4()401()(,)2 2( )x x t a t a td ed f ed at at ξξττϕξξξτξππτ----∞∞--∞-∞=+-⎰⎰⎰(1.7)(1.7)即为定解问题(1.1)—(1.2)的解.在),(t x u 的表达式(1.7)中,函数(;)x t Γ起着一个基本作用. 如果令0≡f ,)()(x x δϕ=,则有(,)()(,)(;).u x t x x t x t δ=*Γ=Γ因此,(;)x t Γ是如下问题的解20, , 0 (1. 8)(,0)(), . (1. 9)t xx u a u x t u x x x δ⎧-=-∞<<∞>⎨=-∞<<∞⎩而(,)x t ξΓ-和(,)x t ξτΓ--分别是下面两问题的解20, , 0 (1. 10)(,0)(), . (1. 11)t xx u a u x t u x x x δξ⎧-=-∞<<∞>⎨=--∞<<∞⎩ 2()(),, 0 (1. 12)(,0)0, . (1. 13)t xx u a u x t x t u x x δξδτ⎧-=---∞<<∞>⎨=-∞<<∞⎩ 由于知道了(;)x t Γ,就可直接写出(1.1)—(1.2)的解(1.7)式. 类似于求解线性方程组0=Ax ,其中A 为n m ⨯矩阵. 如果知道该齐次方程组的一个基解组,则方程的任一解可由基解组的线性组合表出. 因此,),(t x Γ的作用就相当于向量空间中的基,故称),(t x Γ为定解问题(1.1)—(1.2)的基本解(fundamental solution).基本解是线性微分方程的一个很重要的概念,不仅可以表示Cauchy 问题的解,也可用来构造Green 函数表示边值问题的解.基本解有明确的物理解释. 若在初始时刻0t =时在0x =处置放一单位点热源,则此单位点热源在x 轴上产生的温度分布便是(,)x t Γ. 类似地,若在初始时刻0t =时在x ξ=处置放一单位点热源,则此点热源在x 轴上产生的温度分布为(,)x t ξΓ-. 而将初始时刻0t =变为t τ=时,其温度分布就是(,)x t ξτΓ--.注1 在(1.1)—(1.2)解的表达式(1.7)中,如果将其中的第一项和第二项分别记为1(,)u x t 和2(,)u x t ,则1(,)u x t 是相应于(,)0f x t =时齐次方程的解,而2(,)u x t 是相应于0)(=x ϕ时非齐次方程的解.若记1(,)()*Γ(,)(,)u x t x x t M x t ϕϕ==,则由齐次化原理可知20(,)(,)tf u x t M x t d τττ=-⎰.另外,和1(,)u x t 表达式中的卷积形式类似,2(,)u x t 也可表示成某种卷积形式,请同学们试给出这一表示形式. 例1.1 求解如下定解问题20, , 0 (1.14)(,0)(), . (1.15)t xx x u a u bu cu x t u x x x ϕ⎧---=-∞<<∞>⎨=-∞<<∞⎩ 其中,,a b c 均为常数.解 对(1.14)-(1.15)关于x 作Fourier 变换得22ˆ(,)ˆˆˆ(,)(,)(,), 0ˆˆ(,0)()dut a u t bi u t cu t t dtu ωωωωωωωϕω⎧=-++>⎪⎨⎪=⎩解之可得22() ˆˆ(,)().a bi c tut e ωωωϕω---= (1.16)为了求函数22()a bi c teωω---的Fourier 逆变换,利用配方法将其改写为222222224()()42.b a c bi t a t a bi c taaeeeωωω-------=由于222241()(),2x a t a tF eea tωωπ--=利用Fourier 变换的位移性质得000ˆ(())()()()() ,i x F f x e F f f ωωωωωω=-=- 取022biaω=得222222()4221()().2x bibi ixa t a ta aF ee ea tωωπ---=故有2222224()42((,))()b a c bi t a t aaF g x t eeωω----=22(),a bi c teωω---=其中22222244421(,)2b a c x bx taa ta g x t eeea tπ----=22()4.2x bt cta tee a tπ+-=记22()41Γ(,)()2x bt c ta tex t e u t a πt+-=其中)(t u 为单位阶跃函数. 1(;)x t Γ即为定解问题(1.14)—(1.15)的基本解.将(1.16)改写为1ˆˆˆ(,)()(,) ,u t t ωϕωω=Γ.,求Fourier 逆变换得1(,)()Γ(,)u x t x x t ϕ=*1()(,)x t d ϕξξξ∞-∞=Γ-⎰ 22()4() .2x bt cta te ed a tξϕξξπ-+-∞-∞=⎰如果将(1.15)中的齐次方程改为非齐次方程 ,考虑如下定解问题2(,),, 0 (,0)0, . t xx x u a u bu cu f x t x t u x x ⎧=+++-∞<<∞>⎨=-∞<<∞⎩请同学们写出该定解问题的解.例1.2 求解如下定解问题20, , 0(,0)(), .t xx u a u x t u x x x ϕ⎧-=-∞<<∞>⎨=-∞<<∞⎩ 其中0, ()0, .A x x x x x ϕ>⎧=⎨<⎩解 由(1.7)可得该问题的解为22220()()441(,)(),2 2 x x a ta tx A u x t ed ed at at ξξϕξξξππ----∞∞-∞==⎰⎰对积分作变量代换 2x a tξα-=得 02222202020(,) [][]2x x a t x x a t x x a t Au x t e d Aed e d Ae d αααααπααππαπ---∞----∞--==+=+⎰⎰⎰⎰引入下面函数22()xx e d ααπ-Φ=⎰(1.17)该函数称为误差函数. 利用误差函数可得(,)()222x x A A u x t a t-=+Φ. 4.1.2* 二维热传导方程Cauchy 问题为加深对线性微分方程基本解的进一步理解,下面再求解二维热传导方程Cauchy 问题222()(,,), (,)R , 0 (1.18)(,,0)(,), (,)R . (1.19)t xx yy u a u u f x y t x y t u x y x y x y ϕ⎧-+=∈>⎪⎨=∈⎪⎩ 为求解(1.19)—(1.20),先求二维热传导方程的基本解,即如下定解问题的解222()0, (,)R , 0 (1.20)(,,0)()(), (,)R . (1.21)t xx yy u a u u x y t u x y x y x y δδ⎧-+=∈>⎪⎨=∈⎪⎩引入二元函数的Fourier 变换12()12()(,)(,)i x y F f f x y e dxdy ωωωω∞∞-+-∞-∞=⎰⎰和一元函数Fourier 变换的性质相对应,二元函数的Fourier 变换也有类似性质.对(1.20)-(1.21)关于空间变量作Fourier 变换得22ˆ(,)ˆ(,)0, 0ˆ(,0) 1.dut a u t t dtuωωωω⎧+=>⎪⎨⎪=⎩其中2221212(,) , ωωωωωω==+. 解之可得22222212ˆ(,)a ta t a t ut e e e ωωωω---==.故有2212222211222222222()1121221244421(,,)()(,)(2)11=2211=221=.(2)a t i x y a t ix a t iy xy a t a tx y a tu x y t F f eed d eed e e d eea ta t e a t ωωωωωωωωωωωπωωπππππ∞∞-+--∞-∞∞∞---∞-∞--+-==⎰⎰⎰⎰即(1.18)-(1.19)的基本解为222421Γ(,,)().(2)x y a tx y t eu t a πt +-=与(1.7)相对应,(1.20)—(1.21)的解为(,,)(,)*Γ(,,)(,,)*Γ(,,)tu x y t x y x y t f x y x y t d ϕτττ=+-⎰(,)(,,)x y t d d ϕξηξηξη∞∞-∞-∞=Γ--+⎰⎰(,,)(,,).t d f x y t d d τξητξητξη∞∞-∞-∞Γ---⎰⎰⎰作为练习,同学们试用Fourier 变换求解三维热传导方程Cauchy 问题. §4⋅2 波动方程Cauchy 问题4.2.1 一维波动方程Cauchy 问题考虑如下定解问题2(,), , 0 (2.1)(,0)(), (,0)(), . (2.2)tt xx t u a u f x t x t u x x u x x x ϕψ⎧-=-∞<<∞>⎪⎨==-∞<<∞⎪⎩20, , 0 (2.3)(,0)0, (,0)(), . (2.4)tt xx t u a u x t u x u x x x ψ⎧-=-∞<<∞>⎪⎨==-∞<<∞⎪⎩ 若记(2.3)—(2.4)的解为(,)(,)u x t M x t ψ=,则由叠加原理和齐次化原理可得(2.1)—(2.2)的解为0(,)(,)(,)(,)t f u x t M x t M x t M x t d tτϕψττ∂=++-∂⎰ (2.5)因此,只须求解定解问题(2.3)—(2.4).对(2.3)—(2.4)关于空间变量x 作Fourier 变换得2222ˆ(,)ˆ(,)0, 0ˆˆˆ(,0)0, (,0)().t d u t a u t t dt u uωωωωωψω⎧+=>⎪⎨⎪==⎩ 解之可得sin ˆˆ(,)() .a tut a ωωψωω= 记1, 2(;) 0, .x atax t x at ⎧<⎪Γ=⎨⎪≥⎩查Fourier 变换表或直接计算可得sin ˆ((;))()(,)a t F x t t a ωωωωΓ=Γ= 故有ˆˆˆ(,)()(,),ut t ωψωω=Γ 对上式取Fourier 逆变换并利用卷积公式得(,)()*Γ(,)u x t x x t ψ=()(,)x t d ψξξξ∞-∞=Γ-⎰1()2x atx at d aψξξ+-=⎰ . 利用(2.5)便得(2.1)—(2.2)的解为0(,)(,)(,)(,)t f u x t M x t M x t M x t d t τϕψττ∂=++-∂⎰11(())()22x at x at x at x at d d t a aϕξξψξξ++--∂=+∂⎰⎰ ()0()1(,)2tx a t x a t d f d a τττξτξ+---+⎰⎰[]11()()()22x at x atx at x at d a ϕϕψξξ+-=++-+⎰ ()0()1(,)2tx a t x a t d f d a τττξτξ+---+⎰⎰ (2.6)当0f ≡时,(2.6)称为一维波方程Cauchy 问题的达朗贝尔(D ’Alembert )公式.注1 在(2.4)中取()()x x ψδ=,则有(,)(;)u x t x t =Γ,即(;)x t Γ是如下定解问题20, , 0(,0)0, (,0)() .tt xx t u a u x t u x u x x x δ⎧-=-∞<<∞>⎪⎨==-∞<<∞⎪⎩ 的解,称其为一维波动方程的基本解. 利用基本解(;)x t Γ,就可写出(2.1)—(2.2)的解(2.6)式. (;)x t Γ在(2.6)的表达式中也起到一个“基”的作用.4.2.2* 二维和三维波动方程Cauchy 问题下面,首先利用Fourier 变换求解三维波动方程Cauchy 问题,然后用降维法求出二维波动方程Cauchy 问题的解.考虑三维波动方程Cauchy 问题2333(,,,),(,,),0, (2.7)(,,,0)(,,),(,,), (2.8)(,,,0)(,,),(,,). (2.9)tt tu a u f x y z t x y z R t u x y z x y z x y z R u x y z x y z x y z R ϕψ⎧-∆=∈>⎪=∈⎨⎪=∈⎩为求解定解问题(2.7)—(2.9),先求出三维波动方程的基本解,即如下问题的解,23330, (,,), 0 (2.10)(,,,0)0, (,,) (2.11)(,,,0)()()(), (,,). tt t u a u x y z R t u x y z x y z R u x y z x y z x y z R δδδ-∆=∈>=∈=∈ (2.12)⎧⎪⎨⎪⎩记2222123123(,,) , ωωωωωωωω==++. 对定解问题(2.10)—(2.12)关于空间变量 作Fourier 变换得2222ˆ(,)ˆ(,)0, 0ˆˆ(,0)0, (,0) 1.t d u t a u t t dt uu ωωωωω⎧+=>⎪⎨⎪==⎩解之可得sin ||ˆ(,).||a tut a ωωω=故有123312331 ()1233()1233ˆ(,,,)()(,,,)1ˆ =(,)(2)sin 1 =(2)i x y z R i x y z R u x y z t F u x y z t u t e d d d a t e d d d a ωωωωωωωωωωπωωωωπω-++++=⎰⎰⎰⎰⎰⎰为计算上面积分,首先对上面积分作变量代换v A ωT T =,其中123(,,) , v v v v =A 为三阶正交矩阵. 选A 使得将(,,)x y z 变为(0,0,)r ,222 r x y z =++. 根据正交变换的保内积性可得,该变换将123 , x y z ωωωω++分别变为3 , v rv .故有33 1233sin 1(,,,)(2)i rv R a v t u x y z t e dv dv dv a v π=⎰⎰⎰,再利用球坐标变换123cos sin sin sin cos v v v ρθϕρθϕρϕ=⎧⎪=⎨⎪=⎩ 可得22 cos 3000cos 200 201sin (,,,)=sin (2) =sin (2) =sin() ()(2)i r i r i r i r a t u x y z t d d e d a i a t e d ar i a t e e d ar ππρϕπρϕρρρθρρϕϕπρρρπρρπ∞∞∞---⎰⎰⎰⎰⎰22sin() ()81( )()16i r i r i a t i a t i r ir i a t e e d ar e e e e d ar ρρρρρρρρπρπ∞--∞∞---∞=--=---⎰⎰. 注意到()(0)2()2()i i e d F e αραρωρπδωαπδα∞=-∞==-=⎰,221(,,,)( )()1612[()(())()()]161[()()]41().4 i a t i a t i r ir u x y z t e e e e d ar r at r at r at at r ar r at r at ar r at ar ρρρρρππδδδδπδδπδπ∞---∞=---=-⋅++-+----=--+=-⎰记1(,,,)()4x y z t r at arδπΓ=- (,,,)x y z t Γ即为三维波动方程的基本解.因此,当0==ϕf 时,(2.7)—(2.9)的解为33R R (,,,)(,,,) =(,,)(,,,)=(,,)(,,,)()=(,,). 4u x y z t M x y z t x y z x y z t x y z t d d d r at d d d arψψψξηζξηζξηζδψξηζξηζπ=*ΓΓ----⎰⎰⎰⎰⎰⎰其中222()()()r x y z ξηζ=-+-+-.对任一0t >, 记以点(,,)x y z 为心at 为半径的球面为(,,)r S x y z ,即3(,,){ (,,) }r S x y z R r at ξηζ=∈=. 将上面的积分化为累次积分并由δ函数的定义可得(,,)(,,)2(,,)2(,,,)(,,,)(,,)=()()4(,,)=4(,,) =41=(,,4rratS x y z r atS x y z S x y z u x y z t M x y z t r at ds dr ar ds ar ds a t a t ψψξηζδπψξηζπψξηζπψξηζπ∞==-⎰⎰⎰⎰⎰⎰⎰(,,)) . (2.13)at S x y z ds ⎰⎰最后,由叠加原理和齐次化原理便得(2.7)—(2.9)的解为22(,,)(,,)111(,,,)(,,)(,,)44at at S x y z S x y z u x y z t ds ds a t t a t ϕξηζψξηζππ⎛⎫∂=+ ⎪ ⎪∂⎝⎭⎰⎰⎰⎰211(,,,)4r f t d d d a r a ξηζξηζπΩ+-⎰⎰⎰ (2.14) 其中 (,,){ (,,) }a t B x y z r at ξηζΩ==<.(2.14)称为三维波动方程Cauchy 问题的克希霍夫(Kirchhoff )公式.利用Fourier 变换求二维波动方程的基本解比较难. 利用三维空间中已有的结果(2.13),下面用降维法求二维波动方程Cauchy 问题.考虑如下三维波动方程Cauchy 问题2330,(,,),0 (,,,0)0,(,,,0)(,),(,,) tt t u a u x y z R t u x y z u x y z x y x y z R ψ⎧-∆=∈>⎪⎨==∈⎪⎩(2.15)(2.16)对于定解问题(2.15)—(2.16),由于初始数据与z 无关,可推知解u 与z 也无关,故有zz u =0,即定解问题(2.15)-(2.16)其实是一个二维波动方程Cauchy 问题, 由(2.13)可得该问题的解为2(,,)2(,,)1(,,)(,)41=(,) (2.17)2at atS x y z S x y z u x y t ds a t ds a t ψξηπψξηπ+=⎰⎰⎰⎰其中22222(,,){(,,)|()()(),}atS x y z x y z a t z ξηζξηζξ+=-+-+-=≥. 对于上半球面(,,)atS x y z +直接计算得 2222221()() ()()ds d d atd d a t x y ζζξηξηξηξη∂∂=++∂∂=----将上式代入到(2.17)中便得222(,)(,,)(,,)1(,). (2.18)2at B x y u x y t x y t d d a a t rψξηξζπ=Γ=-⎰⎰其中22()()r x y ξη=---,(,){(,)|}at B x y r at ξη=<.和三维情形类似,由(2.18)可得二维波动方程Cauchy 问题2222(,,),(,),0 (2.19)(,,,0)(,),(,), (2.20)(,,,0)(,),(,). (2.21)tt tu a u f x y t x y R t u x y z x y x y R u x y z x y x y R ϕψ⎧-∆=∈>⎪=∈⎨⎪=∈⎩ 的解为222222(,)(,)1(,)1(,)(,,)22at at B x y B x y u x y t d d d d a t a a t r a t r ϕξηψξηξηξηππ∂=++∂--⎰⎰⎰⎰()222(,)1(,,)2()a t tB x y f d d d a a t r τξηττξηπτ---⎰⎰⎰(2.22) (2.22)称为二维波动方程Cauchy 问题的波以松(Poisson )公式.4.2.3 解的物理意义对一维波动方程Cauchy 问题,如果无外力作用,则解由D’Alembert 公式给出,即[]11(,)()()() .22x atx at u x t x at x at d aϕϕψξξ+-=++-+⎰ 将上式改写为(,)()() ,u x t f x at g x at =++-其中011()()() ,22x atf x at x at d a ϕψξξ++=++⎰11()()() .22x at g x at x at d aϕψξξ--=-+⎰ 记1(,)()u x t f x at =+,2(,)()u x t g x at =-,则12(,)(,)(,).u x t u x t u x t =+.首先考虑1(,)() .u x t f x at =+当0t =时1(,0)() .u x f x =在(,)x u 平面上画出函数()f x 的图形,则()f x at +的图形可通过()f x 的图形向左平移at 个单位长度而得. 随着t 的增加,()f x 的图形不断向左平移,移动速度为a ,故称1(,)u x t 为左传播波,a 为波速. 同样道理,2(,)()u x t g x at =-称为右传播波. D’Alembert 公式表明:弦线在t 时刻的振动是初始振动所产生的右传播波和左传播波的叠加.其次,从D’Alember t 公式还可看出:u 在(,)x t 的值(,)u x t 只与x 轴上区间[],x at x at -+上初始值有关,而与其它点的初始值无关. 这是由于波速为a ,在区间[],x at x at -+外的初始扰动在时刻t 还未传播到点x ,故称区间[],x at x at -+为点(,)x t 的依赖区间. 在(,)x t 平面上,过(,)x t 点分别作斜率为1a±的直线,两条直线在x 轴上所截得的区间便是[],x at x at -+(图2.1()a ).给定x 轴上的区间[]12,x x ,过点1(,0)x 作直线1x x at =+,过点2(,0)x 作直线2x x at =-,它们和x 轴构成了一个三角形区域(图2.1()b ).由于该区域内任一点的依赖区间都落在区间[]12,x x 内,因此,解在此三角形区域内的值完全由区间[]12,x x 上的初始值决定,而与此区间外的初始值无关,故称此三角形区域为区间[]12,x x 的决定区域. 同理,过点1(,0)x 作直线1x x at =-,过点2(,0)x 作直线2x x at =+,它们和x 轴构成一个梯形区域(图2.1()c ),该区域称为区间[]12,x x 的影响区域,它表示区间[]12,x x 上初始扰动对弦线振动的作用范围.t (x , t ) t t决定区域 影响区域x x x 0 x at - x a t + 0 1x 2x 0 1x 2x(a ) (b ) (c )图2.1由上面分析可得,波以常速a 沿两族直线x at c ±=向左﹑右两个方向传播,这是波动现象的一个基本特征. 直线x at c ±= 称为一维波动方程的特征线,它们在一维波动问题的研究中起着重要作用.当0f =时,对公式(2.14)和(2.22)进行分析,便可得到和上面类似的结论.对二维波动方程,一点(,,)x y t 的依赖区域是以(,)x y 为心,at 为半径的圆域;而对三维波动方程,一点(,,,)x y z t 的依赖区域是以(,,)x y z 为心,at 为半径的球面,而不是球形区域. 反映在波的传播过程中,平面波有前阵面而无后阵面,正像把一块石子扔在湖中,在湖面上激起层层浪花,这种现象称为波的弥漫现象;而空间波既有前阵面又有后阵面,正像人们听到声音,一会儿就消失了,这种现象称为空间波传播的无后效现象,此即Huygens 原理.§4⋅3 积分变换法举例在前二节中,利用Fourier 变换求出了热传导方程和波动方程Cauchy 问题的解. 下面再进一步举例,说明积分变换法在求解偏微分方程定解问题中的作用.例3.1 求解如下定解问题(,), , 0 (3.1)(,0)(), . (3.2)t x u au f x t x t u x x x ϕ+=-∞<<∞>⎧⎨=-∞<<∞⎩其中a 为实数.解 对(3.1)—(3.2)关于空间变量x 作Fourier 变换得ˆ(,)ˆˆ(,)(,), 0ˆˆ(,0)().du t ai u t f t t dtuωωωωωϕω⎧+=>⎪⎨⎪=⎩ 解之可得()0ˆˆˆ(,)()(,)tait ai t ut e f e d ωτωωϕωωττ---=+⎰ (3.3) 由于(())()ait F x at e ωδω--= ()((()))()ai t F x a t e τωδτω----=故(3.3)可表示为ˆˆˆˆˆ(,)()()()(,)(())()t u t x at f x a t d ωϕωδωωτδτωτ=-+--⎰对上式取Fourier 逆变换得0(,)()()(,)(()tu x t x x at f x x a t d ϕδτδττ=*-+*--⎰()()(,)(())t x at d d f x a t d ϕξδξξτξτδξτξ∞∞-∞-∞=--+---⎰⎰⎰()((),).t x at f x a t d ϕτττ=-+--⎰ 例3.2 求半平面上调和方程边值问题的有界解(,0)(), . (3.5)xx yy u x f x x ⎨=-∞<<∞⎩ 解 对(3.4)—(3.5)关于变量x 作Fourier 变换得222ˆ(,)ˆ()(,)0, 0ˆˆ(,0)().d u y i u y t dy u f ωωωωω⎧+=>⎪⎨⎪=⎩ 解之可得12ˆ(,)y y uy C e C e ωωω-=+ 由于u 有界,故20 .C =结合初始条件可得ˆˆ(,)() y u t f e ωωω-= (3.6) 直接求yeω-的Fourier 逆变换得11()()2y y ix F e x ee d ωωωωπ∞----∞=⎰1cos()y e x d ωωωπ∞-=⎰2201sin()cos()y x x y x e x y ωωωπ∞--=+ 221(,)yg x y x y π==+故(3.6)可表示为ˆˆˆ(,)() g(,)uy f y ωωω= 对上式取Fourier 逆变换得)))(,(*)((),(x y g f y x u ⋅⋅=() g(,)f x y d ξξξ∞-∞=-⎰221().()yf d x y ξξπξ∞-∞=-+⎰ 例3.3* 设有一单位长度均匀杆,侧面绝热,两端温度为零度.若初始温度为sin 2x π,求杆内的温度分布.解 设(,)u x t 为杆内温度分布,则u 满足如下定解问题(0,)(1,)0, 0 (3.8)(,0)sin 2, 0 1. (t xx u t u t t u x x x π==≥=≤≤ 3.9)⎪⎨⎪⎩对(3.7)—(3.9)关于时间变量t 作Laplace 变换,并记(,)u x t 的像函数为(,)u x s 可得222(,)(,)(,0)0(0,)(1,)0.d u x s su x s u x a dx u s u s ⎧--=⎪⎨⎪==⎩即2222(,)1(,)sin 2 (3.10) (0,)(1,)0 (3.11)d u x s s u x s x dx a a u s u s π⎧-=-⎪⎨⎪==⎩(3.10)是常系数二阶线性常微分方程,非齐次项为三角函数. 易得该方程 通解为1222sin 2(,)4s s x x aaxu x s C eC es a ππ-=+++利用边界条件(3.11)得10C =,20,C =故22sin 2(,)4xu x s s a ππ=+取Laplace 逆变换可得224(,)sin 2a tu x t e x ππ-=.例3.4* 求下面半无界弦振动问题有界的解2cos , 0, 0 (3.12) (,0)0, (,0)0, 0 (3.13)(0,)0, 0. tt xx t u a u t x t u x u x x u t t ρω-=>>==≥=≥ (3.14)⎧⎪⎨⎪⎩解 对(3.12)—(3.14)关于时间变量t 作Laplace 变换得222222(,)(,)(0,)0,d u x s s s u x s adx s u s u ρω⎧-=⎪+⎨⎪=⎩有界. 或者2222222(,)(1)(,)()(0,)0,d u x s s su x s dxa a s u s u ρω⎧--=⎪+⎨⎪=⎩有界. 解之可得1222(,)()s s x x aau x s C e C es s ρω-=+++由于u 有界,故20 .C =结合初始条件可得22(,)(1)()s x au x s es s ρω-=-+ (3.15)对(3.15)取Laplace 逆变换可得)()(),(221t s s t x u -⎪⎪⎭⎫ ⎝⎛+=ωρL )()(221t s s e x a s -⎪⎪⎪⎭⎫ ⎝⎛+--ωρL (3.16) 由于)()(221t s s -⎪⎪⎭⎫ ⎝⎛+ωρL =)()1(2221t s s s -⎪⎭⎫ ⎝⎛+-ωωρL =))(())(1(221212t s s t s --ωωρωρ+- L L 2222(1cos )sin 2tt ρρωωωω=-= (3.17) 利用Laplace 变换的延迟性质)()))(()((s f e s t u t f s τττ-=--L其中()u t 为阶跃函数. 取x aτ=得 )()(221t s s e x a s -⎪⎪⎪⎭⎫ ⎝⎛+-ωρL =)()()(221a x t u a x t s s ---⎪⎪⎭⎫ ⎝⎛+ωρL 22()2sin ()2xt x a u t aωρω-=- 22()2sin , 2 0, 0 . x t x a t a x t a ωρω⎧-⎪≥⎪=⎨⎪≤<⎪⎩(3.18)将(3.17)—(3.18)代入到(3.16)中便得222222sin sin () , 22(,)2sin , 0 . 2t x x t t a au x t t x t a ρωωωρωω⎧⎡⎤--≥⎪⎢⎥⎪⎣⎦=⎨⎪≤<⎪⎩ 注1 定解问题(3.7)—(3.9)也可用分离变量法求解. 一般而言,Laplace变换方法的求解过程比较繁琐,而分离变量法已成固定模式,求解过程相对简明.习 题 四1. 用Fourier 变换求解如下定解问题(1) 20, , 0, 2(,0)0, 2.t xx u a u x t A x u x x ⎧-=-∞<<∞>⎪>⎨⎧=⎨⎪<⎩⎩(2) 40, , 0, 1(,0) 0, 1t xx u u x t h x u x x -=-∞<<∞>⎧⎪⎧<⎨⎪=⎨⎪>⎪⎩⎩2*用Fourier 变换求解如下定解问题(1) 20, , 0 , 0(,0) 0, 0.t xx x u a u x t e x u x x -⎧-=-∞<<∞>⎪⎧>⎨=⎨⎪<⎩⎩(2) 2, , 0(,0)0, . t t xx u a u e x t u x x -⎧-=-∞<<∞>⎨=-∞<<∞⎩ 3. 用Fourier 变换求解如下定解问题(1) 2, , 0(,0)sin , .t t x u u xe x t u x x x -⎧+=-∞<<∞>⎨=-∞<<∞⎩(2) 23, , 0(,0)(), . t x u u u x t u x x x ϕ=+-∞<<∞>⎧⎨=-∞<<∞⎩4. 求解如下一维波动方程Cauchy 问题(1) sin , , 0(,0)0, (,0)0, . tt xx t u u t x x t u x u x x -=-∞<<∞>⎧⎨==-∞<<∞⎩(2) 22, , 0 1(,0)sin , (,0), . 1tt xx t u a u tx x t u x x u x x x ⎧-=-∞<<∞>⎪⎨==-∞<<∞⎪+⎩5*求解如下Cauchy 问题(1) 222()0, (,)R , 0(,,0), (,,0)=0, (,)R . tt xx yy t u a u u x y t u x y xy u x y x y ⎧-+=∈>⎪⎨=∈⎪⎩(2) 2222()0, (,)R , 0(,,0)0, (,,0)=, (,)R . tt xx yy t u a u u x y t u x y u x y x y x y ⎧-+=∈>⎪⎨=∈⎪⎩ (3) 2323()0, (,,)R , 0 (,,,0)0, (,,,0)=, (,,)R .tt xx yy zzt u a u u u x y z t u x y z u x y z x z x y z ⎧-++=∈>⎪⎨=∈⎪⎩6. 由三维波动方程Cauchy 问题解的公式,利用降维法求解如下问题20, , 0 (,0)0, (,0)(), .tt xx t u a u x t u x u x x x ψ⎧-=-∞<<∞>⎪⎨==-∞<<∞⎪⎩7. 考虑如下定解问题20, , 0(,0)(), (,0)(), .tt xx t u a u x t u x x u x x x ϕψ⎧-=-∞<<∞>⎪⎨==-∞<<∞⎪⎩ 设()x ϕ和()x ψ为直线R 上奇(偶,周期为T 的)函数,证明该问题的解(,)u x t 关于变量x 也是奇(偶,周期为T 的)函数. 对于一维热传导方程Cauchy 问题,类似结果是否成立?8*设()x ϕ和()x ψ在{0}x x ≥二阶连续可导,(0)(0)0ϕψ==,求解如下波动方程半无界问题20, 0, 0(0,)0, 0(,0)(), (,0)(), 0 . tt xx t u a u x t u t t u x x u x x x ϕψ⎧-=<<∞>⎪=≥⎨⎪==<<∞⎩如将该问题的边界条件换为 (0,)0, 0x u t t =≥,如何求解相应的定解问题?9.考虑如下定解问题000, , 0(), 0, .tt xx t t t u u x t u x u x ϕ==-=-∞<<∞>⎧⎪⎨==-∞<<∞⎪⎩其中初始波形为如下锯齿波1, 12()3, 230, .x x x x x ϕ-<<⎧⎪=-<<⎨⎪⎩其它(1)分别画出1,2t =时刻的(,)u x t 的波形图.(2)如果将初始位移换为1()()()x x x ϕϕϕ=--,分别画出1,2t =时刻的(,)u x t 的波形图.10. 考虑如下定解问题030, , 00, (), . tt xx t t t u u x t u u x x ψ==-=-∞<<∞>⎧⎪⎨==-∞<<∞⎪⎩其中2e , 13()0, .x x x ψ⎧<<⎪=⎨⎪⎩其它 试找出(,)u x t 恒为零的区域,又弦线上10x =-的点在那个时刻开始振动. 11. 考虑如下定解问题2200()0, (,), 00, (,), (,).tt xx yy t t t u u u x y R t u u x y x y R ψ==⎧-+=∈>⎪⎨==∈⎪⎩ 其中, (,)(,)0, .x y x y ψ∈Ω⎧=⎨⎩正值其它 若区域Ω为正方形{ (,) 1 1 , 1 1 }x y x y -<<-<<,试指出(,,10)u x y 恒为零的区域.12. 考虑如下定解问题3300()0, (,,), 00, (,,), (,,).tt xx yy zz t t t u u u u x y z R t u u x y z x y z R ψ==⎧-++=∈>⎪⎨==∈⎪⎩ 若(,,)x y z ψ除在球形域222{ (,,) (1) 1 }x y z x y z -++≤取正值外其它恒为零,试指出(,,,10)u x y z 恒为零的区域.13*求解下面定解问题21200, , 0(,0),0, .tt xx x t t u u u x t u x e u x -=-+=-∞<<∞>⎧⎪⎨==-∞<<∞⎪⎩ 14*考虑下面定解问题20, , 0 (,0)cos , . t xx u a u x t u x x x ⎧-=-∞<<∞>⎨=-∞<<∞⎩求出该定解问题解的有限表达形式.[利用结果2240cos ,0]4b ax ae bxdx ea aπ-∞-=>⎰.15*考虑下面定解问题230, , 0 (,0), . t xx u a u x t u x x x ⎧-=-∞<<∞>⎪⎨=-∞<<∞⎪⎩求出该问题解的有限表达形式.16*利用误差函数求解下面定解问题20, , 0 (,0)(), . t xx u a u x t u x x x ϕ⎧-=-∞<<∞>⎨=-∞<<∞⎩其中, 0(), 0.A x xB x ϕ>⎧=⎨<⎩。
工程数学积分变换答案
工程数学积分变换答案【篇一:复变函数与积分变换是一门内容丰富】建立和发展与解决实际问题的需要联系密切,其理论与方法被广泛应用在自然科学的许多领域,是机械、电子工程、控制工程,理论物理与流体力学,弹性力学等专业理论研究和实际应用中不可缺少的数学工具。
课程包含2部分内容:向量分析与场论,复变函数论与积分变换。
本课程的目的,是使学生掌握向量分析与场论,复变函数论,积分变换的基本理论、基本概念与基本方法,使学生在运用向量分析与场论,复变函数论,积分变换的思想和方法解决实际问题的能力方面得到系统的培养和训练,为在后继专业课程和以后的实际工作打下良好的数学基础向量分析与场论部分第一章向量与向量值函数分析学时:4几何向量,几何向量的加法、数乘、数量积、向量积,向量的混合积与三重向量积,向量值函数的定义,向量值函数的加法、数乘、复合、数量积运算,向量值函数的极限、连续,向量值函数的导数,向量值函数的体积分、曲线积分、曲面积分,高斯公式,斯托克斯公式。
第二章数量场学时:2数量场的等值面,数量场的方向导数、梯度的概念,哈米尔顿算子的用法。
第三章数量场学时:6向量场的向量线,向量场的通量,向量场的散度,向量场的环量,向量场的环量面密度、向量场的旋度,向量场场函数的导数与向量场的散度、旋度及数量场的梯度之间的关系。
第四章三种特殊形式的向量场学时:4保守场,保守场的旋度,保守场的势函数,管形场,管形场的向量势,调和场,调和函数。
复变函数与积分变换部分第一章:复数与平面点集学时:2复数的直角坐标表示法,三角表示法,指数表示法。
复数的模和辐角,复数的四则运算。
平面区域,邻域,聚点,闭集,孤立点,边界点,边界,连通集,区域,单连通区域,多连通区域。
第二章:解析函数学时:6复变函数的概念,复变函数的几何表示。
复变函数的极限,连续性,复变函数可导和解析的概念,复变函数解析的条件,复变初等函数(指数函数,对数函数,幂函数,三角函数)的定义和性质。
数学物理方程第四章 积分变换法(课堂PPT)
❖ 傅里叶变换建立R了信号时域与频域之间的关系,
频率是信号的物理本质之一。
6
❖ 设f(x)为[-π,π]上的有限信号,则f(x)的傅 里叶变换可简化为:
fˆ ( ) π f (x)eix dx π
❖ 对于只在有限区间,例如在上有定义的函数,可 采取延拓的方法,使其成为某种周期函数,而在 上,。然后再对作傅里叶级数展开,其级数和在 区间上代表f(x),由于f(x)在x=0和x=l无定义,因此 可以有无数种延拓方式,因而有无数种展开式, 它们在上均代表.有时,对函数在边界(区间的 端点)上的行为提出限制,即满足一定的边界条 件,这常常就决定了如何延拓。
第四章 积分变换法 傅立叶变换与拉普拉斯变换
数学物理方程
1
1777年以前,人们普遍采用多项式函数P(x)来对信 N 1
号f(x)进行表征:f (x) P(x) anxn。 n0 1777年,数学家Euler在研究天文学时发现某些函
数可以通过余弦函数之和来表达。1807年,法国科学
家傅里叶进一步提出周期为2π的函数f(x)可以表示为
( x ,t 0)
U ' (t; k) k 2a2U (t; k) F(t; k) U (t; k) |t0 0
其中 U (t; k) 为u(x,t)的傅里叶变换。为求解这个非齐次
e 常微分方程,用 k2a2t 遍乘方程各项 18
d [U (t; k)ek2a2t ] F (t; k)ek2a2t dt
19
❖ 交换积分次序
u(x,t) t
1
= 0
f ( , )[2
e e dk] k2a2 (t ) ik (x ) d d
引用积分公式
e2k2 ek dk =
数理方程参考答案4第四章 积分变换法
若 在 点连续,则
1
定义
设函数 f ( x) 在 (−∞, +∞) 上的任意有限区间上满足狄利克雷条件,在 (−∞, +∞) 上绝
对可积,则称广义积分
为
的傅里叶变换,或者称为 定义 称
的像函数。通常记为
,或
。
为
的傅里叶逆变换,或者称为 傅里叶变换及其逆变换的基本性质
的像原函数。记为
.
性质 1(线性性质) 傅里叶变换及其逆变换都是线性变换,即
其中 , 是任意常数。 性质 2(相似性质) 对于任意实常数 ,有 . 性质 3(位移性质)对于任意实常数 ,有 , 性质 4(微分性质)设 , 的傅里叶变换存在,则 . 一般地,若 , ,…, 的傅里叶变换存在,则 . 性质 5(乘多项式性质)设 的傅里叶变换存在,则
2
.
. 性质 6 (积分性质) . 性质 7 (对称性质) . 定义 于所有的 设函数 和 是 上定义的函数。 如果广义积分 对
2 ∂ 2u 2 ∂ u a − = 0 (−∞ < x < +∞, t > 0), ∂t 2 ∂x 2 ∂u u| ψ ( x). ( x), = = t =0 ϕ ∂t t =0
的解为
二维拉普拉斯方程的边值问题
∂ 2u ∂ 2u = 0 ( −∞ < x < +∞, y > 0), ∂x 2 + ∂ y2 u | = f ( x ), x =0 u = 0. |xlim |→+∞ 的解为
2
s2
例3 解
求函数 F ( p ) = 因为
p 的拉普拉斯逆变换 p − 2 p +5
第四章.积分变换法---求解偏微分方程
记作:F [ f ( x)] = f (k ) ,即
F [ f ( x)] = f (k ) = ∫
f f(x): (k ) 的傅里叶逆变换
∞ −∞
f ( x) e −ikx dx
记作: f ( x) = F −1[ f (k )] ,即
1 F [ f (k )] = f ( x) = 2π
−1
∫
∞
9
可以证明: 如果定义在 (−∞, ∞) 的函数在任一有限区间上满足 狄利克莱条件,且绝对可积( ∫ | f ( x) |dx 有界),则在
−∞ ∞
f(x)的连续点处,傅里叶积分存在:
1 f ( x) = 2π ⎡∞ ⎤ ikx −ikξ ∫∞ ⎢−∫∞ f (ξ )e dξ ⎥ e dk − ⎣ ⎦
——对于发生了任意位移x 0 的函数,其傅里叶变换 − ikx 等于 f(x)的傅里叶变换乘以一相位因子 e 0 证明:由定义:
F [ f ( x − x0 )] = ∫ f ( x − x0 ) e −ikx dx
−∞ u = x − x0 ∞
=
∫
∞
−∞
f (u ) e −ik (u + x0 ) du
频率域 波矢域
e − ikx
↔Leabharlann 412.1 傅里叶变换 一、傅里叶级数和复数形式的傅里叶级数 一个以2l为周期的函数f(x),若在区间[-l, l]满足 狄利克莱条件:(1)连续或只有有限个第一类间断 点;(2)只有有限个极值点,则 f(x) 在[-l, l]上可展开 为傅里叶级数
a0 ∞ nπ x nπ x + bn sin ) f ( x) = + ∑ (an cos 2 n =1 l l
积分变换(Fourier)课件与习题
的工程中使用的周期函数都可以用一系列的三角函数的
线性组合来逼近.---- Fourier级数
方波
4个正弦波的逼近
100个正弦波的逼近
4
研究周期函数实际上只须研究其中的一个周 期内的情况即可, 通常研究在闭区间[T/2,T/2]内 函数变化的情况.
T T fT (t )为T 周期函数,在 , 上满足 2 2 Dirichlet条件: fT (t )连续或仅有有限个第一类间断点; fT (t )仅有有限个极值点 则fT (t )可展开为Fourier级数,且在连续点t处成立: a0 fT (t ) an cos nt bn sin nt 2 n1
18
一般地, 对于周期T
1 T2 j n t cn T fT (t )e dt T 2 1 1 j n t e dt T 1 1 1 1 j n t j n j n e e e Tj n Tj n 1 2 sin n 2 sinc( n ) (n 0,1,2, ) T n T
cos nt
e
int
e 2
int
, sin nt
e
int
e 2i
int
6
级数化为: a0 e int e int e int e int an bn 2 n 1 2 2i a0 a n ibn int a n ibn int e e 2 n 1 2 2
1 从 而f (t ) f ( )cos (t )d d 2 1 可得 f (t ) f ( )cos (t )d d , 0 这就是f (t )的Fourier积分公式的三角形式。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
三、典型例题
例1 设 a 是正的实常数,证明
证明
a eiω x dω . 2 a +ω 利用傅里叶变换与逆变换的定义,有 e − a|x| =
π∫
1
∞
−∞
2
F [e − a|x| ] = ∫ e − a|x|e − iω x dx
−∞
∞
和
e − a|x| =
而
∞
1 2π
∫
∞
−∞
F [e − a|x| ]eiω x dω .
0
∞
= 2 Re
{∫ e
∞ 0
− a ( x +iω ) x
dx
}
1 = 2 Re a + iω 2a , = 2 a + ω2
将以上结果代入①式可得
e − a|x| =
例 2 求解热传导方程
π∫
1
∞
−∞
a eiω x dω . 2 a +ω
2
∂u ∂ 2u − a2 = 0 ( −∞ < x < ∞, t > 0) ∂t ∂x 2 的初值问题,已知初始温度为 u |t =0 = sin x 。 解 对方程和初始条件关于 作傅里叶变换,记
上分段连续,则 , 其中 是函数 的 在 上分段
的拉普拉斯变换变换存在, 且对 Re λ > α , 有 增长阶。一般地,如果 , ,…, 连续,则 在
上连续且不超过指数增长,
的拉普拉斯变换变换存在,且有
5
.
性质 5(积分性质)
性质 6(乘多项式性质) 性质 7 (延迟性质) 数 其中 是以下所定义的阶梯函
的解为
,…………………………. ② 其中
3
称为热核,也称为一维热传导方程初值问题的基本解。 定理 如果函数 在 时, 关于 上连续且有界,则由②给出的函数 无穷次连续可微。 是初值问题①
的古典解,且当
无限长杆上的热传导问题 ∂u ∂ 2u a 2 2 f ( x, t ) (−∞ < x < +∞, t > 0), = − ………………………③ ∂x ∂t u | = ϕ ( x). t =0 的解为 u ( x, t ) 定理 1 2a π t
都收敛,则称该积分为 与 的卷积。记为 .
性质 8(卷积性质) 或者 性质 9(像函数的卷积性质) 或者 性质 10(Parseval 等式) .
傅里叶变换的应用 一维热传导方程的初值问题
2 ∂u 2 ∂ u 0 (−∞ < x < +∞, t > 0), −a = …………………………① ∂x 2 ∂t u | = ϕ ( x). t =0
4
拉普拉斯变换及其逆变换的定义 定义 设函数 定义在 上,对于复数 ,如果含参变量 的广义积分
收敛,则称该积分为函数 的像函数。 定理 若函数 使得 数 的增长阶。 如果 原函数。它的形式是 在 ,则
的拉普拉斯变换,记为
,或
。也称 是
上分段连续且不超过指数型增长,即存在常数
和
,
的拉普拉斯变换对于满足 Re λ > α 的所有 都存在。 称常数 为函
其中 , 是任意常数。 性质 2(相似性质) 对于任意实常数 ,有 . 性质 3(位移性质)对于任意实常数 ,有 , 性质 4(微分性质)设 , 的傅里叶变换存在,则 . 一般地,若 , ,…, 的傅里叶变换存在,则 . 性质 5(乘多项式性质)设 的傅里叶变换存在,则
2
.
. 性质 6 (积分性质) . 性质 7 (对称性质) . 定义 于所有的 设函数 和 是 上定义的函数。 如果广义积分 对
2 2
∫
∞
−∞
e
−
s2 4 a 2t
sin( x − s )ds
2
s ∞ − s2 ∞ − 2 4a t 4a t sin cos d cos x sin sds e x s s e = − ∫−∞ ∫ −∞ 2a π t
2
sin x ∞ − 4 a 2t cos sds e = a π t ∫0 sin x 1 − a 2t 4π a 2t e = a πt 2 = e − a t sin x .
p
1
例 4 解常微分方程的初值问题
8
n 2π 2a 2 + T '' t T (t ) = f ( t ), ( ) l2 T ( 0 ) 0, = = T ' ( 0) 0 . 解 对 t 进行拉普拉斯变换 , 设 F (λ ) L = = f ( t ) , T ( λ ) L [T (t )] ,则
ˆ 利用傅里叶变 F [sin x ] = ϕ
换的微分性质,有
7
ˆ du 2 2 ˆ =0, t > 0, +a ω u dt ˆ. ˆ |t =0 = ϕ u
求解以上带参数 的常微分方程,得
对上式作傅里叶逆变换,并运用卷积性质得到
ˆ (ω, t )] u ( x , t ) = F −1 [ u = sin x * F −1[e − a ω t ] = 1 2a π t 1
∫
∞
−∞
ϕ (ξ ) e
在
−
( x −ξ )2 4a t
2
1 dξ + 2a π
∫ dτ ∫
0
t
∞
−∞
f (ξ ,τ ) − 4 a2 (t −τ ) e dξ ………..……④ t −τ 在 上连续且
( x −ξ )2
如果函数
上连续且有界, 是初值问题③的古典解。
有界。则由④给出的函数 一维波方程的初值问题
∫
b
a
K ( s, x ) f
( x )dx = F ( s )
定义了一个从 f ( x ) 到 F (s ) 的变换,称为积分变换, K (s, x ) 为变换的核。 在具体的变换中, K (s, x ) 具有不同的表现形式,对 f ( x ) 有具体的要求。常用的积分变 换有傅里叶变换和拉普拉斯变换。 定理(傅里叶积分定理) 若函数 f ( x) 在 (−∞, +∞) 上的任意有限区间上满足狄利克雷条 件,即满足: (1) 在区间上连续或只有有限个第一类间断点, (2) 在区间上至多有有限个极值点, 且在 (−∞, +∞) 上绝对可积,那么对任意 x ∈ (−∞, +∞) 有
二、知识总结
傅里叶变换及其逆变换的定义 定义 假设 I 是数集(实数或者复数),K (s, x ) 为 I × [a, b] 上的函数,这里 [a, b] 为任意区间。
如果函数 f (x ) 在区间 [a, b] 有定义并且对任意 s ∈ I , K ( s, x ) 和f ( x ) 为 [a, b] 上可积函数, 则 含参变量积分
性质 8(初值定理) 性质 9(终值定理) 性质 10(卷积性质)称 或者 拉普拉斯变换及其逆变换的应用 半无界热传导方程的定解问题
2 ∂u 2 ∂ u − = 0 ( x > 0, t > 0), a ∂t ∂x 2 u | = 0, t =0 u x =0 = f (t ).
解
(方法一) 拉普拉斯变换法。取
,记
,
,对于方程和
边值条件关于 作拉普拉斯变换,得 ˆ − ϕ ( x) − a 2 λu 即 ˆ λ d 2u 1 ˆ] = ˆ + 2 [ϕ ( x ) + f 0, − 2u 2 dx a a
9
ˆ ˆ d 2u = f , 2 dx
把 看作是参数,利用常数变异法, 得到它的通解是 ˆ( x, λ ) u = 1 2a λ x − ∫−∞ ϕ (ξ )e
故对 T ( λ ) 取拉普拉斯逆变换得
T ( t ) = L −1 T ( λ ) nπ a l t = nπ a ∫0 = l
(
)
f ( t ) ∗ sin
nπ a t l nπ a f (τ ) sin (t − τ )dτ . l
例5
求解下列无限长杆上的热传导问题
2 ∂u 2 ∂ u a f ( x, t ) (−∞ < x < +∞, t > 0), = − ∂x 2 ∂t u | = ϕ ( x). t =0
..................................... ①
F [e − a|x| ] = ∫ e − a|x|e − iω x dx
−∞ ∞
=
∫
−∞
e − a|x| cos ω xdx − i ∫ e − a|x| sin ω xdx
−∞
∞
= 2 ∫ e − ax cos ω xdx
λ ( x −ξ )
a
dξ + ∫ ϕ (ξ )e
x a ∞ ^
∞
λ ( x −ξ )
a
dξ
λ ( x −ξ )
a
− 1 x ^ + ∫−∞ f (ξdξ + ∫ f (ξ , λ )e
x
dξ .
由拉普拉斯变换表知 1 − L e λ
第四章 积分变换法
一、重点难点分析
分变换法是求解无界区域上定解问题的有效方法,其基本思想是通过函数的积分变换, 把微分运算转化为代数运算,从而减少偏微分方程中自变量的个数,将线性偏微分方程变为 含有较少变量的线性偏微分方程、常微分方程或代数方程,最终使问题得到简化。 学习本章应以下列两个方面为重点:一是两种积分变换——傅里叶变换和拉普拉斯变换 的定义、性质;二是两种积分变换法在求解无界区域上定解问题的应用。 本章的难点是傅里叶逆变换和拉普拉斯逆变换的计算,需要掌握相关的性质和一些基本 函数的积分变换。