实验3_求解线性方程组迭代法与插值法.pdfVIP

数值分析实验报告三 求解线性方程组的迭代方法和插值法（ 2 学时） 班级专业 统计学 1 班 姓名 何斌 学号 06 日期 2011-4-25 一 实验目的 1．掌握求解线性方程组的简单迭代法； 2. 掌握求解线性方程组的赛德尔迭代法。 3. 掌握不等距节点下的牛顿插值公式以及拉格朗日插值公式。 二 实验内容 1．使用简单迭代法求解方程组（精度要求为 10 6 ）： 20 x 1 2 x 2 3 x3 24 x 8x x 12 1 2 3 2 x1 3x 2 15 x3 30 2 ．使用赛德尔迭代法求解上述方程组（精度要求为 10 6 ）： 3 ．已知函数表： x 1.45 1.36 1.14 y 3.14 4.15 5.65 用拉格朗日插值公式计算 x 1 . 4 以及 y 5 .01 所对应的近似值。 4. 已知函数表： X 93.0 96.2 100.00 104.2 108.7 Y 11.38 12.80 14.70 17.07 19.91 用牛顿插值公式求 x (102 ) 的近似值。 三 实验步骤（算法）与结果 1. 简单迭代法求解方程组（精度要求为 10 6 ）： 1 % 设原方程组的系数矩阵为 A A=[20 2 3;1 8 1;2 -3 15]; b=[24;12;30]; A1=[A,b]; n=length(b); for i=1:n if A1(i,i)0 A1(i,:)=-A1(i,:); end if A1(i,i)==0 error( 系数矩阵的主对角元素不能为零 ) end a(i,:)=A1(i,:)/A1(i,i);%a 是 3 ×4 的矩阵 end if n==3 delta=1e-6;% 精度要求 it=0;% 构造计数器 max=8;% 最大循环次数 x1=0;x2=0;x3=0; % 迭代初值 x_1=0.1;x_2=0.1;x_3=0.1; while (abs(x_1-x1)delta | abs(x_2-x2)delta | abs(x_3-x3)delta) itmax x1=x_1; x2=x_2; x3=x_3; % 迭代公式 x_1=-a(1,2)*x2+(-a(1,3)*x3)+a(1,4); x_2=-a(2,1)*x1+(-a(2,3)*x3)+a(2,4); x_3=-a(3,1)*x1+(-a(3,2)*x2)+a(3,4); it=it+1; end end x=[x_1,x_2,x_3]; for k=1:3