（毕业论文毕业设计）-MATLAB论文_MATLAB在求解差分方程问题中的应用.docVIP

PAGE PAGE 2 《MATLAB语言》课程论文 Matlab在求解差分方程问题中 的应用 姓名： 学号：12010245312 专业：电子信息工程 班级：2010级 指导老师： 学院：物理电气信息学院 完成日期：2011年12月12日 MATLAB在求解差分方程问题中的应用 (姓名： 学号12010245312 2010级) [摘要]将MATLAB应用于差分方程的计算中，充分发挥了,MATLAB强大而方便的功能，通过对公积金问题、有限差分法、线性常系数差分方程组等的试算表面，MATLAB在解决实际的工程和数学问题中，与其它计算机程序设计语言C、FORTRAN相比，具有使用更为简便、语句功能更强，用户界面良好的特点，适合于在工程计算，尤其是物探数据处理领域推广应用。 [关键字]MATLAB语言 差分方程 公积金 线性常系数 一、问题的提出 差分方程是含有未知函数及其导数的方程，满足该方程的函数称为差分方程的解。差分方程是微分方程的离散化。一个微分方程不一定可以解出精确的解，把它变成差分方程，就可以求出近似的解。Matlab是一个为科学和工程计算而专门设计的高级交互式的软件包。它的数值分析、矩阵运算、信号处理和图形显示于一体，构成了一个方便的、界面友好的用户环境。MATLAB在发达国家已广泛应用于科研院所、工程技术等各个部门, 并成为大学生、研究生必备的工具软件，应用MATLAB强大的数学功能通过实例的建模、编程和仿真过程,来简便、直观、高效分析数学问题。 二、差分方程求解公积金问题 本例用Matlab数学工具简化计算的方法用差分方程建模的求解数学问题，来培养和提高运用数学计算的方法解决现实生活中实际问题的能力；实验内容如下： 元旦开始,贷款40万,还款期为15年,采用等额本息发还款(即每月还的本金利息和相等),求解如下问题: 1、15年利率不变，则每月还多少恰好还完？15年总利息多少？ 2、若此人想提前到十年还款完毕，则按1中的方式还一年后需一次性还款多少才 能使他再十年内还清？ 模型建立及理论求解： 上网查阅资料得到公积金贷款的年利率在贷款期限为6~30年时都为3.87%. 由此可得月利率为: 3.87%/12=0.3225% 还款期限总共15年，即： n=15x12=180个月 1、 15 年利率不变，则每月还多少恰好还完？15 年总利息多少？ 设还款第k个月后欠款Ak为元，月还款m元，则由Ak变化到Ak+1除了还款还有利息。设月利率r，得到： 为求m值，设： 易得： 由此可推出Bk的表达式： 由上述推导可得： 由此可得差分方程的解为： 2、若此人想提前到十年还款完毕，则按1中的方式还一年后需一次性还款多少才能使他再十年内还清？ 求解分析可得第k个月后欠款总额为： 由此可得，一年后欠款总额为： 欲十年还完，则必须一次还清在1中接下来五年需要还得欠款。还款六年后即k=6*12=72个月后的总欠款为： 即在一次还清一定数目的欠款之后剩下的九年需要还款总额为: 即需要一次还清的款数为： 分析：用传统计算方法解决时我们需要列出如上传统方程，这样的计算不仅繁琐费时，而且没有图示很难给以直观的感受，现在我们用MATLAB语言来对此例题做以下解析： MATLAB程序如下： r=0.003225; k=180; a0=400000; m=a0*(1+r)^k*r/[(1+r)^k-1]; %差分系数 p=m*180-400000; %总利息 k=12; A12=a0*(1+r)^k-m/r*[(1+r)^k-1]; %12个月后欠款 k=72; A72=a0*(1+r)^k-m/r*[(1+r)^k-1]; %1中方式六年后欠款 Q=A12-A72; %一次性还款数 m p A12 A72 Q 运行结果如下： m = 2.9328e+003 p = 1.2790e+005 A12 = 3.7993e+005 A72 = 2.6710e+005 Q = 1.1283e+005 即： m=2932.76元 总利息为： 需要一次性还款数为： 1、若15年利率不变则每月换2932.76元恰好还完； 支付的总利息为127896.97元； 2、需要在正常还款一年后一次还清112831.84元，才能实现他十年内还清贷款的愿望。 三、有限差分的求解方法 设有一个截面为矩形的无限延伸导体槽，槽的宽为16m,高为10m，槽的盖板（导体）同两侧面绝缘，盖板电位是100V，周边电位为0，如图3所示，试分析槽内的电位分布。 N N N N 图1 方形截面的导体槽 分析：这是一个二维稳定电场的问题，在求解区内（无场源），电位函数满足拉普拉斯方程： 可以应用差分格式： 然后具体给出求