最小二乘法fortran应用。ppt.ppt

最小二乘法 计算地球化学 主要内容 1. 基本原理 2. 程序编写 3. 应用举例 可以证明如果如果A是列满秩的，则该方程存在唯一解，且使得 例 1：用最小二乘法求下列超定方程组的近似解 解： 2. 程序编写 如何拟合线性y=ax+b 如：有6个数据点，(1, 6.9), (2, 9.1), (3, 10.8) (4, 13.2), (5, 14.9), (6, 17.3) 怎么编程求参数 a 和 b？ 3. 应用举例: O2饱和压的拟合与计算 Ps.in 要求： 实验数据共200个，要求从文件（文件名：Ps.in)读入，拟合参数，计算出饱和压与实验数据间的误差，找出最大误差和平均误差，并以文件（文件名：Ps.out）形式输出。 温度、饱和压及相应中间变量都采取双精度实数。 氧气饱合压拟合公式如下： 输出文件Ps.out中部分结果 程 序 作 业 拟合y=ax+b直线中的参数a和b，数据点如下： x, y 1，6.9 2, 9.1 3, 10.8 4, 13.2 5, 14.9 6, 17.3 7