太原理工大学数值计算方法实验报告

1、TAIYUANUNBTRSIIYOFTECHNOLOGY本科实验报告课程名称：计算机数值方法实验项目：方程求根、线性方程组的直接解法、线性方程组的迭代解法、代数插值和最小二乘拟合多项式实验地点：行勉楼专业班级：*学号：*学生姓名：:指导教师：李志，崔冬华2016年4月8日学生姓名实验成绩实验名称实验一方程求根实验内容和要求熟悉使用二分法、迭代法、牛顿法、割线法等方法对给定的方程进行根的求解。选择上述方法中的两种方法求方程：f(x)=x3+4x2-10=0在1,2内的一个实根，且要求满足精度|X*-Xn|<X10-5(1)了解非线性方程求根的常见方法，如二分法、牛顿法、割线法。(2)加深对

2、方程求根方法的认识，掌握算法。(3)会进行误差分析，并能对/、同方法进行比较。实验原理1 .二分法：如果要求已知函数f(x)=0的根(x的解)，那先要找出一个区间a,b,使得f(a)与f(b)异号。根据介值定理，这个区间内一定包含着方程式的根。求该区间的中点m=(a+b)/2,并找出f(m)的值。右f(m)与f(a)正负号相同，则取m,b为新的区间,否则取a,mo重复第3步和第4步，直到得到理想的精确度为止。2 .割线法是利用牛顿迭代法的思想，在根的某个领域内，函数后直至二阶的连续导数，并且不等丁0,则在领域内选取初值xo,x1,迭代均收敛。(1) 在区间m川内输入初值x0,x1.(2) 计算

3、x2。x2=x1f(x1)*(x1x0)/(f(x1)f(x0)(3) x0=x1,x1=x2(4)判断是否达到精度，若是输出x1,若否执行(2)主要仪器设备HP计算机实验记录1.二分法pp:定义控制台应用程序的入口点。gjCNDOWEeyramewe方程组,为！f(x)-x3Tx2-10-0情输入计算区间以空格僵开)121.36523为方程的潜1接任意键继续.续法pp:定义控制台应用程序的入口点。SBC：VVINDOWSsy32cnd.e*e|方程组为f(,)-x"H4/2T0=0请输入两个初始近似值；(以室格隔升)121.26316L33883L366621.36521L3652

4、31.36523为方程的解请按仕意键继续.一123x114012X28241X313高斯分解法：将原方程组化为三角形方阵的方程组：lik=aik/akkaij=aij-lik*akjk=1,2,，n-1i=k+1,k+2,，nj=k+1,k+2,，n+1由回代过程求得原方程组的解：Xn=ann+1/annxk=(akn+1-Eakjxj)/akk(k=n-1,n-2,，2,1)2.LU分解法：将系数矩阵A转化为A=L*U,L为单位卜二角矩阵，U为普通上三角矩阵，然后通过解方程组l*y=b,u*x=y,来求解x.主要仪器设备HP计算机实验记录1.高斯消TC法：#include"&quo

5、t;#include""#include<>doublea56,a056;doublel5,tmp;voidExchange(inti)intj,l,k;doublemax=a0ii,temp;j=i；for(k=i;k<=3;k+)if(a0ki>max)max=a0ki;j=k;for(l=i;l<=4;l+)temp=a0il;a0il=a0皿l;a0jl=temp;for(i=1;i<=3;i+)for(j=1;j<=4;j+)aij=a0ij;voiddisplayA()inti,j;printf("n"

6、;);for(j=1;j<=3;j+)for(i=1;i<=4;i+)printf("%lf",aji);printf("n");voidmain()inti,j,k;for(i=1;i<=3;i+)for(j=1;j<=4;j+)scanf("%lf",&aij);a0ij=aij；displayA();printf("列主元素消元法如下");12314S128241132.0Q0000a.0080SB14.00B0006.0600091.0S09032.0300908.B0S60

7、02.也时的口4.0000001.80000013.H00000列主儿素消元法如F2.腼断电04.0OB909i.900000i3-0H00e1.0000002.0300008.6000D02.aoaaea3.araoaa14.000000112H11-0,0060002.«608004.0000901.03009813.0000000.0600601.0Q00B82.0300808.000B001,6060002.0080003.80000014,000000113H11-0.500网2.阴阴。阴融04.0000001.03000013.00B0B0S.0S006Qi.eanaaa

