Hermite插值的上机实现及应用课程设计

学校代码： 10128学 号： 20122090课程设计说明书题 目：Hermite插值地上机实现及应用学生姓名：学 院：理学院班 级：指导教师：任文秀 曹艳2015年1月16日目录摘要 0第一章Hermite插值地上机实现 0§1.1插值概述 1§1.1.1插值问题地提出 1§1.1.2插值地种类 1§1.2Hermite插值地问题 3§1.2.1Hermite插值地几种形式 3§1.2.2Hermite插值地几个重要定理 10§1.2.3Hermite插值地优点 11§1.3Hermite插值地源程序 11§1.3.1三次Hermite插值地C程序 11§1.3.2二重Hermite插值地matlab程序 12第二章Hermite插值地应用 12§2.1Hermite插值函数地工程应用 12§2.2应用Hermite插值作心电图基线漂移校正 15参考文献 22附录A三次Hermite插值地C程序 23附录B二重Hermite插值地MATLAB程序 26摘要随着计算机技术地普及和应用地日益广泛 ,细分方法在近年来已经成为了计算机辅助设计 (CAD)和计算机图形学 (CG)领域内地一个国际性研究热点 .通过近三十年地发展 ,细分方法日趋完善 ,多数经典地细分方法已经建立起了较为系统地理论知识体系 .1992年Merrien首次提出了 Hermite型地插值细分格式 ,随后Hermite插值型细分方法得到了迅速地发展 ,从一维区间上生成 C1、C2细分曲线地格式到维矩形网格上生成光滑曲面地格式得以在短时间内展现 ,但是对于二维矩形上生成地光滑曲面在直观上与采样函数有不小地差距 .在构造插值时，对所构造地插值，不仅要求差值多项式节点地函数值与被插函数地函数值相同，还要求在节点处地插值函数与被插函数地一阶导数地值也相等对所构造地插值，不仅要求差值多项式节点地函数值与被插函数地函数值相同，还要求在节点处地插值函数与被插函数地一阶导数地值也相等 .关键词Hermite插值；拉格朗日插值； Newton插值；余项；Hermite插值应用第一章Hermite插值地上机实现(n1)!§1.1插值概述§1.1.1插值问题地提出在许多实际问题及科学研究中，因素之间往往存在着函数关系，但这些关系地表达式不一定都知道，通常只是由观察或测试得到一些离散数值，所以只能从这些数据构造函数地近似表达式.有时，虽然给出了解读表达式，不过由于解读表达式过于复杂，使用或计算起来十分麻烦 .这就需要建立某种近似表达，因此引入插值 .§1.1.2插值地种类类型1拉格朗日插值 .定义

1.1若函数

y=f(x)

在若干点

xi地函数值

yi=

fxi

（i=0,1,

,n），则另一个函数

pn(x)：p(

xi)=

yi

,i=0,1,

,n,则称

p(x)为

f(x)

地插值函数，而

f(x)

为被插值函数

.对于

x

[a,b]

，且x xi，用

Pn（x)地值作为

f(x)地近似值或估计值，常称内插法

.对于

x

[a,b]

，用

Pn

(x)地值去估计f(x)地值，又称外插法 .注解1.1拉格朗日插值分为线性插值 L1(x)和n次插值Ln(x).注解1.2拉格朗日插值地余项为Rn(x) f(n1)()W(n1)(x)类型2Newton插值定义1.2任何一个不高于n次多项式，都可以表示成函数1,xx0,(xx0)(xx1),,(xx0)(xx1)(xxn1)地线性组合.既可以把满足插值条件P(xi)yi(i0,1,,n)地n次插值多项式写成如下形式：a0a1(xx0)a2(xx0)(xx1)an(xx0)(xx1)(xxn1)其中，ak为待定系数，这种形式地插值多项式称为牛顿插值多项式，记为Nn(x).注解1.3设x0,x1,...,xn互不相同，则fx关于x0,x1,...,xn地n阶差商为：fx0,x1,...,xnfx1,x2,...,xnfx0,x1,...,xn1.xnx0则一阶差商为fx1fx0.fx0,x1x1x0且二阶差商为fx0,x1,x2fx1,x2fx0,x1.x2x0总结以上可得如下表1-1.表1-1差商表xi f(xi)x0fxoxfx11x2fx2x3fx3

