第一篇:数理方程-分离变量法
第八章
分离变量法
22u2ua0xl,t022tx u(0,t)0,u(l,t)0t0u(x,0)u(x,0)(x),(x)0xlt对于这样的定解问题,我们将介绍分离变量法求解,首先回忆高数中我们如何处理的求解的,高数中处理微分或重积分是把函数分成单元函数
分离变量法的思路:对于二阶线性微分方程变换成单元函数来求解,也就是通过分离变量法把x、t两个变量分开来,即把常微分方程变化为两个偏微分方程来求解。
分离变量法的思想:先求出具有分离形式且满足边界条件的特解,然后由叠加原理做出这些解的线性组合,最后由其余的定解条件确定叠加系数(叠加后这些特解满足边界条件不满足初始条件,再由初始条件确定通解中的未知的数)。
叠加原理:线性偏微分方程的解的线性组合仍是这个方程的解。特点:(1)数学上 解的唯一性来做作保证。(2)物理上 由叠加原理作保证。例:有界弦的自由振动
1.求两端固定的弦的自由振动的规律
22u2ua0xl,t022tx u(0,t)0,u(l,t)0t0u(x,0)u(x,0)(x),(x)0xlt第一步:分离变量(建立常微分方程定解问题)令u(x,t)X(x)T(t)
这个思想可从实际的物理现象可抽象出来,比如我现在说话的声音,它的振幅肯定随时间变化,但到达每个同学的位置不同,振幅又是随位置变化,可把声音分成两部分,一部分认为它随时间变化,一部分随位置变化。
第二步:代入方程
(偏微分就可写成微分的形式,对于u有两个变量,但对于X、T都只有一个变量)
X(x)T(t)a2X(x)T(t)
变形得X(x)T(t)= X(x)a2T(t)左边与t无关,右边与x无关,左右两边相互独立,要想相等,必定等于一个常数。由于x, t 是相互
独立的变量,上式必然等于同一常数。
方程左边为关于x的函数,方程右边为关于t的函数,只有当左右两边都等于常数的时候才成立 令其为(得到的两个常微分方程形式比较标准)
X(x)X(x)0
T(t)a2T(t)0
得到两个常微分方程 第三步:代入边界条件
得到:X(0)T(t)0
X(l)T(t)0,由于是t>0得值,T(t)是一个范围内不固定的值,所以X(0)0
X(l)0
常微分方程含,未知,需要对进行讨论
X(x)X(x)0,X(0)0
X(l)0
特征(固有)值问题:含有待定常数常微分方程子一定条件下的求解问题。特征(固有)函数:和特征值相对应的非零解 第四步:确定特征值并得到它的特征函数 分情况讨论:
1)<0时, 特征方程为R0,特征根为:R- 得通解为X(x)Aex2Bex(A、B为待定系数)
x把定解条件X(0)0
X(l)0代入通解X(x)Ae得到A+B=0
Bex
AelBel0
x于是A=B=0X(x)AeBex即X(x)=0 则u(x,t)X(x)T(t)=0,零解无意义 即<0时,定解问题无解。2)=0时, X(x)X(x)0 有X(x)AxB A=B=0X(x)AexBex即X(x)=0 则u(x,t)X(x)T(t)=0,零解无意义 3)>0时, X(x)X(x)0
令2(为非零实数)
特征方程为R0,特征根为虚数:R-i 通解为X(x)AcosxBsinx(A、B为待定系数)
把定解条件X(0)0,X(l)0代入通解X(x)AcosxBsinx 2X(0)0得到A =0,即X(x)Bsinx X(a)0得到Bsinl0
在B≠0的情况下,有sinl=0,即n为非零实数)
现在就完成了用分离变量法求解X(x)的部分,得到特征值为nn(数为:X(x)Bnsin2n(n=1,2,3,…注意n≠0,若n=0,则=0,0而ln2),所对应的特征函lnx ln2)代入 l下面求解关于t的常微分方程
T(t)T(t)0,将n(2n22Tn(t)aTn(t)0,这种情况的通解与X(x)X(x)0的>0的情况相同。
l2cos即Tn(t)CnnatnatsinDn
(n=1,2,3,…)
ll至此Xn(x)与Tn(t)都求出来了,所以定解问题的n个特解(这n个特解均满足边界条件)为:
un(x,t)Xn(x)Tn(t)=(CncosnatnatnDnsin)sinx
(n=1,2,3,…)lll根据叠加原理,特解的叠加仍是方程的解,所以得到通解
u(x,t)un(x,t))
i1n=(Cncosi1nnatnatnDnsin)sinx(n=1,2,3,…)lllu(x,0)(x)求解)t其中Cn、Dn为待定系数(利用初始条件u(x,0)(x),第五步: 利用本征函数的正交归一性确定待定系数
u(x,t)(Cncosi1nnatnatnDnsin)sinx lll
u(x,t)t0u(x,0)
Cnsini1nnx(x)l(i1nu(x,t)ti1nt0natnatnanatnCnsinccos)sinxt0 lllllnanDnsinx(x)ll(x)与(x)正是傅里叶正弦级数,Cn、Dn是傅里叶系数。
利用三角函数的正交性
l1cos(2n/l)nlxdxdx 00l22lnm1lnmnmsinxxdx[coscosxx]dx0(m≠n)0ll20lllsin2lmnml得到:(x)sinxdxCnsinxxdxCn
00lll2n0l2lm(x)sinxdx 0lll2lm2lm(x)sinxdx(x)sinxdx 同理,Dnnal0lna0l于是得到:Cn回顾整个求解过程,可作出分离变量法流程图
u|t0(x)u|t0(x)uta2uxxu|x0u|xl0X(0)X(L)0分离变量流程图uT(t)X(x)T'/(a2T)X“/XT'a2T0TAexp(a2t)X”X0Xsinx,nlunTn(t)Xn(x)uu(x,t)2.解的性质
uTnXn
un(x,t)=(Cncos于x的函数)natnatnDnsin)sinx---------方程的特解(前面是关于t的函数,后面是关lllun(x,t)=(CncosnnatnatnDnsin)sinx=Ancos(ntn)sinx llllnaD,narctann lCn22其中:AnCnDn,n当xx0时,un(x,t)=Ansin点的振动方程)。
nx0cos(ntn)---------弦上确定的一点以频率n做振动(弦上某lnx----------某一时刻,特解为正弦函数的形式,所有点l当tt0时,un(x,t)=Ancos(nt0n)sin的位置,波动方程(驻波的方程),每个特解代表一个驻波,因此分离变量法又称为驻波法。
标准的驻波方程:y2Acos2xcost
sinn2x的(驻波)波长为nl(n=1,2,3,…)
nl
频率:fnnna 22lna2T la2ln波速:vnfnn3.分离变量法概要:
(1)作分离变量假设,代入方程和边界条件中得到固有值问题(2)确定固有函数和固有值(3)写出定解问题的特解(4)将特解叠加无,给出通解
(5)用初始条件确定通解系数(傅立叶展开)4.回顾整体思路:
2u(x,0)2u2u(x)
定解问题2a初始条件u(x,0)(x), 边界条件u(0,t)0,u(l,t)0 2ttx22u2ua将假设u(x,t)X(x)T(t)代入方程,此偏微分方程得到两个常微分方程t2x2X(x)X(x)0
T(t)a2T(t)0。
将边界条件u(0,t)0,u(l,t)0代入u(x,t)X(x)T(t),得到X(0)0、X(l)0,求解已知定解条件的常微分方程X(x)X(x)0的特征值为nn(n2n),特征方程Xn(x)Bnsinx,llnatnatcossinDn求解T(t)a2T(t)0的特征函数Tn(t)Cn,所以
llnatnatnDnsin)Bnsinx。un(x,t)Xn(x)Tn(t)=(Cncoslll2根据叠加原理,特解的叠加是方程的通解,所以得到:
u(x,t)un(x,t))i1n=
(Cncosi1nnatnatnDnsin)sinxlll,将初始条件u(x,0)(x),u(x,0)(x)代入,求解待定系数Cn、Dn(傅立叶展开)。tx(10x),求弦做
1000分离变量法的适用条件:任何二阶线性(齐次)偏微分方程
例一:设有一根长为10个单位的弦,两端固定,初速度为零,初位移为(x)微小横振动时的位移。
22u4u100x10,t022tx u(0,t)0,u(10,t)0t0x(10x)u(x,0)u(x,0),01000t解:设u(x,t)X(x)T(t),代入
X1T4 X10T4得到:X(x)X(x)0
T(t)10T(t)0
u(0,t)X(0)T(t)0,u(10,t)X(10)T(t)0
X(x)X(x)0,0x10得到本征值问题:,
X(0)0,X(10)0经讨论20时,有非零解,X(x)AcosxBsinx
X(0)A0,X(10)Bsin100,n2n,n=1,2,3,… 10nn22x 得到特征值:
得到特征方程:Xn(x)Bnsin10100于是:T(t)100n2cos10ntDnsin10nt 2T(t)0,其解为Tn(t)Cnun(x,t)Xn(x)Tn(t)
ncos10ntDnsin10nt)x(Cn10nx =(Cncos10ntDnsin10nt)sin10Bnsinu(x,t)un(x,t))=(Cncos10ntDnsin10nt)sini1n1nnx 10将初始条件u(x,0) Cnsin10nt n1x(10x)
1000210x(10x)nsinxdx运用分部积分法求解 010100010110nx(10x)sinxdx
=5000010n为偶数20=44(1cosn)4
n为奇数5n5n44Cnu(x,0)nanDnsinx0,故Dn=0.tlln1所以uun(x,t))=i1n1n(2n1)4sinx cos10(2n1)t105(2n1)4422u2ua0xl,t022txu(l,t)例二: u(0,t)0,0t0xu(x,0)x22lx,u(x,0)00xlt解:设u(x,t)X(x)T(t),代入
X1T2 XaT得到:X(x)X(x)0
T(t)a2T(t)0
u(0,t)X(0)T(t)0u(l,t)0x
X(0)0
u(l,t)X(l)T(t)0xX(l)0
X(x)X(x)0,0xl得到本征值问题:,
X(0)0,X(l)0经讨论0,X(x)AexBex(A、B为待定系数)
x把定解条件X(0)0
X(l)0代入通解X(x)Ae得到A+B=0
Bex
AelBel0
于是A=B=0即X(x)=0 =0时, X(x)X(x)0,有X(x)AxB,A=B=0即X(x)=0 20时,X(x)X(x)0,X(x)AcosxBsinx
X(0)0X(l)0所以n A0
X(x)Bcosl0
(2n1) n=1,2,3,… 2l2(2n1)(2n1)22X(x)Bsinx 写出特征值和特征函数,nn22l4l(2n1)22Tn(t)0 T(t)aT(t)0变为Tn(t)a4l222(2n1)a(2n1)asintDnt,2l2l(2n1)a(2n1)a(2n1)cossintDnt)sinx 所以un(x,t)Xn(x)Tn(t)=(Cn2l2l2lcosTn(t)Cn所以uun(x,t))=(Cncosi1i12nn(2n1)a(2n1)a(2n1)atDnsint)sinx 2l2l2l由初始条件u(x,0)x2lx,u(x,0)0确定Cn、Dn。tu(x,0)Cnsinl(2n1)xx22lx 2l2(2n1)32l22 Cn(x2lx)sinxdx33l02l(2n1)u(x,0)(2n1)a(2n1)Dnsinx0,Dn=0 t2l2luun(x,t))=i1n32l231(2n1)a(2n1)costsinx 32l2li1(2n1)n附录1:二阶常系数微分方程:ypyqy0 特征方程:rprq0 根的三种情况
2r1r2r1r2rri
得到常系数微分方程的通解: yC1er1xC2er2x
附录2:线性方程满足叠加原理。
线性齐次方程(只含未知量的一次项,无零次项)通解为所有线性无关特解的叠加;而线性非齐次方程通解为其特解与相应齐次方程(去掉零次项后的线性方程)通解的叠加。
rxrxyC1eC2xeyex(CcosxCsinx)12附录3:和差化积公式
cos(A-B)=cosAcosB+sinAsinB
第二篇:分离变量法习题
第十章习题解答 求解混合问题
utta2uxx0(0xl,t0)0
u(0,t)0,u(l,t)0,其中(x)v00u(x,0)0,u(x,0)(x)t0xccxc cxl解:用分离变量法:设混合问题的非零解函数为u(x,t)X(x)T(t),则,utt(x,t)X(x)T(t),uxx(x,t)X(x)T(t)
代入混合问题中的微分方程可得:
X(x)T(t)aX(x)T(t)02X(x)X(x)aT(t)T(t)2
由初始条件可得:u(0,t)X(0)T(t)u(l,t)X(l)T(t)0X(0)X(l)0由此可得,X(x)为如下常微分方程边值问题的非零解:
X(x)X(x)0X(0)0,X(l)0(0xl)
若λ<0,则此定解问题的微分方程的通解为 X(x)c1exp(x)c2exp(x),代入边值条件后可得c1c20X(x)0,不符合要求。若λ=0,则此定解问题的微分方程的通解为
X(x)c1c2x,代入边值条件后仍可得c1c20X(x)0,不符合要求。若λ>0,则此定解问题的微分方程的通解为 X(x)c1cos代入边界条件后可得: X(0)c1cos0c2sin0c10X(x)c2sinx,2xc2sinx,X(l)c2sinl0,X(x)0sinnxlnl0,n,l所以可取 X(x)Xn(x)sin
(n1,2,)由T(t)所满足的方程可得:
T(t)a22T(t)0T(t)Tn(t)ancosnatlnatlbnsinnatl,所以,原混合问题的微分方程的满足边界条件的分离变量形式解为 u(x,t)un(x,t)Xn(x)Tn(t)(ancosbnsinnatl)sinnxl,设原混合问题的解函数为 u(x,t)n1(ancosnatlbnsinnatl)sinnxl,则由初始条件可得:0u(x,0)n1ansinnxlan0(n1,2,)
ut(x,t)n1nalbncosnatlsinnxlnxl, (x)ut(x,0)n1natlbnsinbnna2l0(x)sinnxldx,bnna2ccv0sinnxldx2v0lna22(cosn(c)lnxlcosn(c)l)(*)所以,原混合问题的解为 u(x,t)2 求解混合问题
bn1nsinnatlsin,其中的bn由(*)给出。
utta2uxx0(0xl,t0)
u(0,t)E,u(l,t)0
u(x,0)0,u(x,0)0(E为常数)t解:由于边界条件非齐次,需作函数变换如下:设
v(x,t)u(x,t)El(lx)u(x,t)v(x,t)El(lx),则
vxx(x,t)uxx(x,t),vt(x,t)ut(x,t),vtt(x,t)utt(x,t),2vtt(x,t)avxx(x,t)utt(x,t)auxx(x,t)0,v(0,t)u(0,t)
v(x,0)u(x,0)ElEl(l0)u(0,t)E0,v(l,t)u(l,t)00,(lx)El(lx),vt(x,0)ut(x,0)0,所以,u(x,t)是原混合问题的解的充要条件是:v(x,t)是如下混合问题的解:
2vtt(x,t)avxx(x,t)0(0xl,
v(0,t)0,v(l,t)0Ev(x,0)(lx),vt(x,t)0lt0)
(*)
用分离变量法求解此定解问题,由分离变量法的标准步骤可得:
v(x,t)n1(AncosnatlBnsinnatl)sinnxl,代入初始条件可得:,Bn0,An2llEl0(lx)sinnxldx2En(n1,2,)
所以,v(x,t)n12EncosnatlElsinnxl,原混合问题的解函数为u(x,t)3 求解下列阻尼波动问题的解:
(lx)n12Encosnatlsinnxl
utt2huta2uxx0(0xl,t0)
u(0,t)0,ux(l,t)0
u(x,0)(x),u(x,0)(x)t其中,h为正常数,且ha2l。
解:使用分离变量法,设原定解问题的微分方程有如下分离变量形式非零解函数满足边界条件:
u(x,t)X(x)T(t)
则容易算得:uxx(x,t)X(x)T(t),ut(x,t)X(x)T(t),utt(x,t)X(x)T(t),代入方程后化简可得:
T(t)2hT(t)aT(t)2X(x)X(x)
0u(0,t)X(0)T(t)X(0)0,0ux(l,t)X(l)T(t)X(l)0,T(t)2hT(t)aT(t)0
X(x)X(x)0
,X(0)0,X(l)02由X(x)的非零性可得0,此时,X(x)c1cosxc2sinx,X(0)c1cos0c2sin0c10X(x)c2sinx,取c21得:X(x)sin2n1l0n
2l22x,X(l)cos2n1将代入T(t)所满足的方程可得:T(t)2hT(t)aT(t)0
l
22n12ha0nh2l2(2n1)ah
2l222
ha2l(2n1)a2lnh(2n1)a2hi2l(n1,2,)
从而有:
T(t)Tn(t)eht(AncosntBnsinnt),2n1a2l22其中
nh(n1,2,),(1)
设原混合问题的解函数为:
u(x,t)n1eht(AncosntBnsinnt)sin(2n1)2lx,
(x)u(x,0)ln1Ansinl(2n1)2lx,(2n1)xl(1cosdx,0022l2l22l(2n1)xdx(n1,2,)
(2)所以
An(x)sin0l2l而
sin2(2n1)xdx1ut(x,t)n1eht((hAnnBn)cosnt(hBnnAn)sinnt))sin(2n1)x2l
(x)ut(x,0)1n1(hAnnBn)sin(2n1)x2l,Bnn(hAn2ll0(x)sin(2n1)x2ldx)。
(3)
所以,原混合问题的解是u(x,t)n1eht(AncosntBnsinnt)sin(2n1)2lx,其中的 n,An,Bn分别由(1)式、(2)式、(3)式给出。
4 求解混合问题
uxxLCutt(LGRC)utGRu
u(0,t)0,ux(l,t)0GEu(x,0)E,u(x,0)tC(0xl,t0)
其中L、C、G、R为常数,且LG=RC。(提示:作函数变换u(x,t)exp(Rt/L)v(x,t))
解:记a21LC,bGCRL,混合问题的微分方程两边同除LC,方程可化为
a2uxx(x,t)utt(x,t)2but(x,t)b2u(x,t),a22x(u(x,t)exp(bt))t22(u(x,t)exp(bt)),设v(x,t)u(x,t)exp(bt),则有
a2vxx(x,t)vtt(x,t),而且,vx(x,t)ux(x,t)exp(bt),()0,所以
v(0,t)u(0,t)expbtvt(x,t)ut(x,t)exp(bt)bu(x,t)exp(bt),vx(l,t)ux(l,t)expbt()0,vt(x,0)ut(x,0)bu(x,0)0,(0)u(x,0)E, v(x,0)u(x,0)expb所以,若u(x,t)是原混合问题的解函数,则v(x,t)是如下混合问题的解函数:
vtt(x,t)a2vxx(x,t)0
v(0,t)0,vx(x,t)0v(x,0)E,v(x,t)0t(0xl,t0)
用分离变量法求解此混合问题,设方程的分离变量解形式的满足边界条件的非零解为 v(x,t)X(x)T(t),则
vx(x,t)X(x)T(t),vxx(x,t)X(x)T(t),vxx(x,t)X(x)T(t), X(x)X(x)T(t)aT(t)2
由齐次边界条件可得,X(x)为如下定解问题的解:
X(x)X(x)0X(x)c1cosxc2sinx,X(0)0,X(l)0
X(0)0c10,取c21得X(x)sinx,X(l)T(t)aT(t)2(2n1)cosl0n2lnT(t)Tn(t)Ancos(2n1)x2l2(n1,2,),(2n1)at2l
(2n1)at2lBnsin,X(x)Xn(x)sin(n1,2,),设
v(x,t)n1(Ancos(2n1)at2llBnsin(2n1)at2l)sin(2n1)x2l
代入初始条件可得:An2l0v(x,0)sin(2n1)x2ldx4E(2n1),Bn0,所以
v(x,t)(2n1)n14Ecos(2n1)at2lsin(2n1)x2l
所以,原题目所给的混合问题的解函数为:
u(x,t)exp(bt)n14E(2n1)cos(2n1)at2lsin(2n1)x2l。用固有函数法求解
utta2uxxg(const),
u(0,t)0,ux(l,t)0u(x,0)0,u(x,0)0t(0xl,t0)
解:用分离变量法:设原混合问题的微分方程对应的齐次方程有如下分离变量形式的非零解函数:u(x,t)X(x)T(t),利用分离变量法的标准步骤可求得: (2n1)
n,2l2X(x)Xn(x)sin(2n1)x2l(n1,2,)
将f(x,t)g展开成Xn(x)的广义Fourier级数如下:
fn(t)2ll0f(x,t)Xn(x)dx2ll0gsin(2n1)x2ldx4g(2n1),T(t)a2nT(t)fn(t)16gl(2n1)atT(t)T(t)(1cos)n3322l(2n1)aT(0)0,T(0)02[注:方程T(t)aT(t)fn(t)的通解为
Tn(t)Ancos
(2n1)at2lBnsin(2n1)at2l16gl(2n1)a332,代入初始条件即可得此处的结果。] 所以,题目所给的混合问题的解函数为
u(x,t)Tn(t)Xn(x)n1(2n1)16gl3a32(1cos(2n1)at2lt0))sin(2n1)x2l。
ut(x,t)a2uxx(x,t)06.求解混合问题u(0,t)0,ux(l,t)0u(x,0)u(const)0(0xl,。
解:用分离变量法:设混合问题中的微分方程有如下满足边界条件的分离变量形式的非零解函数:u(x,t)X(x)T(t),则
ut(x,t)X(x)T(t),ux(x,t)X(x)T(t),uxx(x,t)X(x)T(t),代入方程后化简再由边界条件可得:
T(t)aT(t)2X(x)X(x)T(t)aT(t)0,22X(x)aX(x)0
u(0,t)X(0)T(t)0X(0)0,ux(l,t)X(l)T(t)0X(l)0,所以,X(x)为如下常微分方程边值问题的非零解函数:
X(x)X(x)0X(0)0,X(l)0
2(0xl)
(2n1)解之得 n,2lX(x)Xn(x)sin(2n1)x2l(n1,2,),2(2n1)a
T(t)na2T(t)0T(t)Tn(t)Anexp(t)。
2l设原问题的解函数为
u(x,t)n1(2n1)x(2n1)a,Anexp(t)sin2l2l2由初始条件可得:
u0u(x,0)An1nsin(2n1)x2l4u0,由此可得:
An2ll0u0sin(2n1)x2ldx(2n1)2(n1,2,),所以,u(x,t)n1(2n1)x(2n1)a exp(t)sin(2n1)2l2l4u0 7 ut(x,t)a2uxx(x,t)0(0xl,7.求解混合问题u(0,t)0,ux(l,t)u(l,t)0u(x,0)(x)t0)
解:用分离变量法:设混合问题中的微分方程有如下满足边界条件的分离变量形式的非零解函数:u(x,t)X(x)T(t),则
ux(x,t)X(x)T(t),uxx(x,t)X(x)T(t),ut(x,t)X(x)T(t),代入方程后化简,并由边界条件可得:
T(t)a2T(t)0,X(x)X(x)0,u(0,t)X(0)T(t)0X(0)0,ux(l,t)u(l,t)(X(l)X(l))T(t)0X(l)X(l)0,所以,X(x)为如下常微分方程边值问题的解函数:
X(x)X(x)0(0xl)
X(0)0,X(l)X(l)0由u(x,t)是非零解可得:0X(x)c1cos
X(0)0c10X(x)sinxxc2sinx
(letc21),X(l)X(l)设
tanlcoslsinl0tanl(n1,2,),则nn
2
n0所以,X(x)Xn(x)sinnx,22((an)t)
T(t)(an)T(t)0T(t)Tn(t)Anexp(n1,2,),设原混合问题的解函数为
u(x,t)An1nexp((an)t)sinnx,2利用Xn(x)的正交性可求得 An(x)sin0lnxdx(n1,2,)。
[注]:可以证明:Xn(x)具有正交性。
l0sinnxdx2 8 ut(x,t)a2uxx(x,t)08.求解混合问题u(0,t),u(l,t)u(x,0)u0(0xl,t0),其中,,,u0为常数。
解:作函数变换 v(x,t)u(x,t)(则
ut(x,t)vt(x,t),lx)u(x,t)v(x,t)(l x),uxx(x,t)vxx(x,t),u(0,t),u(l,t)v(0,t)0,v(l,t)0,u(x,0)u0v(x,0)u0(lx)
所以,u(x,t)是原混合问题的解的充要条件是v(x,t)是如下混合问题的解: 2vt(x,t)avxx(x,t)0(0xl,(*)
v(0,t)0,v(l,t)0v(x,0)u(x)0lt0)
用分离变量法求解(*),由分离变量法的标准步骤可得:
X(x)Xn(x)sinnxl,naT(t)Tn(t)Anexp(t),l2
v(x,t)Tn1n(t)Xn(x)n1nxna,Anexp(t)sinll2代入初始条件可得:u0(l2lx)v(x,0)ln1Ansinnxlnxl
由Xn(x)的正交性可得:An
An0(u0(nlx))sindx,2n((u0)(1)(u0))(n1,2,),2所以,v(x,t)n1nxnan((u0)(1)(u0))exp(t)sinnll2
u(x,t)v(x,t)(lx)。
uxx(x,y)uyy(x,y)0(0xa,9.求解 u(x,0)x(xa),limu(x,y)0yu(0,y)0,u(a,y)0y0)。
解:用分离变量法:设给定的定解问题中的微分方程有如下满足齐次边界条件的分离变量形式非零解:
u(x,y)X(x)Y(y),则
uxx(x,y)X(x)Y(y),uyy(x,y)X(x)Y(y),uxx(x,y)uyy(x,y)X(x)Y(y)X(x)Y(y)0,X(x)X(x)Y(y)Y(y)X(x)X(x)0,Y(y)Y(y)0,
u(0,y)X(0)Y(y)0X(0)0,u(a,y)X(a)Y(y)0X(a)0,所以,X(x)为如下常微分方程边值问题的解函数:
2X(x)X(x)0nn
,X(0)0,X(a)0aX(x)Xn(x)sinnyanxa,从而有:Y(y)Yn(y)Anexp(又由另一个边界条件可得:
nya)Bnexp()(n1,2,)
(limun(x,y)limXn(x)Yn(y)0An0Yn(y)Bnexpyynya),设原定解问题的解函数是u(x,y)n1un(x,y)n1Bnexp(nya)sinnxa,则
u(x,0)x(xa)x(xa)n1Bnsinnxa
Bna2a0x(xa)sinnxandx22aannya333((1)1)n(n1,2,),所以,u(x,y)10.求解边值问题:
4a23n1(1)1n3exp()sinnxa。
uxx(x,y)uyy(x,y)0(0xa,
u(0,y)0,u(a,y)0xxu(x,0)0,u(x,b)sinaa0yb)。
解: 用分离变量法:设给定的定解问题中的微分方程有如下分离变量形式的满足齐次边界条件的非零解:
u(x,y)X(x)Y(y),则有:
uxx(x,y)X(x)Y(y), X(x)X(x)Y(y)Y(y)uyy(x,y)X(x)Y(y),0X(x)X(x)0,Y(y)Y(y)0,u(0,y)X(0)Y(y)0X(0)0,同理 X(a)0,所以,X(x)是如下二阶常微分方程边值问题的解函数:
2X(x)X(x)0nn
,X(0)0,X(a)0aXn(x)sinnyanxa,Y(y)nY(y)0Y(y)Yn(y)Ancoshny,Bnsinha设原定解问题的解为:u(x,y)n1(AncoshnyaBnsinhnya)sinnxa,则
0u(x,0)n1AnsinnxaAn0(n1,2,),xasinxa2au(x,b)nban1aBnsinhnbasinsinnxadx,所以,Bn(sinh)1xa0sinxanxanb2
sinha11(1)n1(1)n(n1)2(n1)2(n2,3,)
axbxxb
B1(sinh)1sinsindx2sinh0aaaaaa21。
所以,原定解问题的解函数为u(x,y)n1Bnsinhnyasinnxa,其中的Bn由以上式子给出。11.求解边值问题
uxx(x,y)uyy(x,y)k(0xa,
u(0,y)0,u(a,y)0u(x,0)0,u(x,b)00yb),提示:令u(x,y)v(x,y)w(x),而w(x)满足条件w(x)k,w(0)w(a)0。解:令w(x,y)k2x(xa),v(x,y)u(x,y)w(x,y),则
vxx(x,y)uxx(x,y)wxx(x,y)uxx(x,y)k,vyy(x,y)uyy(x,y)wyy(x,y)uyy(x,y)
所以,uxx(x,y)uyy(x,y)kvxx(x,y)vyy(x,y)0,u(0,y)0,u(a,y)0v(0,y)0,v(a,y)0,u(x,0)0,u(x,b)0v(x,0)k2x(xa),v(x,b)k2x(xa)
所以,u(x,y)是原定解问题的解的充要条件是v(x,y)是如下定解问题的解: vxx(x,y)vyy(x,y)0(*)v(0,y)0,v(a,y)0,kkv(x,0)x(xa),v(x,b)x(xa)22用分离变量法求解(*),由分离变量法的标准步骤可得:
v(x,y)X(x)Y(y)X(x)X(x)0,n
n,a2Y(y)Y(y)0,Xn(x)sinnxa,nynyYn(y)Anexp()Bnexp()
aa(n1,2,),v(x,y)vn(x,y)Xn(x)Yn(y)设(*)的解函数为v(x,y)n1(Anexp(nyak2)Bnexp(nya))sinnxa
则
v(x,0)n1(AnBn)sinnxa1x(xa),v(x,b)n1(AnDnBnDn)sinnxa,(其中 Dnexp(nba))
若记
Cna2ak20x(xa)sinnxadx2k2aa2n333((1)1),31nb)1CnAnexp(ABCannn则有: ,11ADBDCnbnbnnnnnBexp()exp()1Cnnaa 12 其中,An,Bn,Cn,Dn由以上各式给出。而题目所给的定解问题的解函数为
u(x,y)v(x,y)w(x,y)v(x,y)12.求解边值问题
uxx(x,y)uyy(x,y)0(0xa,
u(x,0)0,u(x,b)0u(0,y)y(yb),u(a,y)00yb)k2x(xa)。
解:用分离变量法求解此定解问题:设u(x,y)X(x)Y(y),由分离变量法的标准过程
nyn可得
n,Yn(y)sinX(x)Y(y)bbX(x)nX(x)0X(x)Xn(x)Anexp(nxb)Bnexp(nxb)(n1,2,)X(x)Y(y)2设原定解问题的解函数为
u(x,y)n1Xn(x)Yn(y)n1(Anexp(nxb)Bnexp(nxb))sinnyb,则由关于x的边界条件可得:y(yb)u(0,y)2bn1(AnBn)sinnyb,AnBnb0y(yb)sinnybdy
0u(a,y)nabn1(Anexp(nabnab1b)Bnexp(nab))sinnyb,Anexp(所以
An
Bn2b)Bnexp(2nabb)0,y(yb)sin)1)12b(exp()1)nyb0dy,nybdy,exp(2na)(exp(2nabb0y(yb)sin所以,u(x,y)所以,……。
13.求解混合问题
(An1nexp(nxb)Bnexp(nxb))sinnyb
3x3at2u(x,t)au(x,t)sinsinxxtt2l2l
u(0,t)0,ux(l,t)0u(x,0)0,u(x,0)0t(0xl,t0)。
解:用分离变量法求解此混合问题:设原给定的混合问题中的微分方程对应的齐次方程有如下分离变量形式的满足边界条件的非零解:
u(x,t)X(x)T(t)ux(x,t)X(x)T(t),uxx(x,t)X(x)T(t),ut(x,t)X(x)T(t),utt(x,t)X(x)T(t),utt(x,t)a2uxx(x,t)0
X(x)X(x)0, 由边界条件可得:u(0,t)X(0)T(t)0X(0)0,ux(l,t)X(l)T(t)0X(l)0,所以,X(x)是如下边值问题的非零解函数:
X(x)X(x)0
X(0)0,X(l)0X(x)X(x)T(t)aT(t)2
(2n1)求解此问题,可当n时,问题有非零解,其解函数集构成一个
2l2一维线性空间,它的一个基向量函数为X(x)Xn(x)sin令
fn(t)2l(2n1)x2l2lsin,dx,l0f(x,t)Xn(x)dx,fn(t)0,l0sin3x2lsin3at(2n1)x2l则
f2(t)sin3at2l(n1,3,4,5,)
令{Tn(t)}为如下初值问题的解函数: T(t)na2T(t)fn(t)
T(0)0,T(0)0(t0),(1)
则Tn(t)0(n1,3,4,5,),对于n=2,可用常数变易法来求:
T(t)2aT(t)0T(t)Acos设(1)的解函数为 T(t)A(t)cos则 T(t)A(t)cos令
A(t)cos3at2lB(t)sin3at2l3at2l3atB(t)sin3a2l2l3at2lBsin3at2l,3at2lB(t)cos3at2l)
(A(t)sin3at2lB(t)sin3at2l0,14 则
T(t)3a2l3a2l(A(t)sin3at2l3atB(t)cos),2lT(t)(A(t)sin3at2lB(t)cos3a2l3at3at3at3a)B(t)sin)(A(t)cos2l2l2l2l3at2l3at B(t)cos)f2(t),2l2
T(t)2a2T(t)f2(t)(A(t)sin3at3at(t)cos(t)sinAB02l2l也就是:
,3a3at3at3at(A(t)sinB(t)cos)sin2l2l2l2l求解此线性方程组得:A(t)22l3asin23at2l,B(t)2l3asin23at2lcos3at2l,3atll
A(t)sintc1,l3a3a3atl B(t)cosc2,l3a所以,(1)的解为:
3atl3at3at3atl
T(t)T2(t) tcosc1cosc2sinsin3a2l3a2l2l2l2由初始条件T(0)0,T(0)0可得:c10,2l22lc2,3a3at2l2所以,T2(t)3asin3at2ll3atcos,所以,题目所给的定解问题的解函数为:
u(x,t)14.求解混合问题
n12l23atl3atXn(x)Tn(t)sintcos(3a)22l3a2l3xsin。2l2x2u(x,t)au(x,t)sin(0xl,xxttl
u(0,t)0,u(l,t)03x2xu(x,0)2sin,u(x,0)sintllt0)。
解:作函数变换v(x,t)u(x,t)w(x),其中w(x)为待定函数,则
vtt(x,t)utt(x,t),vt(x,t)ut(x,t),vxx(x,t)uxx(x,t)w(x),22
vtt(x,t)avxx(x,t)utt(x,t)a(uxx(x,t)w(x))
utt(x,t)auxx(x,t)aw(x),15 设u(x,t)是原定解问题的解函数,2xl取aw(x)sin222xl,则有: 0,即w(x)sinl2a222vtt(x,t)avxx(x,t)utt(x,t)auxx(x,t)aw(x)sin2xl aw(x)0,2而
v(0,t)u(0,t)w(0)000,3xlv(l,t)u(l,t)w(l)0
v(x,0)u(x,0)w(x)2sin2xl2xl,sin2al2
vt(x,0)ut(x,0)sin,所以,v(x,t)为如下定解问题的解函数: v(x,t)a2v(x,t)0ttxx(*)
v(0,t)0,v(l,t)03xlv(x,0)2sinl2a(0xl,2xsin,l2t0),vt(x,0)sin2xl用分离变量法求解此定解问题:由分离变量法的标准过程可得: n
n,l2X(x)Xn(x)sinnatlBnsinnatlnxl,,T(t)Tn(t)Ancos设(*)的解函数为
(n1,2,)
v(x,t)n1un(x,t)n1(AncosnatlBnsinnatl)sinnxl,由初始条件可得:2sin3xl2xlv(x,0)sin2al22n1Ansinnxl
l可得: A10,A2,A32,2aAn0(n4,5,)
natllna
vt(x,t)n1nal(AnsinnatlnxlBncos)sinnxl,sin2xlvt(x,0)n1nalBnsinB2,Bn0(n1,3,4,5,)
2atl2at2x3at3xl所以,v(x,t)(,cossin)sin2cossinl2allll2a2所以,题目所给的定解问题的解函数为u(x,t)v(x,t)w(x)。15. 求解混合问题
2x2sinx(0xl,utt(x,t)auxx(x,t)l
u(0,t)t,u(l,t)sintu(x,0)0,u(x,0)(为常数)tt0)。
[注]:此定解问题中的微分方程非齐次项中的sinx应为sint,才能得到书中答案。
解:先将边界条件齐次化:令v(x,t)u(x,t)((sintt)t),lx则
vtt(x,t)utt(x,t)xlsint,2vxx(x,t)uxx(x,t),若u(x,t)是原定解问题的解函数,则
vtt(x,t)avxx(x,t)utt(x,t)2xl2sintauxx(x,t)
xl22
utt(x,t)auxx(x,t)0lsint0,2tt)t)tt0,v(0,t)u(0,t)((sintt)t)tt0,v(l,t)u(l,t)((sinll
v(x,0)u(x,0)00,vt(x,0)ut(x,0)(xl(cos*0))0,所以,v(x,t)是如下定解问题的解函数:
vtt(x,t)a2vxx(x,t)0
v(0,t)0,v(l,t)0v(x,0)0,v(x,0)0t(0xl,t0)v(x,t)0,所以,原定解问题的解函数为 u(x,t)xl(sintt)t
utt(x,t)a2uxx(x,t)3x2tex16. 求解 ux(0,t)t,ux(l,t)u(l,t)tu(x,0)0,u(x,0)1ext(0xl,t0)。
解:作如下函数变换:v(x,t)u(x,t)t(1ex)u(x,t)ttex,若u(x,t)是原定解问题的解函数,则经验证可得:v(x,t)是如下定解问题的解函数: vtt(x,t)a2vxx(x,t)3x2(1a2)tex
vx(0,t)0,vx(1,t)v(1,t)0v(x,0)0,v(x,0)0t(0x1,t0)
用分离变量法求解此定解问题:设v(x,t)X(x)T(t),T(t)aT(t)2由分离变量法的标准过程可得:
X(x)X(x)X(x)X(x)0,vx(0,t)0,vx(1,t)v(1,t)0X(0)0,X(1)X(1)0 由X(x)所满足的方程可得:X(x)c1cosxc2sinx,由边界条件可得:c20,0,取c11,则得X(x)cos
X(1)X(1)0sincos02所以,nn,X(x)Xn(x)cosnxx,ctg,(n1,2,),其中,n是方程ctg的所有正解。因为
10cosnxdx22100.5(1cos2nx)dx0.5(1sinn),2令
fn(t)1sinn21sin2210f(x,t)cosnxdx
1n0((3x)(1a)te22x)cosnxdx
4sinn(1sinn)3n22(1a)sinn1sinn222tbncnt
则
f(x,t)n1fn(t)cosnx,设原定解问题的解函数为v(x,t)Tn12n(t)cosnx,则
vttavxx2(Tn1n(t)aT(t))cosnx2nn1fn(t)cosnx,22从而有:
Tn(t)anTn(t)fn(t)(n1,2,),由初始条件可得:v(x,0)vt(x,0)0Tn(0)Tn(t)0,所以,Tn(t)为如下初值问题的解函数: 22Tn(t)anTn(t)fn(t)
Tn(0)0,Tn(0)0(t0)
22用常数变易法:Tn(t)anTn(t)0Tn(t)AncosantBnsinant,设此边值问题的解为: Tn(t)An(t)cosantBn(t)sinant,A(t)cosatB(t)sinat0nnnn经简单推导得: ,1A(t)sinatB(t)cosatf(t)nnnnnan1A(t)fn(t)sinantnan解此线性方程级:
1Bn(t)fn(t)cosantan积分并利用初始条件可得:
cn1A(t)((bct)cosatb)sinantnnnn23nanan
,cn1Bn(t)(bncnt)sinant(cosant1)23anan
Tn(t)An(t)cosantBn(t)sinant
1anbn2bncnt1an2(bncosantcnansinant)
an21cosantcnan21tsinatn an所以,u(x,t)Tn1n(t)cosnx,其中的Tn(t)、bn、cn和n均由以上各式给定。[注]课本上的答案为此处的a=1。
ut(x,t)a2uxx(x,t)0(0xl,17. 求解 ux(0,t),ux(l,t)u(x,0)A(A,为常数)t0)。
解:设u(x,t)是原定解问题的解函数,作函数变换v(x,t)u(x,t)x,19 则
vt(x,t)ut(x,t),vx(x,t)ux(x,t),vxx(x,t)uxx(x,t)
vx(0,t)ux(0,t)0,vx(l,t)ux(l,t)0,v(x,0)u(x,0)xAx,所以,v(x,t)是如下定解问题的解函数:
vt(x,t)a2vxx(x,t)0(0xl,t0)
vx(0,t)0,vx(l,t)0
v(x,0)Ax用分离变量法求解此定解问题:设v(x,t)X(x)T(t)为微分方程的满足齐次边界条件的非零解函数,则将v(x,t)代入方程后化简可得:
T(t)aT(t)X(x)X(x)T(t)aT(t)0,2X(x)X(x)0,vx(0,t)0,vx(l,t)0X(0)0,X(l)0,所以,X(x)为如下边值问题的非零解函数:
2nnX(x)X(x)0(0xl)lX(0)0,X(l)lX(x)X(x)cosnxnl(n0,1,2,)
将n代入T(t)的方程可得:
na
T(t)a2nT(t)0T(t)Tn(t)Bnexp(t)lnxna所以,vn(x,t)Tn(t)Xn(x)Bnexp(。t)cosll22(n0,1,2,),设
v(x,t)n0nxna,Bnexp(t)cosll2则由初始条件可得:Axv(x,0)1l2ln0Bncosnxl
可得:
B0
Bn)0l(Ax)dxA12l,(n1,2,),nx2ln(Ax)cosdx(1(1))220lln 20 所以,v(x,t)A
12ln12ln22nxna。(1(1))exp(t)coslln2ut(x,t)a2uxx(x,t)f(x)(0xl,18. 求解 u(0,t)A,u(l,t)B(A,B为常数)u(x,0)g(x)t0)。
解:设F(x)(0xx0f(x)dx)dx,w(x)1a2F(x)(AB)aF(l)al22xA,1a2
v(x,t)u(x,t)w(x)vt(x,t)ut(x,t),vxx(x,t)uxx(x,t)
vt(x,t)a2vxx(x,t)ut(x,t)a2uxx(x,t)f(x)0,1a1a22f(x),v(0,t)u(0,t)w(0)AF(0)(AB)aF(l)al2220A0,v(l,t)u(l,t)w(l)BF(l)(AB)aF(l)al2lA0,v(x,0)u(x,0)w(x)g(x)w(x),所以,v(x,t)是如下定解问题的解函数:
vt(x,t)a2vxx(x,t)0
v(0,t)0,v(l,t)0v(x,0)g(x)w(x)(0xl,t0),用分离变量法可求得:
v(x,t)其中,Ann1nxna,Anexp(t)sinll(g(x)w(x))sin22llnxl20dx(n1,2,)。
所以,u(x,t)n1nxnaAnexp(w(x)。t)sinll21.在扇形区域内求解边值问题
u0(ra,0)
u(r,0)0,u(r,)0。
u(a,)f()解:由极坐标下的Laplace算子表达式可知:
1u1u2
u0rurrruru0。r22rrrr2用分离变量法求解此定解问题:设u(r,)R(r)(),代入以上微分方程化简后可rR(r)rR(r)R(r)2得
()()2:
()()0,rR(r)rR(r)R(r)0
u(r,0)R(r)(0)0(0)0, u(r,)R(r)()0()0,所以,()是如下边值问题的非零解函数:
2nn()()0
(0)0,()0()sinnxn(n1,2,),2n/n/Bnr
rR(r)rR(r)nR(r)0R(r)Rn(r)Anr,n/又显然有:R(0)Bn0,也就是:Rn(r)Anr,所以,un(r,)Rn(r)n()Anrn/sinnsin,n设原定解问题的解函数是 u(r,)n1Anrn/n/,由关于r的边界条件可得:f()u(a,)其
n1Anasinn,中
Anan/20f()sinn2d(n1,2,),n/nr所以,u(r,)f()sindn10asinn。
u0(1r2,0)22 求解边值问题
u(1,)sin,u(2,)0。
u(r,0)0,u(r,)0解:由极坐标下的Laplace算子表达式可知:
1u1u20rurrruru0
ur22rrrr
2用分离变量法求解:设u(r,)R(r)()代入方程中并化简得:
rR(r)rR(r)R(r)2
r2R(r)rR(r)R(r)0,()()()0()
u(r,0)0,u(r,)0(0)0,()0,()()0
(0)0,()02n2nn()()sinnn(n1,2,),将nn2代入R(r)所满足的方程可得:
r2R(r)rR(r)n2R(r)0R(r)Rn(r)AnrnBnrn,n设原定解问题的解函数为 u(r,)Rn1(r)n()(An1nrBnrnn)sinn,nn0u(2,)(An2Bn2)sinnn1由r的边界条件可得:
,sinu(1,)(AnBn)sinnn1容易得到:
AnBn0(n2,3,),11A12A2B10
3,14B11A1B13所以,u(r,)13r43r1sin。2(ra)uxxuyyy23. 求解边值问题 222uraxy,rxy解:作函数变换 v(x,y)u(x,y)112y,24则有:
vxx(x,y)uxx(x,y),vyy(x,y)uyy(x,y)y 此时,有:
vxxvyyuxxuyyyyy0,所以,v(x,y)是如下边值问题的解函数:
222 23 vxxvyy0(ra)
14222vxyy,rxy12ra将此定解问题由直角坐标改为极坐标:
r2vrrrvrv0(ra)
1424v(a,)acossinasin12(xrcos,yrsin),用分离变量法求解此定解问题:设v(r,)R(r)F(),由分离变量法的标准步骤rR(r)rR(r)R(r)2容易得到:
F()F()02,rR(r)rR(r)R(r)0F()F()由v(r,)的实际意义可知:F()是以2为周期的周期函数,R(0) 所以
nn2,F()Fn()AncosnBnsinn(n0,1,2)
22nnn
rR(r)rR(r)nR(r)0R(r)c1rc2r,letRn(r)r,n设
v(r,)Rn0(r)Fn()(An0nncosnBnsinn)r
由关于r的边界条件可得:v(a,)112(An04ncosnBnsinn)a,n而
v(a,)acossin
所以,A013213242asin
12412acos219644a412asin21a,B222196acos4,4a,A224,A4,其余的An、Bn的值均为零。所以,v(r,) u(r,)1324132ar(242124acos212212sin2)1964196rcos4,112rsin。
444ar(124acos22sin2)rcos4u0(ra,0)224.求解边值问题 ur(a,)f()。
u(r,0)0,u(r,)02解:因为其自变量的取值区域是扇形区域,所以可在极坐标系下用分离变量法求解此定 24 解问题,因为,u1rrrur1ur2220,设 u(r,)R(r)(),求出其各阶偏导数并代入方程后化简可得:
rR(r)rR(r)R(r)2
r2R(r)rR(r)R(r)0 ()()()0()(由u(r,)关于的边界条件可得
(0)0,2)0
()()0n4n2所以
(0)0,()0n()sin2n2(n1,2)
r2R(r)rR(r)4n2R(r)0RRn(r)Anr2nBnr2n
u(0,)Rn(0)Rn(r)Anr2n
设原定解问题的解函数为
u(r,)An1nr2nsin2n,则
ur(r,)2nAn1nr2n1sin2n,由边界条件得
f()ur(a,)从而有:
An2na2n12nAn1na2n1sin2n
/20f()sin2nd
(1)
所以,原定解问题的解函数为u(r,)其中的系数由(1)式给出。
An1nr2nsin2n,uxy(ra,0)225.求解边值问题
ur(a,)f()
222u(r,0)0,u(r,)0,rxy2解:设w(x,y)112xy(xy),作函数变换v(x,y)u(x,y)w(x,y),22则
vvxxvyyuxxuyy(wxxwyy)0 在极坐标下:
v(r,)u(r,)w(r,)u(r,)124rsin2,25
vr(r,)ur(r,)
vr(a,)ur(a,)经验算得知:
v(r,0)0,v(r,1616rsin2,asin2,332)0,所以,v(r,)为如下边值问题的解函数:
21v1v(r)20v2rrrr13v(a,)f()asin2r6v(r,0)0,v(r,)02(ra,02)
用分离变量法求解,设v(r,)R(r)()代入方程并化简得:
rR(r)rR(r)R(r)2
r2R(r)rR(r)R(r)0,()()()0()由关于的边界条件可得:(0)0,(2)0,(n1,2,),2由此可得: n4n,n()sin2n222n2n
rR(r)rR(r)4nR(r)0RRn(r)AnrBnr,v(0,)R(0)Rn(r)Anrn2n。
设
v(r,)Rn13(r)n()An1nr2nsin2n,则
f()16asin2vr(a,)22nAn1na2n1sin2n,由可求得: v(r,)An1nr2nsin2na12rsin2,2其中,An2na2n1/20f()sin2nd,124rsin2。
u(r,)v(r,)
第三篇:积分变换与数理方程报告
积分变换与数理方程
班级:电信09103班 学号:200911020309 姓名:何双来
《积分变换与数理方程》学习总结报告
这个学期我们开了《积分变换与数理方程》这门课。这个课是为大三学习《信号与线性系统分析》做准备而开的。
现在,信号与系统的概念已经深入到人们的生活和社会的各个方面。手机、电视机、通信网、计算机网等已经成为人们常用的工具和设备,这些工具和设备都可以看成系统,而各种设备传送的语音、音乐、图像、文字等都可以看成信号。所以《信号与线性系统分析》这门课非常重要,已经成为电子信息类专业的基础必修课。然而,这门课程并不是那么好学,它里面涉及到很多高等数学的知识。要学习这门课程必须有较好的高等数学知识,并且能够运用这些数学知识解决实际问题。除此之外,还要求学生有较强的阅读理解能力,因为本课程的教材里面有很多抽象的概念、定义和公式。总的来说,是运算量大,内容多阅读量大,理解能力要求高。要在一个学期内学好这门课程并不是一件容易的事情。因此,为了减轻大三的时候学习这门课程的负担,我们开设了《积分变换与数理方程》这门课,主要讲授的是《信号与线性系统分析》中三数学变换和其它一些与数学运算有关的知识,目的是在上《信号与线性系统分析》课之前,让学生提前接触这门课程,以减少大三学习这门课程时难度。
经过这个学期对《积分变换与数理方程》这门课程的学习,我学到了很多东西,下面就对我所学到的东西做一个汇总。
一、首先,是信号的概念。信号是信息的一种表示方式,通过信号传递信息。信号有一维信号,也有n维信号,而本课程只讨论一维信号。信号根据不同的分类方式可以分为连续信号和离散信号,也可分为周期信号与非周期信号,又可分为实信号和复信号,还可分为能量信号和功率信号。此外,信号还可以进行某些基本运算,包括加法和乘法运算、反转和平移和尺度变换。
二、在这门课程中我还学习到了一些和信号分析与处理有关的基本的常见的函数。
(1)阶跃函数以及其图像(2)冲击函数及其图像
0,t0def1(t)limn(t)2,t0 (t)limpn(t)
nn1,t0def(t)t o o (t)t 阶跃函数与冲击函数的关系如下:
(t)d(t)dtt
(x)dx (t)
三、此外,我还学到了一些对函数的运算和函数的三大变换。
1、卷积积分。卷积方法在信号与系统理论中占有重要地位。卷积积分的定义如下:
一般而言,如有两个函数f1(t)和f2(t),积分
f(t)=
f1()f2(t)d
(2.3 —
7)
称为f1(t)与f2(t)的卷积积分,简称卷积。式(2.3 — 7)常记作
f(t)=f1(t)*f2(t)=f1()f2(t)d
下面是一些常用函数的卷积积分:
(1)函数与冲击函数的卷积:
f(t)(t)(t)f(t)()f(t)df(t)
f(t)(tt1)(tt1)f(t)f(tt1)
f(tt1)(tt2)(tt1)f(tt2)f(tt1t2)
(tt1)(tt2)(tt2)(tt1)(tt1t2)(2)常用卷积积分:
①f(t)'(t)f'(t)②f(t)(t)f(t)
③f(t)(t)
2、傅立叶变换(1)、正交函数集
如有定义在(t1,t2)区间两个函数1(t)和2(t),若满足
1(t)2(t)dt0
t1t2tf()d ④(t)(t)t(t)
则称1(t)和2(t)在区间(t1,t2)内正交。
若有n个函数1(t),2(t),…,n(t)构成一个函数集,当这些函 在区间(t1,t2)内满足
t2t1i(t)j(t)dt0,ijki0,ij 式中ki为常数,则称此函数集为在区间(t1,t2)的正交函数集。
(2)、周期信号的频谱。如前所述,周期信号可以分解成一系列正弦信号或指数信号之和,即
f(t)A02n12An1cos(ntn)jn 或 f(t)Fennjnt
其中FnAnejn|Fn|e。
(3)、非周期信号的频谱。为了描述非周期信号的频谱特性,引入了频谱密度的概念。令
F(j)limFn1TlimFnTTT 称F(j)为频谱密度函数,简称频谱函数。
对于任意一个非周期信号的时间函数f(t)有
defF(j)limFnTTdeff(t)edjtdt(4.4 — 4)
f(t)12F(j)ejt(4.4 — 5)式(4.4 — 4)称为函数f(t)的傅里叶变换,式(4.4 — 5)称为函数F(j)的傅里叶逆变换。F(j)称为f(t)的频谱密度函数或频谱函数,而f(t)称为F(j)的原函数。f(t)和F(j)的关系可以简记为f(t)F(j)。
(4)、奇异函数的傅里叶变换。
①冲激函数的频谱 F t)](t)t (t)ejtdt1
F(j)(1)1
②冲激函数导数的频谱 F ['(t)]j
F [(n)(t)](j)n ③符号函数的频谱 F [sgn(t)]2j1
④阶跃函数的频谱 F [(t)]()(5)、傅立叶变换的性质。
1()j j①线性,若 f1(t)F1(j),f2(t)F2(j)
则对任意常数a1和a2,有a1f1(t)a2f2(t)a1F1(j)a2F2(j)②对称性,若 f(t)F(j),则 F(jt)2f()③尺度变换,若 f(t)F(j),则对实常数a(a0),有
f(at)Fj|a|a1
④时移特性,若 f(t)F(j),且t0为常数,则有
f(tt0)ejt0F(j)
⑤频移特性,若,f(t)F(j),且0为常数,则有
f(t)ej0tF[j(0)](6)一般周期函数的傅立叶变换。
F [fT(t)]F FnejntnFF [ennjnt]2F(n)
nn
3、拉普拉斯变换
(1)、Fb(s)f(t)e1stdt(5.1 — 4)
stFb(s)eds(5.1 — 5)f(t)2jjj 式(5.1 — 4)和式(5.1 — 5)称为双边拉普拉斯变换对或复傅立叶变换对。式中复变函数Fb(s)称为f(t)的双边拉普拉斯变换(或象函数),时间函数f(t)称为Fb(s)的双边拉普拉斯逆变换(或原函数)。(2)、单边拉普拉斯变换
F(s)L [f(t)]defdef0f(t)e0,stdt
t0stf(t)=L 1[F(s)]12jjjF(s)eds,t0
其变换与逆变换的关系也简记作f(t)F(s)。(3)、拉普拉斯变换的性质
①线性,若 f1(t)F1(s),Re[s]1
f2(t)F2(s),Re[s]2
且有常数a1,a2,则 a1f1(t)a2f2(t)a1F1(s)a2F2(s),Re[s]max(1,2)
②尺度变换,若 f(t)F(s),Re[s]0 则 L [f(at)]0f(x)e(sa)xdxa1sFaa
③时移特性,若 f(t)F(s),Re[s]0 且有正实常数t0,则
f(tt0)(tt0)estF(s),Re[s]0
0 ④复频移特性,若 f(t)F(s),Re[s]0 且有复常数saaja,则
f(t)estF(ssa),Re[s]0a
a(4)、几种常用函数的拉普拉斯变换
①L ['(t)]s ② L [(t)]1 ③L [(t)] ⑤L [sin(t)]
3、z变换
(1)如果有离散序列f(k)(k0,1,2,),z为复变量,则函数
1s ④L [b0et]b0s
s22 ⑥L [cos(t)]ss22 ⑦L [sinh(t)]s22
F(z)kf(k)zk(6.1 — 7)
F(z)k0f(k)zk(6.1 — 8)
式(6.1 — 7)称为序列f(k)的双边z变换,式(6.1 — 8)称为序列f(k)的单边z变换。
(1)几种常用函数的z变换
①Z [(k)]1 ②Z [(km)]zm ③Z [(k)]④Z [(km)]zz1zmzz1zza
⑤Z [ak] 以上就是我这个学期所学到的内容。通过这个学期的学习,我对《信号与线性系统分析》这门课程中涉及到的数学运算进行了初步的学习。这将为我大三的时候进一步学习《信号与线性系统分析》打下坚实的基础。到时候,我一定能把《信号与线性系统分析》这门课学好。
第四篇:关于欧拉方程变量代换后系数递推关系的一点总结
关于欧拉方程变量代换后系数递推关系的一点总结
光信1104 李号
我们知道,对于欧拉方程anxyn(n)an1xn1y(n1)a1xya0yf(x)'(a2,a3,,an不全为0可以通过变量代换xet或tlnx化简。本文主要介绍如何用低阶导数来表示高阶导数以及线性表示时的系数递推关系。
先用一个例子来说明我们要探讨的问题。
dydydy'2'''''
已知:xe,。,,求xy,xy,xy(此处均为对x的导数)23dtdtdtt23
显然,由xet可知tlnx,则
y'dtdx1xdydt
dydx2dydtdtdx1xdydtxy'y''dydxddx1x321dy1dy1dydt1dydydydy2''()()2()xy2222dxdxdxxdtdtxdtdxdtdtxxdtdtdydx322ddyd222y'''()ddxxdydt(4)22[12(dydtdydt22dydt)]2x3(dydt22dy1dy1)2(3)2dtxxxdtdtdydt22dy132
(dydt332)xy3'''dydt23332dydt
同理可求出xy4dydt446dydt3311dydt26dydt
我们把系数提出,如下排列: n=1 1 n=2 1-1 n=3 1-3 2 n=4 1-6 11-6 为了方便讨论,我们作出以下两点规定: i)m用“Bn”表示第n排第m列的数(显然nm);
ii)(1)n1(n1)!(1n)!即(-1)n!(-n)!
n由上文中的迭代求导不难得出下面三点规律:
1i)
Bn1;
nn1ii)
Bn(1n)Bn1;mmm1iii)BnBn1(1n)Bn1nm1 该规律可用数学归纳法归纳得出,限于篇幅,此处省去不证。显然,只要找出了Bnm的通式,就可以表达出xny(n)。
n(n)1nxyBdydtnnB2ndn1ydtn1Bnndydt。
1n为求Bnm,我们有两条路出发。一是由“Bn1”着手,另一个是由Bn着手。
1注意到Bnm既与Bnm1有关,也与Bnm有关,我们选择从Bnn着手。后面我们会看到,带入数1字“1”计算会因为讨论n的取值而使得表达式无法统一。
1
由Bnn(1n)Bnn可知: 1Bn(1n)Bn1(1n)(2n)Bn1(1n)(2n)(1)B1(1n)!nn1n11
Bnn(1n)!
显然我们可以把平行于主对角线的数看成一组数列。取m=n-1,则:
Bnn1Bn1(1n)Bn1n2n1n2Bn1(1n)Bn2(1n)(2n)Bn3(1n)(2n)(3)B2(1n)(2n)(3)(2)B(1n)!1nnn1n3212
(1n)!2n(1n)!3n(1n)!2(1n)!1
(1n)!(i211i)
取m=n-2,则
Bnn2Bn1(1n)Bn1
n3n421n2n3Bn1(1n)Bn2(1n)(2n)Bn3(1n)(2n)(4)B3(1n)(2n)(3)B3n2(1n)!1nn11ii21(1n)!n212nj1i21i(1n)!331ii21(1n)!221i
i21n(1n)![(j311j1i)]
i21n故由数学归纳法可求出Bmn(1n)!a1nm1{[1a1a2nm1a21a111a21()]}
a3nm11现在我们再从“Bn1”入手,看看会有什么情况。Bn1
BnBn1(1n)Bn1Bn11n
12(1n)(2n)(2)B2B2 22121(123n1)BnBn1(1n)Bn1 332n(1n)2
(1n)Bn1(2n)Bn2(3)B3B3 2223n(1n)22(n1)(2n)224(14)22B3
3可以看出求Bn2是需要B22的值,求Bn3时需要B33的值而且还要用立方和与平方和公式。当m较大时,需要Bm的值以及m次方和与m-1次方和的公式。更重要的是,具体化m后,表达式无法统一!因此可以看出,成列分布并不是该数组的真正特性,而平行于主对角线的分布才能使该数组统一。
原式可化为D(D-1)(D-2)(D-3)…(D-i+1)
D^n=d^ny/dx^n m
第五篇:“诉访分离”终结信访不信法
“诉访分离”终结信访不信法
2014-03-20 02:30:50 新京报
近年来,随着越来越多的社会矛盾以案件形式进入司法领域,出现了诉讼与信访交织、法内处理与法外解决并存的现状,导致少数群众“信访不信法”。
近日,中共中央办公厅、国务院办公厅印发了《关于依法处理涉法涉诉信访问题的意见》,提出实行诉讼与信访分离制度;强调建立涉法涉诉信访依法终结制度;坚决杜绝一切“拦卡堵截”正常上访人员的错误做法。
新京报讯 今后,涉诉涉法信访事项的上访群众不需要到信访部门上访了,直接到相关的政法部门反映问题即可。
中办国办近日印发的《关于依法处理涉法涉诉信访问题的意见》提出:实行诉讼与信访分离制度;建立涉法涉诉信访依法终结制度;对于涉法涉诉信访中的错案、瑕疵案应依法纠正。
解决涉法涉诉信访问题入法治轨道
去年11月,党的十八届三中全会通过的《全面深化改革若干重大问题的决定》首次提出改革信访工作制度,把涉法涉诉信访纳入法治轨道解决。《意见》正是根据《决定》的精神细化了规定。
此次改革涉法涉诉信访工作机制、依法处理涉法涉诉信访问题的总体思路是:改变经常性集中交办、过分依靠行政推动、通过信访启动法律程序的工作方式,把解决涉法涉诉信访问题纳入法治轨道,由政法机关依法按程序处理,依法纠正执法差错,保护合法信访、制止违法闹访,努力实现案结事了、息诉息访。
在实行诉讼与信访分离制度方面,《意见》提及,把涉及民商事、行政、刑事等诉讼权利救济的信访事项从普通信访体制中分离出来,由政法机关依法处理。
值得关注的是,在处理涉诉涉法信访案件的程序上,《意见》强调:对于已经进入法律程序处理的案件,各级政法机关应当依法按程序在法定时限内公正办结。对于经复议、审理、复核,确属错案、瑕疵案的,各级政法机关应依法纠正错误、补正瑕疵。
杜绝违法限制上访人员人身自由行为
本轮改革的一大核心是建立涉法涉诉信访依法终结制度。
《意见》称,对涉法涉诉信访事项,已经穷尽法律程序的,依法做出的判决、裁定为终结决定。对在申诉时限内反复缠访缠诉,经过案件审查、评查等方式,并经中央或省级政法机关审核,认定其反映问题已经得到公正处理的,除有法律规定的情形外,依法不再启动复查程序。
此外,对于“拦卡堵截”正常上访人员的做法,《意见》再次强调,坚决杜绝一切“拦卡堵截”正常上访人员的错误做法,坚决杜绝违法限制或变相限制上访人员人身自由的行为。
亮点1
实行诉讼与信访分离制度
意见:实行诉讼与信访分离制度。各级信访部门对到本部门上访的涉诉信访群众,应当引导其到政法机关反映问题;对按规定受理的涉及公安机关、司法行政机关的涉法涉诉信访事项,收到的群众涉法涉诉信件,应当转同级政法机关依法处理。
【解读】
诉访分离倒逼司法公正
中央政法委有关负责人:诉讼与信访分离的好处是信访群众直接到政法机关反映诉求,可以少走弯路,便于政法机关及时处理。
北京大学法学院教授傅郁林:司法是社会公正的最后一道防线。以往以信访这个行政方式来处理涉法涉诉这类司法问题,容易形成民众“信访不信法”的思维,对民众的法律意识、整个国家的司法权威会带来负面影响。
意见提出将诉讼与信访分离,实际上就是要解决这一问题。在以往涉诉涉法信访事项中,包括刑事、民商事、行政案件等多种类型,此次信访制度改革意在让民众意识到法治的重要性,让涉法信访回归司法程序,这一方面有助于更为专业的处理信访案件,同时对司法部门提升办案质量也起到推动作用。
亮点2
错案瑕疵案应依法纠正
意见:严格落实依法按程序办理制度。各级政法机关对于已经进入法律程序处理的案件,应当依法按程序在法定时限内公正办结。对经复议、审理、复核,确属错案、瑕疵案的,依法纠正错误、补正瑕疵。
【解读】
“纠错”能堵住程序漏洞
最高法立案一庭副庭长高莎薇:目前涉法涉诉信访事项有三种情形,第一种是确实出现法官枉法裁判需要再审改判的案件,第二种是有审判瑕疵的案件,第三种是当事人无理缠访的案件。
第一种和第三种是少数,第二种情形较多。对于第一种和第二种情形,意见明确了在法律程序内的案件法官要在法定时限内公正办结时,这对法官在审判此类案件时如何运用自由裁量权起到了监督作用。对于需要司法救济再审的案件,经过多部门审核后,确实发现错案和瑕疵案的,意见明确要纠正错误、补正瑕疵。这对涉法涉诉信访事项在程序上如何保证司法公正起到了非常好的作用,堵住了可能出现的程序漏洞。
亮点3
穷尽程序的判决为终结决定
意见:建立涉法涉诉信访依法终结制度。对涉法涉诉信访事项,已经穷尽法律程序的,依法做出的判决、裁定为终结决定。对在申诉时限内反复缠访缠诉,经过案件审查、评查等方式,并经中央或省级政法机关审核,认定其反映问题已经得到公正处理的,除有法律规定的情形外,依法不再启动复查程序。
【解读】
节约访民成本减少资源浪费
中央政法委有关负责人:以往一些涉诉涉法信访事项终而不结、无限申诉,反复启动法律处理程序。建立涉法涉诉信访依法终结制度,一方面减少了上访群众的负担,另一方面也可减少行政资源和司法资源的浪费。
吉林市保民律师事务所律师修保:以往很多涉诉涉法上访群众直接到信访部门上访,有的直接到北京的国家信访局上访,上访后信访部门一般会将信访事项转交给政法机关处理,这期间无论是时间成本、上访生活成本都是比较大的。建立涉法涉诉信访依法终结制度后,省去了一个环节,对于访民来说可以有效节约成本。
亮点4
办案质量将终身负责
意见:要完善执法司法责任制,严格落实办案质量终身负责制,健全执法过错发现、调查、问责机制,严格倒查执法办案中存在问题的原因和责任,严肃查处错案背后的执法不公、不廉等问题;要把加强执法公开、扩大群众参与、接受群众监督作为依法处理涉法涉诉信访问题的重要内容,以公开确保公正、促进息诉。
【解读】
信访问题突出倒查领导责任
中央政法委有关负责人:此次涉法涉诉信访改革一个重要内容是对于涉诉涉法信访事项,负责处理的政法机关要公正办案。如何在制度上予以强化?就是要严格落实办案质量终身负责制。
中央政法委要求各级政法机关在对信访事项办案时严查案件中存在的司法不公问题,对于群众诉求不及时受理、不按期办结、有错不纠的,要依纪依法追究办案人员和相关领导责任。对于信访问题突出的地方,还要倒查政法单位领导班子责任。同时,中央和省级政法机关将设立审核程序,对反映问题已经得到公正处理的,不再启动复查程序,对未得到公正处理的,要倒查追责。
■ 专家说法
“信访制度改革提升民众法律意识”
北京大学法学院教授傅郁林认为,此次涉诉涉法信访制度改革,实质上是涉诉涉法信访程序和强化纠错方式上的重大变革。
以往,涉诉涉法信访上访群众有冤错案件、瑕疵案件不知道去哪里上访,出现了信访部门、政法部门多头上访的现象。上访群众和办理信访事项的政法机关都浪费了很大的成本。上访群众会浪费时间、精力、金钱等成本,政法机关会在处理时增加时间、政务成本,信访事项的办理也会因此出现效率不高的问题。诉讼与信访分离后,无论对上访群众还是政法机关,都会减少程序,节约成本,办理的效率也会提升。
群众之所以上访,是因为所经历的案件大多数存在司法不公的问题,因此如何在改革过程中强化纠错方式尤为重要。此轮改革通过办案质量终身负责、上级政法部门审查监督、倒查追责机制强化了纠错方式,这对于维护司法公正将起到重要作用。
此外,涉诉涉法信访制度改革实际上是行政信访方式逐步退出司法领域,有利于维护司法权威,提升民众的法律意识。
涉法涉诉信访终结
涉法涉诉信访反映
向政法机关而不是向党政信访部门反映问题。
向有管辖权的政法单位反映问题,而不是多头访、越级访。
提醒:属于哪一级管的,到哪一级申诉;属于哪一个部门办的,到哪一个部门申诉。
审查、甄别
符合法律规定的信访事项进入复议、复核、再审程序处理。
不符合法律规定的信访事项,或是正在法律程序办理中,当事人直接上访的,政法机关依法不予受理。
提醒:反复缠访甚至违法闹访的,将受到依法处理。
依法按程序办理
对已经进入法律程序处理的涉法涉诉信访问题,政法机关依法按程序在法定时限内公正办结。
依法终结
已经穷尽法律程序的,依法做出的判决、裁定为终结决定。对于反复缠访缠诉的,经过案件审查、评查,由中央或省级政法机关审核,认定其反映问题已得公正处理的,除有法律规定情形外,依法不再启动复查程序。有关部门重点做好对信访人解释、疏导工作。
涉法涉诉信访三情形
第一种是确实出现法官枉法裁判需要再审改判的案件,占比少
第二种是有审判瑕疵的案件,占比较多
第三种是当事人无理缠访的案件,占比少
涉法涉诉信访占比下降
去年,中央政法委在涉法涉诉信访改革方面分批部署各省市试点,并于去年10月全面推开。经过1年时间的改革试点,涉法涉诉信访占到整个信访的比例下降明显。
要求一
畅通信访渠道。坚决杜绝一切“拦卡堵截”正常上访人员的错误做法;坚决杜绝违法限制或变相限制上访人员人身自由的行为。
要求二
提高基层化解能力。坚持分级分类处理,加大分流疏导力度,劝导当事人依法按程序反映问题,减少越级上访。
要求三
防止案件积压。政法机关要规范案件流程,加快案件流转,确保依法处理涉法涉诉信访问题有序高效进行。
要求四
严肃处理违法上访行为。坚持教育与处罚并重。对采取极端方式闹访、借上访之名煽动闹事的,依法严肃处理。
各地公安机关不得限制正常信访活动。同时对少数人在信访活动中以信访活动为名实施的违法犯罪行为,要区分情形,依法予以处置。处罚只是一种手段,不是最终目的,从根本上还是为了保护信访群众的合法权益,维护正常的信访秩序和社会秩序。——公安部信访办副主任周新 据新华社