8、2.0000008.Q00B000.S600S0B.88BO0B2.5000007.5000002.«600004.905的。i.uaaaao13.0030006.0600001.009002.0000008.080800S.OS00S0B.8SB00B2.5309007.500000113H21=0.00Q00S2.必咐附MU04.05。酗1.S0OBEH13.000000阴.阴咐阴阴301.009032.0990008,0000006.6006600.(3000082.50110007.500000.000000xta=2.0000Hx3=3.000000Pressdnykeyt

9、ocontinue-1.0300002.0060903.0000001.6399002.030000-5.0090Q0输出L的矩阵为±Mwm0.0000001.0060302.000009B.KHMa1.000000所求方程组解因、=：|.口回回0回乩=2.080800,x3=3.00B000PFes£a”keytocontinue10x1X22x37.2Xi10x22x38.3X1X25X34.2pleaseinput,>2,x2：000nxirnx2(nlxSLnl00.O0G0SQ0.0000089.00030016.200600.8308000.8499092

10、0.971U8Q1.0700B81.1&08903d.8570001.1571001.24S2084l.eB53501.1B5340i.2gZ8Z5i.esEBfe1.1950991.294L3B61.6963391.1989371.298B3?71.0994421.1994421,2?933581.69981111.2?91.0799361.1959361.2?24101.09V9?91.1999791.2?99?5111.69999311121.09997811”99B1.2?131.8999991.1999971.2”9H141.1WH801.1.300900151.100600

11、1.260000i.ssaoao161.1006001.2OB0OB1.3»H03B171.1000601.2000001.300900习使用拉格朗日插值法或牛顿插值法求解方法。2.了解最小二乘法的多项式拟合的具体计算方法并且注意克服正规方程组的病态。给定数据点(Xi,yi)如下：Xi0yi1(i)使用拉格朗日插值法或牛顿插值法,求f的近似值.(2)用最小二乘法拟合数据的(n次)多项式，求f的近似值.(3)对比、分析上两结果实验原理设函数在区间a,b上n+1互异节点xo,xi，,xn上的函数值分别为yo,yi,yn,求n次插值多项式B(x)，满足条件Pn(Xj)=yj,j=0,1，,

12、n令Ln(x)=yolo(x)+yili(x)+ynln(x)=Eyili(x)其中(x'lKx),，ln(x)为以xo,xi,xn为节点的n次插值基函数，则Ln(x)是一次数不超过n的多项式，且满足Ln(xj)=yj,L=0,i，,n再由插值多项式的唯一性，得Pn(x)三Ln(x)主要仪器设备HP计算机实验记录(写出实验内容中的程序代码和运行结果)(可分栏或加页)拉格朗日插值法：#include""intmain()doublem=,a=,l=0;inti,j;doublex6=,;doubley6=,;for(i=0;i<=5;i+)for(j=0;j&l

13、t;=5;j+)if(i=j)continue;m=m*(a-xj)/(xi-xj);l+=yi*m;m=1;printf("结果为%lf",l);return0;bressanykbytocointinue最小二乘法：#include""#include""intmain()doublex7=0,y7=1,a0,a1,sum1=0,sum2=0,sum3=0,sum4=0,sum5=0,l,r;intm=6,i,k;for(i=0;i<7;i+)sum1+=xi;sum2+=xi*xi;sum3+=yi;sum4+=xi*yi

14、;sum5+=yi*yi;l=sum1/(m+1);a1=(sum4-l*sum3)/(sum2-l*sum1);a0=(sum3-sum1*a1)/(m+1);doubles=sum3*a0+sum4*a1;r=sum5-s;printf("y=a0+a1*x'n");printf("a0=%fa1=%ftn",a0,a1,r);doubleq=,p;p=a0+a1*q;printf("y=%fn",p);return0;y=a0+ai*xa0=0.8782Glal=1-978261L=2.571652Pressanykeytocontinue心得体会拉格朗日插值的优点是插值多项式特别容易建立，缺点是增加节点是原有多项式不能利用，必须