一阶差商 二阶差商 三阶差商 n阶差商f[x0,x1]f[x1,x2] f[x0,x1,x2]f[x2,x3] f[x1,x2,x3] f[x0,x1,x2,x3]xn fxn fxn1,xn fxn2,xn1,xn fxn3,...,xn fx0,x1,...,xn类型3分段插值定义1.3对给定区间 a,b做划分a x0 x1 ... xn b在每个小区间 [xi,xi1]上作fx以xi,xi1为节点地线性插值，记这个插值px pix，pixxxi1fxixxifxi1，(xixxi1)xixi1xi1xi把每一个区间地线性插值函数连接起来，得到fx地以ax0x1...xnb为剖分节点地分段性函数 px.注解1.4分段插值地基本思想将被插值函数 f x地插值节点由小到大排序，然后在每对相邻地两个节点为端点地区间上用n次多项式去近似

f x

.类型4Hermite插值定义 1.4Hermite插值是利用未知函数

f x在插值节点上地函数值及导数值来构造插值多项式；分为带导数地插值与不带导数地插值二类 .类型5三次样条插值样条插值是一种改进地分段插值 .定义1.5函数Sx a,b，且在每个小区间 xi,xi1上是三次多项式，其中a x0

x1

...

xn

b是给定节点，则称

Sx

是节点

x0,x1,...,xn上地三次样条函数

.若在节点xi上给定函数值yi fxi,i 0,1,..,n，并且Sxi yi,i 0,1,..,n，则称Sx为三次样条插值函数 .注解1.5本文着重介绍 Hermite插值§1.2Hermite插值地问题§1.2.1Hermite插值地几种形式类型一Hermite插值地一般形式求一个次数不大于 n+r+1地代数多项式 H(x),满足H(xi)f(xi),i=0,1,2,...,n.（1.1）H(xi) f(xi), i 0,1,2, ,r(r n)称以上地插值问题为 Hermite插值问题注解1.6Hermite插值多项式地推导（即建立 Hermite插值多项式地方法）令n rH(x) hk(x)f(xk)k 0

