数学物理方法大作业

合集下载

数学物理方法习题解答(完整版)

数学物理方法习题解答(完整版)

数学物理方法习题解答一、复变函数部分习题解答第一章习题解答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 ⎧++≠⎪=+⎨+⎪⎩。

数学物理方法习题及解答

数学物理方法习题及解答

2. 试解方程:()0,044>=+a a z44424400000,0,1,2,3,,,,i k iiz a a e z aek aez i i ππππωωωωω+=-=====--若令则1.计算:(1)iii i 524321-+-+ (2)y =(3)求复数2⎝⎭的实部u 和虚部v 、模r 与幅角θ(1) 原式=()()()123425310810529162525255i i i i i i +⋅+-⋅+-++=+=-+--(2) 332()102052(0,1,2,3,4)k i e k ππ+==原式(3)2223221cos sin cos sin ,3333212u v 1,2k ,k 0,1,2,23i i i e r ππππππθπ⎛⎫==+=+==-+ ⎪⎝⎭⎝⎭=-===+=±±原式所以:,3.试证下列函数在z 平面上解析,并分别求其导数.(1)()()y i y y ie y y y x e x x sin cos sin cos ++-3.()()()()()()()()cos sin ,cos sin ,cos sin cos ,sin sin cos ,cos sin sin sin ,cos sin cos ,,,x x x x x x x x u e x y y y v e y y x y ue x y y y e y x ue x y y y y y ve y y x y e y y x ve y y y x y yu v u v x y y x u v z f z u iv z u f z =-=+∂=-+∂∂=---∂∂=++∂∂=-+∂∂∂∂∂==-∂∂∂∂=+∂'=∂证明:所以:。

由于在平面上可微所以在平面上解析。

