方程求解省公开课金奖全国赛课一等奖微课获奖PPT课件.pptxVIP

数学建模与数学试验数学实验二方程求解(1)云南大学数学系任课教师：朱娟萍1/22定义：一些由实际问题列出方程中经常包含三角函数、指数函数等，它与n(≥2)代数方程一起统称为非线性方程(组)，记作引言求高次代数方程根是一个基本、古老数学问题例：exsin(x)+5x3=02/220.252.5-0.51.52.5引言计算机技术和数学软件飞速发展，为一些古老数学问题处理提供了极大方便。如利用MATLAB图形功效就能帮助我们判断方程有没有根，确定根近似位置，比如：3/22对于普通非线性方程，怎样得到比较准确近似解?引言引例:放射性废物处理问题【问题背景】一段时间,美国原子能委员会是按以下方式处理浓缩放射性废物.他们将废物装入密封性能很好圆桶中,然后扔到水深300英尺海里.这种做法是否会造成放射性污染,很自然地引发了生态学家及社会各界关注.原子能委员会一再确保,圆桶非常坚固,决不会破漏,这种做法是绝对安全.然而一些工程师们却对此表示怀疑,他们认为圆桶在海底相撞时有可能发生破裂.由此双方展开了一场笔墨官司.终究谁意见正确呢?只能让事实说话了!4/22引例:放射性废物处理问题问题分析问题关键在于圆桶到底能承受多大速度碰撞?圆桶和海底碰撞时速度有多大?工程师们进行了大量破坏性试验,发觉圆桶在直线速度为40ft/s冲撞下会发生破裂,剩下问题就是计算圆桶沉入300ft深海底时,其末速度终究有多大?5/220y问题假设1.使用55加仑圆桶;(1加仑=3.7854升)2.装满放射性废物时圆桶重量为W=527.436磅(1磅=0.4526千克)3.在海水中圆桶受到浮力B=470.327磅4.圆桶下沉时受到海水阻力D=CvC为常数,经测算得:C=0.08.5.建立坐标系,取垂直向下为坐标方向y,海平面为坐标原点.6/22依据牛顿第二定律,圆桶下沉时应满足微分方程建立模型7/22建立模型为了求出圆桶与海底碰撞速度v(t),需要求出圆桶下沉到海底300英尺时时间t,再计算v(t)，要做到这一点是十分困难.若将速度v看成是海水深度y函数,即由复合函数求导法知8/22建立模型借助数值方法求出v(300)=45.1ft/s，显然大于40ft/s结论：放射性废物不能随意放入公海！9/22模型准备模型假设模型组成模型检验模型分析模型求解模型应用数学建模普通步骤方程惯用求解方法图形放大法简单迭代法加速迭代法10/221．图形放大法方程f(x)=01）建立坐标系，画曲线f(x)2）观察曲线f(x)与x轴相交交点3）将其中一个交点进行局部放大4）该交点横坐标值就是方程根例:求方程x5+2x2+4=0一个根该方程有几个根？欲寻找其中一个实根，而且到达一定精度11/221．图形放大法画方程曲线图x=-6:0.01:6;y=x.^5+2.*x.^2+4;y1=x;plot(x,y,x,y1)或ezplot(‘f(x)’,[a,b])由此判断：方程一个根在区间[-2，2]内，所以将区间[-6，6]缩小至[-2，2]，再观察！逐次缩小区间，观察一个根在-1.55-1.5之间12/2213/22Z2.简单迭代算法引例：3x?ex=01）该方程有多少个根？怎样判断？2）怎样进行迭代求解？方程变形：x=ex/301/31/30.46520.46520.530814/222.简单迭代算法方程3x?ex=0迭代求解表：x=ex/315/22y=xy=ex/30.6162.简单迭代算法迭代过程如图所表示16/22Z2.简单迭代算法方程：f(x)=0经过简单变形：x=j(x)或j(x)=f(x)+xx被称为不动点迭代过程以下：xn+1=j(xn)，n=0，1，…x0定义为迭代初值1.表示式x=j(x)是否唯一？2.迭代产生序列是否一定会收敛？3.迭代与初始值x0是否相关？17/222.简单迭代算法例：用迭代方法求解方程x3-x2-x-1=0解：第一步结构迭代函数：x=j(x) 用MATLAB编程x=1;y=1;z=1;(初始点）fork=1:20x=x^3-x^2-1;%j1(x)y=(y^2+y+1)^(1/3);%j2(y)z=1+1/z+1/z^2;%j3(z)endX，y，z第二步迭代设定初值x0=1，xn+1=j(xn)，n=0，，…18/22jjjj(x)(x)(x)(x)序号序