基于Matlab的数值实验设计和实现___课程设计.ppt

实验目的和任务：理解非线性方程与方程组数值求解思想，设计并掌握MATLAB中非线性方程及方程组的数值解法，理解迭代法的收敛性与收敛速度. 要求精度为 ，给定方程与方程组分别为： ； （1-1） ， . （1-2） 1.1 线性迭代和斯蒂芬森加速迭代 对于方程（1-1），分别用线性迭代和斯蒂芬森加速迭代法计算如下： 线性迭代法(不动点迭代法)： function[xstar,index,it]=iterate(phi,x,ep ,it_max) if nargin4 it_max=100;end if nargin3 ep=1e-6;end index=0;k=1; while k=it_max x1=x;x=feval(phi,x); if abs(x-x1)=ep index=1;break; end k=k+1; end xstar=x;it=k; 函数phi： function f=phi(x) f=(x^2+2-exp(x))/3; [xstar,index,it]=iterate(@phi,0.5,1e- 5) 结果：xstar = 0.2575 index = 1 it = 9  （2）斯蒂芬森加速迭代法 function[xstar,index,it]=steffensen(phi, x,ep,it_max) if nargin4 it_max=100;end if nargin3 ep=1e-6;end index=0;k=1; while k=it_max x1=x;y=feval(phi,x); z=feval(phi,y), x=x-(y-x)^2/(z-2*y+x) if abs(x-x1)=ep index=1;break; end k=k+1; end xstar=x;it=k; 函数phi： function f=phi(x) f=(x^2+2-exp(x))/3; [xstar,index,it]=steffensen(@phi,0.5,1e-6) 结果：xstar = 0.2575 index = 1 it = 3  1.2 牛顿法及数值比较  function[xstar,index,it]=Newton(fun2,x ,ep,it_max) if nargin4 it_max=100;end if nargin3 ep=1e-6;end index=0;k=1; while k=it_max x1=x;f=feval(fun2,x); if abs(f(2))ep break; end x=x-f(1)/f(2); if abs(x-x1)=ep index=1;break; end k=k+1; end xstar=x;it=k; 函数fun2： function f=fun2(x) f=[x-(x^2+2-exp(x))/3,1+x]; [xstar,index,it]=Newton(@fun2,0.5,1e-6) 结果：xstar = 0.2575 index = 1 it = 4 1.3 MATLAB数值计算 function y=fun3(x) y=x^2-3*x+2-exp(x); [x,fval, iterations]=fzero(fun3,1) 结果：x = 0.2575 iterations=5 function F=fun4(x) F=[3*x(1)^2-x(2)^2;3*x(1)*x(2)^2-x(1)^3-1]; [x,fval,iterations]=fsolve(@fun4,[1;1]) 结果：x =[0.5000,0.8660]’ iterations = 1 2.1 LU分解和列主元高斯消去法 ②列主元高斯消去法 输入：[x,det,index,j]=Gauss_col(A1,b1) 得到方程的解为:x 1= 1592.6 x2=-631.9 x3=-493.6 得到系数矩阵行列式的值为det =-0.0305 ③小结 从代数上看,直接分解法与高斯消去