hk(x)f(xk)（1.2）k 0其中hk(x)(k 0,1, ,n)和hk(x)(k 0,1, ,r)都是n r 1次待定多项式，并且它们满足以下条件：hk(xi1iki,k0,1,,n)ik0（1.3）hk(xi)0,k0,1,,n;i0,1,,rhk(xi1iki,k0,1,,r)ik0（1.4）hk(xi)0,k0,1,,r;i0,1,,n显然满足条件式（1.3）,（1.4）地多项式（1.2）地次数不大于nr1次，且满足插值条件式1.1）.形式一求解hk(x)（不带导数地Hermite插值）由条件式（1.3）知xi(i 0,1, ,r;i k)是hk(x)地二重零点.且由条件式（ 1.3）知xi(i r 1,r 2, ,n;i k)是hk(x)地零点.当0kr时hk(x)具有如下形式：hk(x)2x222xr1(xxn)(AxB)(xx0)xk1xxk1(xxr)xr(xxi)2n(AxB)xxi（1.5）i0ir1ik其中，A,B是待定系数.由条件式（1.3）知hk(xk) 1,hk(xk) 0即rxi)2n(AxkB)(xk(xkxi)1i0ir1ikrxi)2nrrxi)2nA(xk(xkxi)2(AxkB)(xkxj)(xk(xkxi)i0ir1j0i0ir1ikijiknrxi)2n(AxkB)(xk(xkxi)0jr1i0ir1ikij由上述两式解得r1n12xjjr1xkxjAj0xkrxi)2n.(xk(xkxi)i0ir1ikB

1Axk.rn(xkxi)2(xkxi)i0ir1ik将A,B代入式（1.5），得hk(x) {1 (x xk)[lkn(xk) lkr(xk)]}lkn(x)lkr(x)（1.6）k 0,1, ,r其中，lkn(x)nxxi.i0xkxiikrxxi.lkr(x)i0xkxiikn1lkn(xk)xk.i0xiikr1lkr(xk)xk.i0xiik当r 1 k n时，hk(x)具有如下形式rxi)2nhk(x)C(x(xxi).（1.7）i0ir1ik由条件式（1.3）知hk(xk) 1C1.rn(xkxi)2(xkxi)i0ir1k将C代入式（1.7），得h(x)wr(x)l(x),kr1,r2,,n（1.8）kw(x)knrk其中，rwr(x)(xxi).i0rwr(xk)(xkxi).i0lkn(x)nxxi.i0xkxiik综合式（1.1）、（1.2）可以得到hk(x)(k0,1,n)，即式（1.6）、（1.8）形式二求解hk(x)(即带导数地Hermite插值）由条件式（1.4）知xi(i0,1,,r;ik)是hk(x)地二重零点，且由条件式（1.4）知xi(ik,r1,r2,,n)是hk(x)地零点.当0kr时，hk(x)具有如下形式：nrhk(x)D(xxi)(xxi)（1.9）i0i0由条件式（1.4）知hk(xk)1ikD1nnrrnr(xkxi)(xkxi)(xkxi)(xkxi)j0i0i0j0i0i0ijikjkikij将D代入式（1.9），得hk(x)(xxk)lkn(x)lkr(x),k0,1,,r.(1.10)其中，lkn(x)nxxi.i0xkxiikrxxi.lkr(x)i0xkxiik由式（1.2）,（1.6）,（1.8）,（1.10）所表示地多项式称为Hermite插值多项式，其中由式（1.6）,（1.8），（1.10）所表示地多项式称为Hermite插值基函数.Hermite插值多项式地余项为R2n1(x)=f(2n1)()W（2n1）(x).（2n2)!类型二二重Hermite插值多项式一般地Hermite插值为m=2地情况，即给定地插值节点nxi i 0均为二重节点，更具体些f(x) C2 a,b ，n及插值节点 xi i 0，若有H2n1(x) P2n1满足H2n1(xi)f(xi).H2'n1(xi)f'(xi),i0,1,n.就称H2n1(x)为fx关于节点xin地二重Hermite插值多项式.i0类型三三次Hermite插值设yfx是区间[a,b]上地实函数,x0,x1是[a,b]上相异两点,且xx,01yfx在xi上地函数值和一阶导数值分别为yfxi0,1和mfxii0,1,求iii三次多项式H3x,使其满足：H3(xi)yiH3'(xi)(i0,1).miH3(x)称为三次埃尔M特插值多项式.注解1.7误差估计定理1.1设f(x)在包含x0、x1地区间[a,b]内存在四阶导数，则当R31f(4)22(x)f(x)H3(x)()(xx0)（xx1）（4!设M4maxf(4)(x)x0xx1则当

x∈[a,b]时有余项(a,b)且与x有关）xx0,x1时,余项有如下估计式R3(x)M4h4.384类型四两点三次Hermite插值设f(x)在节点x0、x处地函数值为y0、y，在节点x0、x处地一阶导数值为y0'、y1'，两111个节点最高可以用3次Hermite多项式H3(x)，作为插值函数H3(x)应满足插值条件H3(x0)y0H3(x1)y1.H3(x0)y0H3(x1)y1.H3(x)应用四个插值基函数表示，设H3(x)地插值基函数为hi(x)0,1,2,3，H3(x)a0h0(x)a1h1(x)a2h2(x)a3h3(x)如果希望插值系数与 Lagrange插值一样简单，那么重新假设H3(x) y00(x) y11(x) y00(x) y11(x)H3(x) y00(x) y11(x) y00(x) y11(x)其中0(x0) 1 0(x1) 0 0(x0) 0 0(x1) 01(x0) 0 1(x1) 1 1(x0) 0 1(x1) 00(x0) 0 0(x1) 0 0(x0) 1 0(x1) 01(x0) 0 1(x1) 0 1(x0) 0 1(x1) 1可知x1是 0(x)地二重零点，即可假设0(x)(xx)2(axb).1由0(x0) 1 0(x0) 0可得a23.x1)(x0b12x0x1)2(x0x1)3(x0a0(x)(xx1)(axb)(xx1)2{2x3122x03}(x0x1)(x0x1)(x0x1)(xx1)212x02x(x0x1)2x0x1x0x1(12xx0xx1)2x1x0)(x1x0(1 2l1(x))l02(x)......Lagrange插值基函数如下式所示20(x)(12l1(x))l02(x)12xx0xx1x1x0x0x1类似可得1(x)(12l0(x))l12(x)12xx1x0x120(x)(xx0)l02(x)xx0xx1x0x121(x)(xx1)l12(x)xx1xx0x1x0将以上结果代入H3(x)y00(x)y11(x)y00(x)y11(x)得两个节点地三次Hermite插值公式H3(x)y00(x)y11(x)y00(x)y11(x)y0(12l1(x))l02(x)y1(12l0(x))l12(x)y0(xx0)l02(x)y1(xx1)l12(x)xx0xx12xx1xx122y0y1y0y1xx012x0x112xx0x1xx1x0x1x0x0x1x0x1注解1.8二点三次 Hermite插值地余项R(x)=f4()(xx0)2(xx1)2x0x134!§1.2.2Hermite插值地几个重要定理定理1.2误差定理若fC2n2(a,b)，则fx关于a,b上节点{xi}n地二重Hermite插值多项式误差为R2n1(x)f(x)H2n1(x)f(2n2)()wn2(x)(2n2)!定理1.3唯一性定理Hermite插值问题地表达式H(xi)f(xi),i0,1,2,,n.H(xi)f(xi),i0,1,2,,r(rn).地解H(x)存在而且唯一.定理1.4Hermite插值余项定理Hermite插值公式地余项为f(x)H(x)f(nr2)()wn(x)wr(x).(nr2)!其中,是插值区间a,b内地某一点.§1.2.3Hermite插值地优点分段线性插值地算法简单，计算量小，然而从整体上看，逼近函数不够光滑，在节点处，逼近函数地左右导数不相等 .Hermite插值地逼近函数与被逼近函数不仅在插值节点上取相同地函数值，而且逼近函数与被逼近函数在插值节点上去相同地若干阶导数值 .Hermite插值法结合了函数地导数值，使得插值地精度更为提高.Hermite插值具有少节点得到高次插值多项式地特点.Hermite插值插值多项式灵活多样.Hermite插值在节点一定地条件下，可以多种构造插值条件.§1.3Hermite插值地源程序§1.3.1三次Hermite插值地C程序例1.1已知函数y=1/(1+x2)在区间[0,3]上取等距插值节点，求区间[0,3]上地分段三次埃尔M特插值函数，并利用它求出f(1.5)地近似值（0.3075）表1-2例题数据表xi012yi10.50.2yi'0-0.5-0.16注解1.9本例题程序流程图及 C程序详见附录A.1.3.2二重Hermite插值地matlab程序注解1.10程序及程序演示详见附录 B第二章

Hermite

插值地应用§2.1Hermite插值函数地工程应用对于同一个问题运用不同地方法或许都能得到相同地结果，但是每一个方法都有其得天独厚地优势以及劣势 .特别是在现在这个现代化地信息时代，计算已经变得越来越重要，对计算结果地要求也十分苛刻 .插值方法在实际问题中有着广泛地应用它能使一个有着大量数据地问题变得简单明了、易于观察，因此，地位自然不喻 .Hermite插值为使插值函数能更好地和原来地函数重合，不但要求二者在节点上函数值相等，而且还要求相切，对应地导数值也相等，甚至要求高阶导数也相等 ,凭借其精度高，计算严谨被大多数人应用了起来 .算例分析在土方工程中，土地最大干密度与最优含水量是确保路基压实质量地两个关键指标 ,利用埃尔M特插值函数求得地干密度、含水量能更好地逼近实验得到地 pd一 曲线，求解精度较高

.通过绘制干密度与含水量地相关曲线，即

pd

一

曲线，求得最大干密度与最优含水量地方法为图解法.图解法因简便直观而在实际工作中被广泛采用，但图解法随意性大，易产生人为误差.目前，数解法主要有两类：一是利用曲线拟合法求解，二是利用代数插值求解 .用上述方法分别对实验地工程实例进行了求解，发现所得结果 地差值较大，其中最大干密度差值达 0.01～0.06g/cm3，最佳含水量差值达0.5％～1.4％.在本研究中利用埃尔M特插值问题，试图更加精确地求解最大干密度与最优含水量.例2.1某公路工程路基填七地一组室内标准击实实验结果见表2-1，由表2-1可知，其最火干密度应在含水量11.6％附近.表2-1室内标准击实实验结果实验序号12345含水量%5.87.411.615.517.6干密度pd1.771.801.851.821.78gcm3)根据图解法将最大干密度定为1.85g/cm3对应地最优含水最为l1.6％而根据pd曲线图，,,.最优含水量在12％附近更为恰当.下面利用埃尔M特插值函数求解最大干密度与最优含水量.取0,1,2分别为7.4、l1.6、15.5，对应地f(i)分别为1.80、1.85、1.82，得到f[0,1]=0．01l905,f[0,1,2]=-0.0024194.步骤一建立干密度、含水量地埃尔M特插值函数.利用式f[0,1,2]A(0)(1)(2)建立干密度、含水量地埃尔M特插值函数为H()=1.8+0.0l1905(-7.4)-0.0024194(-7.4)(-11.6)+A(-7.4)(-11.6)(-15.5)，利用式子Af(1)f[0,1](10)f[0,1,2].(10)(12)可得A=-0.06105f'(1)+0.00010644.步骤二求解最大干密度与最优含水量.取3=5.8％，对应地pd3f(3)=1.77g/cm3,根据插值条件H(3)f(3),代入式A=-0.06105f'( 1)+0.00010644,令H'( ) 0,得2 10.379 24.168 0,解此方程得最优含水量为 12.3％.得最大干密度为 1．85g／cm3.步骤三误差分析：由表 2-1中实验数据可得f[ 0,1,3, 4] 6.184105，由f(4)()(0)(1)2(2)R()4!和f(4)()1,2,3]f[,4!可得R()f[0,1,2,3](0)(1)2(2)=6.184105(12.37.4)(12.311.6)2(12.315.5)=4.751 105,根据误差分析可知，此法求解最大干密度与最优含水量地精度较高，能更好地逼近实验中得到地pd曲线.模糊矩阵综合评价得：DWR[0.28820.22420.07860.13320.18540.0959]T100000000100.480.5200=001000000.380.611111=[0.3787,0.1336,0.1368,0.2996,0.4313]以上计算结果表明，

I级水地隶属度为

0.3787，II

级水地隶属度

0.1336，III

级水地隶属度为0.1368，IV

级水地隶属度为

0.2996，V

级水地隶属度为

0.4313.由于

V级水地隶属度最大，因此鉴湖水体综合评价等级应为

V级.总结应用模糊数学原理综合评判鉴湖水质等级，比采用单因子极值评价更为合理 .评判结果表明，鉴湖所在地区由于经济社会地快速发展，已经造成了严重地水体污染，因此水质等级很快由Ⅲ类变成V类.水体地污染引起地一系列问题应该引起足够地重视，如果这样发展下去，鉴湖将失去它原来地价值，因而政府应该采取措施，防止和减轻水污染，努力提高鉴湖水质等级，从而使之能发挥更好地作用 .§2.2应用Hermite插值作心电图基线漂移校正消除心电图地基线漂移是个重要向题.采用分段三次Hermite插值来作基线漂移校正,提出了当心率变化引起插值区间信号长度变化时,插值墓函数地线性变换规则.由此可以保持拟合地高精度,又减少计算量.有可能用于实时心电监护.如果监护仪中地CPU能力有限,本文还提出了一种计算Hermite插值函数硬件电路,使每一点地计算时间缩短为12微秒.心电图(ECG)信号地计算机处理历来国内外十分重视.国内外其临床应用主要分为二大类：一是ECG计算机辅助诊断,主要用于医院地心电分析中心,常为离线分析,使用地计算机也多为中小型机甚至大型机；二是作ECG实时监护,主要用于临床危重病人、手术病人地监护,强调实时性要求,计算机多是由ECG等集成片构成,计算能力与存贮容量均受到限制.尽管ECG计算机分析已有二十多年地历史,国内外已做了大量地工作,但是仍然存在不少困难问题未予妥善解决.例如：消除ECG基线漂移是实时监护中地一个重要而又困难地问题.导致ECG基线漂移地主要因素有：电极地极化电位地变化,心电放大器地直流偏置漂移,人体由于呼吸或其它肌肉、体位地缓慢移动等.尽管可以努力消除产生基线漂移地原因例如努力使病人静卧不动，改善电极材料与导电膏地性能，改善心电放大器地特性等，但基线漂移仍然是不可避免地，因而会造成诊断疾病地困难

.消除基线漂移地困难在于基飘地频率很低，其范围为

0.05Hz

至

1Hz，主要分量在

0.1Hz

左右，如图

2.1所示，而

ST段地频率成分也很低，其最大分量在

0.6Hz-0.7Hz

左右，它们地频谱非常接近

.所以若使用高通频率滤波地方法以消除基飘，即使采用线性相位滤波器，仍会引起

ST段地严重失真，而

ST段在临床上有重要地价值

.图2.1基线漂移与 ST段地频谱目前解决基线漂移地方法，除了高通滤波外，常采用某种数学函数校正法，如分段直线校正，三次样条函数校正，二次函数校正及三次函数校正法 .在每个心电周期中选取 1-2个零电位点作为插值结点，俩结点之间地心电漂移，以消除基飘 .若采用直线进行逼近，是为直线校正法，这种方法计算量小，可实时实现，对慢变化地基线漂移效果尚好，对变化较快地基飘误差就严重.应用三次样条函数插值，可以获得较高地精度，本次报告就三次样条函数插值进行谈论 .今设二个相邻结点为 t0和t1，并已知这二个结点地函数值和一阶导数值为：y0,y0',y1,y1'，则三次Hermite插值函数为：s(t) s(t0)1(t) s'(t) 2(t) s(t1) 3(t) s'(t1) 4(t).满足下述条件：s'(t0) y0,s(t0) y0',s(t1) y1,s'(t1) y1'.上述式子中1(t)[2(tt0)(t1t0)](t1t)2(t1t0)32(t)(tt1)2(tt0)(t1t0)2t0tt1(2.1)(tt0)23(t)[(t1t)(t1t0)](t1t0)34(t)(tt0)2(tt1)(t1t0)2是为插值基函数 .由于实际心电信号地心率是不断随机变化地，所以不能按照等间隔计算，即(t1 t0)值将随心率地变化而变化地 .由于这个变化，将使得上述 4个插值基函数随之改变 .因而必须重新计算新地插值基函数，因此用一种简化插值基函数地计算方法，令样周期，k=0,1...m,t=kT,则可将式 2.1插值基函数写成离散形式：

t1t0，T为采T(mk)2(2km)1(k)m3(km)2(kT)2(k)m2k2(3m2k)3(k)m3(2.2)k2(km)T4(k)m2如若将m 'k k'm代入式（2.2），可得插值基函数为：(mmk')2(2mk'm)')m'm'1(km3(m' k')2(2k' m')/(m')3(m'k'm)2(m'k'T)'mm2(K)m2(k' m')2(k'T)(m)/(m')2m'm( '(k')2(3m'(m'')m4(k

' 2 m 'k)(3m 2 'k)m32k')/(m')3' 2 m 'k)( 'k m)Tm2(k')2(k'm')T(m')/(m')2(2.3)m比较式（2.2）与（2.3），可得1(k)1(k'),2(k)2(k')mm'3(k)3(k')4(k)4(k')mm'由此可见，当(t1t0)变化时，插值基函数1(k)、（3k)地幅值不变，只是时间轴发生线性变化.而2(k)、4(k)地幅值也将发生线性变化.因而可得变换公式如下.mk'k'm2(k')(m')2(k')m4(k')(m')4(k')m这样地变换可以使计算简化，图所示为各插值基函数随着，地变化而变化地图形 .图2.2当(t1 t0)变化时地插值基函数确定插值结点，即选择心电信号地零电位点 .一般常可选 TP段，为此可先估计 T波地终了点T1.根据临床心电图学 Bazett（巴泽特）公式：QT 0.39RR式中QT表示地是 Q波起点Q1至T1地间期.RR是二R波地间隔.若确定了Q1则由QT值可得T1点.所以可先检测 R波峰值 ，再往后定 H点.该H点约在T1之后30至70ms处，便可作插值结点图 2.3.这样可以吧连续二个心拍地 H点作为二个插值节点，进行三次 Hermite计算，然后作基线漂移校正 .具体步骤如下 .图2.3 确定插值结点Step1确定信号长度m.如下确定了几个典型地m值，如表2-2所示.表2-2几个典型地m值mH2'H1心率（次/分)压缩方式d修正系数m'm375401.00033545170.8932955070.7872556040.6802157030.5731758520.46713511510.360Step2

压缩方式是指由于

m压缩后，基函数地点数也需作相应减少

.Step3

计算插值地边界条件

.实验是用

8组不同心率地心电信号，迭加上不同频率地基线漂移（

0.1Hz,0.2Hz10.3Hz)来进行基漂校正

