网站维护一般都是维护什么,外包公司哪家好,东莞营销型网站建设找火速,济南高端网站设计策划一、实验目的1#xff0e; 了解一般非线性方程的求根是比较复杂的事情#xff1a;要讨论(或知道)它有无实根#xff0c;有多少实根#xff1b;知道求近似根常用的几种方法#xff0c;每种方法的特点是什么。2#xff0e; 用通过二分法(区间半分法)、不动点(也Picard)迭代…一、实验目的1 了解一般非线性方程的求根是比较复杂的事情要讨论(或知道)它有无实根有多少实根知道求近似根常用的几种方法每种方法的特点是什么。2 用通过二分法(区间半分法)、不动点(也Picard)迭代法及Newton迭代(切线)法求其它非线性方程的根并尽可能估计误差。二、实验原理三、实验程序四、实验内容1. 用二分法求方程x3-3x-10在的所有根.要求每个根的误差小于0.001.提示与要求: (1) 利用精度找到迭代次数;(2) 由f(x)3(x2-1)可取隔根区间[-2,-1].[-1,1].[1,2]);(3) 用程序求各隔根区间内的根.2. 用不动点迭代求: (1)x32x210x-200的所有根.或: (2)9x2-sinx-10在[0,1]上的一个根.3. 用Newton迭代法求解下列之一,准确到10-5:(1) x3-x-10的所有根;(2) ex2-x2cosx-60位于[0,2]上的根.五、实验程序• 二分法文件代码function y f(x)yx^3-3*x-1;end程序代码function Bipart(a0,b0,tol)%a0为左区间b0为右区间tol为区间误差限aa0;bb0;mceil(log((b-a)/tol)/log(2));for k1:mp(ab)/2;if f(p)*f(b)0ap;elsebp;endenddisp([经过二分法求得的跟为x,num2str((ab)/2,%.6f)])disp([共经过,num2str(k),次计算])命令窗口Bipart(-2,-1, 0.001)Bipart(-1, 1, 0.001)Bipart (1,2, 0.001)运行结果• 不动点法文件代为function y f(x)y9*x^2-sin(x)-1;end程序代码function Budongdian(x0,tol,m)%x0为初始值tol为误差容限m为最大迭代次数syms xF(x)sqrt(sin(x)1)/3;mm;text;x[];x(1)x0;for k1:mx(k1)F(x(k));if (abs(x(k1)-x(k))tol)1text迭代成功;disp(text);disp([经过不动点迭代法求得的跟为x,num2str(x(k1),%.7f)]);disp([共经过,num2str(k),次计算])breakendendif isempty(text)1disp(Method failed)end命令窗口Budongdian(1, 0.00001, 100)运行结果• 牛顿迭代法文件代码function y f(x)yexp(x)2^(-x)2*cos(x)-6;end程序代码function Newton_gen(x0,tol,m)%x0为初始值tol为误差容限m为最大迭代次数syms xF(x)x-f(x)/diff(f(x));mm;text;x[];x(1)x0;for k1:mx(k1)F(x(k));if (abs(x(k1)-x(k))tol)1text迭代成功;disp(text);disp([经过Newton迭代法求得的跟为x,num2str(x(k1),%.7f)]);disp([共经过,num2str(k),次计算])breakendendif isempty(text)1disp(Method failed)end命令窗口Newton_ .gen(2, 0.00001, 100)运行结果另解• 二分法定义函数function Bipart_2(a0,b0,tol,Tol)%a0为左区间b0为右区间tol为区间误差限Tol为f误差限aa0;bb0;mceil(log((b-a)/tol)/log(2));for k1:mp(ab)/2;if f(p)*f(b)0ap;elsebp;endif abs(f((ab)/2))break;endenddisp([经过二分法求得的跟为x,num2str((ab)/2,%.6f)])disp([共经过,num2str(k),次计算])命令窗口Bipart. .2(-1, 1,0.001,0.1)Bipart_ .2(-1, 1, 0.001, 0.000001)运行结果:• 不动点法定义函数function [x,k]budong(fun,x0,tol,m)for k1:mxfun(x0);if abs(x-x0)break;endx0x;endxvpa(x,8);function tfun(x1)syms x;f9*x.*x-sin(x)-1;ssubs(diff(f,x),x,x1);xx1;f9*x.*x-sin(x)-1;tx-f/s;命令行窗口输入[x,k]budong(fun,0.5,1e-5,100)运行结果
