三次样条插值函数在调洪演算中的应用
详细讲解三次样条插值法及其实现方法61页PPT
39、没有不老的誓言，没有不变的承 诺，踏 上旅途 ，义无 反顾。 40、对时间的价值没有没有深切认识 的人， 决不会 坚韧勤 勉。
66、节制使快乐增加并使享受加强。 ——德 谟克利 特 67、今天应做的事没有做，明天再早也 是耽误 了。——裴斯 泰洛齐 68、决定一个人的一生，以及整个命运 的，只 是一瞬 之间。 ——歌 德 69、懒人无法享受休息之乐。——拉布 克 70、浪费时间是一桩大罪过。——卢梭
三次NURBS曲线的插值与应用_陈绍平
第 20 卷 第 5 期 2001 年 9 月
机械科学与技术 M ECHA N ICA L SCIEN CE A ND T ECHNO L OG Y
文章编号: 1003-8728( 2001) 05-0692-02
SVeop lt.e2m0b erNo2.0501
陈绍平
三次 NU RBS 曲线的插值与应用
陈绍平
( 武汉理工大学, 武汉 430063)
摘 要: N U RBS 方法已经 被广泛地使用 在 CA D 中。本 文 将 通 过 三 次 N U RBS 曲 线 的 矩 阵 表 示 给 出 三 次 N U R BS 曲线 插值 方 法, 即反 求三 次 NU RBS 曲 线 控 制 顶点的 算法, 并给出 数字算 例和三 次 N U RBS 曲 线 插值方 法在船舶型线设计中的一个应用。
2
,
2
n32 = - ( 3n11 + n33) , n34 = 0
n41 = -
n11, n43 = -
(
1 3
n33
+
n44 +
( ¨i+ 3 ) 2 ¨i2+ 3 ¨3i+ 2
)
,
n42 =
n11 -
n43 -
n44, n44 =
( ¨i+ 3) 2 ¨i3+ 3 ¨2i+ 3
三次样条插值ppt
把以上各式由后向前代入,可得
Nn (x) f (x0) f [x0, x1](x x0) f [x0, x1, xn](x x0) (x xn1)
Rn (x) f (x) Nn (x) f [x, x0, x1, xn ](x x0) (x xn)
yi
n1 ( x) ( x xi )n' 1 ( xi )
（2）插值误差估计
定理2 设 f (n) (x) 在[a,b] 上连续，f (n1) (x)在 (a,b) 内存在, 节点 a x0 x1 xn b，Pn (x) 是拉格朗日插值多项 式，则对任意 x [a,b] , 插值余项
x4 f ( x4 ) f [x3, x4 ] f [x2 , x3 , x4 ] f [x1, x2, x3, x4 ] f [x0, x1, x2, x3, x4 ]
(2) Newton插值公式
由差约定义 x [a,b]
f (x) f (x0 ) f [x, x0 ](x x0 )
f [x, x0 ] f [x0, x1] f [x, x0, x1](x x1)
xn1] f [x1, x2 , x0 xn
xn ] n 阶差商
差商表
xk
f
(xk )
一阶 差商
二阶差商
三阶差商 四阶差商
x0 f (x0 )
x1 f (x1) f [x0, x1]
x2 f (x2 ) f [x1, x2 ] f [x0 , x1, x2 ]
x3 f (x3 ) f [x2, x3] f [x1, x2 , x3 ] f [x0, x1, x2, x3]
三次样条插值方法的应用
CENTRAL SOUTH UNIVERSITY数值分析实验报告三次样条插值方法的应用一、问题背景分段低次插值函数往往具有很好的收敛性，计算过程简单，稳定性好，并且易于在在电子计算机上实现，但其光滑性较差，对于像高速飞机的机翼形线船体放样等型值线往往要求具有二阶光滑度，即有二阶连续导数，早期工程师制图时，把富有弹性的细长木条（即所谓的样条）用压铁固定在样点上，在其他地方让他自由弯曲，然后沿木条画下曲线，称为样条曲线。
样条曲线实际上是由分段三次曲线并接而成，在连接点即样点上要求二阶导数连续，从数学上加以概括就得到数学样条这一概念。
下面我们讨论最常用的三次样条函数及其应用。
二、数学模型样条函数可以给出光滑的插值曲线（面），因此在数值逼近、常微分方程和偏微分方程的数值解及科学和工程的计算中起着重要的作用。
设区间[]b ,a 上给定有关划分b x x n =<<<=Λ10x a ，S 为[]b ,a 上满足下面条件的函数。
● )(b a C S ,2∈；● S 在每个子区间[]1,+i i x x 上是三次多项式。
则称S 为关于划分的三次样条函数。
常用的三次样条函数的边界条件有三种类型：● Ⅰ型 ()()n n n f x S f x S ''0'',==。
● Ⅱ型 ()()n n n f x S f x S ''''0'''',==，其特殊情况为()()0''''==n n x S x S 。
● Ⅲ型 ()()Λ3,2,1,0,0==j x S x S n j j ，此条件称为周期样条函数。
鉴于Ⅱ型三次样条插值函数在实际应用中的重要地位，在此主要对它进行详细介绍。
三、算法及流程按照传统的编程方法，可将公式直接转换为MATLAB可是别的语言即可；另一种是运用矩阵运算，发挥MATLAB在矩阵运算上的优势。
样条函数及三次样条插值PPT课件
(x)
lim
x xk
Sk 1( x)
lim
x
x
k
Sk (x)
lim
x
x
k
Sk1( x)
k 1,2,,n 1
------(4)
lim
x
x
k
Sk( x)
lim
x
x
k
Sk1( x)
共4n 2个条件
5
Sk (x)是[xk , xk 1 ]上的三次样条插值多项式,应有4个待定的系数 即要确定S(x)必须确定4n个待定的系数 少两个条件 并且我们不能只对插值函数在中间节点的状态进行限制 也要对插值多项式在两端点的状态加以要求 也就是所谓的边界条件:
例. 使用不同的插值方法于函数
y
1
1 x2
x [5,5]
最后，介绍一个有用的结论
定理 . 设f (x) C 2[a,b], S(x)是以xk (k 0,1,, n)
为节点, 满足任意边界条件的三次样条插值函数，
设hi
xi 1
xi
,
h
max
0in1
hi
,
min
0in1
hi
,
则当 h
c 时
S(x)和S(x)在[a,b]上一致收敛到f (x)和f (x)
------(6)
13
由(11)式,可知
S0( x0
)
6( x0
x1 h03
2 x0
) ( y1
y0 )
6 x0
2 x0 h02
4 x1
m0
6 x0
4 x0 h02
2 x1
m1
6 h02
(
MATLAB大作业 给定一个时间序列,使用三次样条插值方法进行均匀内插
MATLAB作业给定一个时间序列，使用三次样条插值方法进行均匀内插（题目的相关说明：按题目要求编写一个MATLAB程序函数，并把自己编制程序所得的结果与MATLAB库函数分析结果进行对比。
）理论基础：时间序列的概念：时间序列是一种定量预测方法，又称简单外延法，时间序列分析是根据系统观测得到的时间序列数据，通过曲线拟合和参数估计来建立数学模型的理论与方法，时间序列分析可分为以下三种情况（1）把一个时间序列的值变动为N 个组成部分，通常可以分为四种 a、倾向变动，又称长期趋势变动 b、循环变动，又称周期变动 c、季节变动，即每年有规则的反复进行变动 d、不规则变动，即随机变动。
然后把这四个综合到一起得出预测的结果。
虽然分成这四部分，但这四部分之间的相互关系是怎么样的呢，目前一般采用相乘的关系，其实各个部分都是在其他部分作用的基础上进行作用的，所以采用相乘是有一定依据的，此种方法适合于短期预测和库存预测（2）把预测对象、预测目标和对预测的影响因素都看成为具有时序的，为时间的函数，而时间序列法就是研究预测对象自身变化过程及发展趋势，如果未来预测是线性的，其数学模型为YT+L=aT+bTL,YT+L为未来预测值，aT为截距，bT为斜率，L为由T到需要预测的单位时间数（如5年、10年等）（3）根据预测对象与影响因素之间的关系及影响程度来推算未来，与目标的相关因素很多，只能选择那些因果关系较强的为预测影响的因素，此时间序列法用于短期预测比较有效，若要用于长期预测，还需要结合其他方法才行。
三次样条插值的实际应用：在制造船体和汽车外形等工艺中传统的设计方法是，首先由设计人员按外形要求，给出外形曲线的一组离散点值，施工人员准备好有弹性的样条（一般用竹条或有弹性的钢条）和压铁，将压铁放在点的位置上，调整竹条的形状，使其自然光滑，这时竹条表示一条插值曲线，我们称为样条函数。
从数学上看，这一条近似于分段的三次多项式，在节点处具有一阶和二阶连续微商。
三次样条插值函数在城区河道水面线计算中的应用
三次样条插值函数在城区河道水面线计算中的应用
随着我国城镇化的不断发展，河道曲线是居民社区的重要组成部分。
为确定城区河道的水面线，需要手动测量实测点之间的距离，并且根据测量结果拟合表达河道水面线的曲线。
然而，由于河道曲线的复杂性，传统拟合方法很难从误差数据中拟合出精确的曲线。
因此，已成为重要研究课题，如何从实测数据中准确拟合城区河道水面线，甚至以及如何提高城区河道水面线拟合的精度。
三次样条插值函数是一种经典的函数拟合方法，它按照实测点的位置，根据三次函数的性质，经过一些误差的最小化，用二次函数的系数来估计曲线的拐点。
它可以有效地表示给定的实测数据，并可以用于曲线拟合和数据拟合等应用中。
因此，三次样条插值函数在城区河道水面线计算中得到了广泛应用。
首先，三次样条插值函数可以有效地反映河道形状的复杂性。
这种插值函数可以最大程度地拟合河道水面线，从而得到更准确的曲线。
此外，由于曲线
拟合后求解出的拐点点是以矩阵形式存在的，可以更加容易的对曲线的起伏进行分析和改善。
这样可以有效提高河道水面线计算的准确度。
另外，三次样条插值函数还可以用于河道水面线设计。
研究者首先根据实际河道情况，利用三次样条插值函数拟合出一个基础曲线，然后增加实测点以改善曲线的拐点，或者增加实测点的数量，来对河道水面线进行修正，达到期望的设计要求。
综上所述，三次样条插值函数可以有效地表示河道的复杂性，可以拟合出更准确的曲线，同时也可以用于曲线设计，具有很强的实际意义。
因此，三次样条插值函数可以在城区河道水面线计算中得到有效应用，对于城区河道水面线计算和设计具有重要的意义。
数值计算方法三次样条插值
4.4 三次样条插值
令
A1
j1 (u )
(1
2
u
x hj
j 1
)(
u
xj hj
)2
A2
j (u )
(1
2
u
x hj
j
)(
u
x hj
j 1
)2
B1
j1 (u )
(u
u x j 1 )(
xj hj
)2
B2
j (u )
(u
x
j )(
u
x hj
j
)2
分段三次Hermite插值算法
I2(x)
I
n
(
x
)
x ( x0 , x1)
x ( x1, x2 ) ...... x ( xn1, xn )
其I中 j xxj1 xxj j yj1xxj xxjj 11yj yj1(xxj1)(yj yj1)/(xj xj1)
缺点：I(x)连续，但不光滑，精度较低,仅在 hm 1jan{xhj xj xj1}足够小才能较好。的逼近
ss((xxn0
) )
f f
( x0 ) (xn )
m0 mn
三次样条插值在工程拟合中的应用
三次样条插值在工程拟合中的应用摘要: 介绍了工程实验、勘测、设计中常见的列表函数之数值插值方法、程序实现及工程应用, 应用此法可方便地将任何列表函数计算到工程设计、施工所需要的精确程度, 给出了各参数随主要参数变化而变化的光滑曲线, 并将其应用推广到一般情况.关键词: 列表函数; 数值拟合; 三次样条插值; MA TLAB 程序设计与应用在实际工程中, 广泛存在这样的问题: 根据设计要求和具体的工程条件, 在初始设计阶段会勘测得到若干组该工程的控制参数, 但这些参数之间彼此离散、不够密集, 利用它们来施工则不能满足施工的精度要求. 为了解决这一问题, 需要对已知的参数数据进行分析处理, 进行必要的插值、拟合, 以达到施工所需要的数据精度.本文以工程实例为基础, 对实际工程中插值方法的选取、插值的实现和插值曲线的拟合加以讨论, 提出能得到较合乎实际的插值方法, 给出一般工程人员就能实现的计算方法以及能得到光滑曲线的拟合方法.1 工程应用实例表1 所示的为某双曲拱坝体形原始参数[ 1对于这一类工程列表参数有一个显著的特点:尽管不同工程的参数多寡不同, 但都是由n 行k列的离散的列表数据给出, 虽然同一行代表某工程特定位置的几个参数(或高程参数, 或上游半径参数⋯) , 但相邻两行由于位置距离太大, 两行各参数之间究竟存在什么数值关系, 对工程设计、施工有何影响, 这是工程技术人员需要弄清楚的[ 2 ].以双曲拱坝为例, 它沿整个高程的变化是一个连续光滑的空间曲面. 从施工需要来看, 这些数据太稀疏, 难以满足设计、施工放样与钢筋配置等要求, 如果照此施工, 则有可能达不到工程精度、降低工程效率; 从计算机图形模拟来看, 要生成这个曲面仅由这一列表函数是得不到光滑曲面的, 是不可取的. 所以, 为使计算精确, 满足工程施工过程中任何断面位置、任意水平位置、任意高程位置所必需的施工数据与设计图纸, 保证工程施工的高品质,就要求作精确的数据处理.进一步分析可知, 在这些参数表中, 各行的参数都随某一主要参数的变化而变化, 如上游半径参数随高程的变化而变化⋯, 它们的这种函数关系,在数值分析中有许多的方法可以求得. 但是哪种方法能更好、更合乎实际地给出平滑曲线呢? 下面所选的插值方法能够较好地满足这一要求.2 插值方法的选择在数值分析中, 这种插值过程可具体使用线性( 1inear ) 插值、三次样条( sp line ) 插值、立方(cub ic) 插值等方法, 在曲线插值法中最常用的是线性插值法, 它是估计两个主干点之间数值的最简单、最易实现的方法, 但采用线性插值法会有以下缺点:一是使得曲线不能显示连接主干点间的凸状弧线;二是使得从曲线导出远期曲线时会形成人为的“尖头”(sp ikes) [ 2 ].因此, 通常采用样条法来构造曲线. 样条法是用一平滑曲线来对各主干点进行拟合的方法. 它是通过构造多项式(一个或一组不同阶多项式) 来形成一条把所有主干点连接起来的平滑曲线. 一般常常选择三次曲线(根据三次插值样条函数所得的曲线) 进行拟合.通常, 在[a, b ]上的以x i ( i=0, 1, 2, ⋯, n) 为节点的三次插值样条函数[ 3 ] 定义如下: 给定区间[a, b ]的一个划分$: a= x 0< x 1< x 2< ⋯< x n = b和区间[a, b ]上的一个函数f (x ) , 若函数S (x ) 满足下列条件:(1) 一致通过n+ 1 个插值点(x i, y i) , 即S (x i) = f (x i) = y i ( i= 0, 1, 2, ⋯, n) ;(2) 二阶连续, 即S (x ) ∈C2 [a, b ];(3) 三次分段, 即在每一个小区间[ x i- 1, x i ]( i= 1, 2, ⋯, n) 上均为三次多项式.则称S (x ) 为函数f (x ) 的三次插值样条函数. 在构造三次插值样条函数时, 为确定S (x ) 应根据n+ 1个插值条件, 3n- 3 个连续条件以及给定的边界条件, 再利用节点处的一阶导数或二阶导数就可构造出三次插值样条函数. 在构造曲线过程中, 关键是估计三次多项式函数和确定样条函数形式.从以上理论分析可知, 三次活动曲线具有优良的数学特征, 而且用三次曲线去拟合时, 其结果要比线性插值估计更接近于工程实际情况[ 4 ]. 三次曲线法又可分为三次样条插值法和立方插值法. 在数值分析中有许多的方法, 限于篇幅, 本文仅以工程上用得较多的、具有优良效果的三次样条插值为例介绍插值方法.3 插值计算原理三次样条函数的数学原理及其子程序, 可见于多种数学著作[ 5 ]与算法手册. 这里作简单介绍.由于拱坝或其他工程曲面都是连续而光滑的空间曲面, 它的断面高程自坝底至坝顶均满足a= j 1< j 2< ⋯< j n= b,且每一位置(高程) 都对应有一组几何参数: y 1, y 2,⋯, y n. 如上游半径、下游半径、拱厚等(见表1 所列) , 因此对于一组高程插值点j 1= t1< t2< t3< ⋯< tm ≤j n ,可用三次自然样条函数S (x ) 求解它们在各插值点的函数值及其一阶导数S ′(x ) 和二阶导数S ″(x ).三次样条函数S (x ) 是用分段三次多项式逼近函数y = f (x ) , 且满足S (x ) 为区间[a, b ]上曲线y= f (x ) 的三次样条插值函数的三个条件.经两次积分, 可得三次样条插值函数S (x ) 的表达式为利用函数S (x ) 在样点x i 处具有连续二阶导数的条件, 再根据三次自然样条插值法, 增加自然边界条件得到如下方程组:解上述方程组, 求得M i ( i= 0, 1, 2, ⋯, n) 代入S (x ) 公式, 即可得每个子区间[ x i- 1, x i ] ( i= 1, 2,⋯, n) 上的三次样条函数.根据上述原理, 对工程原始列表数进行插值计算, 即可满足多种施工要4 插值方法的实现由以上可以看出, 三次样条插值的关键是寻找插值函数, 但插值函数寻找相当复杂, 对于一般的工程人员很难完成, 那么怎样才能使三次样条插值这一优秀的插值方法被人们所掌握呢?M athworks公司推出了功能强大的数学计算软件MA T2LAB[ 6 ] , 它不但使源程序编写简单、源程序代码简短(因为现成的三次样条插值函数可供使用) , 而且可以利用其强大的作图功能方便地拟合出光滑曲线. 因此, 本文选用MA TLAB 语言作为计算语言MA TLAB 程序设计原理:在以上参数表中, 各行的各参数都随高程这一主要参数的变化而变化, 根据它们变化的这种函数关系, 以高程为插值的已知节点(其中已知节点个数n = 6) , 为使插值结果一致通过这些节点, 以1. 36为步长调用插值函数进行插值.MA TLAB 程序设计算法:( 1) 写入原始参数矩阵, 以同一组参数为行,以同一种参数为列;(2) 产生插值的精度矩阵, 在最小值与最大值之间以1. 36 为步长, 产生矩阵;(3) 调用MA TLAB 中的三次样条插值函数,产生插值结果矩阵, 以对每一种参数的插值结果为行产生矩阵, 再转置.MA TLAB 程序设计:x 0= [470∶1. 36∶504 ];ou t= [x 0; sp line (x (1∶6) , x (7∶12) , x 0) ; sp line (x (1∶6) , x (13∶18) , x 0) ; sp line (x (1∶6) , x (19∶24) , x 0)sp line (x (1∶6) , x (25∶30) , x 0) ; sp line (x (1∶6) , x (31∶36) , x 0) ; sp line (x (1∶6) , x (37∶42) , x 0) ]′运算数据分析:(1) 这组运算数据一致通过已知节点, 而且偏差较小、数学处理和程序设计都大大简化(与文献[1 ]相比).(2) 经过以上的运算, 可以使原来仅有的6 组数据变为26 组, 而且还可以根据工程人员的需要对上述程序步长进行修改, 就可任意提高精度, 从而使工程人员能够更好地了解各种参数在各点的数据, 使工程精度大大提高5 插值曲线拟合当然, 无论以多么小的数为步长、无论给出多少组数据, 这些参数还是一些离散的数据, 在有些情况下, 工程人员要了解某些数据随某一主要参数的变化而变化的连续曲线, 这时, 可以在数据插值的基础上, 发挥MA TLAB 在图形处理上的强大功能, 对以上插值所得的数据进行曲线拟合, 以便更好地了解各参数随某一主要参数变化而变化的趋势.在以上插值数据的基础上,在上面程序的尾部编写MA TLAB 作图程序, 作图程序如下, 运行后得到图1 所示插值拟合曲线.p lo t (x 0, ou t (27∶52) ,‘- ’)ho ld onp lo t (x 0, ou t (53∶78) ,‘- + ’)p lo t (x 0, ou t (79∶104) ,‘∶’)p lo t (x 0, ou t (105∶130) ,‘- - ’)p lo t (x 0, ou t (131∶156) ,‘- 3 ’)p lo t (x 0, ou t (157∶182) ,‘- . ’)legend (‘上游半径’,‘下游半径’,‘拱厚’,‘半中心角’,‘圆心距’,‘淤沙高程’)ho ld offgrid on从图1 中, 可以看到各参数随高程的变化而变化的曲线, 从而更好地去了解各参数的变化规律,实现对工程各参数的整体把握, 这是一般数值处理方法所无法实现的.6 小结以上仅为三次样条插值及其实现方法的一个实例, 本文在插值方法的选择上选取了能够得到平滑曲线的、具有优良数学特征的三次样条插值法;在插值的实现上选取了具有强大计算功能的数学软件MA TLAB, 它能够以较少的编码, 较简单的语句实现这一复杂的计算, 并能得到较合理的结论; 在曲线的拟合上我们在插值的基础上同样选取具有强大图形处理功能的MA TLAB 软件, 从而形成较准确、较平滑、较合实际的曲线. 总之, 以上所提供的方法是三次样条插值和MA TLAB 科学计算语言在工程中应用的一个实例, 它能使计算较简便, 又能很好地满足光滑性要求, 使曲线也不失真.实现了工程数学、计算数学、程序设计的结合与简化.三次样条插值不仅在工程方面, 而且在测绘、勘察、预测等方面都有着十分广泛的应用参考文献:[ 1 ]彭荣利, 靳萍, 欧阳建国. 工程列表函数的数值拟合与应用[J ]. 武汉大学学报(工学版) , 2002, 35 (4) : 42～45.[ 2 ]王瑞华. 水利工程数据插值计算及图形处理[J ]. 农田水利与小水电, 1994 (8) : 15～19.[ 3 ]鞠时光, 郭伟刚. 实用三次样条插值函数[J ]. 小型微型计算机系统, 1992, 13 (9) : 20～23.[ 4 ]谢赤, 钟钻. 插值法在零息收益曲线构造中的实证研究[J ]. 数量经济技术经济研究, 2002 (4) : 31～34.[ 5 ]曾绍标, 韩秀芹. 工程数学基础[M ]. 北京: 科学出版社, 2001.[ 6 ]王沫然. MA TLAB 与科学计算(第2 版) [M ]. 北京: 电子工业出版社, 2003。
三次样条插值算法详解
三次样条插值算法要求数据点数量较多，且在某些情况下可能存在数值不稳定性，如数据 点过多或数据点分布不均等情况。此外，该算法对于离散数据点的拟合效果可能不如其他 插值方法。
对未来研究的展望
01
02
03
改进算法稳定性
针对数值不稳定性问题， 未来研究可以探索改进算 法的数值稳定性，提高算 法的鲁棒性。
3
数据转换
对数据进行必要的转换，如标准化、归一化等， 以适应算法需求。
构建插值函数
确定插值节点
根据数据点确定插值节点，确保插值函数在节点处连续且光滑。
构造插值多项式
根据节点和数据点，构造三次多项式作为插值函数。
确定边界条件
根据实际情况确定插值函数的边界条件，如周期性、对称性等。
求解插值函数
求解线性方程组
06
结论
三次样条插值算法总结
适用性
三次样条插值算法适用于各种连续、光滑、可微的分段函数插值问题，尤其在处理具有复 杂变化趋势的数据时表现出色。
优点
该算法能够保证插值函数在分段连接处连续且具有二阶导数，从而在插值过程中保持数据 的平滑性和连续性。此外，三次样条插值算法具有简单、易实现的特点，且计算效率较高 。
根据数据点的数量和分布，合理分段，确保 拟合的精度和连续性。
求解线性方程组
使用高效的方法求解线性方程组，如高斯消 元法或迭代法。
结果输出
输出拟合得到的插值函数，以及相关的误差 分析和图表。
03
三次样条插值算法步骤
数据准备
1 2
数据收集
收集需要插值的原始数据点，确保数据准确可靠。
数据清洗
对数据进行预处理，如去除异常值、缺失值处理 等。
第五章(3)三次样条插值
6( xi xi 1 2 x ) ( yi 1 yi ) 3 hi 1
而
2 4 6 S ( xi 0) mi 1 mi 2 ( yi yi 1 ) hi hi hi 4 2 6 S ( xi 0) mi m i 1 2 ( yi 1 yi ) hi 1 hi 1 hi 1
n
当n 时，Ln ( x )只在 | x | 3.63 内收敛，而在该区间外 是发散的。
从图中可以看出，在 0 附近插值效果是好的，即余项较 小，另一种现象是插值多项式随节点增多而振动。这种插值 多项式当节点增加时反而不能更好地接近被插值函数的现象， 称为龙格现象。
上述现象告诉我们用高次插值多项式是不 妥当的，从数值计算上可解释为高次插值多项 式的计算会带来舍入误差的增大，从而引起计 算失真。因此，实践上作插值时一般只用一次、 二次最多用三次插值多项式。
式中x [ xi 1 , xi ] (i 1,2,, n)
第(2)步
为了确定mi，需要用到S ( x )的二阶导数在节点连续 的条件， S ( x )在[ xi 1 , xi ]和[ xi , xi 1 ]上的二阶导数分别为
Si( x ) 6 x 2 xi 1 4 xi 6 x 4 xi 1 2 xi mi 1 mi 2 2 hi hi ( x [ xi 1 , xi ])
若记hi xi xi 1，则上式可写为
( x x i ) 2 hi 2( x x i 1 ) ( x x i 1 ) 2 hi 2( x i x ) Si ( x) y i 1 yi 3 3 hi hi ( x x i ) 2 ( x x i 1 ) ( x x i 1 ) 2 ( x x i ) m i 1 mi 2 2 hi hi
三次样条插值在实际问题中的应用和实现
曲线 ， 这样的曲线通常是光滑的。 ２ 水库是实际生活 中用来 防洪兴 利 的主要 建筑 物。 当 、 洪水来临时 ， 库调度人 员应根据 当前时段入库 的水量迅 速 水 作 出正确 的水库调度方案 。在水库调洪过程 中， 始终遵循水 库水量平衡方程 ：
三 次 样 条 插 值 在 实 际 问 题 中 的 应 用 和 实 现
杨 虎 山
（ 州 师 范 学 院 ， 西 忻 州 ０ ４０ ） 忻 山 ３０ ０
摘
要： 文章通过 对 三 次样 条插 值 函数 的研 究 ， 并借 助 可视 化 编程 工具 Ｖｓａ Ｂｓ ｉ ｌ ａｉ 言 ｕ ｃ语
的强 大功 能 ， 程解 决 了水库调 度 和风 洞 实验 中的一 些 实际 问题 ， 而为 它们 日后 的应 用提 供 编 从
维普资讯
第２ ４卷
第 ２期
忻 州 师 范 学 院 学 报
Ｊ ＯＵＲ ＮＡＬ ＯＦ ＸＩ ＨＯＵ Ｔ ＡＣＨＥ ＵＮ ＶＥＲ Ｉ Ｙ ＮＺ Ｅ ＲＳ Ｉ ＳＴ
Ｖｏ ． ４ Ｎｏ． １２ ２ Ａｐ ．２０ ｒ ０８
２０ ０ ８年 ４月
件 的具 体路 径 位 于 ： 工 程 ＼ 件 ＼ ｃ ｓｆＦｅｇ ｄＣ ｎｒｌ ． 部 Ｍｉ ｏ ｌ ａ ｏｔ ０ ｏ ｒ ｔ ｘ ｏ６
２ 问题 解 决 的 思路 和 方 法
其 中ｇ 。Ｑ ，。ｑ—— 时段 始 、 的入库 流 量和 出库 流 Ｑ ， ｑ ，２ 末
三次样条插值法
三次样条插值法根据李庆阳的《数值分析》这本教材中的讲解编写的程序，使用的是第一边界条件，用追赶法求解了M矩阵。
为了调用方便，我将整个函数所有的信息构造成一个结构体，输入插值点的坐标和数量，定义边界条件后，将这个结构体的指针作为参数传给Spline3（）函数，就完成了函数计算，计算结果也存储在该结构体中。
程序如下：/*=================================================================== ====*///=====================三次样条插值的函数S(x)实现=============================// 创建人：汪雅楠// 北京交通大学 QQ312818820// 说明：根据研究生教材《数值分析》（李庆杨）第51页~第56页编写/* 初始条件： 1. 已知两端的一阶导数值2. 已知两端的二阶导数值3. 周期样条函数### 此函数选择1条件 ###函数建立： 1. 设定样条函数S(x)的一阶导数为变量ki,用分段三次Hermitte差值2. 设定样条函数S(x)的二阶导数为变量Ki,用分段积分### 此函数选择2方法 ###矩阵求解：追赶法求解严格三对角占优矩阵{M},根据教材第195页编写*//*=================================================================== ====*/#include <stdio.h>///////////////////////////////////////////////////////////////////// ///////////#define MAXNUM 50 //定义样条数据区间个数最多为50个typedef struct SPLINE //定义样条结构体，用于存储一条样条的所有信息{ //初始化数据输入float x[MAXNUM+1]; //存储样条上的点的x坐标，最多51个点float y[MAXNUM+1]; //存储样条上的点的y坐标，最多51个点unsigned int point_num; //存储样条上的实际的点的个数float begin_k1; //开始点的一阶导数信息float end_k1; //终止点的一阶导数信息//float begin_k2; //开始点的二阶导数信息//float end_k2; //终止点的二阶导数信息//计算所得的样条函数S(x)float k1[MAXNUM+1]; //所有点的一阶导数信息float k2[MAXNUM+1]; //所有点的二阶导数信息//51个点之间有50个段，func[]存储每段的函数系数float a3[MAXNUM],a1[MAXNUM];float b3[MAXNUM],b1[MAXNUM];//分段函数的形式为Si(x) = a3[i] * {x(i+1) - x}^3 + a1[i] * {x(i+1) - x} +// b3[i] * {x - x(i)}^3 + b1[i] * {x - x(i)}//xi为x[i]的值，xi_1为x[i+1]的值}SPLINE,*pSPLINE;typedef int RESULT; //返回函数执行的结果状态，下面为具体的返回选项#ifndef TRUE#define TRUE 1#endif#ifndef FALSE#define FALSE -1#endif#ifndef NULL#define NULL 0#endif#ifndef ERR#define ERR -2#endif///////////////////////////////////////////////////////////////////// //////////////*=================================================================== ============*** 函数名称：Spline3()*** 功能说明：完成三次样条差值，其中使用追赶法求解M矩阵*** 入口参数：(pSPLINE)pLine 样条结构体指针pLine中的x[],y[],num,begin_k1,end_k1*** 出口参数：(pSPLINE)pLine 样条结构体指针pLine中的函数参数*** 返回参数：返回程序执行结果的状态TRUE or FALSE===================================================================== ===========*/RESULT Spline3(pSPLINE pLine){float H[MAXNUM] = {0}; //小区间的步长float Fi[MAXNUM] = {0}; //中间量float U[MAXNUM+1] = {0}; //中间量float A[MAXNUM+1] = {0}; //中间量float D[MAXNUM+1] = {0}; //中间量float M[MAXNUM+1] = {0}; //M矩阵float B[MAXNUM+1] = {0}; //追赶法中间量float Y[MAXNUM+1] = {0}; //追赶法中间变量int i = 0;////////////////////////////////////////计算中间参数if((pLine->point_num < 3) || (pLine->point_num > MAXNUM + 1)){return ERR; //输入数据点个数太少或太多}for(i = 0;i <= pLine->point_num - 2;i++){ //求H[i]H[i] = pLine->x[i+1] - pLine->x[i];Fi[i] = (pLine->y[i+1] - pLine->y[i]) / H[i]; //求F[x(i),x(i+1)] }for(i = 1;i <= pLine->point_num - 2;i++){ //求U[i]和A[i]和D[i]U[i] = H[i-1] / (H[i-1] + H[i]);A[i] = H[i] / (H[i-1] + H[i]);D[i] = 6 * (Fi[i] - Fi[i-1]) / (H[i-1] + H[i]);}//若边界条件为1号条件，则U[i] = 1;A[0] = 1;D[0] = 6 * (Fi[0] - pLine->begin_k1) / H[0];D[i] = 6 * (pLine->end_k1 - Fi[i-1]) / H[i-1];//若边界条件为2号条件，则//U[i] = 0;//A[0] = 0;//D[0] = 2 * begin_k2;//D[i] = 2 * end_k2;/////////////////////////////////////////追赶法求解M矩阵B[0] = A[0] / 2;for(i = 1;i <= pLine->point_num - 2;i++){B[i] = A[i] / (2 - U[i] * B[i-1]);}Y[0] = D[0] / 2;for(i = 1;i <= pLine->point_num - 1;i++){Y[i] = (D[i] - U[i] * Y[i-1]) / (2 - U[i] * B[i-1]);}M[pLine->point_num - 1] = Y[pLine->point_num - 1];for(i = pLine->point_num - 1;i > 0;i--){M[i-1] = Y[i-1] - B[i-1] * M[i];}//////////////////////////////////////////计算方程组最终结果for(i = 0;i <= pLine->point_num - 2;i++){pLine->a3[i] = M[i] / (6 * H[i]);pLine->a1[i] = (pLine->y[i] - M[i] * H[i] * H[i] / 6) / H[i];pLine->b3[i] = M[i+1] / (6 * H[i]);pLine->b1[i] = (pLine->y[i+1] - M[i+1] * H[i] * H[i] / 6) /H[i]; }return TRUE;}///////////////////////////////////////////////////////////////////// /////////////SPLINE line1;pSPLINE pLine1 = &line1;///////////////////////////////////////////////////////////////////// /////////////main(){line1.x[0] = 27.7;line1.x[1] = 28;line1.x[2] = 29;line1.x[3] = 30;line1.y[0] = 4.1;line1.y[1] = 4.3;line1.y[2] = 4.1;line1.y[3] = 3.0;line1.point_num = 4;line1.begin_k1 = 3.0;line1.end_k1 = -4.0;Spline3(pLine1);return 0;}///////////////////////////////////////////////////////////////////// /////////////。
三次样条插值方法在工程实践中的应用
h[i],x[i],(y[i]-M[i]*h[i]*h[i]/6)/h[i],x[i-1]);
}
四、实验过程中需要记录的数据：
在程序中，主要记录的数据有hi，ui，vi，gi以及最重要的是在区间[xi-1,xi]的三次插值函数的二阶导数Mi，最后输出在各个小区间内的三次样条插值函数即可。
当x属于〔100.00，150.00〕时
s(x)=-0.000024(150.000000-x)3-0.000014(x-100.000000)3+0.932286(150.000000-x)
+1.033857(x-100.000000)
当x属于〔150.00，200.00〕时
s(x)=-0.000014(200.000000-x)3-0.000024(x-150.000000)3+1.033857(200.000000-x)+
南昌航空大学
数学与信息科学学院
实 验 报 告
课程名称：计 算 方 法
实验名称：三次样条插值方法在工程实践中的应用
实验类型：验证性√综合性□设计性□
实验室名称：D504
班级学号：08061115
学生姓名：杨朝峰
任课教师〔教师签名〕：
成绩：
实验日期：2021-11-13
公路平面曲线的设计
一、实验目的
由实验一，我们已经对插值的思想和原理有了更深入的理解。从中可以知道，分段低次插值虽然解决了高次插值的振荡现象和数值不稳定现象，使得插值多项式具有一致收敛性，保证了插值函数整体的连续性，但在函数插值节点处不能很好地保证光滑性要求，这在某些要求光滑性的工程应用中是不能接受的。如飞机的机翼一般要求使用流线形设计，以减少空气阻力。因此，在分段插值的根底上，引进了一种新的插值方法，在保证原方法的收敛性和稳定性的同时，又使得函数具有较高的光滑性。
三次样条插值法
鉴于高次插值不收敛又不稳定的特点，低次插值既具有收敛 性又具有稳定性，因此低次值更具有实用价值，但是低次插 值的光滑性较差，比如分段线性插值多项式在插值区间中仅 具有连续性，在插值节点处有棱角，一阶导数不存在；分段 三次Hermite插值多项式在插值区间中仅具有一阶导数即一 阶光滑性但不具备二阶光滑性，不能满足某些实际应用比如 汽车、轮船、飞机等的外形中流线形设计。样条是在二十世 纪初期经常用于图样设计的一种富有弹性的细长条，多个样 条互相弯曲连接后沿其边缘画出的曲线就是三次样条曲线。 后来数学上对其进行了抽象，定义了m次样条函数，并成为 数值逼近的重要研究分枝，进一步扩大了样条函数的应用范 围。
立即可得下式：
21
n m1 n mn1 2mn gn
其中：
h0 hn1 n , n 1 n h0 hn1 h0 hn1
y1 y0 yn yn1 g n 3 n n h0 hn 1
联合基本方程得一个广义三对角或周期三对角方程组：
S ( xn )
6 4 2 ( yn yn1 ) mn 1 mn M n 2 hn 1 hn 1 hn1
19
稍加整理得
y1 y0 h0 2m0 m1 3 M 0 g0 h0 2 yn yn 1 hn1 mn1 2mn 3 M n gn hn 1 2
hi i hi hi 1 hi 1 i 1 i hi hi 1
称为三转 角法基本 方程组
yi yi 1 yi 1 yi gi 3(i i ) hi 1 hi
共个n 1个方程, n 1个未知量
以上推导还没有考虑边界 条件！针对不同类型的三 次样条问题，就可以导出 不同的方程组！
(建筑工程管理]三次样条插值在工程拟合中的应用
![(建筑工程管理]三次样条插值在工程拟合中的应用](https://img.taocdn.com/s3/m/8434864e360cba1aa811dae4.png)
（建筑工程管理）三次样条插值在工程拟合中的应用三次样条插值于工程拟合中的应用摘要:介绍了工程实验、勘测、设计中常见的列表函数之数值插值方法、程序实现及工程应用,应用此法可方便地将任何列表函数计算到工程设计、施工所需要的精确程度,给出了各参数随主要参数变化而变化的光滑曲线,且将其应用推广到壹般情况.关键词:列表函数;数值拟合;三次样条插值;MATLAB程序设计和应用于实际工程中,广泛存于这样的问题:根据设计要求和具体的工程条件,于初始设计阶段会勘测得到若干组该工程的控制参数,但这些参数之间彼此离散、不够密集,利用它们来施工则不能满足施工的精度要求.为了解决这壹问题,需要对已知的参数数据进行分析处理,进行必要的插值、拟合,以达到施工所需要的数据精度.本文以工程实例为基础,对实际工程中插值方法的选取、插值的实现和插值曲线的拟合加以讨论,提出能得到较合乎实际的插值方法,给出壹般工程人员就能实现的计算方法以及能得到光滑曲线的拟合方法.1工程应用实例表1所示的为某双曲拱坝体形原始参数[1对于这壹类工程列表参数有壹个显著的特点:尽管不同工程的参数多寡不同,但均是由n行k列的离散的列表数据给出,虽然同壹行代表某工程特定位置的几个参数(或高程参数,或上游半径参数⋯),但相邻俩行由于位置距离太大,俩行各参数之间究竟存于什么数值关系,对工程设计、施工有何影响,这是工程技术人员需要弄清楚的[2].以双曲拱坝为例,它沿整个高程的变化是壹个连续光滑的空间曲面.从施工需要来见,这些数据太稀疏,难以满足设计、施工放样和钢筋配置等要求,如果照此施工,则有可能达不到工程精度、降低工程效率;从计算机图形模拟来见,要生成这个曲面仅由这壹列表函数是得不到光滑曲面的,是不可取的.所以,为使计算精确,满足工程施工过程中任何断面位置、任意水平位置、任意高程位置所必需的施工数据和设计图纸,保证工程施工的高品质,就要求作精确的数据处理.进壹步分析可知,于这些参数表中,各行的参数均随某壹主要参数的变化而变化,如上游半径参数随高程的变化而变化⋯,它们的这种函数关系,于数值分析中有许多的方法能够求得.可是哪种方法能更好、更合乎实际地给出平滑曲线呢?下面所选的插值方法能够较好地满足这壹要求.2插值方法的选择于数值分析中,这种插值过程可具体使用线性(1inear)插值、三次样条(spline)插值、立方(cubic)插值等方法,于曲线插值法中最常用的是线性插值法,它是估计俩个主干点之间数值的最简单、最易实现的方法,但采用线性插值法会有以下缺点:壹是使得曲线不能显示连接主干点间的凸状弧线;二是使得从曲线导出远期曲线时会形成人为的“尖头”(spikes)[2].因此,通常采用样条法来构造曲线.样条法是用壹平滑曲线来对各主干点进行拟合的方法.它是通过构造多项式(壹个或壹组不同阶多项式)来形成壹条把所有主干点连接起来的平滑曲线.壹般常常选择三次曲线(根据三次插值样条函数所得的曲线)进行拟合.通常,于[a,b]上的以xi(i=0,1,2,⋯,n)为节点的三次插值样条函数[3]定义如下:给定区间[a,b]的壹个划分$:a=x0<x1<x2<⋯<xn=b和区间[a,b]上的壹个函数f(x),若函数S(x)满足下列条件:(1)壹致通过n+1个插值点(xi,yi),即S(xi)=f(xi)=yi(i=0,1,2,⋯,n);(2)二阶连续,即S(x)∈C2[a,b];(3)三次分段,即于每壹个小区间[xi-1,xi](i=1,2,⋯,n)上均为三次多项式.则称S(x)为函数f(x)的三次插值样条函数.于构造三次插值样条函数时,为确定S(x)应根据n+1个插值条件,3n-3个连续条件以及给定的边界条件,再利用节点处的壹阶导数或二阶导数就可构造出三次插值样条函数.于构造曲线过程中,关键是估计三次多项式函数和确定样条函数形式.从之上理论分析可知,三次活动曲线具有优良的数学特征,而且用三次曲线去拟合时,其结果要比线性插值估计更接近于工程实际情况[4].三次曲线法又可分为三次样条插值法和立方插值法.于数值分析中有许多的方法,限于篇幅,本文仅以工程上用得较多的、具有优良效果的三次样条插值为例介绍插值方法.3插值计算原理三次样条函数的数学原理及其子程序,可见于多种数学著作[5]和算法手册.这里作简单介绍.由于拱坝或其他工程曲面均是连续而光滑的空间曲面,它的断面高程自坝底至坝顶均满足a=j1<j2<⋯<jn=b,且每壹位置(高程)均对应有壹组几何参数:y1,y2,⋯,yn.如上游半径、下游半径、拱厚等(见表1所列),因此对于壹组高程插值点j1=t1<t2<t3<⋯<tm≤jn,可用三次自然样条函数S(x)求解它们于各插值点的函数值及其壹阶导数S′(x)和二阶导数S″(x).三次样条函数S(x)是用分段三次多项式逼近函数y=f(x),且满足S(x)为区间[a,b]上曲线y=f(x)的三次样条插值函数的三个条件.经俩次积分,可得三次样条插值函数S(x)的表达式为利用函数S(x)于样点xi处具有连续二阶导数的条件,再根据三次自然样条插值法,增加自然边界条件得到如下方程组:解上述方程组,求得Mi(i=0,1,2,⋯,n)代入S(x)公式,即可得每个子区间[xi-1,xi](i=1,2,⋯,n)上的三次样条函数.根据上述原理,对工程原始列表数进行插值计算,即可满足多种施工要4插值方法的实现由之上能够见出,三次样条插值的关键是寻找插值函数,但插值函数寻找相当复杂,对于壹般的工程人员很难完成,那么怎样才能使三次样条插值这壹优秀的插值方法被人们所掌握呢?Mathworks公司推出了功能强大的数学计算软件MAT2LAB[6],它不但使源程序编写简单、源程序代码简短(因为现成的三次样条插值函数可供使用),而且能够利用其强大的作图功能方便地拟合出光滑曲线.因此,本文选用MATLAB语言作为计算语言MATLAB程序设计原理:于之上参数表中,各行的各参数均随高程这壹主要参数的变化而变化,根据它们变化的这种函数关系,以高程为插值的已知节点(其中已知节点个数n=6),为使插值结果壹致通过这些节点,以1.36为步长调用插值函数进行插值.MATLAB程序设计算法:(1)写入原始参数矩阵,以同壹组参数为行,以同壹种参数为列;(2)产生插值的精度矩阵,于最小值和最大值之间以1.36为步长,产生矩阵;(3)调用MATLAB中的三次样条插值函数,产生插值结果矩阵,以对每壹种参数的插值结果为行产生矩阵,再转置.MATLAB程序设计:x0=[470∶1.36∶504];out=[x0;spline(x(1∶6),x(7∶12),x0);spline(x(1∶6),x(13∶18),x0);spline(x(1∶6),x(19∶24),x0) spline(x(1∶6),x(25∶30),x0);spline(x(1∶6),x(31∶36),x0);spline(x(1∶6),x(37∶42),x0)]′运算数据分析:(1)这组运算数据壹致通过已知节点,而且偏差较小、数学处理和程序设计均大大简化(和文献[1]相比).(2)经过之上的运算,能够使原来仅有的6组数据变为26组,而且仍能够根据工程人员的需要对上述程序步长进行修改,就可任意提高精度,从而使工程人员能够更好地了解各种参数于各点的数据,使工程精度大大提高5插值曲线拟合当然,无论以多么小的数为步长、无论给出多少组数据,这些参数仍是壹些离散的数据,于有些情况下,工程人员要了解某些数据随某壹主要参数的变化而变化的连续曲线,这时,能够于数据插值的基础上,发挥MATLAB于图形处理上的强大功能,对之上插值所得的数据进行曲线拟合,以便更好地了解各参数随某壹主要参数变化而变化的趋势.于之上插值数据的基础上,于上面程序的尾部编写MATLAB作图程序,作图程序如下,运行后得到图1所示插值拟合曲线.plot(x0,out(27∶52),‘-’) holdonplot(x0,out(53∶78),‘-+’)plot(x0,out(79∶104),‘∶’)plot(x0,out(105∶130),‘--’)plot(x0,out(131∶156),‘-3’)plot(x0,out(157∶182),‘-.’)legend(‘上游半径’,‘下游半径’,‘拱厚’,‘半中心角’,‘圆心距’,‘淤沙高程’)holdoffgridon从图1中,能够见到各参数随高程的变化而变化的曲线,从而更好地去了解各参数的变化规律,实现对工程各参数的整体把握,这是壹般数值处理方法所无法实现的.6小结之上仅为三次样条插值及其实现方法的壹个实例,本文于插值方法的选择上选取了能够得到平滑曲线的、具有优良数学特征的三次样条插值法;于插值的实现上选取了具有强大计算功能的数学软件MATLAB,它能够以较少的编码,较简单的语句实现这壹复杂的计算,且能得到较合理的结论;于曲线的拟合上我们于插值的基础上同样选取具有强大图形处理功能的MATLAB软件,从而形成较准确、较平滑、较合实际的曲线.总之,之上所提供的方法是三次样条插值和MATLAB科学计算语言于工程中应用的壹个实例,它能使计算较简便,又能很好地满足光滑性要求,使曲线也不失真.实现了工程数学、计算数学、程序设计的结合和简化.三次样条插值不仅于工程方面,而且于测绘、勘察、预测等方面均有着十分广泛的应用参考文献:[1]彭荣利,靳萍,欧阳建国.工程列表函数的数值拟合和应用[J].武汉大学学报(工学版),2002,35(4):42～45.[2]王瑞华.水利工程数据插值计算及图形处理[J].农田水利和小水电,1994(8):15～19.[3]鞠时光,郭伟刚.实用三次样条插值函数[J].小型微型计算机系统,1992,13(9):20～23.[4]谢赤,钟钻.插值法于零息收益曲线构造中的实证研究[J].数量经济技术经济研究,2002(4):31～34.[5]曾绍标,韩秀芹.工程数学基础[M].:科学出版社,2001.[6]王沫然.MATLAB和科学计算(第2版)[M].:电子工业出版社,2003。
三次样条插值
三次样条插值0 引⾔三次样条插值以构造简单，使⽤⽅便，拟合准确，具有“保凸”的重要性质等特点成为了常⽤的插值⽅法。
⼀般三次样条插值解算过程中通过追赶法求解三弯矩阵，但使⽤计算机求解时会表现出解的精度不⾼的问题，导致其计算结果⽆法应⽤到⼯程实践之中。
因此需要找出⼀种提⾼解精度的⽅法。
1 基本概念三次样条函数的定义：在区间内对于给定的函数值，其中，如果函数满⾜条件：（1）在每个⼦区间，上都是不⾼于三次的多项式；（2）、、在上都连续；（3），。
则称为函数关于节点的三次样条函数。
想要求解三次样条插值函数，只需在每个⼦区间上确定⼀个三次多项式共有4个系数，确定它们需要 4n 个条件，因此要完全确定共需 4n 个条件。
由所满⾜的条件（1）、（2）、（3），可确定个条件，仍然缺少两个条件。
这两个条件通常由实际问题对三次样条插值函数在端点的状态要求给出，也称之为边界条件，常见的边界条件有：1）夹持边界条件（Clamped Spline）：给定两端点的⼀阶导数值，即，；2）⾃然边界条件（Natural Spline）：使两端点的⼆阶导数值为零，即；3）⾮扭结边界条件（Not-A-Knot Spline）：强制第⼀个插值点的三阶导数值等于第⼆个点的三阶导数值，最后第⼀个点的三阶导数值等于最后第⼆个点的三阶导数值，即，。
2 计算⽅法设三次样条函数，（0），，，由三次样条函数定义（1）（2）（3）可得：，（1）如下构造式（1）矩阵：（2）由式（1）可知：，，，，（3）1）在夹持边界条件时，，，，；，，，；2）在⾃然边界条件时，，，，；，，，；3）在⾮扭结边界条件时，，，，；，，，；由n个未知数的⾮齐次⽅程组有惟⼀解的充分必要条件是，可知矩阵⽅程（2）在以上三种情况下都有惟⼀解。
对矩阵⽅程（2）采⽤⾼斯列主元消去法即可求解得出。
最后，代⼊式（0）可以得出：，，，，3 应⽤算例有点集，在⾮扭结边界条件下进⾏插值。
同时使⽤Matlab R2010a和⽂章所述⽅法进⾏插值计算，对⽐计算结果。
三次样条插值方法在太湖水质评价中的应用
三次样条插值方法在太湖水质评价中的应用毛新伟;陆铭锋;高琦;林其彬【摘要】The cubic spline interpolation method was used to linearly interpolate the indexes of water quality,eutrophication, and algae density in the Taihu Lake. The spatial variation map of the concentration of main water quality indexes in the Taihu Lake was drawn, which directly reflected the overall situation of water quality, and well compensated for the Lack of monitoring sites. Although the evaluation results had some deviations in certain areas where water quality changed greatly, the evaluation method can provide a reference for water quality evaluation of large lakes.%采用三次样条播值方法对全太湖代表性水质指标和反映富营养化、藻类密度的指标进行线性插值,绘制太湖主要水质指标浓度空间变化图,直观反映了太湖水质总体状况和空间分布特征,较好地解决了监测站点密度不够的问题.虽然评价结果在水质变化较大的局部湖区存在一定误差,但该方法对于大面积湖泊水质评价工作具有借鉴价值.【期刊名称】《水资源保护》【年(卷),期】2011(027)004【总页数】5页(P58-61,78)【关键词】太湖;水质评价;水质指标;三次样条插值【作者】毛新伟;陆铭锋;高琦;林其彬【作者单位】太湖流域管理局水文水资源监测局,江苏,无锡,214024;太湖流域管理局水文水资源监测局,江苏,无锡,214024;太湖流域管理局水文水资源监测局,江苏,无锡,214024;无锡滨湖都市建设投资有限公司,江苏,无锡,214000【正文语种】中文【中图分类】X8241 问题的提出太湖是我国第三大淡水湖泊,水域广阔,水面面积2338km2。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
样 条 插 值 函 数 进 行 入 库 洪水 的调 洪 演 算 。在 分 析 三次 样 条 插 值 函数 与 调 洪 演 算 结 合 的 基 础 上 ，通 过 ｄ ｅ ｌ ｐ ｈ ｉ 平 台编 写
２ ． Ｇｕ ａ ｎ ｇ ｘ ｉ Ｎａ ｎ ｎ ｉ ｎ ｇ Ｄｅ ｓ ｉ ｇ ｎ Ｉ ｎ ｓ ｔ ｉ ｔ ｕ ｔ ｅ ｏ ｆ Ｗａ ｔ ｅ ｒ Ｃ ｏ ｎ ｓ ｅ ｒ ｖ ａ ｎ ｃ ｙ＆ Ｈｙ ｄ ｒ ｏ ｐ ｏ ｗｅ ｒ ， Ｎａ ｎ ｎ ｉ ｎ ｇ ５ ３ ０ ０ ０ ４ ， Ｇｕ ａ ｎ ｇ ｘ ｉ ． Ｃ ｈ ｉ ｎ ａ ；
（ １ ． Ｃｏ ｌ ｌ ｅ ｇ ｅ ｏ ｆ Ｃ ｉ ｖ ｉ ｌ Ａｒ ｃ ｈ ｉ ｔ ｅ ｃ ｔ ｕ ｒ ａ ｌ Ｅ ｎ ｇ ｉ ｎ ｅ ｅ ｉｎ ｒ ｇ ， Ｇｕ ａ ｎ ｇ ｘ ｉ Ｕｎ ｉ ｖ ｅ ｒ ｓ ｉ ｔ ｙ ， Ｎａ ｎ ｎ ｉ ｎ ｇ ５ ３ ０ ０ ０ ４ ， Ｇｕ ａ ｎ ｇ ｘ ｉ ， Ｃ ｈ ｉ ｎ ａ ；
Ｋｅ ｙ Ｗｏ ｒ ｄ ｓ ： ｃ ｕ ｂ ｉ ｃ ｓ ｐ ｌ ｉ ｎ ｅ ｉ ｎ ｔ ｅ ｒ ｐ ｏ ｌ ａ ｔ ｉ ｏ ｎ ； ｌ ｆ ｏ ｏ ｄ ｒ ｏ ｕ ｔ ｉ ｎ ｇ ； ａ ｐ ｐ ｌ ｉ ｃ ａ ｔ ｉ ｏ ｎ ｓ ｔ ｕ ｄ ｙ
水 力 发 电
ቤተ መጻሕፍቲ ባይዱ
第３ ９ 卷第 １ ０ 期
２ ０ １ ３ 年１ ０ 月
三 － 灭样 条 搐 值 函 数 在 调 洪 ， ， 演 ， 篡 甲 的 应 用
莫 崇 勋 ，魏 炜 ，夏 秉龙 ，王 大彬 。 ，孙桂 凯 ，刘 俐 ，姜 庆玲
（ １ ． 广 西 大 学土 木 建 筑 工 程 学 院 ，广 西 南 宁 ５ ３ ０ ０ ０ ４ ；２ ． 广 西 南 宁水 利 电力 设 计 院 ， 广 西 南宁 ５ ３ ０ ０ ０ ４ ；３ ． 山东 垦 利 水 利 局 ， 山东 东 营 ２ ５ ７ ０ ０ ０ ）
３ ． Ｋｅ ｎ ｌ ｉ Ａ ｄ ｍｉ ｎ ｉ ｓ ｔ ｒ ａ ｔ ｉ ｏ ｎ ｏ ｆ Ｗａ ｔ ｅ ｒ Ｒｅ ｓ ｏ ｕ ｒ ｃ ｅ ｓ ， Ｄｏ ｎ ｇ ｙ ｉ ｎ ｇ ２ ５ ７ ０ ０ ０ ， Ｓ ｈ ａ ｎ ｄ ｏ ｎ ，Ｃｈ ｇ ｉ ｎ ａ ）
Ａｂ ｓ ｔ ｒ ａ ｃ ｔ Ｉ ｎ ｏ ｒ ｄ ｅ ｒ ｔ ｏ ｂ ｅ ｔ ｔ ｅ ｒ ｉ ｆ ｔ ｔ ｈ ｅ ｃ ｕ ｖｅ ｒ ｓ ｏ ｆ ｗａ ｔ ｅ ｒ ｌ ｅ ｖ ｅ ｌ － ｓ ｔ ｏ ｒ ａ ｇ ｅ ａ ｎ ｄ ｗ ａ ｔ ｅ ｒ ｌ ｅ ｖ ｅ ｌ － ｄ ｉ ｓ ｃ ｈ ａ ｒ ｇ ｅ ｃ ａ ｐ ａ ｃ ｉ ｔ ｙ ｉ ｎ ｌ ｆ ｏ ｏ ｄ ｒ ｏ ｕ ｔ ｉ ｎ ｇ ， ａ ｎ ｄ ｉ ｍｐ ｒ ｏ ｖ ｅ ｔ ｈ ｅ ａ ｃ ｃ ｕ ｒ ａ ｃ ｙ ｏ ｆ ｌ ｆ ｏ ｏ ｄ ｒ ｏ ｕ ｔ ｉ ｎ ｇ ， ｔ ｈ ｅ ａ ｐ ｐ ｌ ｉ ｃ ａ ｔ ｉ ｏ ｎ ｏ ｆ ｃ ｕ ｂ ｉ ｃ ｓ ｐ ｌ ｉ ｎ ｅ ｉ ｎ ｔ ｅ ｒ ｐ ｏ ｌ ａ ｔ ｉ ｏ ｎ ｆ ｕ ｎ ｃ ｔ ｉ ｏ ｎ ｉ ｎ ｌ ｆ ｏ ｏ ｄ ｒ ｏ ｕ ｔ ｉ ｎ ｇ ｏ ｆ ｉ ｎ ｌ ｆ ｏ ｗ ｌｏ ｆ ｏ ｄ ｉ ｓ ｓ ｔ ｕ ｄ ｉ ｅ ｄ ｈ ｅ ｒ ｅ ｉ ｎ ． Ｂ ａ ｓ ｅ ｄ ｏ ｎ ｔ ｈ ｅ ａ ｎ ａ ｌ ｙ ｓ ｉ ｓ ｏ ｆ ｔ ｈ ｅ ｃ ｏ ｍｂ ｉ ｎ ａ ｔ ｉ ｏ ｎ ｏ ｆ ｃ ｕ ｂ ｉ ｃ ｓ ｐ ｌ ｉ ｎ ｅ ｉ ｎ ｔ ｅ ｐｏ ｒ ｌ ａ ｔ ｉ ｏ ｎ ｆ ｕ ｎ ｃ ｔ ｉ ｏ ｎ ａ ｎ ｄ ｌ ｆ ｏ ｏ ｄ ｒ ｏ ｕ ｔ ｉ ｎ ，ｔ ｇ ｈ ｅ ｌ ｆ ｏ ｏ ｄ ｒ ｏ ｕ ｔ ｉ ｎ ｇ ｐ ｒ ｏ ｃ ｅ ｄ ｕ ｒ ｅ ｓ ａ ｒ ｅ
了 调 洪 演算 程 序 。实 例计 算 结 果 证 明 ，三 次样 条插 值 函数 在 调 洪 演算 中是 可 行 的 。
关 键 词 ：三 次样 条 插 值 ；调洪 演 算 ；应 用 研究
Ａｐｐ ｌ ｉ ｃ ａ ｔ ｉ ｏ ｎ ｏ ｆ Ｃｕ ｂ ｉ ｃ Ｓ ｐ ｌ ｉ ｎ ｅ Ｉ ｎ ｔ ｅ ｒ ｐ ｏ ｌ ａ ｔ ｉ ｏ ｎ Ｆｕ ｎ ｃ ｔ ｉ ｏ ｎ ｉ ｎ Ｆ ｌ ｏ ｏ ｄ Ｒｏ ｕ ｉ ｔ ｎ ｇ
ｗ ｉｔ ｒ ｔ ｅ ｎ ｂ ｙ Ｄｅ ｌ ｐ ｈ ｉ ｐ ｌ ａ ｆ ｔ ｏ ｒ ｍ． Ｔ ｈ ｅ ｃ ａ ｓ ｅ ｃ ａ ｌ ｃ ｕ ｌ ａ ｔ ｉ ｏ ｎ ｓ ｈ ｏ ｗｓ ｔ ｈ ａ ｔ ｔ ｈ ｅ ｃ ｕ ｂ ｉ ｃ ｓ ｐ ｌ ｉ ｎ ｅ ｉ ｎ ｔ ｅ ｒ ｐ ｏ ｌ ａ ｔ ｉ ｏ ｎ ｆ ｕ ｎ ｃ ｔ ｉ ｏ ｎ ｃ ａ ｎ ｂ ｅ ｕ ｓ ｅ ｄ ｉ ｎ ｌ ｆ ｏ ｏ ｄ ｒ ｏ ｕ ｔ ｉ ｎ ｇ ．
ＭＯ Ｃ ｈ ｏ ｎ ｇ ｘ ｕ ｎ ， ＷＥ Ｉ Ｗｅ ｉ ， Ｘ Ｉ Ａ Ｂ ｉ ｎ ｇ ｌ ｏ ｎ ｇ ， ＷＡ Ｎ Ｇ Ｄ ａ ｂ ｉ ｎ 。 ， Ｓ Ｕ Ｎ Ｇ ｕ ｉ ｋ ａ ｉ ， Ｌ Ｉ Ｕ Ｌ ｉ ‘ ， Ｊ Ｉ Ａ Ｎ Ｇ Ｑ ｉ ｎ ｇ ｌ ｉ ｎ ｇ