三转角方程推导的提示.ppt

代入 的分段表示式，得到： 三转角方程推导的提示 问题：求三次样条函数插值多项式 ， 使其满足： 1） 2） 在所有的内节点有二阶连续导数 要求：结合三类边界条件，写出以 为待定参数的求解方法。 方法：利用埃尔米特插值公式，其中用到 的一阶导数现在未知，故将其作为 参数。 设 从而可得： 在 的表达式为： 由 的表达式可求得 ，进而求出 得出 及 ，由 在内节点连 续，可得 个方程，从而分别结合边界 条件，写出三转角方程。 三弯矩算法 待定系数法要解一个 阶的线性方程组 本方法只需求解一个不超过 阶的线性方 程组，而且有明确的力学含义。 对于 上的节点 ，记 注意到 为折线，若记 在节点 处的待定值为： 三弯矩算法 则 可用分段线性插值表示为： 对于 上的节点 ，记 注意到 为折线，若记 在节点 处的待定值为： 则 可用分段线性插值表示为： 这隐含了 满足内节点处的二阶导数连续 条件。对 作二次不定积分，可得： 其中 为待定参数，利用插值条件： ，可得： 这隐含了 满足内节点处的二阶导数连续 条件。对 作二次不定积分，可得： 其中 为待定参数，利用插值条件： ，可得： 故有： 故有： 从而有： 于是由 ，得： 即： 从而有： 于是由 ，得： 即： 若记： 则上式可简记为： 这就是待定值 满足的线性方程组， 因为其中含有三个弯矩值 ，故 称之为三弯矩方程。 若记： 则上式可简记为： 这就是待定值 满足的线性方程组， 因为其中含有三个弯矩值 ，故 称之为三弯矩方程。 三弯矩方程中共有 个未知数，但只 有 个方程，另外两个由边界条件得到。 对于第一种边界条件，即： 则有： 故有： 对于第二种边界条件，直接可得： 对于第三种边界条件，有： 故有： 对于第二种边界条件，直接可得： 对于第三种边界条件，有： 这样，有： 即： 令： 这样，有： 即： 令： 得： 把三种边界条件所得到的关系分别与 前面的方程联立，得到对应的线性方程组 如下： I: II: III: 把上面三个方程写成统一形式： 注意到 且 ，其中的系数矩阵 是三对角矩阵或仅比三对角矩阵多两个 元素的严格对角占优矩阵，因此各方程组 的解存在唯一，都可用追赶法求解，而且 是数值稳定的。把求出的弯矩值 代 入 的表达式，即得到三次样条插值多 项式： 三弯矩算法 对于区间 的划分 ， 设已知函数 在插值节点的值 ， (1) 计算参数 (2) 计算与边界条件有关的参数 I: II: III: (3) 求解与边界条件对应的三弯矩方程，得 到弯矩值 ； (4) 把 代入公式，得到三次样条插 值多项式 。 实例 设 给定边界条件 试求三次样条函数 解：先求出三弯矩方程的参数： 于是，三弯矩方程组为： 求出的解为：