實(shí)驗(yàn)3_求解線性方程組迭代法與插值法

1、數(shù)值分析實(shí)驗(yàn)報(bào)告三求解線性方程組的迭代方法和插值法（2學(xué)時(shí)）班級(jí)專業(yè) 信科3 姓名 梁嘉城 學(xué)號(hào)201130760314日期一 實(shí)驗(yàn)?zāi)康?掌握求解線性方程組的簡(jiǎn)單迭代法；2.掌握求解線性方程組的賽德爾迭代法。3. 掌握不等距節(jié)點(diǎn)下的牛頓插值公式以及拉格朗日插值公式。二 實(shí)驗(yàn)內(nèi)容1使用簡(jiǎn)單迭代法求解方程組（精度要求為）：2使用賽德爾迭代法求解上述方程組（精度要求為）：3已知函數(shù)表：x1.451.361.14y3.144.155.65用拉格朗日插值公式計(jì)算所對(duì)應(yīng)的近似值。4. 已知函數(shù)表：X93.096.2100.00104.2108.7Y11.3812.8014.7017.0719.91用牛頓插

2、值公式求的近似值。三 實(shí)驗(yàn)步驟（算法）與結(jié)果1#include<stdio.h>main()float a34=20,2,3,24, 1,8,1,12, 2,-3,15,30;for(int i=0;i<=2;i+)for(int j=0;j<=2;j+)aij=(-1)*aij;a00=20;a11=8;a22=15;float x=0,y=0,z=0;float X,Y,Z;for(int q=0;q<=1000;q+)X=(y*a01+z*a02+a03)/a00;Y=(x*a10+z*a12+a13)/a11;Z=(x*a20+y*a21+a23)/a22

3、;x=X;y=Y;z=Z;printf("方程組的解是X=%9.6f,Y=%9.6f,Z=%9.6fn",X,Y,Z);2#include<stdio.h>main()float a34=20,2,3,24, 1,8,1,12, 2,-3,15,30;for(int i=0;i<=2;i+)for(int j=0;j<=2;j+)aij=(-1)*aij;a00=20;a11=8;a22=15;float x=0,y=0,z=0;for(int q=0;q<=1000;q+)x=(y*a01+z*a02+a03)/a00;y=(x*a10+z*

4、a12+a13)/a11;z=(x*a20+y*a21+a23)/a22;printf("方程組的解是X=%9.6f,Y=%9.6f,Z=%9.6fn",x,y,z);3.#include<stdio.h>main()float x3=1.14,1.36,1.45;float y3=5.65,4.15,3.14;float Y;Y=(1.4-x2)*y1/(x1-x2)+(1.4-x1)*y2/(x2-x1);float X;X=(5.01-y1)*x0/(y0-y1)+(5.01-y0)*x1/(y1-y0);printf("由拉格朗日插值公式得當(dāng)X

5、=1.4時(shí),Y=%f,當(dāng)Y=5.01時(shí),X=%fn",Y,X);4.#include<stdio.h>main()float x5=93.0,96.2,100.00,104.2,108.7;float y5=11.38,12.80,14.70,17.07,19.91;float dy1,dy2,dy3,dy4;float ddy1,ddy2,ddy3;float dddy1,dddy2;float ddddy;dy1=(y0-y1)/(x0-x1);dy2=(y1-y2)/(x1-x2);dy3=(y2-y3)/(x2-x3);dy4=(y3-y4)/(x3-x4);ddy1=(dy1-dy2)/(x0-x2);ddy2=(dy2-dy3)/(x1-x3);ddy3=(dy3-dy4)/(x2-x4);dddy1=(ddy1-ddy2)/(x0-x3);dddy2=(ddy2-ddy3)/(x1-x4);ddddy=(dddy1-dddy2)/(x0-x4);float Y;Y=y3+(102-x3)*dy3+(102-x3)*(102-x2)*ddy2+(1002-x3)*(102-x2)*(102-x1)*dddy1;printf("由牛頓