三次样条插值法《数值分析》上机实验作业.docxVIP

昆明理工大学研究生《数值分析》上机实验作业 姓名: 学号: 专业: 课题名称 课题七三次样条插值法 1问题提出 设已知数据如下: 0.2 0.4 0.6 0.8 1.0 0.9798652 0.9177710 0.8080348 0.6386092 0.3843735 求f(x)的三次样条插值函数S(X)。 2、要求 满足自然边界条件s(0.2) = s(1.0) = 0 ； 满足第一类边界条件 s(0. 2) = 0. 20271，s(1.0) = 1.55741 o 打印输出用追赶法解出的弯曲向量 (M。，Mi, M2 , M3, MJ和 S(0.2 0.1i) (i=0,1,2,3,4,5,6,7,8 )的值。并画出 y 二 S(x)的图形。 班级、姓名、学号 三、目的和意义 由于航空、造船等工程设计的需要而发展起来所谓样条插值方法， 既保留了分段低次插值多项式的各种优点， 又提高了插值函数的光滑 性，而且具有较好的稳定性。今天，样条插值方法已成为数值逼近的 一个极其重要的分支，在许多领域里得到越来越广泛地应用。其中， 尤以三次样条插值函数应用最为广泛，如在高速飞机的机翼形体和船 体放样等方面的应用，同时在计算机作图方面更是大有作为。 它能够 解决一些既有二阶光滑度，又有二阶连续导数的方程，具有良好的收敛性和稳定性。 通过本次实验进一步了解三次样条插值函数，并通过求解三 弯矩方程组得出曲线函数组； 通过MATLA编程实现求三次样条插值函数的算法，分别考虑 不同的边界条件，同时用追赶法解出弯曲向量和 S(0.2 + 0.1) (i=0,1,2,3,4,5,6,7,8)的值。 四、计算公式 首先我们利用S(x)的二阶导数值S(Xj)二M(j二0,1,2，…，n)表 达S(x),因为在区间 凶“了刊］上S(x) =S(x)是不高于三次的多项式， 其二阶导数S(x)必是线性函数，所以可表示为： S (x)二 Mj3 X M j 1 - , x ? [x j, Xj 1 ]对 S( S (x)二 Mj h 6 利用 S (Xj )二 yj, S (xj，1) = yj .1，可定出积分 )X )Xj 1 一 Xhj (Xj 1 - X )3 …(x - Xj )3 , M j( S (x)二 M j —— M j 1 - (yj - hj 6hj 6hj 6 (yj 1 6M j 1 J、X (yj 1 6 hj ，j 0,1,2, 这里M(j 这里M(j 二 0,1, ,n)是未知的。为了确定M(j 二 0,1 厂,n)，对S(x) 求导得 S(x)… S(x)…Mj(j^ j 2hj (x - Xj )2 - yi M」1= F 由此可得 s (Xj 0) = yj 1 一* 一 匕 M j 一 也 M j j hj 6 j 3 j 类似的可求出a；x）在区间［Xj-1,Xj］上的表达式，进而得 S(Xj - 0) = y：y2 ^Mjj hj」 6 利用 S (xj - 0) = S (xj 0)可得 Pj Pj M j j 2M j M j d = d j j二1,…，n 一 1,（三弯矩方程） 其中“J 其中“J， hj hj j hj dj 二 6 dj 二 6 f Xj , Xj 1 J- fhXjJ，xJ^6fXj_1,Xj,Xj.11. j =1, ,n-1, hj_1 hj 其中有（n T ）个未知数M°,M1,...Mn，而方程只有（n-1 ）个, 当满足第一种边界条件时，可得另两个方程 2M 0 2M 0 M 一 f X0,x1 h。 —仃 - f Xnj,Xn hn 4 如果令■ 0 = 1, d 如果令■ 0 = 1, do 6 f〔X。，Xi L fo ,=1, ho 汁二 fnXn -1 , xn 1，将 上述方程综合后的一下矩阵形式: Mo | M1 入一12| 入一1 2 |M -Mn n —1 4_1 -dn 可以证明此方程组满足追赶法的条件，我们用追赶法可得 （n 1）值，将其带入公式即得s（x）。 对第二种边界条件，直接的端点方程M。二M厂f 对第二种边界条件，直接的端点方程 M。二M厂f；并且令 ■ 0 = 0,d° =2f°；dn =2fn , 则又得三弯矩方程同理即可求得解。 五、结构程序设计 1.满足自然边界条件时 自定义函数：followup.m %追赶法求m %A为线性方程组的系数矩阵 %b为常数向量 fun ctio n m=followup(A,b) n=ran k(A); fori=1: n if A(i,i)==0 disp(error:对角元素中有数据为 0); return; end end d=on es( n,1); a=on es( n-1,