(1.3)--2023数值实验——练习2.pptVIP

1x0=0;Tol=1/2*10^(-8);N=10^4;k=1;x=g(x0);whileabs(x-x0)TolkNk=k+1;x0=x;x=g(x0);endx,k2、迭代法算法:不动点迭代给定初始近似值x0，求x=g(x)的解.输入:初始近似值x0;容许误差TOL;最大迭代次数N输出:近似解x或失败信息.Step1k=1;Step2x=g(x0);Step3Whileabs(x-x0)TOLk?NStep4k=k+1;Step5x0=x;Step6x=g(x0);Step7Outputthesolutionofequation:xanditerativetimes第二章：非线性方程的数值解法

2functiony=g(x)y=(x+1).^(1/3);迭代函数的定义functiony=newton(x)y=x-f(x)./(3*x.^2-1);Newton迭代函数functiony=f(x)y=x.^3-x-1;由函数定义的方程为：x0=0;Tol=1/2*10^(-8);N=10^4;k=1;x=newton(x0);whileabs(x-x0)TolkNk=k+1;x0=x;x=newton(x0);endx,k%Newton法x0=0;x1=1;Tol=1/2*10^(-8);N=10^4;k=1;x=x0-f(x0)*(x1-x0)/(f(x1)-f(x0));whileabs(x-x0)TolkNk=k+1;x0=x;x=x0-f(x0)*(x1-x0)/(f(x1)-f(x0));endx,k%割线法

第三章：线性方程组的数值求解一、三角形方程组的解法?考虑下三角形方程组

将计算出来，并存入算法3.1.1下三角形方程组Ly=b的前代法：forend

?考虑上三角形方程组将计算出来，并存入

算法3.1.2上三角形方程组Ux=y的回代法：forend两种算法的工作量(加减乘除运算次数之和)均为

?三角分解法的基本思想（MATLAB命令lu(A)）:记方程组可化为下面两个易求解的三角方程组设已知方程组系数矩阵的三角分解为其中,为下三角矩阵,为上三角矩阵.

forforfor三角分解（Doolittle分解）算法工作量约为：,L和U都存于Aendendend

9程序1：编写用一般迭代法（不动点迭代）求解的程序；请分别编写以下程序题目1：求方程在中的根。程序2：编写用牛顿迭代法求解的程序；程序3：编写用割线法求解的程序。

10程序4：编写回代法的求解题目2的程序；程序5：编写前代法求解题目3的程序；题目2：求解线性方程组程序7：编写三角分解法求解题目4的程序。程序6：编写高斯消去法求解题目4的程序；题目4：求解线性方程组题目3：求解线性方程组