南开大学计算物理讲义2.ppt

a. 自由边界(Natural) 首尾两端没有受到任何让它们弯曲的力，即 。具体表示为 和 (3) 则要求解的方程组可写为： (10) (11) b. 固定边界(Clamped): 首尾两端点的微分值是被指定的，这里分别定为A和B。则可以推出 则要求解的方程组可写为： (10) (2) (11) 1.3 算法总结 假定有n+1个数据节点 a. 计算步长 (i = 0, 1, …, n-1) b. 将数据节点和首位端点条件带入矩阵方程 c. 解矩阵方程，求得二次微分值 。该矩阵为三对角矩阵，具体求法参见三对角矩阵的求解。 d. 计算样条曲线的系数： 下图可以看出不同的端点边界对样条曲线的影响 3. 例子 以y=sin(x)为例, x步长为1，x取值范围是[0,10]。对它使用三次样条插值，插值前后对比如下： * One way to achieve the Lagrange interpolation efficiently is by performing a sequence of linear interpolations. This scheme was first developed by Aitken(1932). We can fir