数值计算实验报告气温变化问题

《数值计算实验报告气温变化问题》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、为了适应公司新战略的发展，保障停车场安保新项目的正常、顺利开展，特制定安保从业人员的业务技能及个人素质的培训计划数值计算实验报告气温变化问题　　数值计算实践实验报告　　学号　　姓名指导教师　　成绩　　XX－XX学年第1学期　　一.方程求根　　?一?实验目的?利用Newton法,割线法,抛物线法分别对三次方程、三角函数方程、四次方程和五次方程进行求根，并对这三种求根方法进行收敛性大小的比较。????????二?问题描述???依据学号的后三位乘以倒数第二位加1再除以1000为方程的根，以学号的后三位数分别作为方程的三次项,二次项,一次项的系数,根据所给的根

2、以及三个系数确定常数项，得到三次方程；以sinx+4x2+x+a0=0的形式构造三角函数方程；将三次多项式再乘以(x-p*)2得到对应的五次多项式(p*为已经确定的方程的根,显然,得到的五次方程有重根)；用三角函数方程同样乘以(x-p*)得到一个新的方程；用Newton法,割线法,抛物线法对上述方程进行求根。????算法介绍??目的-通过该培训员工可对保安行业有初步了解，并感受到安保行业的发展的巨大潜力，可提升其的专业水平，并确保其在这个行业的安全感。为了适应公司新战略的发展，保障停车场安保新项目的正常、顺利开展，特制定安保从业人员的业务技能及个人素质

3、的培训计划　　Newton法：将非线性方程逐步化为线性方程来求解，设fx=0有近似根xk，将f(x)在xkf??(?)*处作一阶Taylor展开,即f(x*)?f(xk)?f?(xk)(x*?xk)?(x?xk)2，x在xk和x之间，　　2!　　f(x)　　)(xx*?将(x*?xk)2看成高阶小量,则有:f(xk)?f?(xk*?xk)?0,则xk?k。　　f(xk)　　f(x)?f(x)f(xk)(xk?xk?1)　　割线法：切线斜率=割线斜率?f?(xk)?kk?1?xk?1?xk?　　f(xk)?f(xk?1)xk?xk?1　　????????

4、抛物线法：?如果已知f(x)=0的三个近似根，以这三个点为插值节点构造二次插值多项　　式P2(x),然后以P2(x)的一个零点作为新的近似根。这样的迭代过程称为抛物线方法。设已知的三点为xk?2,xk?1,xk则插值多项式为　　p2(x)有两个零点xk?1?xk?　　其中，ω=fxk,xk?1+fxk,xk?2kk?1程序?Newton法：　　functionx=newton(fname,dfname,x0,e)e=1e-6;x=x0;x0=x+2*e;whileabs(x0-x)>ex0=x;　　x=x0-feval(fname,x0)/feval(

5、dfname,x0);end　　?目的-通过该培训员工可对保安行业有初步了解，并感受到安保行业的发展的巨大潜力，可提升其的专业水平，并确保其在这个行业的安全感。为了适应公司新战略的发展，保障停车场安保新项目的正常、顺利开展，特制定安保从业人员的业务技能及个人素质的培训计划　　割线法：　　functionx=secant(fname,x0,x1,e)e=1e-6;x=x1;x1=x0;whileabs(x1-x)>ex0=x1;x1=x;　　x=x1-(feval(fname,x1)*(x1-x0))/(feval(fname,x1)-feval(fna

6、me,x0));disp(x);end　　抛物线法：?　　functionx=paowx(fname,x0,x1,x2,e)e=1e-6;　　x=x2;x2=x1;x1=x0;whileabs(x2-x)>ex0=x1;x1=x2;x2=x;h1=x1-x0;h2=x2-x1;　　q1=(feval(fname,x1)-feval(fname,x0))/h1;q2=(feval(fname,x2)-feval(fname,x1))/h2;w=q2+(q2-q1)*(x2-x1)/(x2-x0);ifw>0,　　x=x2-2*feval(fname,x2

7、)/(w+sqrt(w^2-4*feval(fname,x2)*(q2-q1)/(x2-x0)));elsew数值计算实验报告气温变化问题)???三角函数方程为：????????+5???，四次方程为：???????????+5??2???+，五次方程为：??5???4+??3+??2+???。Newton法：目的-通过该培训员工可对保安行业有初步了解，并感受到安保行业的发展的巨大潜力，可提升其的专业水平，并确保其在这个行业的安全感。为了适应公司新战略的发展，保障停车场安保新项目的正常、顺利开展，特制定安保从业人员的业务技能及个人素质的培训计划　　fu

8、n=inline('x^3+5*');>>dfun=inline('3*x^2+5');>>n