常微分方程初值问题问题 式3中的x为什么不需要加绝对值号 错了吗

2007年3月阴山学刊;Mar.2007;V01.21;NO.1;第21卷第1期;ⅥNSHANACADEMIC;JOURNAL;关于微分方程解的几点注释;(包头师范学院数学科学学院,内蒙古包头01403;摘要:本文通过几个典型例子,指出了在微分方程的求;关键词:微分方程;求解过程;解题技巧;中图分类号:0175.1文献标识码:A文章编号:;微分方程是数学专业
2007年3月阴山学刊
Mar.2007
V01.21
第21卷第1期
ⅥNSHANACADEMIC
JOURNAL
关于微分方程解的几点注释
(包头师范学院数学科学学院,内蒙古包头014030)
摘要:本文通过几个典型例子,指出了在微分方程的求解过程中容易发生的一些错误,并给出了正确的求解方法。
关键词:微分方程;求解过程;解题技巧
中图分类号:0175.1文献标识码:A文章编号:1004―1869(2007)01--0016―02
微分方程是数学专业的一门专业基础课,学生在学习微分方程的过程中,除了要掌握基本概念和相关定律外,熟练掌握解题的方法和技巧也是不可缺少的一个环节。笔者从多年的教学过程中发现,在求解微分方程的过程中,
有一些容易忽视的细节,往往为初学者所忽略,结果是不
旷饯,如果在求解这个微分方程时不加绝对值符号,分离
变量后,两边积分,得:产:产
解出lny=lnx+C.。
即Y=elnHq=e。1X=C2X(其中C!=e。1)这里c2>0,实际上当C2<0时,Y=C2X也是此方程的解。因此,正确的解法是,
能正确地求解方程”“。我们通过以下几个例题加以说明,
以便与读者共享。1
第一类问题(丢解问题)
在求解一阶微分方程的过程中,常常容易发生丢解现
象,这类问题常常发生在方程两端同除以同一个因子,或
分离变量后,两边积分得:庠:庠,
者用分离变量法解方程的过程中。考察下面的两个例子:
解得:lnIyj=lnlxl+C.Y由上式得:酬=e叫十q=eqH,于是
Y=±eqx=CX(其中C=±e‘。)
例l:求解方程一dy+(3-2xy)dx=0
如果用jcz同除方程两端,可将原方程变为
一dy一型:一三,这是一个一阶线性非齐次方程,用常dxx
例:求解方程x兰~一√x2+y2=例3:求解方程x孚~,Y一0巧t五0
数变易法容易求出通解为Y=CX‘+一。但是,工s0也是这个方程的解,它不包含在通解当币,由于方程两端同除,这个因子,所以z兰0这个解很容易被丢掉。
r――――●r_――――_
解:用x同除方程两端,当x>O时,得到
例2:求解方程x√l―Y2dx+y41一x2dy=0
一,,I一/一
专;当x<o时得到
用分离变量法将原方程变为丽十丽一u
一手一㈤+当。令“:上,将原方程化为
厂―――=-厂―――_
,积分得出通解为、/l―x2+、/1一y2=c,其中0<c<2。
J.Y=1
(一】<x<1),(一1<y<1),
夕=一1x=一1
(一1<J<1)(一l<y<1)
丽2±了,两边积分,对于这两种情形可以分别得
厂―=_――_
y+√x2+y二=蹦
都是原方程的解,由于在分离变量过程中,需要假定l―x2≠0,l―Y二≠0,所以很容易将这四个解丢掉。
另外,有些微分方程在求解过程中,要用到积分公式
y+、/x。+Y。=一蹦畔o)。于是,我们看到这个方程
的积分曲线分成两个曲线族,一个位于',轴左侧,另一个位
于y轴右侧。
厂i―――了
兰=inH+c,这时千万不能忘记对x加绝对值符号,否…
2第二类问题(解的值域问题)
在用初等积分法求解微分方程时,有时会遇到对数运
则也会出现丢解情况。例如微分方程烹2=的通解是
收稿日期:2006―10一20
作者简介:陈向华(1954一),女,广西桂平人,副教授,研究方向:微分方程的应用和数学教育。
运算、三角函数运算、取绝对值运算以及开方运算等,这时要做认真、具体的分析,否则就会出现计算错误。
例如在解微分方程时经常用到积分公式
3第三类问题(通解中的任意常数问题)
众所周知,微分方程通解中含有一些任意常数,有
些时候,微分方程通解中的任意常数的取值是受到一定限制的。那么,只有准确的理解微分方程通解中的任意常数的取值限制,才能给出微分方程解的正确结果,我们通过下面例子说明。
例6:求解方程dy/dx+3y=-0
J一2In[xl+c,而在积分时就会遇到这样的问题:应该取J坐:ln.Y+f,还是』坐=In(一一)+c,或者是J一2InI.Y卜c’?为此,做如下分析:我们知道,微分方
程的解的定义总是离不开一定的区间,因此在考虑积分
这是一个一阶线性齐次方程,分离变量后,妙,
l’=一3救解此齐次方程,两边积分得:
J一的时候也要考虑积分是在自变量工的哪个区间上进
行。考察下面的例子:
IJ’I=一3x+cl,(cl∈尺)
取指数形式得到H=c2e。。,这里c2=e(1>0。
去掉绝对值得到:J一_+c2e一3、=∞一3。p删。
J?(一1):l。首
由上述求解过程得到的解J’=cal’,对常数c有一个限制f为,这是由求解过程中的数学运算决定的,因为在求解过程0eNy同除方程两端是假是》为的,但是注意
例4:求解初值问题÷dv:坠,
一.I’一
1一、.一
先将方程分离变量得到
、?d1.=二―二一d、|,两边积分得
二(-、2+,-2):p。
d、I+(.。由题设可知所求的初值问题
,d.r,
到r=0也是此齐次方程的解,所以C为这个限制可以去
掉。因此,该齐次方程的通解是y=ce。、,(c∈R
应考虑x<O,因此,应选取公式J__2Jn(叫J+c。于是
通过此例说明,对于常数的取值限制,是求解过程中某些数学运算的结果,在方程通解表达式中,问题本身的性质和意义可以将这个限制去掉。
所求通解为=(.、I。+.、‘。)=In(一一)+(.,代入初始条件得c=l。最后得到满足初始条件的解为
二(.、I一+、’一)=Ill(一.、I)+1。
又如微分方程2.’’等+2.x3’二=_、e。的通解为
例5:求解方程(.、’+3)d.v+cot.、.妙=0。用分离变
量法将原方程化为I/O斗3)d.)一---taILvdx,在积分时遇到这样的问题f
.1‘二=e一、?(÷.r二+f),由通解表达式可以看出如果要使解的
定义域非空,常数c必须满足(。s0。因而,通解中的常数f的取值限制,对讨论方程的通解起着非常重要的作用,它还可以描绘出解的性态。考察下面例子。
例7:讨论方程dy/dx=1120,2一1)的通解。将方程分离变量得:2dl?/0,2一1)=dx,
dx=一incosx,还是』tanx
=一ln(一co时),或者是』tattrdx=--InlcosxJ?同理,
在考虑积分f
dx的时候也要考虑积分是在自变量x
的哪个区间上进行。
如果是在区间(0,=)上考虑上述方程的解,那么应当是』tanx
dx=一In(cosx),于是得到Inl3十l’FIn
两边积分得到:InlO,一1)/0,+1)I--v+ci,即10’一1)/
O件1)1=e、“‘‘=c2e。,或0,一1)/0,+1)=±c2e。=ce.v,所以原方程的通解为
y=-(1+ce')/(1一ce。1
对于每个确定的负值c,(?)式给出一个定义于(一
(cosx)+c,从而得到13十l’I=ccosr(f>0),1i113+r=ccosx(f
不等于o)。但是因为.、’=一3也是方程的一个解,所以方
程的通解是.、’=一3+(’COS.Y((‘∈R)。
反之,如果是在区间(=,丌)上考虑上述方程的解,那么应当是ftall
oq+哟上的解,它整个位于条形区域.1<J<l内。对于每
一个正值c,(+)式给出了两个解,一个是定义于区间(一oo,lnc)上,而位于半平面y>l内的解;另一个是定义于区间(.Inc,+固上,而位于半平面l<一l内的解。当c=0和c为时,给出两个特解,y=I和,.--l。
由以上两个例子可以看出,确定常数(.的取值范围,需要注意两点,第一是在进行积分运算时,应注意到常数c的来源;第二是注意通解表达式对于常数c的取值限制。
综上所述,我们看到,学习过程中除了理解基本概念,掌握基本运算和基本解题方法之外,还应该特别注意解题中的细节,强调严谨的分析和认真的演算,只有
AdA"~In(COS.Y),运算结果表明,最
后得到的通解也是1’=-3+(。COS.Y((。∈R)。现在的问
题是,在这里写成ftan,vdv=一ln(cosx)行不行?正确
的答案是:如果要求的是微分方程的通解,并且没有指明是在哪个区间上求微分方程通解,这样写是可以的;如果题目已经指明在区间【0.:)求解,应写成
ftan.rd_v=一In(cos.v);若在区间(=,7r)上求解,应写
成ffanY小=一1n(~cos.Y);如果问题是求解微分方程
通解,并且没有指明是在哪个区间上求微分方程通解,
这样才能准确、无误。(下转53页)
就应该写成J
r、d.v=一lnlcos.YI。
由于一NO:具有拉电子共轭效应(一C),所以,I、Ⅲ
Ⅲ不稳定,Ⅱ稳定。故一NOz为第二类定位基。
一NO::具有一C和一I效应,故为致钝基团。
2.2关于取代基致活和致钝问题
第一类定位基(邻、列-位定位基)使苯环再发生亲电取代反应活性增加(致活,卤素除外),第二类定位基(问位定位基)使苯环再发生亲电取代反应活性降低(致钝)。原有取代基决定苯环再发生亲电取代反应活性高低的原因是:原有取代基的电子效应(含共轭效应和诱导效应)使苯环电子密度增加,则再发生亲电取代反应活性增加,该取代基为
致活基团:反之为致钝基团。
以一OH,一CI,一N02为例说明:
苯环定位规律是有机化学芳香烃一章的重点内容。对于它的理论解释是本章节的难点。笔者基于20多年有机化学教学经验,提出直接用中间体正离子稳定性解释定位规律,比用共振论法更简单更易馇。
【参考文献】
(1】袁履冰.有机化学【M】.北京:高等教育出版
社.1999.164一171.
一OH:具有给电子共轭效应(+c)和拉电子诱导效应(一I),由于+c>一I,故一OH为致活基团。
一cl:因为一cl与苯环共轭,其提供的P轨道为3P轨道,而Cl的3P轨道与苯环c的2P轨道大小差异较大,故共轭效果不好,所以cl给电子共轭效应(+C)较小,总的结果是:+C<一I,所以一cl为第一类定位基,但却是致钝基
【2】邢其毅,徐瑞秋,等.基础有机化学Dq.第2版.北京:高
等教育出版社.1993.272―280.
【3】荣国斌,苏克曼.大学有机化学基础瞰】.上海:华东理工大学出版社.2000.187―192.
【4】高鸿宾.有机化学瞰】.第3版.北京:高等教育出版
社.19919.259―267.
ANovelSimpleExplanationabouttheOrientationintheBenzeneRing
LIANG
Cheng-gang
(FacultyofChemistry,BaotouTeachersCollege;Baotou014030)
.Abstract:NownearlyalltheexplanationsappearedinorganicbooksRingisbased
theruleoforientationinthebenzene
thetheoryof
resonance.
Aftermanyyearsthinkandreseat’chtheauthorproposedbenzenering.Tireexplanationwastestified
simpleexplanationabouttheorientationinthe
begood
Afterfiveyearsapplicationinteaching.
Keywords:1)e,nzenering:l'u|e“orientation;theoretiealsolution
(上接17页)【参考文献】
【2】李苹荪,高等数学微积分习题详解【M】,武汉:武汉大学出版社,1994年。
【115E高雄,周之铭,朱思铭,等.常微分方程[MI.第2版.北
京;高等教育出版社,1982年.
【3】薛嘉庆,高等数学题库精编嗍,沈阳:东北大学
出版社,2001年。
AFewAnnotationabouttheSolutionofDifferentialEquation
CHENXiang-hua
(FacultyofMathematics,BaotouTeachersCollege,Baotou014030)
Abstract:Inthispaper,passthrough
takenplaceinprocessofsolution,and
fewtypicalexample,anumberofsomemistakenhasbeenpointoutitmethodofsolve
eolTeel
problemhasbeengiven.
Keywords:differentialequation;processofsolution;techniqueofsolve
problem.
关于微分方程解的几点注释
作者:作者单位:刊名:英文刊名:年,卷(期):被引用次数:
陈向华, CHEN Xiang-hua
包头师范学院,数学科学学院,内蒙古,包头,014030阴山学刊(自然科学版)YINSHAN ACADEMIC JOURNAL)0次
参考文献(3条)
1.王高雄.周之铭.朱思铭 常微分方程 19822.李苹荪 高等数学微积分习题详解 19943.薛嘉庆 高等数学题库精编 2001
相似文献(10条)
1.学位论文 蒋青松 计算机代数求解微分方程的方法研究及其机器实现 2006
本课题研究用计算机代数求解微分方程的理论与方法及其机器实现。研究了基于吴特征列方法的双曲正切函数展开法,给出了求解比较复杂的偏微分方程的通用求解过程;在Adomian分解原理的基础上,给出了利用Adomian分解原理求解一类奇异初值问题的算法及其Maple环境下的程序,解决了用手工过程求解无法实现的求解问题;研究了微分代数方程的幂级数展开法和Pade级数近似解,给出了相应的机器求解过程,并给出了Pade级数的精确解和近似数值解在同一坐标系下的图形,给出了二者的逼近比较;作出了将符号计算方法和数值计算方法结合起来求解微分方程的研究工作,这是求解比较复杂的偏微分方程的新途径。在求解过程中,运用了计算机代数的有关知识并且借助了计算机代数系统中的Maple和Matlab软件.求解过程相对手工过程而言,具有极大优势。
2.学位论文 王小兵 热机电耦合压电智能薄板结构不确定性及其H&,∞&振动控制研究 2007
本文以‘热机电耦合压电智能薄板结构’作为研究对象,对其不确定性和鲁棒H∞振动控制进行了研究,主要工作包含如下内容:
1、以每层材料各不相同的层叠板为研究对象,在热传导、对流以及热辐射边界条件共同作用下,采用4节点矩形单元,利用伽辽金法建立了该结构瞬态温度场的有限元微分方程,进而推导出了温度场对任意设计参数的灵敏度表达式,发现了灵敏度表达式中存在的公共项具有惯性效应。
2、对于线性瞬态温度场,当所有结构参数均具有随机性时,为了避开矩阵指数的随机性分析难题,通过方程解耦、特征对随机性分析以及一次二阶矩法,给出了瞬态随机温度场数字特征的求解过程。以三层复合板结构为例进行了求解,并将所得结果与Monte Carlo数值模拟法的结果进行比较,证明了求解过程的实用性。
3、对于非线性瞬态温度场,当所有结构参数均具有随机性时,为了节省随机温度场:Monte Carlo数值模拟的时间,在引入随机因子后,提出了一类关于热传导矩阵和热容矩阵的近似处理方法,并给出了误差的定量计算方法,即通过灵敏度和一阶泰勒展开线性化的方法去近似求解由于近似处理所造成的温度场数字特征的误差。算例表明大量的Monte Carlo数值模拟时间得以节省,且所得的近似误差具有一定的精度。
4、将3中获得的思路运用到了平面桁架结构的随机性分析中,对结构的质量矩阵和刚度矩阵的随机性提出了两种近似处理方法,通过对误差的定性分析表明第一种方法能获得结构动力特性均方根上限,第二种方法能获得结构动力特性随机性的近似解,这两种方法都能显著节俭Monte Carlo数值模拟法的计算量。算例表明所给方法不但节省了数值模拟的时间,且具有较好的精度。
5、对于3和4中提及的随机因子,这里首次从数学证明的角度对其进行了论述:将随机物理参数结构系统的质量矩阵和刚度矩阵均以随机因子的齐次形式表示,通过对QR变换法求解随机矩阵特征值过程的论述,证明了系统中各阶固有频率具有相同的随机因子,并导出了求解固有频率随机变量概率密度的公式;在此基础上,论述了系统正则振型矩阵中各元素亦具有相同的随机因子,并导出了求解正则振型随机因子概率密度的公式。算例表明了所得结论的正确性。
6、对于非线性瞬态温度场,为了避开对复杂分布进行假设检验时面临的统计量寻找和分位数确定等经典难题,通过综合运用差分法、Monte
Carlo数值模拟法、Q―Q图法、Box-Cox幂变换法,正态假设检验以及参数区间估计,提出了一整套近似求解温度响应主体区间的方法和表达式,并通过算例表明所提出的分析和求解策略的有效性与合理性。
7、针对热机电耦合压电智能薄板结构,提出了一种新型的包含4个位移节点、2个电势节点和8个温度节点的混合有限元模型,并基于虚功原理推倒了热机电耦合有限元微分方程,通过与已有文献的比较,表明所得有限元方程的正确性;同时也给出了这类有限元模型具体的编程实现方法。以智能悬臂薄板为例,对位移响应、温度响应和输出电压进行了数值仿真,凸现了在小挠度变形下,变化幅值较大的温度对输出电压具有重大影响。
8、对于7中所得有限元微分方程,当所有参数皆具有随机性时,通过各响应量对各个随机参数灵敏度的解析求解和一次二阶矩法,给出了依次求解结构固有频率、固有振型、温变、位移以及输出电压等随机响应数字特征的方法和步骤。并通过算例表明处理方法可行且有较好精度,但求解过程十分费时费力。
9、针对8中存在的求解困难问题,研究了综合应用差分法和基于灵敏度求解的矩法对结构位移和温度响应的数字特征进行了求解;算例表明改进后的方法运算快捷且具有一定精度。
10、为了实现对小区间参数不确定性智能薄板结构的鲁棒H∞振动控制,这里先研究了‘如何由原小区间参数不确定性系统的不确定性求解平衡降阶后系统的不确定性’这一难题。针对平衡降阶过程,提出了一种近似求解方法,即:通过正交变换获得实对称矩阵的特征对,在求解矩阵指数灵敏度和实对称矩阵特征对灵敏度的基础上,获得了平衡变换矩阵的灵敏度,在探讨了排序矩阵的选取后,进而利用一阶泰勒展开式和区间运算法则求得了降阶后的不确定性系统。在对原不确定性系统和降阶后的不确定性系统分别施加相同的输入后,通过二者输出区间的比较对所提方法的正确性进行了验证。
11、针对10中获得的不确定性系统,提出了一类维数较低的不确定性矩阵的分解形式;利用H∞控制理论对降阶后的增广被控对象进行了鲁棒控制器求解,并给出了该控制器对原不确定性系统的鲁棒性验证过程。算例表明lO中的研究成果起到了连接降阶与鲁棒控制的桥梁作用。
12、为了实现对小区间参数不确定性智能薄板结构的鲁棒H∞振动控制,在完成10和11针对一般系统的研究内容后,还需研究确定性智能薄板结构的振动控制问题。为此,这里利用平衡降阶法对其进行了降阶处理;论述了增广被控对象的构建以及利用Hoo控制理论对振动进行抑制的过程。并以悬臂薄板为例,其控制效果表明所给处理过程的合理性,也反映出该类结构在同时进行弯曲抑制和拉伸抑制时,压电致动层所具有的控制作用与干扰作用的双面性。
13、通过综合运用10~12中的研究成果,实现了对小区间参数不确定性智能薄板结构的鲁棒H∞振动控制。
3.期刊论文 展丙军 常微分方程求解过程的简化及特解的公式化 -通化师范学院学报)
通过对常微分方程的常规解法的进一步探讨,总结出使求解过程简化的具体作法,同时实现了几类微分方程求解的公式化或半公式化,提高了求解速度和准确性.
研究一类二阶实常系数非齐次微分方程y″+py′+q=(a0+a1x)eαxsinβx的解法,应用叠加原理和Euler公式,将其化为二阶线性非齐次方程,并利用对应的特征方程给出了这一类方程特解的一般公式,简化这一类微分方程的求解过程.
5.期刊论文 汤光宋.徐利姣.TANG Guang-song.XU Li-jiao 三类一阶微分方程的求解公式 -沙洋师范高等专科学校学报)
文献对三类一阶微分方程的求解,采用先找积分因子,再利用积分因子转化为全微分方程,然后按全微分方程的求解方法求解,其过程较繁复.本文借用变量替换法,化为变量可分离的方程,直接给出通解的积分形式,推广了方程的可积类型,并使求解过程大为简化.
6.学位论文 高灵芝 可压流体圆柱绕流计算模拟 2006
本文对可压流体圆柱绕流计算模拟进行了研究。文章主要内容为:对区域进行网格剖分、研究各个结点处可压流体的速度、压强和密度的变化情况、列出相应的微分方程、用差分法把微分方程离散为线性方程组、用Jacobi迭代法进行方程组的求解过程,最后用C程序来实现迭代的求解过程。主要工作可分以下几步进行:首先,对具体的问题进行分析,建立数学模型(N---S方程),对此方程进行化解、整理,用相应的微分方程组表示出来。其次,对二维求解域(圆域)进行网格剖分,在此网格图中用向后有限差分法离散微分方程,得到有关速度、压强和密度各个变量的方程组。再次,用Jacobi迭代法求解线性方程组,讨论Jacobi迭代的敛散性。最后,用C程序来实现Jacobi迭代的求解过程,通过运行程序得出相应结果,再把结果用图表表示出来。
7.会议论文 陈永进.官会林.李琴 用类比法求解极坐标下的平衡微分方程 2007
在数理学科中得到广泛的应用的类比法是一种从特殊到特殊和富有创造性的逻辑推理过程。本文应用类比法,根据直角坐标下求解平衡微分方程的思路,对极坐标下的平衡微分方程进行了求解。整个求解过程目的明确,思路清晰,简明严谨。对弹性力学平面问题的解的研究、教学与实践具有参考价值。
8.学位论文 王振 微分和差分方程的构造性变换与机械化求解 2008
本文以数学机械化思想和AC=BD模式为指导,以构造性的变换及符号计算为辅助工具,主要研究了数学物理中微分方程的构造性变换与机械化求解问题,包括微分、微分差分方程组的精确求解算法,渐进求解算法,线性变系数微分方程间的代数变换和初等函数等价问题的算法。
第一章介绍了数学物理机械化及其相关学科的发展,围绕微分方程的求解理论和计算机代数的关系,简述了相关方面的国内外研究和发展概况,在本章最后介绍了本文的选题和主要工作。
第二章考虑了微分方程的AC=BD理论,介绍了C-D对理论的基本内容和思想,AC=BD理论侧重于对微分方程变换的机械化构造,可以把复杂问题转化为简单问题,把变系数问题转化为常系数问题,把非线性问题转化为线性问题来解决。
第三章从一个新的角度来考虑非线性发展方程的精确求解理论,类似于代数方程的理论,微分方程也可以看作在特定环上的求解。通过构造微分多项式环及其相应的导子,在这个环上求解满足给定方程的元素。可以不断的扩大构造的微分多项式环及其导子,以扩大给定方程特解的形式。并且给出了多项式环上的多项式次数的确定方法。这种方法确定的次数不仅和给定的方程有关,还和未定元的导数(相当于辅助方程)有关。并且给出了计算非线性发展方程精确解的两个机械化算法,改进了射影Riccati方法,提出了多行波变换的辅助方程法,这个方法可以获得更多类型的精确解,包括多孤立子解,周期和孤子相互作用解,双周期解和孤子作用解,并将其推广到变系数的情形,以2+1维KdV方程为例,模拟了解的性态。给出了n+1维Klein-Gordon方程,n+1维Klein-Gordon-Zakharov方程的双行波形式的解。这种模式的优点主要有:第一,可以统一很多方法,并在此基础上提出新方法。第二,解的表达式不再是假设的,而是按照微分多项式环的基展开,减少了原来算法中的假设.第三,可以给出更有效的展开阶次的确定方法。
第四章首先将第三章中确定展开阶次方法推广到微分差分方程,通过对非线性微分差分方程的差分项和微分项分别处理,即构造差分方程来约化差分项,构造微分方程来约化微分项。通过将所得多项式化为代数方程组的求解。研究了两种辅助方程,第一种为Riccati方程,另一种为椭圆方程。这个方法比原来的方法提高了效率,因为最终所解方程化为代数方程,而原来的解法有时候会有非代数的方程。并且该方法更有统一性,离散的Riccati方程法将非线性微分差分方程的孤立波解,周期解,和有理解的求解统一起来了。只需计算一次即可得到三种类型的解。而椭圆方程法则将几种主要的椭圆函数的求解统一起来了。接下来通过构造方程的不同差分项的递推关系构造了微分差分方程的射影Riccati方程法,对其解分别求出其差分项的递推关系,构造了一个机械化算法,它可以得到广泛的精确解。
第五章讨论了几种求解微分方程的近似解的方法在微分差分方程中的应用,有同伦分析法,Adomian分解法,这种思想对Lyapunov小参数展开法,摄动方法等都同样适用,在第一节以离散的KdV方程为例说明了同伦分析法的应用,将结果与精确解相对比,它们吻合良好。在第二节我们以Hybrid方程为例说明了Adomian分解法在微分差分方程的应用,并且求解了它的冲击波解和扭结孤立波解。对连续情形的求解过程进行分析说明了同伦分析中的几个假设量之间可以相互确定的,这样就减少了假设量的个数,使求解过程更加有章可循。
第六章利用符号积分中的Order函数的性质,设计了一个算法确定两个三阶变系数常微分方程之间的所有自变量之间的有理代数变换(因变量的变换不做限制),若存在就给显式的求出来,如不能求解可以证明所给方程间不存在有理变换。特别的,我们可以给出变系数线性常微分方程通过变换变为常系数线性微分方程。
有时计算机输出的结果异常庞大和复杂,而且有些结果是可以化简的,计算机并没有自动将这些结果化简为我们想要的相对简单的形式,针对这个问题,我们设计了一个算法判断两个孤立子解是否是等价的,亦即是否是同一个解。另外,我们给出了无平方分解的另外一种算法,不同于已有的先求最低次数的无平方因子,反过来,尽可能的从具有较高的次数的无平方因子开始求.这样就可以一次消掉次数较高的因子,这种算法对次数较高的多项式效率更高。
9.期刊论文 熊灿.谢建新.XIONG Can.XIE Jian-xin 二阶常系数微分方程解法的简化 -南昌工程学院学报)
给出了二阶常系数齐次线性微分方程通解的三角函数形式或双曲函数形式,同时得出了利用位移定理,结合待定系数法解几类特殊的二阶常系数非齐次线性微分方程的方法,简化了此类微分方程的求解过程.
10.学位论文 张立达 非均匀变截面梁(柱)自由振动的差分法求解 2010
变截面梁在工程中得到广泛的应用。由于其控制方程为变系数微分方程,一般很难得到解析解,可以运用数值计算方法求解其近似解。本文基于Euler-Bernoulli梁的理论并采用中心差分法分析了置于弹性地基上的变截面梁在轴向荷载作用下的自由振动响应,获得了固有频率和振型,讨论了各种参数对频率的影响。主要内容包括以下几个方面:
采用中心差分格式将任意四阶变系数常微分方程两点边值词题离散为系数矩阵为五对角的线性代数方程组。然后,利用追赶法推导出了该代数方程组的递推求解过程。利用该递推求解过程编程时只需要5个一维数组即可存储代数方程组的系数矩阵。具有便于编程、存储空间小、计算速度快等优点。它可以用于求解任意四阶变系数常微分方程的两点边值问题。
2.基于Euler-Bernoulli梁理论并利用上述中心差分法的递推求解过程,研究了弹性地基上非均匀变截面梁在轴向压力作用下的自由振动响应。在分析中考虑了梁的横截面尺寸和弹性模量沿轴线连续变化的情形。作为算例,分别计算了两端夹紧和两端简支边界条件下弹性地基上的圆形非均匀变截面梁在轴向压力作用下的自由振动响应。在退化为无弹性支撑和无轴向压力作用均匀梁的特殊情形,所得无量纲固有频率和振型与解析解相比十分接近。讨论了截面变化参数、材料弹性模量参数、弹性地基参数、轴向载荷参数以及边界条件对梁的自由振动固有频率的影响。结果表明,固有频率随轴向压力的增大而减小,随地基参数的增大而增大。
三亿文库包含各类专业文献、各类资格考试、文学作品欣赏、行业资料、幼儿教育、小学教育、专业论文、中学教育、外语学习资料、关于微分方程解的几点注释_图文13等内容。 
 偏微分方程的MATLAB解法_数学_自然科学_专业资料。偏微分方程的解法 引 言 偏微分方程定解问题有着广泛的应用背景。人们用偏微分方 程来描述、解释或者预见各种...  各类微分方程的解法大全_数学_自然科学_专业资料。南京林业大学 各类微分方程的解法 1.可分离变量的微分方程解法 一般形式:g(y)dy=f(x)dx 直接解得∫g(y)dy...  关于常系数非齐次线性微分方程通解的注释 [J]. 大学数学,-145. [9]张云艳.常系数非齐次线性微分方程的几个解法[J].黄山学院学报,-9...  t . 3 ?1 就此,我们已经用微分方程求解的理论解释了差分方程求解的特征根方 法,并说明了两者之间的关系. 第二部分、非齐次常系数线性差分方程(组)暨欧拉-...  微分方程的通解和特解如何理解_数学_自然科学_专业资料。微分方程的通解和特解如何理解,右端等于0,说明它是齐次方程;右端不等于0,说明它是非齐次方程。...  常微分方程解法_理学_高等教育_教育专区。常见的有解析解的常微分方程 1、可分离变量方程: f1 ( x) g1 ( y ) dx + f 2 ( x) g 2 ( y ) dy = ...  第7章 微分方程课本基础几点_理学_高等教育_教育专区。自己考研总结的 1、 2、 微分方程的解是 0 阶的函数,而不是具体的值。 若微分方程的解函数中含有的...  微分方程的符号解与数值解_数学_自然科学_专业资料。实验报告实院(系) 班级 专业 实验 名称 所用 软件 实验 目的 课程名称:数学模型 验 报学号 姓名 告(六)...  ? ? 对于求解微分方程 u? ? u ? t 与三阶 Taylor 级数法的计算格式的形式完全相同。 14.对 Heun 二阶方法(2.10)式作出如图 2.3 那样的几何解释。 15...}

我要回帖

更多关于 常微分方程边值问题 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信