数值计算大作业——刘概念.docVIP

课 程 设 计 课程名称： 数值分析 设计题目： 数值计算大作业 学 号： S315070064 姓 名： 刘峰 完成时间： 2015年10月25日  题目一、非线性方程求根 1.题目 假设人口随时间和当时人口数目成比例连续增长，在此假设下人口在短期内的增长建立数学模型。 （1）如果令表示在时刻的人口数目，表示固定的人口出生率，则人口数目满足微分方程，此方程的解为； （2）如果允许移民移入且速率为恒定的，则微分方程变成， 此方程的解为； 假设某地区初始有1000000人，在第一年有435000人移入，又假设在第一年年底该地区人口数量1564000人，试通过下面的方程确定人口出生率，精确到；且通过这个数值来预测第二年年末的人口数，假设移民速度保持不变。 2.数学原理 采用牛顿迭代法，牛顿迭代法的数学原理是，对于方程，如果是线性函数，则它的求根是很容易的，牛顿迭代法实质上是一种线性化方法，其基本思想是将非线性方程逐步归结为某种线性方程来求解。 设已知方程有近似根(假定)，将函数在点进行泰勒展开，有 于是方程可近似地表示为 这是个线性方程，记其根为，则的计算公式为 fplot(1000*exp(x)+(435*x)*(exp(x)-1)-1564,[0 3]);grid 大概估计出初始值x=0.5 function [p1,err,k,y]=newton(f,df,p0,delta,max1) % f是非线性系数 % df是f的微商 % p0是初始值 % dalta是给定允许误差 % max1是迭代的最大次数 % p1是牛顿法求得的方程近似解 % err是p0误差估计 % k是迭代次数 p0,feval(f,p0) for k=1:max1 p1=p0-feval(f,p0)/feval(df,p0); err=abs(p1-p0); p0=p1; p1,err,k,y=feval(f,p1) if(errdelta)|(y==0), break,end p1,err,k,y=feval(f,p1) end function y=f(x) y=1000000*exp(x)+435000*(exp(x)-1)/x-1564000; function y=df(x) y=1000000*exp(x)+435000*(exp(x)/x-(exp(x)-1)/x^2); 4.结果分析与讨论 在MATLAB中的command window输入 newton(f,df,1.2,10^(-4),10) 运行后得出结果 p0 =0.5000 p1 =0.1679 err =0.3321 k =1 y =9.2415e+004 p1 =0.1031 err =0.0648 k =2 y =2.7701e+003 p1 =0.1010 err =0.0021 k =3 y =2.6953 p1 =0.1010 err =2.0129e-006 k =4 y = 2.5576e-006 ans =0.1010 运算后的结果为，通过这个数值来预测第二年年末的人口数， t=2时候对于 实践表明，当初始值难以确定时，迭代法就不一定收敛了，因此要根据问题实际背景或者二分法先得一个较好的初始值，然后再进行迭代；再者迭代函数选择不合适的话，采用不动点迭代法也有可能出现不收敛的情况；因此我采用的是牛顿法。 题目二：线性方程组求解 1.题目 假设一个物体可以位于个等距点的任意位置，当物体在位置时，它只能等可能的移动到或者，而不能直接移动到其他任何位置，概率表示物体从位置开始在到达右端点之前到达左端点的概率，显然，且有 既有下面方程组： 取对方程组进行求解（迭代法或者直接法）。 2.数学原理 在解微分方程的边值问题、热传导方程以及船体数学放样中建立的三次样条函数等工程技术问题时，经常遇到下面形式的线性方程组： = 方程简记，该线性方程称为三对角线方程组，其系数矩阵A满足条件 所以为弱对角阵可以采用追赶法进行计算，利用三对角矩阵的LU分解建立计算量更少的线性方程组求解公式。将系数矩阵A进行克劳特分解，即A分解为下三角矩阵和单位上三角矩阵的乘积； A== 其中，，为待定系数，直接利用矩阵乘法公式可得 ，， ，， ， 于是推得计算，，的公式 ，； ，，； ，； 由此计算出L和U中的全部元素，完成了系数矩阵A的克劳特分解。求解线性方程组等价于求解和。 因而得到解三对角线性方程组的追赶法公式 计算的递推公式： 解 解 我们将计算系数和称为追的过程，将计算方程组的解称为赶的过程。整个过程为追赶法的思想。 3.程序设计 function