数值计算方法实验报告.doc.docVIP

《数值计算方法》实验报告 班 级 数学132班 学 号 201300144402 姓 名 袁媛 2016年 1月3日 实验报告一 实验名称 解线性方程组的直接法 2.实验题目 用追赶法求解下列方程组 3.实验目的 熟练运用已经学过的方法计算方程组，巩固已经学到的解决方程组的方法，培养使用计算机进行科学计算和解决问题的能力，熟悉了解这样的系数矩阵，能运用追赶法进行方程组的求解。 4.基础理论 设A有如下形式的分解 其中和为待定常数，则有 由可得如下计算公式： 即在A满足条件的情况下，可以把和完全确定出来，从而实现上面 给定形式的LU分解，且等于。 这样，求解三对角阵方程组Ax=f就等价于求解两个三角形方程组 从而得到公式： （1）计算和的递推公式 （2）求解 （3）求解 通常把计算和的过程称为追的过程，而把计算方程组的解的过程称为赶的过程，这一方法称为解三角方程组的追赶法。 5.实验环境 C语言，C++, Windows7 6.实验过程 第一步：将这种形式的系数矩阵进行分解，由; 第二步：由得到y; 第三步：求解得到x即可。 7.结果分析 8.附录 程序清单 function x=Hungarian a=[0 -1 -2 -3];c=[-1 -2 -3];b=[2 3 4 5];d=[6 1 -2 1]; n=length(b); u0=0;y0=0;a(1)=0; L(1)=b(1)-a(1)*u0; y(1)=(d(1)-y0*a(1))/L(1); u(1)=c(1)/L(1); for i=2:(n-1) L(i)=b(i)-a(i)*u(i-1); y(i)=(d(i)-y(i-1)*a(i))/L(i); u(i)=c(i)/L(i); end L(n)=b(n)-a(n)*u(n-1); y(n)=(d(n)-y(n-1)*a(n))/L(n); %“赶”的过程 x(n)=y(n); for i=(n-1):-1:1 x(i)=y(i)-u(i)*x(i+1); end 实验报告二 1.实验名称 Lagrange 插值法与最小二乘拟合法 2.实验题目 在某化学反应里，测得某物质的浓度y（单位：%）随时间t（单位：min）的变化数据如表5—7所列。 表5-7 变化数据 t y t y t y t y 1 4.00 5 9.22 9 10.00 13 10.50 2 6.40 6 9.50 10 10.20 14 10.55 3 8.00 7 9.70 11 10.32 15 10.58 4 8.80 8 9.86 12 10.42 16 10.60 理论上已知y与t间的关系为 其中和为待定系数。上式两端取对数可得。，，，， 根据所测数据，利用最小二乘直线拟合法先确定系数A和B，进而给出y与t间的关系并绘图展示。 3.实验目的 在求解的过程中，对数据进行拟合法,可以熟悉最小二乘拟合法的理论和方法；更加理解最小二乘拟合法的原理和求解过程，达到熟能生巧的程度，牢记该方法，学以致用。 4.基础理论 最小二乘拟合法 5.实验环境 Visual C++语言 6.实验过程 第一步：做变量替换z=lny, x=1/t， 第二步：记A=lna,B=b， 第三步：解方程z=A+Bx.， 第四步：由A,B确定a,b， 第五步：列出y与t间的关系 7.结果分析 由运行的结果可以知道A=-1.0567，B=2.4270, 8.附录 程序清单 #includestdio.h #includemath.h void main() { FILE *f; ?int?n,i;? ?float?tx,ty,x[20],y[20],sum_x=0,sum_y=0,sum_x2=0,sum_xy=0,D,A,B,a;?? ??f=fopen(.\\FittingData.txt,r);? ?fscanf(f,%d,n); ??for(i=0;in;i++)? ?{?? ?fscanf(f,%f%f,tx,ty);?? ?x[i]=1/tx;?y[i]=log(ty);? ?}? ?fclose(f);? ?for(i=0;in;i++)? ?{?? ?sum_x=sum_x