数学实验之2.非线性方程(组)求解(1).pptVIP

计算结果 线性连接图 震荡发散情形 收敛情形 a=0.5;x1=[]; x1(1)=0.5; for i=2:20 x1(i)=a*x1(i-1)*(1-x1(i-1)); end n=1:20; subplot(2,2,1),plot(n,x1),title(a=0.5,x0=0.5) MATLAB程序 Died4.m 2）蛛网图 在直角坐标系中，首先画出直线y= x和曲线y= f(x)，其中f(x)为迭代曲线。迭代过程如下图。 ① 从直线y= x上点An(xn,xn)到曲线y= f(x)上点Bn(xn,xn+1)；（垂直） ② 从Bn(xn,xn+1)到An(xn+1,xn+1); (水平） ③ 重复。 a=2.9, x0=0.2 An Bn An y = x y = f(x) died8.m * * 重庆大学信息与计算科学系 报告人: 何仁斌 非线性方程(组)求解(1) 引言及引例 非线性方程的常用求解方法 使用MATLAB的命令进行方程求解 范例:波音公司最优定价策略问题 实 验目的与内容 求高次代数方程的根是一个基本的、古老的数学问题。 高阶方程能否求解？ 定义：一些由实际问题列出的方程中常常包含三角函数、指数函数等，它与n(≥2)代数方程一起统称为非线性方程(组)，记作 引言 例：exsin(x)+5x3=0 一元二次方程 高阶多项式方程 计算机技术和数学软件的飞速发展，为一些古老数学问题的解决提供了极大的方便。如利用MATLAB的图形功能就能帮助我们判断方程有没有根，确定根的近似位置，例如： -0.5 1.5 2.5 0.25 2.5 对于一般的非线性方程，怎样得到比较精确的近似解，正是该试验要研究或解决的问题。 下面给出一个实际问题，它们都归结为求非线性方程的近似解。 放射性废物的处理问题 问题背景 一段时间, 美国原子能委员会是按以下方式处理浓缩放射性废物的. 他们将废物装入密封性能很好的圆桶中, 然后扔到水深300英尺的海里. 这种做法是否会造成放射性污染, 很自然地引起了生态学家及社会各界的关注. 原子能委员会一再保证, 圆桶非常坚固, 决不会破漏, 这种做法是绝对安全的. 然而一些工程师们却对此表示怀疑, 他们认为圆桶在海底相撞时有可能发生破裂. 由此双方展开了一场笔墨官司. 究竟谁的意见正确呢? 只能让事实说话了! 引例 问题的关键在于圆桶到底能承受多大速度的碰撞? 圆桶和海底碰撞时的速度有多大? 工程师们进行了大量破坏性的实验, 发现圆桶在直线速度为40 ft/s 的冲撞下会发生破裂, 剩下的问题就是计算圆桶沉入300 ft 深的海底时, 其末速度究竟有多大? 问题分析 1. 使用55加仑的圆桶; ( 1加仑 = 3.7854升 ) 2. 装满放射性废物时的圆桶重量为 W = 527.436磅 (1 磅 = 0.4526公斤 ) 3. 在海水中圆桶受到的浮力 B = 470.327磅 4. 圆桶下沉时受到海水的阻力 D = C v C 为常数, 经测算得: C = 0.08. 5. 建立坐标系, 取垂直向下为坐标方向 y , 海平面为坐标原点. y 0 假设: 根据牛顿第二定律, 圆桶下沉时应满足微分方程: 为了求出圆桶与海底的碰撞速度v(t), 需要求出圆桶下沉到海底300英尺时的时间 t, 再计算v(t)，要做到这一点是十分困难的. 若将速度v 看成是海水深度y 的函数, 即 由复合函数的求导法知 如何从上述方程中解出函数 v = v(y)？只能借助数值方法求出v(300)=45.1ft/s，显然大于40ft/s。 结论：工程师们的怀疑有道理。 数学建模的一般步骤 模型准备 模型假设 模型构成 模型检验 模型分析 模型求解 模型应用 引言及引例 非线性方程的常用求解方法——迭代 使用MATLAB的命令进行方程求解 范例:波音公司最优定价策略问题 实 验目的与内容 求解方法 图形放大法 简单迭代算法 数值求解法 非线性方程 三种形式的图形描述 f(x) = 0 软件求解 方程 f(x)=0 1）建立坐标系，画曲线f(x)； 2）观察曲线f(x)与x轴相交的交点； 3）将其中一个交点进行局部放大； 4）该交点的横坐标值就是方程的根。 1．图形放大法 例: 求