《求线性差分方程的通项之数学实验.doc》由会员分享，可在线阅读，更多相关《求线性差分方程的通项之数学实验.doc（4页珍藏版）》请在金锄头文库上搜索。
1、钦州学院数学与计算机科学学院数 学 实 验 报 告实验完成日期 2010 年 11 月 8 日 ， 第 11 周 ， 星期 一 成绩等级（五级分制） 评阅教师 评阅日期 年 月 日数学实验报告填写要求：思路清晰，中间结果和最终结果真实；字迹工整，报告完整。实验题目及内容 实验题目：借助计算机，求解下列线性差分方程（即求出数列的通项公式）。 内容：为了能直观了解数列的特性，结合运用Matlab软件计算出数列的通项。问题描述（用自己组织的相关数学语言重述现实问题；注意对约定的条件作说明）数列满足递推关系，称这样的递推关系为二阶线性差分方程。根据观察，可以猜测具有指数形式，可设，得到，从而有，运用M
2、atlab求方程得、的解。 和都是差分方程的解，都是数列的通项，但数列不会有两个通项。猜测与的线性组合仍是差分方程的解。设，代入差分方程进行检验，猜测是否成立。若成立，求出差分方程的解。 根据初始条件，,可能确定常数、的值。由此得出，数列的通项公式。模型建立或思路分析（建立合理，可解释的数学模型，通过公式、表格或图形直观明确地描述模型的结构；无法通过建立模型解决的，给出解题的思路及办法。）设，二阶线性差分方程转换成 ，所以有，利用Matlab解得差分方程的解 ，的值，求得的解也是数列的通项，但数列不会有两个通项。由于与满足线性关系，所以可设，代入差分方程进行检验。根据初始条件、及、的值可算出常
3、数,由此可求得数列的通项。实验结果（通过数学表达式、 列表或图形图像的方式显示实验结果。） 、满足线性组合，所以因为初始条件，从而有因此，数列的通项公式为：结果分析及结论（对实验结果进行定量分析、合理性分析或误差分析； 对所讨论的问题重新认识或提出相关类似问题的拓延； 给出自己的意见和合理建议。）（1）解题过程中，线性方程的基本知识要了解，还要大胆假设与猜测具有的指数形式，如果不理解题目，则想不出解题思路。结合解题步骤中的思路，以及借助计算机的Matlab软件，求出数列的通项公式为：（2）类似的问题可以了解Fabonacci数列的特性，利用求通项的方法，可以计算出Fabonacci数列的通项为
4、。求解方法或解题步骤（针对所建模型或解题思路，给出具体的求解方法或解题步骤。对通过编程解决的问题，画出流程图，给出细节部分的算法，给出相关软件的代码；其他方法解决的，给出详细的解题步骤。）（1）数列满足递推关系，称这样的递推关系为二阶线性差分方程。根据观察，可以猜测具有指数形式。不妨设为进行尝试。将代入差分方程：得到 消去因子有 运用Matlab软件：solve(x2-2*x-2,x)ans = 1+3(1/2) 1-3(1/2)所以，解得 ，由此可知这两个都是差分方程的解。 （2）和都是差分方程的解，都是数列的通项，但这是不怎么可能，因为数列不会有两个通项。猜测与的线性组合仍是差分方程的解。设，代入差分方程进行检验，猜测确实成立！因此，差分方程的解为： （3）根据初始条件，,有 可以确定常数,的值。运用Matlab软件输入：syms x yx,y=solve(1+3(1/2)*x+(1-3(1/2)*y-3,(1+3(1/2)2)*x+(1-3(1/2)2)*y-8)得出 x =1/3*3(1/2)+1/2y =-1/3*3(1/2)+1/2因此，求解得 ， 所以，数列的通项公式为：1