()()()cos sin cos cos sin sin .x x x x vi e x y y y e y i e y y x y e y x x∂+=-++++∂由下列条件求解析函数()iv u z f += (),1,22i i f xy y x u +-=+-=解:()()()()()()()222222222212,2,212,2,,,2112,22111,0,1,1,,221112.222u v x y v xy y x x y v u v y x y x x x x x c x y x f z x y xy i xy y x c f i i x y c c f z x y xy i xy x y ϕϕϕϕ∂∂==+∴=++∂∂∂∂∂''=+=-=-+∴=-=-+∂∂∂⎛⎫=-+++-+ ⎪⎝⎭=-+==+==⎛⎫=-++-++ ⎪⎝⎭而即所以由知带入上式,则则解析函数2. ()21,3,,.ii i i i i e ++试求()()(((()()()2(2)Ln 144(2)4ln32Ln32ln32ln1222Ln 21cos ln sin ,0,1,2,3cos(ln 3)sin(ln 3),0,1,2,i i k k i ii i k i i k i i k i k i k i ii ii eeeei k e e e e i k i eeeππππππππππππ⎛⎫⎛⎫+ ⎪⎪-+++⎝⎭⎝⎭-++-+-⎛⎫⎛⎫++-+ ⎪⎪⎝⎭⎝⎭+====+=±±====+=±±=== 解:()222,0,1,2,cos1sin1.k i i k e e e e i π⎛⎫ ⎪⎝⎭+=±±=⋅=+3. 计算 2,:122c dzc z z z =++⎰()2222220110,1,1,11,220,022z z z z i z i z c z z z c z z ++=++=+==-+=≤++≠=++解:时,而在内,故在内解析,故原式 1.计算221(1),21c z z dz c z z -+=-⎰: ()2221(2),21cz z dz c z z -+=-⎰:(1)212(21)=4 z i z z i ππ==-+解:原式 (2)2112(21)=2(41)6z z i z z i z i πππ=='=-+-=解:原式. 计算2sin()114,(1):1,(2):1,(3): 2.122c z dz c z c z c z z π+=-==-⎰其中1sin (1)sin 442.112c z z z z i i z z πππ=-⎡⎤-⎢⎥===⎢⎥+-⎢⎥⎣⎦⎰解:(1)原式1sin (1)sin 442.11c z z z z i i z z πππ=⎡⎤+⎢⎥===⎢⎥-+⎢⎥⎣⎦⎰(2)原式 12(3):2,1,11,.c z z z c c ===-以分别以为中心,为半径,做圆1222sinsin44.11c c z zdz dz i i i z z ππ=+=+=--⎰⎰原式 3、将下列函数按()1-z 的幂级数展开,并指明收敛范围。

数学物理方法大作业

数学物理方法大作业

基于分离变量法的波导中的电磁波研究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

数学物理方法大作业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一、单项选择题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 -=的实部,则这个解析函数为。

数学物理方法习题集

数学物理方法习题集

数学物理方法习题集第一章 复数与复变函数习题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 -≤+≤+。

数学物理方法试题汇总

数学物理方法试题汇总

12届真题1. 求下列各小题(2*5=10分):(1)用几何图形表示0arg(1)4z π<-<; (2)给出序列(1/)sin 6n n z i n π=+的聚点; (3)在复数域中求解方程cos 4z =的解;(4)给出二阶偏微分方程的基本类型;(5)给出解析函数所满足的柯西-黎曼方程。

2.按给定路径计算下列积分(5*2=10分):(1)320Re izdz +⎰,积分路径为线段[0,3]和[3,3+2i]组成的折线;(2)11,==⎰积分路径由z=1出发的。

3.利用留数定理计算下列积分(5*2=10分):(1)241x dx x +∞-∞+⎰; (2)3||1zz e dz z =⎰。

4.求常微分方程20w z w ''-=在0z =邻域内的两个级数解(15分)。

5.求下列线性非奇次偏微分方程的通解:2222u u xy y x y∂∂-=-∂∂(15分)。

6.利用分离变量法求解:(20分)2222000(),|0, |0,0, 0.x x l t t u u x l x t x u u u u t ====⎧∂∂-=-⎪∂∂⎪⎪==⎨⎪∂⎪==∂⎪⎩7.用拉普拉斯变换方法求解半无解问题(20分)220, 0,0,(0,)1, lim (,) 0, (,0)|0, 0.x u u x t t x u t u x t t u x x κ→∞⎧∂∂-=>>⎪∂∂⎪⎪=>⎨⎪=>⎪⎪⎩有界,2005级一、填空(请写在答题纸上,每题6分,共计48分)1. 三维泊松方程是______________________________2. 边界为Γ的区域Ω上函数u 的第二类边界条件为___________________。

3. 极坐标下的二维拉普拉斯方程为__________________________。

4. 定解问题2002||0tt xx t t t u u x u x u ===-∞<<+∞⎧⎪⎨==⎪⎩, ,的解__________________________。

数学物理方法习题解答

数学物理方法习题解答

第八章习题P201:1,2,5,6,11,12,13,16,17,201.长为l 的弦,两端固定,弦中张力为T ,在距一端为0x 的一点以力0F 把弦拉开,然后突然撇除这力,求解弦的振动。

解:此题的定解问题为200000000,(0),(0,)(,)0,,(0),(,0)(),(),0.tt xx t t u a u x l u t u l t F l x x x x T l u x F x l x x x l T l u =⎧-=<<⎪==⎪⎪-⎧⎪<<⎪⎪⎨=⎨⎪⎪⎪-<<⎪⎩⎪⎪=⎩)4()3()2()1(令(,)()()u x t X x T t =代入泛定方程(1)中得X T X aTλ''''==- 可得20T a T X X λλ''⎧+=⎨''+=⎩ (0)()0X X l ==求解关于x 本征值问题,得到本征值和本征函数()2/n l λπ= (1,2,3,n =⋅⋅⋅⋅⋅⋅()sinn X x C x lπ= 将本征值代入关于t 的常微分方程,得到22220a n T T lπ''+= 其解为 ()cossin n n n n a n aT t A x B t l lππ=+ 1(,)()()cos sin sin n n n n a n a n u x t X x T t A t B t x l l l πππ∞=⎛⎫∴==+ ⎪⎝⎭∑将u 的级数解代入初始条件(4)得到001|sin cos sin t t n n t n n a n a n a n a n u A x B t xl l l l l πππππ∞===⎛⎫=-+ ⎪⎝⎭∑1sin 0nn n a n B x l lππ∞===∑ 0n B ∴=则1(,)cossin n n n a n u x t A t x l lππ∞=∴=∑ 根据初始条件(3)有0001000,(0),(,0)sin (),(),n n F l x x x x n T lu x A x F x l l x x x l T l π∞=-⎧<<⎪⎪==⎨⎪-<<⎪⎩∑02()sin l n n A d l l πϕξξξ=⎰ 000000022sin ()sin x l x F l x F x n n d l d l T l l l T l l ππξξξξξξ-=+-⎰⎰ 02000022222sin cos cos x lx F l x F x l n l n n l n l T l n l n l l T l n l ππξππξξξπππ⎧⎡⎤-⎪=--⎨⎢⎥⎣⎦⎪⎩020022sin cos lx F x l n n n T l n l l l ππξπξξπ⎫⎪⎡⎤--⎬⎢⎥⎣⎦⎪⎭000000000220()2sin cos cos cos xF l x l n x n x n x F x n x n l T n l l l T n l πππππππ⎧-⎪⎡⎤⎡⎤=---⎨⎢⎥⎢⎥⎣⎦⎣⎦⎪⎩0000022cos sin cos F x l n x n x n x n n T n l l l ππππππ⎫⎡⎤---+⎬⎢⎥⎣⎦⎭ 002221sin F l n x T n lππ=∴ 00221121(,)cos sin sin cos sin n n n F l n x n a n n a n u x t A t x t x l l T n l l l ππππππ∞∞==∴==∑∑2.求解细杆热传导问题,杆长l ,两端保持为零度,初始温度分布20/)(l x l bx u t -==。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

基于分离变量法的波导中的电磁波研究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=σ;波导内的电磁场为时谐场,角频率为ω。

2.2矩形波导中的电磁波现在我们求解矩形波导中的电磁波解。

选一直角坐标系,如图1所示。

取波导内壁面为0=x 和a ,0=y 和b ;z 轴沿电磁波传播方向。

在一定频率下,管内电磁波是方程(5)的解。

次解在管壁上还满足边界条件0=⨯E n ]4[,即电磁场在管壁上的切向 分量为零。

由于电磁波沿z 轴方向传播,它应有传播因子t i ik z e ω-。

因此,我们把电场E 取为E (x,y,z)= E (x,y)z ik ze . (16)将(16)代入(5)式得()()()0,,222222=-+⎪⎪⎭⎫ ⎝⎛∂∂+∂∂y x E k k y x E y x z ]4[ (17) 用直角坐标分离变量]12[,设()y x u ,为电磁场的任一直角分量,它满足方程(17)。

设()y x u ,)()(y Y x X = (18) 从而方程(17)可以分解为两个方程0222=+X k dxX d x (19) 0222=+Y k dy Y d y (20) 2222k k k k z y x =++ (21) 求解方程(19)式和方程(20),得()y x u ,的特解()()y k D y k C x k D x k C y x u y y x x sin cos sin cos ),(2211++= (22)其中,1C ,1D ,2C ,2D 为任意实数,当()y x u ,具体表示E 某特定分量时,考虑边界条件0=⨯E n 及0=∂∂nE n,可以得到对这些常数的一些限制条件。

2.2.1矩形波导TM 波矩形波导的横截面如图1所示,波导内传播TM 波时,有0=z H 。

波导内的电磁场由z E 确定。

在给定波导管中,z E 满足下面的波动方程和边界条件]1[:02222222=+∂∂+∂∂+∂∂z z z z E k z E y E x E (23) 0,00====ax zx zE E (24)0,00====by zy zE E (25)由均匀波导中,设z z z e y x E z y x E γ-=),(),,( (26) 将(26)式代入方程(23)中,得0),(22222=⎥⎦⎤⎢⎣⎡+∂∂+∂∂y x E k y x z c (27) 其中222k k c +=γ为截止波数。

用分离变量法]12[求解。

其解为:()XY y x E z =, (28) 将式(28)代入方程(27)中,然后方程两边除以XY ,于是方程可分裂为两个常微分方程0222=+X k dxX d x (29) 0222=+Y k dyY d y (30) 且222c y x k k k =+方程(29)的通解为()x k D x k C X x x sin cos 11+=,并考虑边界条件(24)得am k x π=其中⋅⋅⋅=3,2,1m 故x a m C X πsin 1= (31) 同理,得方程(30)的通解为:()y k D y k C Y y y sin cos 22+=考虑边界条件式(25)得bn k y π=其中⋅⋅⋅=3,2,1n 故)sin(2y b n C Y π= 所以,得到矩形波导中TM 波的纵向场分量 z m z z z e y bn x a m E e y x E z y x E γγππ--==)sin()sin(),(),,( (32) 式中,21C C E m =由激励场源强度决定。

由式222c y x k k k =+得截止波数222222⎪⎭⎫⎝⎛+⎪⎭⎫ ⎝⎛=+=+=b n a m k k k k yx c ππγ (33)利用式(12*)~(15*)可求得TM 波的其它横向场分量]1[z m cx e y b n x a m E a m k z y x E γπππγ--=)sin()cos()(),,(2(34) z m cy e y bn x a m E b n k z y x E γπππγ--=)cos()sin()(),,(2(35) z m cx e y b n x a m E b n k i z y x H γπππωε-=)cos()sin()(),,(2(36) z m cy e y b n x a m E a m k i z y x H γπππωε--=)sin()cos()(),,(2(37) 2.2.2 矩形波导TE 波对于TE 波,因为TE 波在传播方向上0=z E 。

故由式(12+)~(15+)可得TE 波的纵向场分量与横向场分量关系]1[为:xH k H zcx ∂∂-=2γ (38)yH k H zcy ∂∂-=2γ (39)y H k i E zcx ∂∂-=2ωμ (40) xH k i E zc y ∂∂=2ωμ (41) 故波导内的电磁场由z H 分量确定,在给定的矩形波导中,z H 满足下面的波动方程和边界条件:022=+∇z z H k H (42)0,00=∂∂=∂∂==ax z x zxH xH (43)0,00=∂∂=∂∂==by z y zyH yH (44)用分离变量法求解]12[,可以得到TE 波的纵向场分量)cos()cos(),(y bn x a m H y x H m z ππ= ⋅⋅⋅=3,2,1,0,n m (45) 式中m H 由激励场源决定。

将式(45)代入式(38)~(41)中,得:z m cx e y b n x a m H b n k i z y x E γπππωμ-=)sin()cos()(),,(2(46) z m cy e y b n x a m H a m k i z y x E γπππωμ--=)cos()sin()(),,(2(47) z m c x e y bn x a m H am k z y x H γπππγ-=)cos()sin()(),,(2(48) z m cy e y bn x a m H b n k z y x H γπππγ-=)sin()cos()(),,(2(49) 下面以TE 10模为例研究其分布其场分量为cos z ik zz x H e a π=sin z ik zz x ik a x H e aππ-=sinz ik z y i axE e aωμππ=由matlab 作图程序如下,取a=5cm ,b=3cm ,令z 为1cm w=10^11;u=4*pi*10^(-7); a=.05; b=.03;e=10^7/(4*pi*9*10^16); kz=sqrt(w^2*u*e-(pi/a)^2); x=(0:0.001:0.05); y=(0:0.001:0.03); [X,Y]=meshgrid(x,y);hz=cos(pi*X/a).*cos(kz*0.01);hx=sin(kz)*0.01*kz*(a/pi)*sin(pi*X/a); ey=-(w*a*u/pi)*sin(pi*X/a)*sin(kz*0.01); subplot(1,3,1) mesh(X,Y ,hz) subplot(1,3,2) mesh(X,Y ,hx) subplot(1,3,3) mesh(X,Y ,ey)下图波形依次为Hz 波,Hx 波,Ey 波基于Microsoft Visual Studio 2010编程实现施密特正交化设x 为向量空间m R 的子空间,且{21,,x x }为其基底,因此n m n X ≥=,dim 。

相关文档
最新文档