数学实验课件--方程求解.pptVIP

fsolve()函数的第三个输入是options,它是一个结构型数据,可以通过函数optimset ()进行设定.当不进行设定时采用缺省设置. fsolve函数还可以有后面的参数设定,这一功能在有的时候非常有用. fsolve()语句的用法 例如,求解函数sin(ax)-x=0的最小正解. 分析:有已知的知识我们知道,该方程没有解析形式的解,也就是说很难得到这个解和a之间的具体关系.数值的方法可以对a先取确定的值,这样变成一个一元方程,容易进行求解,当我们变化a时,就得到很多这样的解,通过这样的方法,我们可以得到解和a之间的一些大致关系. 1) 掌握方程求解的三种解法：解析 法、数值解法以及图形表示解的方法; 2) 学会使用MATLAB软件求解析解、数值解 和图形解； 3） 通过范例学习怎样建立方程模型和 分析问题的思想； 实验目的 一，建立方程 例子【问题背景】一段时间, 美国原子能委员会是 按以下方式处理浓缩放射性废物的. 他们将废物装入 密封性能很好的圆桶中, 然后扔到水深300英尺的海 里. 这种做法是否会造成放射性污染, 很自然地引起 了生态学家及社会各界的关注. 原子能委员会一再保 证, 圆桶非常坚固, 决不会破漏, 这种做法是绝对安 全的. 然而一些工程师们却对此表示怀疑, 他们认为 圆桶在海底相撞时有可能发生破裂. 由此双方展开了 一场笔墨官司. 究竟谁的意见正确呢? 只能让事实 说话了! 数学建模的一般步骤 模型准备 模型假设 模型构成 模型检验 模型分析 模型求解 模型应用 二，方程求解 1，解析方法 2，图形放大法 3，迭代方法 4，区间方法 方程求解之解析方法 主要针对一些比较简单的方程以及方程组，比如多项式方程等。同学们以前对方程的求解也是针对这样一些方程进行的。该方法的优点是可以利用纸笔得到简单有效并且精确的解；缺点是可以求解的方程数量太少。Matlab和Maple提供了求方程解析解的函数，可以说对数学演算提供了不少方便。 方程求解之图形放大法 图形最大的有点就是直观，试想如果我们有了函数准确的图形，那么曲线和x轴的交点就是我们要求的方程的解。因此我们可以利用图形工具得到方程的解。当然，计算机上的图形不可能等同于函数的真实图形，因为计算机上的图形是曲线上部分点的轨迹而不是全部，因此通过图形不可能得到方程的精确解，甚至它只是一个比较粗略的解，当然，通过对图形的放大可以得到更精确一些的解。同时，这种方法也不适应大量的数据处理。 方程 f(x)=0 1）建立坐标系，画曲线f(x)； 2）观察曲线f(x)与x轴相交的交点； 3）将其中一个交点进行局部放大； 4）该交点的横坐标值就是方程的根。 2．图形放大法 例: 求方程 x5 +2x2 + 4 = 0 的一个根. 画方程曲线图（tuxfd.m） 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]，再观察！ 该方程有几个根？欲寻找其中一个实根，并且达到一定的精度。 2．图形放大法 逐次缩小区间，观察一个根在-1.55~-1.5之间。 2．图形放大法 方程求解之迭代法 迭代法的理论以及方法的出现，对方程求解有着里程碑式的意义。其基本思想如下: 需要求解的方程： f (x) = 0 （1） 经过某种变形得：x = j (x) （2） 从而求解方程（1）转化成为求解（2）得不动点。（满足条件x*=j(x*)的点x*称为不动点） 为得到方程的不动点，可以构造迭代过程如下： xn+1 =j (xn)，n =0，1，… x0 定义为迭代初值。 解： 第一步 构造迭代函数： x=j (x) 例：用迭代方法求解方程 x3 -x2 -x-1 = 0。 第二步 迭代 设定初值 x0=1， xn+1 =j (xn)，n =0，1，… 用 MATLAB 编程（died2.m文件) x=1;y=1;z=1;(初始点） for k=1:20 x=x^3-x^2-1; % j1 (x) y=(y^2+y+1)^(1/3); % j2 (y) z=