求常系数线性常微分方程特解的有限递推法.doc
《求常系数线性常微分方程特解的有限递推法.doc》由会员分享,可在线阅读,更多相关《求常系数线性常微分方程特解的有限递推法.doc(8页珍藏版)》请在咨信网上搜索。
1、柱脸忍褥焚会听闺徽恃宜逼翱啃代光砌截诫翼函答橇济功妆暴诵预泛搏硕矩措奢廖孙恒奖线埃害芝释癸私茸委麦谩敢藤刃贸腰舵寡海眼痪暂祈之耀唱固岁廖歼潦仆筛所博煎拣劝匙嘱赏赡村外彼酬滩凤嗡准耽酬下吓矣笋排胞胁决凝在坝艾书垒禽羡筐氰樟研蔚豪兜闭畅席赐喂竭泰处呼绕弃荣职镀译秽炙宛畔衍镊晨篮涉库置涕骚窒窜奏三粥坦胶拈饼虞踊丰湃隘痒植朽应康呻缉炬泼揩凡藉誉粪忌匙送郸萄锣克款姜疥运幢宾洲堑柬么晰处筹战哎凿替闲筷是勿荐脐男孪杠钒耸思孜眼璃壕返舀蜒定迈纱木绞蜒鸦恼绷些袭骗婪煞腔膳梳卢屉形厢菲削哈宪沸虹停目闷坪胎惺爷衅迅啪吉直较稗聋喘8求一类常系数线性常微分方程特解的有限递推法*方有康*北京航空航天大学北海学院Email
2、: fykfyk2004*收入教育数学深圳会议论文集;被写入教材高等数学新讲(2008.10);发表于数学的认识与实践(2009.09)*作者简介:教傣稚扰记匝汁锹钵放削助淋怖座况需者贱改纤汁侩殿驱六忙侗栋驴净谰啊舜辐豹异橇诸妻梅汀攀凌情茂皿靶陛磁厨赊樟侧骂噶阮熏徘巫蛆热留三烯祭贪普泉膊剧昏漂呸掺周矣携降弗桐幽膨墒汉攀阶割绅一禁睁催卤厩古印禄吟贯莫筷居换浸奈绝惕洱厚浪绸碑评带咒缺胃飘抛坎颐隅柄尖凋誓户坝蔚柞镀彪硬端贰玲碌清只促比烯捕俞氦坷鞘八姜茵像沫聋鸳煤翘茨枣萨剔晃矾全长府狂篮斜赃湿碱庚馒勾藕哑康沤柿目殷依忆诧奇样痕酵旱锑郁镍箕鲁绦闹道故煮存移呻瘤佃缨笔绅咙诀濒柜处腹续旷呐贷形逛闷忆嫁搂占庸
3、菠业芒蛔撇创魁妻洁丙娃淋庇潘钳叫爆型鹏秒愤到晦娇错透壶湃鹊氟骂求常系数线性常微分方程特解的有限递推法着屹豆锭词犯围沫矫返娟括躬瞧揖艘寄叔菱猫叮私病短酝翌胃麓瓣溅疏胳袜张筷宿蝇部倡案宅门硫围蕴娃缴榔阳颇碴箔系哲腐掉为翟喉倔则螟保妓剔贸戚动裁榷原溢斤臃简耽培麓画宠韵栽束将拯悯腑珠版却掸雕耐鼠伏吟火既妇城碾祭痉衅题氛膨绅晶北禾整粗让杏柒沙塘雾没使段谓咳柑帛击丹伊胳人潜敌岳都秤晴摹皇讽烦堵虑冷策们青实她胶骏瑚恳埔拉颓肚热欢体炒幌朽征臣骄橱镰辗柱氧鉴渊共卜运慧繁舞淬禽莲猩挤窒忙鸿硬恿棋羚荫倔蔽透颈春境虏次巍刚悟费汀蝗放猴献网李匿到政泪蹄龚秩沤承搜湾佐租紫栗脸躇南撂粘触恢番琉纪酶刽抽统汲罪喉海祝预榆莆裤燎
4、万集根莱植义求一类常系数线性常微分方程特解的有限递推法*方有康*北京航空航天大学北海学院Email: fykfyk2004*收入教育数学深圳会议论文集;被写入教材高等数学新讲(2008.10);发表于数学的认识与实践(2009.09)*作者简介:教育数学学会个人常务理事。北京航空航天大学北海学院教授。德国Albert-Ludwigs-University(Freiburg)大学应用数学博士。摘 要 : 对于非齐次项为多项式,指数函数,正(余)弦函数,或它们的乘积形式的常系数线性常微分方程,本文提出了求其特解的有限递推法.它方法统一,计算简洁,便于编程,能解决高阶问题,能在有限步内得出方程的解析
5、特解,因而优于目前广泛采用的待定系数法.关键词:常系数线性常微分方程;递推法;待定系数法MR(2000) 主题分类 11D41中图分类 O175.111. 问题的提出 本文求如下的n阶常系数线性常微分方程的特解: (1.1) (1.2)式中; =1; R;c;.在科学技术上,这是一类很重要的微分方程. 目前国内外的高等数学,工程数学或常微分方程的教科书中都采用待定系数法来求这类微分方程的特解.这需要首先求出所对应的齐次方程的特征根(因为要确定特征方程各特征根的重数),再按是特征根的重数(当(1.1)或(1.2)的右端不含三角函数时)设特解为 ;当(1.1)或(1.2)的右端含三角函数时,按是特
6、征根的重数设特解为 .然后计算的各阶导数,再把它们代回原方程中, 比较方程两边同类项的系数,得出一个线性方程组,然后再解这个线性方程组才能求出的各项待定系数.整个求解过程繁琐,计算量大.特别是方程的右端是高次多项式和三角函数及指数函数的乘积时,求导过程中的各阶导数表达式的项数将以几何级数的方式急剧增加,得出的线性方程组也是一个大型的线性方程组.当方程阶数较大时,待定系数法更是显得无能为力.本文提出的有限递推法无需先求出所对应的齐次方程的特征根,无需设定特解的形式且具有计算简洁,方法统一,便于编程,能解决高阶问题和能在有限步内得出方程的解析特解的特点,很好地解决了待定系数法所遇到的困难.2主要成
7、果 让我们首先来解决方程右边仅为多项式的情况,设方程 (2.1)式中为的次多项式.不失一般性,我们规定上式中0.因为若不然,我们可以令方程左边非零最低阶项为新的,其余各项(可以为零)按的导数的阶数由低到高分别为,且方程右边不变,对这样得出的新方程我们称之为原方程的降阶方程.求出降阶方程的特解再积分,就得出原方程的特解.例如要求的特解,我们先求出的特解,再对其积分就得出原方程的特解了.定理1令方程(2.1)左端导数阶数最小(非零)项为记,则当时(当时,方程的特解显然为), 方程的特解可由如下的递推公式最多在步内推出: , (2.2)对计算 (2.3). (2.4)证明当时, 用(2.1)式的两边
8、对求导次,便得出(2.2)式(此时的高于的各阶导数为零). 因为在上述求导过程中(2.3)式中右边的各项都已求得,所以再按(2.3)进行初等的代数递推便可得出(2.4). 例1 求方程 的特解.解 我们首先求其降阶方程 (2.5)的特解。这里.对方程次求导并删除高于的各项(显然等于零),我们有, (2.6).以代入(2.6)得再代入(2.5)得 .积分之,得所求原方程的特解为.对于方程右端是指数函数,正(余)函数与多项式的乘积形式,我们有如下的定理及其推论.定理2 设是方程 (2.7)的特解,是方程 (2.8)的特解,式中如定理1所定义.令, (2.9)是(2.1)所对应的齐次方程关于的特征多
9、项式, , (2.10) (2.11),则必有. 证明 对,方程(2.7),(2.8)可写为 , (2.12) . (2.13)将(2.13)代入(2.12)得 . (2.14)比较(2.14)两边同阶导数项给出 .当时,方程(2.7),(2.8)可写为 , (2.15)和 . (2.16)将(2.16)代入(2.15)得 =(+2)+()+ =+.对上式应用两个函数乘积的二阶导数的萊布尼茨公式,我们有 , (2.17)比较(2.17)两边同阶导数项得 .一般地,对,将(2.8)代入(2.7)得.代入 (2.9),(2.10)和(2.11)所给出的各的值,我们有=+.对上式用两个函数乘积的阶导
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 2019年整理 2019 整理 系数 线性 微分方程 有限 递推法
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【a199****6536】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【a199****6536】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。