计算方法大作业—质心型拉格朗日插值算法与Vandermonde矩阵逆矩阵算法.pdfVIP

东南大学 仪器科学与工程学院 戴天宇 计算方法算法设计实践 戴天宇 1 东南大学 仪器科学与工程学院 戴天宇 一、 拉格朗日插值算法的改进 1、 质心型拉格朗日插值算法原理 2、 质心型拉格朗日插值算法讨论 3、 质心型拉格朗日插值算法Matlab 实现 4、 数值算例 二、 Vandermonde 矩阵逆矩阵的快速求法 1-1、算法一原理阐述 1-2、算法一Matlab 实现 1-3、算法一数值算例以及优缺点 2-1、算法二原理阐述 2-2、算法二Matlab 实现 2-3、算法二数值算例以及优缺点 3、结论 2 东南大学 仪器科学与工程学院 戴天宇 一、 拉格朗日插值算法改进： 1、原理阐释： ( ) ∏ ( ) = =1 − 3 东南大学 仪器科学与工程学院 戴天宇 2、问题： 1）、上述2 种方法的工作量如何多少（如何尽可能的降低工作量）？ 解答： 方式一： 首先，计算 Φn （x ）需要n 次乘法 其次，每次计算wj 需要n 次乘法与加法 再其次，乘n 次，而后连加n 次 故，所需工作量为： （连乘n 次+n 次加法）*连加n 次+乘法n 次+连乘n 次=2n^2+2n 次 即，2n^2 量级的工作量 方式二： 首先，每次计算wj 需要n 次乘法与n 次加法 其次，将wj 保存后，算分子与分母的加法各需n 次，分子还需n 次乘法 故，所需工作量为： （连乘 n 次+n 次加法）*连加 n 次+乘法 n 次+连加2n 次=2n^2+3n 次 即，2n^2 量级的工作量 2 ）、上述两种方法的优缺点是什么？ 解答： 方式一： 优点是不需在添加节点时重新计算每个插值多项式，降低了一个量级的计算量， 方式二： 在方式一优点的基础上，只需要在计算wj 的时候用到连乘，省去了计算Φn （x ）的步骤， 实现了wj 的复用。并且实际上乘法比加法要更为复杂，所以只需要进行一次乘法的此方式 更加优越。 4 东南大学 仪器科学与工程学院 戴天宇 3 ）、如果是均匀节点，wj 应当如何算？ 解答： 由于是均匀节点，故步长均一致，此时的计算公式大幅简化。 令Δ为相邻节点差值，即为步长，则： 1 Wj = ∏ ( ) ≠ − 1 = ( ) ( ) ( ) ( ) ( ) − 1 ∗ − 2 ∗ … ∗ ∗ − ∗ −2 … ∗ j − n 1 = ( ) ( ) ( ) ( ) ( ) −1 − 1 ∗ − 2 ∗ … ∗ 1 ∗ −1 ∗ −2 … ∗ j − n ∗