.图

2.4所示为其中一例，心率为

105次/分.迭加三种不同频率地正弦波作为基线漂移

.图中C1

为原始信号

,C2为由插值函数计算所得地基漂，

C3为经校正后地心电信号

.图2.4实时基漂校正结果.HR=105次/分.C1原始ECG，C2由插值拟合地基漂， C3校正后地ECG.而实际临床情况，心率一般均在 60次/分以上，基漂频率为 0.17Hz至0.33Hz之间，所以基漂校正地仿真结果误差都在 1.0%以下，可以满足临床要求 .ST段地计算也是令人满意地 .硬件电路实现虽然上述插值方法经过改进与简化，计算量已有很大减少，但对小型实时心电监护仪来说，CPU还可能不能承担

.因此又用专用硬件电路实现了上述地插值计算，并且还构成了一个

ST段检测仪

.

其框图如图

2.5所示

.

其中插值基函数电路是将

Hermite

插值基函数1(k),

2(k),

3(k),

4(k)，其中（

k=12m,m=375).计算并量化后写入

EPROM片.再在乘法控制线地控制下可依次读出

.插值条件寄存电路则由由

CPU送入地每段插值结点地边界条件

y0,y1,y0',y1',它们也可在乘法控制线地控制下依次读出

.这样每当由插值基函数电路端口读出函数值时，乘法控制线变回产生含有4个负脉冲地脉冲序列，乘法电路就依次产生4个对应地乘积y01(k),y0'2(k),y13(k)和y1'4(k),这四个乘积经累加电路累加后送至输出端口，完成一次基漂拟合值地计算 .由此连续运行 k=0至m，即可完成一个周期地基漂校正 .这个电路具有高速性能，插值基函数地确定、乘法运算、累加、翻转、技数、清零、重复等操作均是在乘法控制线地控制下同步进行地，有一部分操作室并行进行地，最大限度地减少了运算时间，提高了运算速度，可在12微秒内完成一个点地插值计算，时钟脉冲频率为 100MHz.且整个电路地成本也很低 .此监测仪对于心率在 40次/分以上，基线漂移频率在 0.4Hz以下地ECG基线漂移能相当好地进行校正.对ST段地检测，在一般情况下也能满足临床要求 .图2.5插值计算硬件电路框图参考文献文畅平.人民黄河.湖南：邵阳学院，2006.李庆阳，王能超，易大义.数值分析[M].北京：清华大学出版社，2008.白峰杉.数值计算引论.北京：高等教育出版社，2004.李庆阳.计算科学方法基础.北京：清华大学出版社，2006.冯康等.数值计算方法.北京：国防工业大学，1978.张雪敏.MATLAB基础及应用.北京：中国电力出版社，2009.附录A三次Hermite插值地C程序流程图开始输入xi,yi,xy=0,j=02.C程序代码#include<stdio.h>#include<math.h>floatf0(floatx){return((x-1)*(x-1)