数学物理方法大作业1
数学物理方法习题解答(完整版)
数学物理方法习题解答一、复变函数部分习题解答第一章习题解答1、证明Re z 在z 平面上处处不可导。
证明：令Re z u iv =+。
Re z x =，,0u x v ∴==。
1ux∂=∂，0v y ∂=∂，u v x y ∂∂≠∂∂。
于是u 与v 在z 平面上处处不满足C －R 条件， 所以Re z 在z 平面上处处不可导。
2、试证()2f z z=仅在原点有导数。
证明：令()f z u iv =+。
()22222,0f z z x y u x y v ==+ ∴ =+=。
2,2u u x y x y ∂∂= =∂∂。
v vx y∂∂ ==0 ∂∂。
所以除原点以外，,u v 不满足C －R 条件。
而,,u u v vx y x y∂∂∂∂ , ∂∂∂∂在原点连续，且满足C －R 条件，所以()f z 在原点可微。
()0000x x y y u v v u f i i x x y y ====⎛⎫∂∂∂∂⎛⎫'=+=-= ⎪ ⎪∂∂∂∂⎝⎭⎝⎭。
或：()()()2*000lim lim lim 0z z x y z f z x i y z∆→∆→∆=∆=∆'==∆=∆-∆=∆。
22***0*00limlim lim()0z z z z z z zzz z z z z z z z z=∆→∆→∆→+∆+∆+∆∆==+−−→∆∆∆。
【当0,i z z re θ≠∆=，*2i z e z θ-∆=∆与趋向有关，则上式中**1z zz z∆∆==∆∆】3、设333322()z 0()z=00x y i x y f z x y ⎧+++≠⎪=+⎨⎪⎩，证明()z f 在原点满足C －R 条件，但不可微。
证明：令()()(),,f z u x y iv x y =+，则()33222222,=00x y x y u x y x y x y ⎧-+≠⎪=+⎨+⎪⎩， 33222222(,)=00x y x y v x y x y x y ⎧++≠⎪=+⎨+⎪⎩。
数学物理方法大作业
基于分离变量法的波导中的电磁波研究1 空间当中的电磁波在迅变情况下，电磁场以波动形式存在，电磁场的基本方程是麦克斯韦方程组，对于在0==J σ情况下的迅变场，麦克斯韦方程组为]4[⎪⎪⎪⎭⎪⎪⎪⎬⎫=⋅∇=⋅∇∂∂=⨯∇∂∂-=⨯∇00B D t D H t B E (1)为了便于求解，通常将（1）式化为⎪⎪⎭⎪⎪⎬⎫=∂∂-∇=∂∂-∇010122222222t BcB t E c E (2) 必须指出的是，（2）式中第一式E 的三个分量X E ，y E ，z E 虽然是三个独立方程，但是其解却是相互关联的，因为（1）式到（2）式麦克斯韦方程变为二阶的麦克斯韦方程，故解的范围变大了。
为了使波动方程（2）的解是原方程（2）的解，必须是波动方程的解满足条件 0=⋅∇E 。
求解方程（1），即为求解⎪⎪⎪⎭⎪⎪⎪⎬⎫∂∂-=⨯∇=⋅∇=∂∂-∇t BE E t Ec E 0012222(3)(3)式在给定的边界条件下，可以求得定解. 对于定态电磁波，场量可以表示为t i e z y x E E ω-=),,( (4)考虑(4)式，(3)式可表示如下：⎪⎪⎭⎪⎪⎬⎫⨯∇-==⋅∇=+∇E iB E E k E ω0022(5)设电磁波为时谐波，并考虑到关系H B μ=，由（5）式可得到z y x ,,三个分量的6个标量方程：x y xH i E yE ωμγ-=+∂∂ (6) y x zH i E xE ωμγ-=-∂∂-(7) z xy H i yE xE ωμ-=∂∂-∂∂ (8) x y zE i H yH ωεγ=+∂∂ (9) y x zE i H xH ωεγ=-∂∂-(10) z xy E i yH xH ωε=∂∂-∂∂ (11) 以上6个方程经过简单运算，可以将横向场分量y x y x H H E E ,,,用两个纵向场分量z z H E ,来表示，即：)(12yE i x H k H zz cx ∂∂-∂∂-=ωεγ(12) )(12x E i y H k H zz cy ∂∂+∂∂-=ωεγ (13) )(12y H i x E k E z z cx ∂∂+∂∂-=ωμγ (14) )(12x H i y E k E z z cy ∂∂-∂∂-=ωμγ (15) 式中222k k c +=γεμω=kTM 波的纵向场分量与横向场分量关系[]1为：yE k i H zc x ∂∂=2ωε (12*) x E k i H zcy ∂∂-=2ωε (13*) xE k E zcx ∂∂-=2γ (14*)y E k E zcy ∂∂-=2γ (15*)TE 波的纵向场分量与横向场分量关系为[]1：xH k H zcx ∂∂-=2γ (12+)yH k H zc y ∂∂-=2γ （13+）yH k i E zc x ∂∂-=2ωμ （14+） x H k i E zcy ∂∂=2ωμ （15+） 2 波导内的电磁场 2.1波导的几个假设这里所讨论的波导,有以下假设:波导的横截面沿z 方向是均匀的,即波导内的电场与磁场只与坐标y x ,有关,与z 无关;构成波导壁的导体是理想导体,即∞=σ;波导内的介质各向同性,并且0=σ;波导内的电磁场为时谐场,角频率为ω。
数学物理方法习题解答(完整版)
数学物理方法习题解答一、复变函数部分习题解答第一章习题解答1、证明Re z 在z 平面上处处不可导。
证明：令Re z u iv =+。
Re z x =，,0u x v ∴==。
1ux∂=∂，0v y ∂=∂，u v x y ∂∂≠∂∂。
于是u 与v 在z 平面上处处不满足C －R 条件， 所以Re z 在z 平面上处处不可导。
2、试证()2f z z=仅在原点有导数。
证明：令()f z u iv =+。
()22222,0f z z x y u x y v ==+ ∴ =+=。
2,2u u x y x y ∂∂= =∂∂。
v vx y∂∂ ==0 ∂∂。
所以除原点以外，,u v 不满足C －R 条件。
而,,u u v vx y x y∂∂∂∂ , ∂∂∂∂在原点连续，且满足C －R 条件，所以()f z 在原点可微。
()000000x x y y u v v u f i i x x y y ====⎛⎫∂∂∂∂⎛⎫'=+=-= ⎪ ⎪∂∂∂∂⎝⎭⎝⎭。
或：()()()2*000lim lim lim 0z z x y z f z x i y z∆→∆→∆=∆=∆'==∆=∆-∆=∆。
22***0*00limlim lim()0z z z z z z zzz z z z z z zz z=∆→∆→∆→+∆+∆+∆∆==+−−→∆∆∆。
【当0,i z z re θ≠∆=，*2i z e z θ-∆=∆与趋向有关，则上式中**1z zz z∆∆==∆∆】3、设333322()z 0()z=00x y i x y f z x y ⎧+++≠⎪=+⎨⎪⎩，证明()z f 在原点满足C －R 条件，但不可微。
证明：令()()(),,f z u x y iv x y =+，则()33222222,=00x y x y u x y x y x y ⎧-+≠⎪=+⎨+⎪⎩， 33222222(,)=00x y x y v x y x y x y ⎧++≠⎪=+⎨+⎪⎩。
数学物理方法大作业1
目录一．实际现象的描述 3二．问题的求解4（一）求弦振动泛定方程 4（二）解弦振动方程 (6)Ⅰ.达朗贝尔法求“无限和半无限的”弦振动函数 (6)Ⅱ.分离变量法求两端固定弦振动方程 (7)三．各种情形下的弦振动求解及图像 (9)四．总结21一·实际现象的描述演奏者在演奏弦乐器（如二胡、提琴）时，用弓在弦上来回拉动，并通过另一只手指在按不同弦的不同地位的协调作用，奏出各种不同的美妙的音乐。
演奏者所用的乐器不同，奏出音乐的悦耳度也就不同。
演奏者虽然用弓所接触的只是弦的很小一段，似乎应该只引起这个小段的振动，而事实上，振动总是传播到整根弦。
这振动是怎样传播的呢？如何利用数学方法来求解这种物理问题？如何通过直观的方程来说明不同乐器演奏出的音乐效果不同的原因？可否利用matlab来将这种振动直观表示出来？通过对于弦振动方程的学习，及对matlab的初步了解，我对于不同定解问题下弦振动方程的求解做了初级小结。
也尝试利用matlab 直观表述不同定解条件下的弦振动动态图像。
二·问题的求解（一）求弦振动泛定方程在求解时，我们不妨认为弦是柔软的，就是说在放松的条件下，把弦完成任意的形状，它都保持静止。
由于弦乐器所用的弦往往是很轻的，它的重量只有力的几万分之一。
跟拉力相比，弦的重量完全可以略去，这样，真实的弦就抽象为“没有重量”的弦。
把没有重量的弦绷紧，它在不振动时是一根直线，就取这直线作为x轴。
把弦上各点的横向位移记作u。
这样，横向位移u是x和t的函数，记作u（x，t）。
要求解弦振动，首先应找出u所遵从的方程。
把弦细分为许多极小的小段，拿区间(x,x+dx)上的小段B为代表加以研究。
B既然没有重量而且是柔软的，它就只受到邻段A和C的拉力和。
弦的每小段都没有纵向（即x方向）的运动，所以作用于B的纵向合力应为零。
弦的横向加速度记作。
按照，小段B的纵向和横向运动分别为式中时弦的线密度，即单位长度的质量。
数学物理方法习题及解答1
数学物理方法习题及解答1试题1一、单项选择题1.复通区域柯西定理（）（A ）0)(=?dz z f l（B ）0)(1=∑?=n i l idz z f （C ）0)()(1=+∑??=ni l lidz z f dz z f （l 是逆时针方向，i l 也是逆时针方向）(D)0)()(1=+∑??=ni l lidz z f dz z f （l 是逆时针方向，i l 是顺时针方向）2.周期偶函数:,cos)(10为其中k k k a lxk a a x f ∑∞=+=π：（）（A ）?=lk d l k f l a 0cos )(1ξπξξ （B ）?-=ll k d l k f l a ξπξξcos )(1（C ） ?=lk k d l k f l a 0cos )(1ξπξξδ （D ）?lkk d lk f l a 0cos)(2ξπξξδ 3.柯西公式为：（）（A ）ξξξπd z f i n z f l ?-=)(2!)( (B) ξξξπd z f i z f l ?-=)(21)( (C) ξξξπd z f i z f l n ?-=)()(21)( (D) ξξξπd z f i n z f l n ?-=)()(2!)( 4.在00=z 的邻域上把()=z f 2zz )(sin 展开为（）（A ）+-+-!6!4!21642z z z(B) +-+-!7!5!31642z z z (C) +-+-6421642z z z(D) +-+-!7!5!31864z z z5.求()z z f sin 1=在z 0=πn 的留数为（）（A ）!1n （B ）n （C ）n )1(- （D ）16.以下那一个是第一类边界条件（）（A ）)(),(t f t x u ax == （B ）)(,()t f t x u ax n == （C ）)()(t f H u ax n u =+= （D ）lx ttlx xu Mg t x u ==-=),(7.下列公式正确的为：（A ）)()()(0x f dx x x f t =-?+∞∞-δ （B ）0)()(0=-?+∞∞-dx x x f t δ （C ）∞=-?+∞∞-dx x x f t )()(0δ （D ）)()()(0t t f dx x x f =-?+∞∞-δ8.勒让德方程为（A ）0)1(2)1(222=++--y l l dx dy x dx yd x（B ）0]1)1([2)1(22222=--++--y x m l l dx dy x dx y d x（C ）0)(22222=-++y dx dy x dx ym x d x（D ）0)(22222=+-+y dxdy x dx y m x d x9.m 阶贝塞尔方程为：（A ）0)(22222=--+R m x dx dR x dx R d x （B ）0)(22222=-++R m x dx dR x dx R d x （C ）0)(22222=+-+R m x dxdR x dx R d x （D ）0)(2222=-++R m x dxdR x dx R d x 上 10Z 0是方程W ‘’+P （Z ）W ‘+Q （Z ）W=0的正则奇点，用级数解法求解时，这个方程的“判定方程“为（A ）0)1(21=++---q sp s s （B ）0)1(21=++--q sp s s （C ）0)1(11=++---q sp s s （D ）0)1(22=++---q sp s s二、填空题1、已知解析函数22),()(y x y x u z f -=的实部，则这个解析函数为。
数理方法答案
i
1
+
1 1!3i
(z
− i)
−
2 2!9i2
(z
− i)2
+ ... +
1 11 k ! 3 3
−
1
1 3
−
2
...
1 3
−
k
+ 1 i−k
(z
− i)k
+ ...
R
=
lim
k →∞
ak ak +1
= lim k →∞
1 k!
1 3
1 3
( ) ( ) 解 z 3 = (x + iy)3 = x3 − 3xy 2 + i 3x2 y − y3 = ρ 3ei3ϕ = ρ 3 (cos 3ϕ + i sin 3ϕ )
其中
ρ = x2 + y2
ϕ
=
arctg
y x
(6) e1−i
解 e1+i = e × ei = e(cos1 + i sin1)
e x cosy
∴ v = −e x cos y + c
则f(z) = e x sin y − ie x cos y + ic = e x [sin y + i(− cos y)]+ ic
= e x( − i)[cos y + i sin y]+ ic = −ie x+iy + ic = −ie z + ic
故必须求出 ∂ 2u ∂x 2
= ... = 2F '+4x2 F"
∂2u ∂y 2
数学物理方法习题集
数学物理方法习题集第一章 复数与复变函数习题1，计算：（1），1)(1i ---。
（2），iii i 524321-+-+。
（3），5(1)(2)(3)i i i ---。
（4），4(1)i -。
（5），bi a +。
2，求下列复数的实部u 与虚部v ，模r 与幅角θ：（1），ii i i 524321----。
（2），1(2n+, 4,3,2=n 。
（3），i +1。
（4），3)i -。
（5），231i -。
3，设211i z +=，i z -=32，试用三角形表示21z z 及21z z 。
4，若21=+Z z θcos ，证明21=+m m zz θm cos 。
5，求下列复数z 的主幅角z arg ：（1），iz 312+-=。
（2），6)z i =-。
6，用指数形式证明：（1），(1)2i i -+=+。
（2），i ii2125+=+。
（3），7(1)8(1)i i -+=-+。
（4），1011(12(1)--=-。
7，试解方程44(0)z a a +=>。
8，证明：（1），1212Re()Re()Re()z z z z +=+ ；一般1212Re()Re()Re()z z z z ≠。
（2），1212Im()Im()Im()z z z z +=+ ；一般1212Im()Im()Im()z z z z ≠。
（3），2121z z z z = ；一般2121z z z z +≠+。
9，证明：（1），2121z z z z +=±。
（2），2121z z z z ⋅=。
（3），1122(z zz z = （02≠z ）。
（4），121212122Re()2Re()z z z z z z z z +==。
（5），()z z ≤Re ，()z z ≤Im 。
（6），2121212z z z z z z ≤+。
（7），222121212()()z z z z z z -≤+≤+。
数学物理方法课后答案 (1)
充分性。设任给ε > 0，存在N(ε ) > 0，使当n>N时，zn+ p − z0 < ε成立。由
xn+ p − xn ≤ （xn+ p − xn）（2 yn+ p − yn）2 = zn+ p − zn < ε
yn+ p − yn ≤ （xn+ p − xn）（2 yn+ p − yn）2 = zn+ p − zn < ε
②
将①式与②式相除，易见 c 3 = 1,即 c = 1，由此得证。
8.试利用 Re z = x ≤ x2 + y2 = z 证明 z1 + z2 ≥ z1 + z2 , z1 − z2 ≥ z1 − z2
证 将第一个不等式两边平方，则不等式右边的式子为
z1 + z 2 2 = ( z1 + z 2 )( z1 + z 2 )∗ = z1 z1∗ + z 2 z 2∗ + z1 z 2∗ + z1∗ z 2
4x ≥ 0
x≥0
3 见课上例题
4. 求在ω = 1 下，直线 Re z = c （常数）映射为什么图形？
z
μ ν 解：在ω = 1 的映射下, 与 满足什么方程？ z
右半平面（包括 y 轴）
由c = Re z = Re 1 w
= Re 1 μ + iν
=
μ μ2 +ν 2
∴c(μ 2 +ν 2 ) = μ
25(x2 − 6x + 9) + 25 y2 = 625 −150x + 9x2
16x2 + 25 y2 = (20)2 ，点集为椭圆： ( x )2 + ( y )2 = 1 54
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
目录一．实际现象的描述 3二．问题的求解4（一）求弦振动泛定方程 4（二）解弦振动方程 (6)Ⅰ.达朗贝尔法求“无限和半无限的”弦振动函数 (6)Ⅱ.分离变量法求两端固定弦振动方程 (7)三．各种情形下的弦振动求解及图像 (9)四．总结21一·实际现象的描述演奏者在演奏弦乐器（如二胡、提琴）时，用弓在弦上来回拉动，并通过另一只手指在按不同弦的不同地位的协调作用，奏出各种不同的美妙的音乐。
演奏者所用的乐器不同，奏出音乐的悦耳度也就不同。
演奏者虽然用弓所接触的只是弦的很小一段，似乎应该只引起这个小段的振动，而事实上，振动总是传播到整根弦。
这振动是怎样传播的呢？如何利用数学方法来求解这种物理问题？如何通过直观的方程来说明不同乐器演奏出的音乐效果不同的原因？可否利用matlab来将这种振动直观表示出来？通过对于弦振动方程的学习，及对matlab的初步了解，我对于不同定解问题下弦振动方程的求解做了初级小结。
也尝试利用matlab 直观表述不同定解条件下的弦振动动态图像。
二·问题的求解（一）求弦振动泛定方程在求解时，我们不妨认为弦是柔软的，就是说在放松的条件下，把弦完成任意的形状，它都保持静止。
由于弦乐器所用的弦往往是很轻的，它的重量只有力的几万分之一。
跟拉力相比，弦的重量完全可以略去，这样，真实的弦就抽象为“没有重量”的弦。
把没有重量的弦绷紧，它在不振动时是一根直线，就取这直线作为x轴。
把弦上各点的横向位移记作u。
这样，横向位移u是x和t的函数，记作u（x，t）。
要求解弦振动，首先应找出u所遵从的方程。
把弦细分为许多极小的小段，拿区间(x,x+dx)上的小段B为代表加以研究。
B既然没有重量而且是柔软的，它就只受到邻段A和C的拉力和。
弦的每小段都没有纵向（即x方向）的运动，所以作用于B的纵向合力应为零。
弦的横向加速度记作。
按照，小段B的纵向和横向运动分别为式中时弦的线密度，即单位长度的质量。
ds为小段B的弧长。
因考虑的振动为小围振动，这时、为小量，如果忽略、以上的高阶小量，则，，，，,又，。
这样，（1）和（2）简化为因此，弦中力不随x而变，它在整根弦中取同一数值。
另一方面，在振动过程中的每个时刻都有长度ds=dx，即长度ds不随时间而变，所以作用于B段的力也不随时间而变。
弦中力即跟x无关，又跟t无关，只能是常数，记为T。
则（4）式为由于dx取得很小，。
这样，B段的运动方程就成为（5）由于B是作为代表任选的，所以方程（5）适用于弦上各处，是弦做微小横振动的运动方程，简称弦振动方程。
由此就求得了自由振动状态下的弦振动方程为若为受迫振动，则方程为（二）解弦振动方程Ⅰ .达朗贝尔法求“无限和半无限的”弦振动函数弦振动方程为：即 (6)先求其通解：依据方程（6）的形式作代换，即在此代换之下，由此，方程（6）可化为（7）先对η积分，得（8）其中f是任意函数，再对ξ积分，就得到通解u dξ+f= （9）通解中的与可用定解条件确定。
因我们在此求解的为“无限长或半无限长的”弦，因而此种情况下就不存在边界条件，设初始条件是（10）将一般解（9）带入初始条件，得即由此解得以此带回（9）式即得满足初始条件（10）的特解dξ.即所谓的达朗贝尔公式。
Ⅱ.分离变量法求两端固定弦振动方程定解问题为：泛定方程（11）边界条件 (12)初始条件 (0) (13) 解：令 u(x,t)=X(x)T(t)带入泛定方程及边界条件得X (14)(15)因T(t)不恒等于零，故而（15）式即为 X(0)=0,X(l)=0 (16) 用遍除（14）式各项即得因此式左边是时间t的函数，与坐标x无关；右边是坐标x的函数，跟时间t无关。
若两边相等，则两边比为一常数。
将此常数记为-，即由此化为(19)先求解X，将（1）当λ时，方程（17）的解为X(x)=积分常数由条件（18）确定，即由此解出，进而有X，所以此解为无意义之解，故排除λ的可能。
（2）当λ时，方程（17）的解为X(x)=积分常数由条件（18）确定，即由此解出，进而有X，所以此解为无意义之解，故也排除λ的可能。
（3）当λ时，方程（17）的解为X(x)=积分常数由条件（18）确定，即由于，为使函数有意义，只能，进而解得由此可得 X(x)=（21）其中为任意常数将式（20）代入方程（19）得这个方程的解为其中A,B为任意常数将（21）和（22）代入u(x,t)=X(x)T(t) 即得分离变数的解进而可得弦振动方程的解为其中系数情况一：初速度不为零，初位移为零设初速度为ψ其解析解为其中系数取a=2，l=1，则弦振动动画的源程序为：function u(x,t)N=50;t=0:0.005:2.0; x=0:0.001:1;ww=u1fun1(N,0);h= plot(x,ww,'linewidth',3);axis([ 0, 1, -1, 1])sy=[ ];for n=2:length(t)ww=u1fun1(N,t(n));set(h,'ydata',ww);drawnow;sy=[sy,sum(ww)];endfunction wtx=u1fun1(N,t)x=0:0.001:1; a=2; wtx=0;for k=1:NBk=5/(k*k*pi*pi)*(cos(2*k*pi/5)-cos(4*k*pi/5)); wtx=wtx+Bk*sin(2*k*pi*t)*sin(k*pi*x);end此种情况与可用差分法求解clearN=4010; dx=0.0024;dt=0.0005; c=dt*dt/dx/dx;u(1:420,1)=0; x=linspace(0,1,420)';u(2:280,1)=0.05/280*(2:280)';u(281:419,1)=0.05/(419-280)*(419-(281:419)');u(2:419,2)=u(2:419,1)+c/2*(u(3:420,1)-2*u(2:419,1)+u(1:418,1));h=plot(x,u(:,1),'linewidth',3)axis([0,1,-0.05,0.05]);set(h,'EraseMode','xor','MarkerSize',18)for k=2:Nset(h,'XData',x,'YData',u(:,2));drawnow;u(2:419,3)=2*u(2:419,2)-u(2:419,1)+c*(u(3:420,2)-2*u(2:419,2)+u(1:418,2)); u(2:419,1)=u(2:419,2);u(2:419,2)=u(2:419,3);end此种情况下n=1，2，3,4,5,6时的本征振动随时间分布图cleara=1;l=1;x=0:0.05:1;t=0:0.001:3;u=0;n=1;[X,T]=meshgrid(x,t)R=5/(n*n*pi*pi)*(cos(2*n*pi/5)-cos(4*n*pi/5)) u=u+(R*sin(n*pi*a*T/l)).*sin(n*pi*X/l);figure(1)axis([0 1 -0.05 0.05])mesh(X,T,u)title('n=1本征振动随时间分布图')xlabel('x')xlabel('t')xlabel('u')情况二：初速度为零，初位移不为零设初位移为取a=1，l=1得其中系数此种情况下弦振动动画的源程序为：function FN=50t=0:0.005:2.0;x=0:0.001:1;ww=wfun(N,0);ymax=max(abs(ww));h= plot(x,ww,’linewidth’,3);axis([ 0, 1, -ymax, ymax])sy=[ ];for n=2:length(t)ww=wfun(N,t(n));set(h,’ydata’,ww);drawnow;sy=[sy,sum(ww)];endfunction wtx=wfun(N,t)x=0:0.001:1; a=1; wtx=0;for I=1:Nif I~=5wtx=wtx+0.05*( (sin(pi*(5-I)*4/5)-sin(pi*(5-I)*2/5)).../(5-I)/pi-(sin(pi*(5+I)*4/5)-sin(pi*...(5+I)*2/5))/(5+I)/pi )*cos(I*pi*a*t).*sin(I*pi*x);elsewtx=wtx+0.05/5*cos(I*pi*a*t).*sin(I*pi*x);endend情况三：有驱动力的弦振动方程定解问题为因有，所以此类问题可以用傅里叶级数法求解，如果初始条件为零，则同样也可用冲量法求解傅里叶级数法求解过程为：因边界条件为下的本征函数是（n=0,1,…）。
所以可以把所求的解展开为傅里叶余弦级数将这个级数代入泛定方程得将傅里叶级数展开等式两边比较系数得将的傅里叶余弦级数代入初始条件，得由傅里叶级数基本函数族的正交性可得的常微分方程在以上初始条件下的解为由此可得，所求的解是初始条件为零时的冲量法求解过程为：应用冲量定理法，先求解因边界条件为下的本征函数是（n=0,1,…）。
所以可以把所求的解展开为傅里叶余弦级数将此代入泛定方程得由此分离出的常微分方程为这个常微分方程的解为这样解v的傅里叶级数是将上式代入初始条件得比较等式两边系数得至此求得：因故而可得当l=1; a=1，在驱动力作用下的弦振动图像为w=0.5; al=pi*a/l; A=2;x=0:0.05:1;m=moviein(301);for k=0:300u=A/al/(w^2 - al^2)*(w*sin(al*k*0.1)...-A*al*sin(w*k*0.1))*cos(pi*x./l);plot(x,u);axis([0 1 -0.5 0.5])m(:,k+1)=getframe;end总结通过这次对于弦振动解法的小结，让我更深刻的了解了实际问题中为什么不同的乐器会奏出不同的音乐，主要是因为求解出的振动函数展开为傅里叶级数所包含的项不同，外加弹奏着所使用的驱动力不同，也会使振动效果不同，自然悦耳度也就不同。
通过利用matlab对于不同情况下的弦振动的动态表述，让我对于matlab有了更深的认识，由于时间有限，外加对于matlab的掌握程度比较低，也让我的操作受了很大的限制。
在今后的学习中更因该加强这方面的操作，真正将利用软件和原始计算求解能够完美结合，以此来加深自己对于数理方程的更深刻了解。