变量分裂法(通用3篇)
变量分裂法 篇1
法律战, 是指依据国内法、国际法和国际惯例, 通过各种渠道所进行的有利于己而不利于敌的法律对抗活动。开展法律战, 有利于夺取法理优势, 争取国际社会的理解和支持, 打击和瓦解敌方意志, 从而更好地凝聚民心、军心, 为己方争取更大的军事行动空间, 赢得战略上的主动。2005年3月, 全国人大审议通过了《反分裂国家法》, 这正是我们做好反“台独”法律斗争准备的一项重要举措。《反分裂国家法》为我们遏制“台独”势力、促进祖国统一、妥善解决台湾问题提供了有力的法律武器, 我们应当深入研究并切实运用好这一法律武器, 打好反“台独”斗争的法律战。
一、制定《反分裂国家法》的法律依据
法律战是以法为兵、先法后兵, 是运用法律这一武器来争取法理上的主动、占据法律上的优势, 从而达到征服军心民意、掌握军事行动主动权、赢得国际社会理解和支持的目的。这就要求作为武器的法律之剑本身要过硬, 也就是说这些法律的制定要有着合理合法的依据。《反分裂国家法》的制定就有着充分的法理依据、明确的宪法依据和严格的立法程序。
(一) 制定《反分裂国家法》有着充分的法理依据
分裂国家的活动历来被各国视为国家和民族利益的大敌, 任何主权国家都不会容忍、更不会允许任何分裂国家的言行。从世界各国来看, 无论是大陆法系国家还是英美法系国家, 通过立法来解决国家统一问题都是一种通行的国际惯例。世界上任何一个主权国家都有义务也有权利采取法律手段来制止分裂、维护统一。近几年来, “台独”势力分裂国家的活动愈演愈烈, 妄图通过所谓“公投”、“宪改”搞“法理台独”, 改变大陆和台湾同属一个中国的事实, 把台湾从中国分裂出去。为了遏制“台独”分裂国家的活动, 维护祖国和平统一的前景, 全国人大通过立法形式维护国家主权和领土完整、促进祖国和平统一, 是完全必要的, 《反分裂国家法》的制定有着充分的法理依据。
(二) 制定《反分裂国家法》有着明确的宪法依据
我国宪法明确规定:“台湾是中华人民共和国的神圣领土的一部分。完成统一祖国的大业是包括台湾同胞在内的全中国人民的神圣职责。”这是我们制定《反分裂国家法》的宪法依据。《反分裂国家法》第一条就开宗明义地规定:“为了反对和遏制‘台独’分裂势力分裂国家, 促进祖国和平统一, 维护台湾海峡地区和平稳定, 维护国家主权和领土完整, 维护中华民族的根本利益, 根据宪法, 制定本法。”台湾是中国的一部分, 维护国家主权和领土完整是包括台湾同胞在内的全中国人民的共同义务, 对于祖国的和平统一大业, 我们寄希望于台湾人民, 并尽最大努力增进台湾民众的福祉。同时, 我们也绝不允许“台独”分裂势力以任何名义、任何方式把台湾从中国分裂出去, 对此我们将会采取一切必要措施维护国家的主权和领土完整。
(三) 制定《反分裂国家法》有着严格的立法程序
《反分裂国家法》是由全国人民代表大会制定的有关国家主权和领土完整的法律, 整个立法过程是十分慎重并且完全符合立法程序的。在立法准备阶段, 进行了全面细致的调查研究, 广泛征求社会各界的意见, 听取专家学者的意见, 然后提请全国人民代表大会常务委员会进行审议, 审议后根据委员们的意见做了修改, 形成一个草案, 再拿到全国人民代表大会全体会议上进行审议, 又根据全国人大代表的审议意见进行了必要的修改, 最后形成这部由全国人民代表大会表决通过的法律。
二、制定《反分裂国家法》的战略意义
(一) 遏制“台独”分裂势力
近年以来, 台湾当局不断加紧推行“台独”分裂活动, “台独”势力的所作所为, 一步步地把两岸关系推向了危险的边缘, 严重威胁着中国的主权和领土完整, 严重破坏和平统一前景, 严重损害中华民族的根本利益。《反分裂国家法》的出台, 不仅适时、必要, 而且必将对“台独”势力形成有效遏制。这从台湾当局对《反分裂国家法》的反应中就已经可见一斑了。全国人大将审议《反分裂国家法》草案的消息甫一公布, 台湾当局即有反弹, 或称要制订“反并吞法”, 或派出“宣达团”四处游说, 并极尽“妖魔化”这部法律之能事, 称其为“战争法”、“战争动员令”, 企图混淆视听, 误导民众, 这些都恰恰说明《反分裂国家法》击中了“台独”势力的要害、打到了“台独”势力的痛处。《反分裂国家法》的出台, 拉开了运用法律战手段打击“台独”势力的序幕, 使我们不再仅仅依靠舆论斗争手段、政治斗争手段揭露“台独”阴谋, 还能够积极运用法律斗争手段沉重打击“台独”分裂势力, 增强遏制“台独”分裂活动的威慑力和主动性, 从而更为有效地遏制“台独”分裂势力。
(二) 钳制外部干涉势力
台湾问题是历史遗留问题。由于种种原因, 尤其是外部势力的干涉, 台湾问题变得更为复杂, 两岸迄今尚未统一。而且, 在外部势力的纵容与支持下, “台独”分裂势力更是日益恶性膨胀。《反分裂国家法》的制定, 向全世界宣示了中国政府和中国人民反对“台独”分裂、促进国家统一的共同意志和坚定决心, 以国内法方式建构了国际社会普遍认知的一个中国政策的合法性, 占据了不容他国势力干涉我国内政的法理主动和法律优势, 为反对、阻止外部势力支持“台独”分裂势力提供了有力的法律武器, 特别是针对美国的“与台湾关系法”等, 实行“以法制法”, 从政治到法律层面排除外部势力插手干涉台湾问题。
(三) 赢得各界广泛支持
《反分裂国家法》本着以人为本、以民为重的精神, 从维护台湾同胞利益和福祉的角度出发, 积极谋求台湾地区的繁荣稳定, 积极推动海峡两岸的和平统一与共同发展, 必将赢得广大台湾民众的理解与支持。《反分裂国家法》旨在反对和遏制“台独”分裂势力分裂祖国, 维护国家主权和领土完整, 维护中华民族的根本利益, 必将赢得包括海内外炎黄子孙在内的绝大多数中国人的衷心拥护。《反分裂国家法》是一部反分裂、求稳定、求和平、求统一的法律, 对于维护台海乃至整个亚太地区的和平与稳定有着积极而深远的影响, 必将赢得国际社会的广泛赞同。《反分裂国家法》以法律的形式向国际社会宣示中国大陆与台湾的关系, 用法律手段处理一国的内政问题, 更容易得到国际社会的广泛认可和理解支持。
三、发挥《反分裂国家法》的战略作用
(一) 以法制“法”
所谓以法制“法”, 就是运用《反分裂国家法》这一法律武器与“台独”势力的“法理台独”作斗争、与外部势力的法律干涉作斗争。
“台独”分裂势力近年来不断谋求和推动“法理台独”, 妄图通过“公民投票”、“宪政改造”等方式, 利用所谓“宪法”和“法律”形式, 为实现“台独”分裂势力分裂国家的目标提供所谓“法律”支撑, 改变大陆和台湾同属一个中国的事实, 把台湾从中国分裂出去。对此, 我们应当积极运用法律武器予以坚决回击。台湾是中国领土不可分割一部分的地位从来都是明确的, 两岸同属一个中国的事实从来就未曾发生过改变, 这早已深得包括美国在内的国际主流社会的承认。《反分裂国家法》把我国一贯坚持和主张的、以及国际主流社会一贯认可和遵守的一个中国原则法律化, 明确规定“世界上只有一个中国, 大陆和台湾同属一个中国, 中国的主权和领土完整不容分割”, “台湾是中国的一部分”。这就使得一个中国原则从政策层面上升到了法律层面, 从而更具强制力和约束力, 换句话说, 这就意味着任何一部违反一个中国原则的法律都将是无效的。“台独”分裂势力妄图利用所谓的“宪法”和“法律”把台湾从中国分裂出去, 是行不通的。我们应当充分运用《反分裂国家法》这一法律武器, 有效开展反“台独”的法律对抗活动, 揭示所谓“法理台独”的荒谬性、欺骗性和非法性。
《反分裂国家法》在将一个中国原则法律化的基础上, 进一步清晰明确地规定“台湾问题是中国内战的遗留问题”, “解决台湾问题, 实现祖国统一, 是中国的内部事务, 不受任何外国势力的干涉”。这为我们排除外部势力干涉台湾问题提供了有力的法律支持。“外部势力的干涉”, 既包括政治干涉、经济干涉、军事干涉, 也包括法律干涉。比如, 美国单方面制定的与台湾关系法, 这无疑是一部地地道道的干涉他国内政的法律, 无论是从《反分裂国家法》的角度, 还是从国际法、中美外交关系的角度来看, 这部法律都是非法的。在反对和阻止外部势力插手干预台湾问题的斗争中, 我们应当积极运用各种法律武器, 扫除外部干涉势力恶意设置的阻挠我国统一大业的种种法律障碍。
(二) 以法遏“独”
所谓以法遏“独”, 是指运用《反分裂国家法》这一法律武器遏制“台独”势力的分裂活动。
《反分裂国家法》明确规定:“‘台独’分裂势力以任何名义、任何方式造成台湾从中国分裂出去的事实, 或者发生将会导致台湾从中国分裂出去的重大事变, 或者和平统一的可能性完全丧失, 国家得采取非和平方式及其他必要措施, 捍卫国家主权和领土完整。”这实际上是明确而清晰地为“台独”势力设下了不可逾越的法律红线, 警告“台独”分子不要心存幻想, 务必悬崖勒马, 以免把台湾民众带入万劫不复的险境。《反分裂国家法》划出“台独”红线, 可以让台湾民众清楚地认识到大陆反“台独”的决心, 正确地判断民进党“台独”政策的危害性, 可以让岛内政党、民众清楚知道如何反对“台独”、维护两岸和平, 从而堵死“台独”之路。在“台独”的模糊空间被完全抹煞之后, 台湾当局过去靠着一再试探大陆底线而步步推进的危险伎俩, 也将无用武之地。这些都将使得“台独”势力的活动空间日益狭小, 使得我们在对台法律斗争中占得先机, 从而有效地遏制住“台独”势力的分裂活动, 把对台斗争的主动权牢牢掌握在我们手中。另外, 用“非和平方式”取代过去的“使用武力”, 不但有利于增强岛内民众和国际社会的认知, 更增加武力以外的其他选项, 加大对台威慑的力度与空间, 增加回旋空间。同时, 对分裂国家行为的界定是通则性的, 避免划线过度明确而束缚手脚, 从而在解释和操作方面预留空间, 掌握两岸关系的主动权。以法律形式划出“台独”势力不可碰触的三条“红线”, 无疑加大了这一措施的强制力与可信度, 也增强了对“台独”势力的威慑力, 为台海地区的和平稳定构筑了一道坚强的、“台独”势力不可逾越的法律防线。
(三) 以法促“统”
所谓以法促“统”, 就是运用《反分裂国家法》这一法律武器积极宣传我们以最大的诚意、尽最大的努力争取和平统一的一贯主张, 积极争取台湾民众和国际社会对我们的理解与支持。
《反分裂国家法》以大量篇幅重点突出了以和平方式实现国家统一的主题, 主要表现在:一是以法律的形式确立实现国家统一的方式、前提, 即“坚持一个中国原则, 是实现国家和平统一的基础。以和平方式实现国家统一, 最符合台湾海峡两岸同胞的根本利益”。二是以法律的形式规定大陆在“维护台海地区和平稳定、发展两岸关系”方面采取的措施, 即“五个鼓励和推动”。三是以法律的形式规定实现和平统一的途径是“台湾海峡两岸平等的协商和谈判”, 两岸协商和谈判涉及六大议题, 即正式结束两岸敌对状态、发展两岸关系的规划、和平统一的步骤和安排、台湾当局的政治地位、台湾地区在国际上与其地位相适应的活动空间、与实现和平统一有关的其他任何问题。“五个鼓励和推动”及“平等协商和谈判”六大议题的法制化, 充分体现了“没有人比我们更希望通过和平的方式实现祖国统一”的诚意, 体现了“以人为本”的思想, 体现了“尽最大努力维护广大台湾同胞福祉”的善意。六大议题不仅与我们一贯的“在一个中国原则下, 什么都可以谈”的态度一致, 同时也使得议题更加具体化、扩大化。对于法律的这些规定, 我们要加大宣传力度, 消除因台湾当局对民众的误导而造成的民众对《反分裂国家法》的误识, 争取更多的台湾同胞对祖国和平统一大业的有力支持, 争取国际社会的广泛理解与支持, 促进两岸的和平统一。
变量分裂法 篇2
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。
变量分裂法 篇3
求二维随机变量函数的密度函数是概率论中的一个重要内容, 由于变量分布的差异性, 如何求又是一个难题, 一般没有统一的公式可循。一般教材中介绍了常用的方法, 即先求分布函数, 然后对分布函数求导就得密度函数, 但计算比较麻烦, 学生掌握困难很大。变量变换法和增补变量法相对于常用方法而言, 计算过程更加简捷。其中变量变换法许多学者都有研究, 而增补变量法甚少提及, 总结多年教学经验发现, 学生难以掌握这部分内容的精髓, 运用起来容易犯错, 特别表现在确定被积函数的积分区域上, 本文针对这个问题理清了一个简单通用的确定方法。
二增补变量法
增补变量法实质上是变量变换法的一种应用:为了求出二维连续随机变量 (X, Y) 的函数Z=g (X, Y) 的密度函数, 增补一个新的随机变量V=h (X, Y) 。先用变量变换法求出 (Z, V) 的联合密度函数pZV (z, v) , 再对pZV (z, v) 关于v积分, 从而得到关于Z的边际密度函数pZ (z) 。
问题: (1) 如何增补新随机变量? (2) 如何确定p (z, v) 中v的积分区域?
第二个问题, (Z, V) 的联合密度函数pZV (z, v) 的解析式中只含有z和v两个量, 对v积分时, 积分上下限一定是关于z的表达式。在二维连续随机变量 (X, Y) 的联合密度函数pXY (x, y) 的非零区域里, 画出函数Z=g (X, Y) 的曲线, 用曲线与非零区域相交点确定积分的上下限。
解:本例题作为示范题, 先记V=Y求出Z=X+Y的密度函数, 再换V=X用同样的方法求出Z=X+Y的密度函数, 如果两次求出的密度函数相同, 即验证方法的正确性。
由例题可以看出, 只要掌握了增补变量和确定积分区域的技巧, 增补变量法是一个极易掌握而且便于计算的方法。
增补变量法将比较难求的多维连续随机变量函数的密度函数, 转化为求变换后的两个变量的联合密度函数, 然后利用联合密度函数与边际密度函数之间的关系, 积分求出要求的变量的密度函数。相对于常用方法, 可简化运算。
摘要:二维连续随机变量函数的密度函数的计算是概率论教学中的一个重点, 更是一个难点, 其中增补变量法是一个简洁明了易掌握的方法, 但学生不能准确确定联合密度函数的积分区域, 本文针对这个问题给出了确定积分区域的方法。
关键词:连续随机变量,密度函数,增补变量法
参考文献
[1]茆诗松、程依明、濮晓龙编著.概率论与数理统计教程 (第二版) [M].北京:高等教育出版社, 2011:169~171
[2]王凡彬、严雳.多维随机变量函数密度的求解方法[J].内江师范学院学报, 2012 (10) :17~19
[3]余本国.一般二维连续型随机变量函数分布的讨论[J].华北工学院学报, 2004 (2) :94~96
[4]李思齐、李昌兴、柳晓燕.二维连续型随机变量函数的分布密度的计算[J].大学数学, 2011 (5) :162~166