数值分析实验报告三.docVIP

数学与信息工程学院 实 验 报 告 课程名称： 计算方法 实 验 室： 7404 实验台号： 班 级： 姓 名： 实验日期： 2014 年 5　月 21 日 实验名称 非线性方程求解及方程组直接解法 实验目的 和 要 求 1.二分法的Matlab实现； 2.牛顿法的Matlab实现； 3.牛顿下山法、割线法、艾特金加速法、重根迭代法、非线性方程 组牛顿法中任选其一。 4.列主元消元法的Matlab实现； 5.二选一：（1）平方根法的Matlab实现；（2）追赶法的Matlab实现。 实验内容和步骤： 实验内容： 1、Compare the number of computation for finding the root of with accuracy . (1)、Use the bisection method starting with the interval [0,1]. (2)、Use the iteration method ,the initial value . 2、The equation has two roots near 0.1. Determine them by means of Newton’s method. (with accuracy ) 3、用迭代法求方程附近的一个根。方程写成下 列等价形式，并建立相应的迭代公式。 ，迭代公式 4、用列主元消法解线性方程组Ax=b。 5、用平方根法解线性方程组Ax=b。 （1）（2） 实验步骤： 实验编程 运行结果 二分法 实验编程 1. 新建M文件newtonqx.m输入： function [k,x,wuca,yx]=erfen(a,b,abtol) a(1)=a; b(1)=b; ya=fun(a(1)); yb=fun(b(1)); %程序中调用的fun.m 为函数 if ya* yb0, disp(注意：ya*yb0,请重新调整区间端点a和b.), return end max1=-1+ceil((log(b-a)- log(abtol))/ log(2)); % ceil(u)是大于u的最小取整数 for k=1: max1+1 a;ya=fun(a); b;yb=fun(b); x=(a+b)/2; yx=fun(x); wuca=abs(b-a)/2; k=k-1; [k,a,b,x,wuca,ya,yb,yx] if yx==0 a=x; b=x; elseif yb*yx0 b=x;yb=yx; else a=x; ya=yx; end if b-a abtol , return, end end k=max1; x; wuca; yx=fun(x); 2. 建立名为fun.m的M文件 function fun=fun(x) fun=exp(x)+10*x-2; 3.在MATLAB工作窗口输入： x=-4:0.1:4; y=exp(x) +10*x-20; plot(x,y) grid [k,x,wuca,yx]=erfen (﹣1,1,10^-5) 运行结果 ans = 0 -1.0000 1.0000 0 1.0000 -11.6321 10.7183 -1.0000 ans = 1.0000 0 1.0000 0.5000 0.5000 -1.0000 10.7183 4.6487 ans = 2.0000 0 0.5000 0.2500 0.2500 -1.0000 4.6487 1.7840 ans = 3.0000 0 0.2500 0.1250 0.1250 -1.0000 1.7840 0.3831 ans = 4.0000 0 0.1250 0.0625 0.0625 -1.0000 0.3831 -0.3105 ans = 5.0000 0.0625 0.1250 0.0938