数值计算方法某实验报告材料

《数值计算方法某实验报告材料》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实用标准文案《数值计算方法》实验报告班级数学132班学号201300144402姓名袁媛2016年1月3日精彩文档实用标准文案实验报告一1.实验名称解线性方程组的直接法2.实验题目用追赶法求解下列方程组3.实验目的熟练运用已经学过的方法计算方程组，巩固已经学到的解决方程组的方法，培养使用计算机进行科学计算和解决问题的能力，熟悉了解这样的系数矩阵，能运用追赶法进行方程组的求解。4.基础理论设A有如下形式的分解其中和为待定常数，则有由可得如下计算公式：精彩文档实用标准文案即在A满足条件的情况下，可以把和完全确定出来，从而实现上面给

2、定形式的LU分解，且等于。这样，求解三对角阵方程组Ax=f就等价于求解两个三角形方程组从而得到公式：（1）计算和的递推公式（2）求解（3）求解通常把计算和的过程称为追的过程，而把计算方程组的解的过程称为赶的过程，这一方法称为解三角方程组的追赶法。5.实验环境C语言，C++,Windows7精彩文档实用标准文案6.实验过程第一步：将这种形式的系数矩阵进行分解，由;第二步：由得到y;第三步：求解得到x即可。7.结果分析8.附录程序清单functionx=Hungariana=[0-1-2-3];c=[-1-2-3];b=[2345

3、];d=[61-21];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);fori=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);endL(n)=b(n)-a(n)*u(n-1);y(n)=(d(n)-y(n-1)*a(n))/L(n);%“赶”的过程x(n)=y(n);fori=(n-

4、1):-1:1x(i)=y(i)-u(i)*x(i+1);end实验报告二1.实验名称Lagrange插值法与最小二乘拟合法2.实验题目在某化学反应里，测得某物质的浓度y（单位：%）随时间t（单位：min）的变化数据如表5—7所列。表5-7变化数据tytytyty14.0059.22910.001310.5026.4069.501010.201410.5538.0079.701110.321510.5848.8089.861210.421610.60理论上已知y与t间的关系为精彩文档实用标准文案其中和为待定系数。上式两端取对数

5、可得。做变量替换，，并记，，则有根据所测数据，利用最小二乘直线拟合法先确定系数A和B，进而给出y与t间的关系并绘图展示。3.实验目的在求解的过程中，对数据进行拟合法,可以熟悉最小二乘拟合法的理论和方法；更加理解最小二乘拟合法的原理和求解过程，达到熟能生巧的程度，牢记该方法，学以致用。4.基础理论最小二乘拟合法5.实验环境VisualC++语言6.实验过程第一步：做变量替换z=lny,x=1/t，第二步：记A=lna,B=b，第三步：解方程z=A+Bx.，第四步：由A,B确定a,b，第五步：列出y与t间的关系7.结果分析精彩文档

6、实用标准文案由运行的结果可以知道A=-1.0567，B=2.4270,8.附录程序清单#include#includevoidmain(){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;i

7、f",&tx,&ty);   x[i]=1/tx; y[i]=log(ty);  }  fclose(f); 精彩文档实用标准文案 for(i=0;i

8、D;     a=exp(A);  printf("A=%7.4f B=%6.4fy=%7.4f*e^(%6.4f/t)",A,B,a,B); } 实验报告三1.实验名称复化求积法2.实验题目测得飞机在高度h时的上升速度的数据如表6-4所列。表6-4的数据024685