文档介绍：数值计算方法实验报告
《数值计算方法》实验报
告
行分解，由 A LU;
第二步：由 Ly f 得到 y; 第三步：求解 Ux y 得到

x 即可。
7. 结果分析
附录程序清单
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
实验报告二
实验名称
Lagrange 插值法与最小二乘拟合法
实验题目
在某化学反应里，测得某物质的浓度 y（单位： %）随时间 t （单位：
min）的变化数据如表 5—7 所列。
表 5-7
变化数据
t
y
t
y
t
y
t
y
1

5

9

13

2

6

10

14

3

7

11

15

4

8

12

16

理论上已知 y 与 t 间的关系为
y
ae b / t
其中 a
0 和 b
0 为待定系数。上式两端取对数可得
ln y lna
b / t 。做变量替换 z
lny，x 1/ t ，并记 A
ln a ，B
b ，
则有
z
A Bx
根据所测数据，利用最小二乘直线拟合法先确定系数 A和 B，进而给
出 y 与 t 间的关系并绘图展示。
实验目的
在求解的过程中，对数据进行拟合法 , 可以熟悉最小二乘拟合法
的理论和方法；更加理解最小二乘拟合法的原理和求解过程，达到熟
能生巧的程度，牢记该方法，学以致用。
基础理论
最小二乘拟合法
实验环境
Visual C++ 语言
实验过程
第一步：做变量替换 z=lny, x=1/t ，
第二步：记 A=lna,B=b ，
第三步：解方程 z=A+Bx.，
第四步：由 A,B 确定 a,b ，
第五步：列出 y 与 t 间的关系
结果分析
由运行的结果可以知道 A=- ，B=, y
附录程序清单
#include<>
#include<>
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(".\\","r");
fscanf(f,"%d",&n);
for(i=0;i<n;i++)
{