Fortran使用第七章.ppt

电容器电压随时间变化曲线 讨论 两点三次Hermite插值多项式 分析 分析 实验数据和拟合多项式曲线 dimension x(7),y(7),a(3) data x/1,2,3,4,6,7,8/ data y/2,3,6,7,5,3,2/ N=7 M=3 call LSM(x,y,a,N,M) write(*,*) write(*,20) (I,a(I),I=1,M) 20 format(1X,‘a(,I2, )=,D15.6) write(*,*) end 最小二乘法求解正规方程组系数c,d ! subroutine LSM(x,y,a,N,M) dimension x(N),y(N),a(M),c(M,M),d(M) do 15 i=1,m do 10 j=1,m c(i,j)=0.0 d(i)=0.0 do 40 k=1,n c(i,j)=c(i,j)+x(k)**(i+j-2) 40 d(i)=d(i)+y(k)*x(k)**(i-1) 10 continue 15 continue write(*,*) c write(*,*) d call agaus(c,d,m,a) return end 计算程序 例14：在某电路实验中，测得电压 与电流 的一组数据如下： 试用最小二乘法求最佳数据拟合函数。 解：将数据经过绘制，它近似为一条指数曲线 ( 为待定常数) 两边取常用对数得 即 因此取指数曲线： 其中 测量数据 得正规方程组： ∴ 测量数据和拟合曲线 作 业 EX7-12: 编程计算 个数据点 的 次多项式 的最佳拟合曲线。 ） （即采用最小二乘法求 EX7-13: 利用EX7-12题程序，计算例13中 时的多项式 拟合曲线。 作 业 EX7-14: 设 间存在函数关系，对它们进行观测所 试用最小二乘法寻求经验公式，以拟合以上数据。 和 获得数据如下： EX7-15: 在某次实验中，需要观察水分的渗透速度，测得 与水的重量 的数据如下： 时间 设已知 与 之间的关系为 ，试用最小二乘法确定 参数 本章小结 一、拉格朗日插值法 点插值公式 判断 值在序列 中的位置。 注意 三点公式及计算程序 二、牛顿插值法 定义 阶差商的差商为 阶差商 构造并计算差商表 牛顿插值多项式 EX7-7: 求三次样条插值函数 ，已知 的值如下： 边界条件为 ， 。 作 业 §7.5 数值微分 一、插值型求导公式 当函数 以表格形式给出： 但 的解析表达式并不知道时，要求 在节点 处的 导数值 数值微分问题 插值型求导公式 的近似值： 作为 由拉格朗日插值原理，可建立 次插值多项式 取 的值作为 的近似值，从而得到数值公式 插值型求导公式 即使 与 的值相差不多，但导数的近似值 与 注意 的真值可能有较大差别，在使用以上求导公式 导数 时应特别注意误差分析。 由拉格朗日插值公式可知，用 近似 所产生的截断 故上式对应的导数误差为： 误差为： 其中 。 由于 是 的未知函数，无法对它的第二项 是无法估计的。 讨论 作出精确的估计，因此对任意给定的点 ，误差 但如果限定求某个节点 处的导数值，由于 ，等号右边第二项变为零，此时余项公式： 下面仅考察节点处的导数值。为简化讨论，假定所给的 高阶插值的不稳定性，实际应用时 节点是等距的。由于 的两点、三点和五点插值型求导公式。 多采用 1. 两点公式 上的函数值 ,作线性插值函数 设已给出两个节点 将上式两边对 求导，并记 ，有 即 截断误差 2. 三点公式 设已给出三个节点 上的函数值 令 ，则上式可表示为 。作二次插值函数 分别为 求导（实际是对 求导）得 将上式两端对 取 ，得三点求导公式： 用拉格朗日插值多项式 作为 的近似函数， 说明 还可以建立高阶数值微分方程公式 如对 再进行求导，可得 中的 3.五点公式 设已给出五个节点 上的函数值 同三点公式类似，可以导出五点公式： 相邻的节点，若一侧的节点不是两个（即一个或没有）， 依次为 五个相邻节点的选法，一般是在所考察的节点两侧各取 说明 则用另一侧节点补足。 解： 例10：利用 的数据表，按三点公式求节点处的 导数值。 节点及其对应的导数值的计算结果列表如下： 各点数值导数都准确到小数点后