数值分析综述报告.pptx

淮 阴 工 学 院;;则称该近似值具有位有效数字（significant digit）。 此时，该近似值的相对误差为;k k;5;otherwise x(:,1)=normrnd(0,9,5,1); k=1 while k=50*5 x(:,k+1)=M2*x(:,k)+f2; if norm(x(:,k+1)-x(:,k))=R k=k+1; else X=x(:,k+1) disp(Causs-Seidel迭代法迭代次数为：) IterN=k break end end end 第四章 非线性代数方程（组）的数值解法： 一、二分法： 首先要确定适当的包含根的区间，这可以依据闭区间上连续函数的介值定理来确定，例 如该方程： f ? x? ? x3 ? x2 ? 0.8x ? 0.75 ? 0;;8;? ?;10;11;8285950923799759686915051372346805*t^2 + 8477820147617239232740036583612328*t + 512657632901869980094703083834211328)/22484983899204823688364170343605760;; ;;n;17;;;20;Runge—Kutta 法是按选取区间上函数变化率数目的多少和截断误差的阶数来区分的一系列 方法。 例题：给定如下常微分方程 Cauchy 问题： ;z(k) = YY(:,k) + h*K3(k); K4(k) = F(tt(k)+h, z(k)); YY(k+1) = YY(k) + 6\h*(K1(k) + 2*K2(k) + 2*K3(k) + K4(k)); end hold on plot(tt, YY, LineWidth, 3) title(An Example of Applying The Fourth-Order Rong - Kutta Method) legend(Euler Curve, The Curve Obtained By Order 4 R-K,The Curve Obtained By ode45) xlabel(Time) ylabel(Orbit) %输出三种方法最后的数值解 C=cell(2,3); C={Euler法,4阶R-K法,ode45法;Y(end),YY(end),z(end)} 运行结果： 最终得到的三种方法的数值解分别为： C =