CN113342881B - 核电厂仪控系统测试曲线构建方法、装置、设备、介质 - Google Patents

Abstract

本发明提供了一种核电厂仪控系统测试曲线构建方法、装置、设备及介质，方法包括：获取原始输入曲线，根据曲线时间长度t及原始采样周期T，计算原始采样点数N以及坐标；设定目标输入曲线的曲线参数；根据原始输入曲线的曲线斜率对原始输入曲线进行分段，得到m个曲线段；根据插值精度K对每个曲线段进行分段，得到多个曲线分段；对每个曲线分段进行K阶牛顿插值表达式的计算，得到各曲线分段的插值表达式；根据每个曲线分段的时间长度、目标采样点数以及插值表达式计算插值；将已完成插值的曲线分段重新拼接，获得目标曲线。本发明提高统计测试的完备性、准确性、有效性，为核电厂仪控系统统计测试提供技术支撑。

Description

核电厂仪控系统测试曲线构建方法、装置、设备、介质

技术领域

本发明涉及核电领域，具体而言，涉及一种核电厂仪控系统测试曲线构建方法、装置、设备、介质。

背景技术

核电厂仪控系统是保证核电厂安全运行和防止放射性泄露的重要系统，它根据电厂运行参数的变化情况产生与保护设备(或动作)有关的信号，通过驱动反应堆紧急停堆装置和专设安全设施动作，防止反应堆状态超过规定的安全极限或减轻超过安全极限所造成的损坏。仪控系统的可靠性对核电厂安全性和经济性有至关重要的影响，对仪控系统的可靠性测试及评估是确保核电厂安全运行的重要措施。

统计测试是系统可靠性评估的重要方法，通过仿真核电厂的工况运行曲线作为输入，测试仪控系统的保护功能。统计测试的核心是测试曲线的构建，测试曲线反映了核电厂运行时(尤其是事故发生时)仪控系统所面临的输入信号变化，测试曲线的构建直接决定了统计测试及可靠性评估结果的信心。测试曲线构建越科学，与真实工况越接近，测试结果越有效，可靠性评估信心越高。

现有的统计测试输入曲线示意如图1所示。同一时间轴上，多个参数将按照对应工况进行变化，真实工况下，该曲线变化是连续的，仪控系统通过传感器及IO板卡周期性采样。在统计测试时，也将通过模拟传感器及IO板卡周期性采样获得离散点逐个输入。

核电厂仪控系统的工作特点包括：

1)可靠性要求：系统的可靠运行决定了保护系统和设备的有效性，保证核电厂运行的安全性和经济性；

2)参数曲线变化：系统的输入参数值在某一工况下表现为连续变化的曲线，系统通过对参数的周期性采样获取参数值进行逻辑判断，当工况发生后，曲线会根据参数之间的相互影响产生剧烈且不规则的变化，进而出现很多变化率较大的峰值数据(峰值指曲线斜率从小于0到大于0或从大于0到小于0的快速变化)；

3)保护逻辑结构复杂，体现在：

a)常规阈值判断逻辑：根据参数的当前采样输入值与阈值进行比较判断保护动作输出；

b)趋势预测判断逻辑：根据参数变化趋势的预测值与阈值进行比较判断保护动作输出。预测值是根据参数的当前时刻采样输入值与前一时刻相邻采样输入值的差值，计算出该参数当前的变化趋势，并在趋势方向上生成一个预测值，差值越大(趋势越大)，在趋势方向上生成的预测值越大，即触发保护动作的提前量越大(典型保护功能：超温超功率保护、蒸汽压力高保护)。

实际值经过趋势预测算法后为预测值，预测值相比实际值的变化，当变化趋势较小时，预测值放大量相对较小，当变化趋势较大时，预测值放大量相对较大。由实际值构成的曲线称为实际曲线，由各实际值对应的预测值构成的曲线称为趋势预测曲线，趋势预测曲线示意如图2所示。

介于上述核电厂仪控系统参数曲线的复杂特点，在进行统计测试时，如何构建真实有效的测试曲线，是确保测试有效性，进而确保系统有效性的关键。因此，提高曲线构建真实度，是提高统计测试质量，提高可靠性评估信心，进而保证核电厂安全的重要方法。

当前工程中测试曲线仅通过核电厂工况仿真模型计算生成各参数的变化曲线。仿真模型当前的计算能力最高实现周期为50ms的采样精度，并未采用补偿措施来提高精度。然而因现有计算模型的精度不足，采用现有计算模型生成的曲线产生的失真较大。模型生成的测试曲线数据采样周期远大于仪控系统实际数据采样周期，使得数据失真，测试数据表现为阶跃式变化，相邻数据的差值远大于真实数据。主要体现在：

(1)、覆盖率低：对于常规阈值保护逻辑，数据采样点明显减少，对参数运行值的覆盖率降低，尤其对于阈值边界的覆盖和测试不充分，降低了测试的充分性和完备性；

(2)、趋势预测错误：对于趋势预测判断逻辑，采样周期较大的相邻输入值差值将远大于采样周期较小的相邻输入值差值，使得原本为连续变化的信号变为阶跃式信号，导致预测的变化趋势远大于实际趋势，得出错误的预测值，导致反应堆保护系统输出错误的保护动作，使得测试缺失了准确性和有效性；且对同一条实际曲线，在不同采样周期下，会生成不同精度的趋势预测曲线，若因采样周期过大导致趋势预测曲线过大，则会失真于真实系统的趋势预测曲线。测试曲线失真导致错误阈值判断的示意如图3所示。现有模型生成的测试曲线采样周期为50ms，而真实系统采样周期为15ms，50ms相邻采样值变化率远大于的15ms相邻采样值变化率，趋势预测曲线失真。

(3)峰值数据缺失：对于变化剧烈，如斜率变化较大的曲线拐点，部分重要的峰值数据会因采样精度不足而缺失。

发明内容

有鉴于此，本发明的目的在于提出一种核电厂仪控系统测试曲线构建方法、装置、设备、介质，以改善上述问题。

本发明采用了如下方案：

一种核电厂仪控系统测试曲线构建方法，其包括：

获取原始输入曲线，根据曲线时间长度t及原始采样周期T，计算原始采样点数N以及各个原始采样点的坐标；

设定目标输入曲线的曲线参数；所述曲线参数包括目标采样周期T’、目标采样点数N’以及插值精度K；

根据原始输入曲线的曲线斜率对原始输入曲线进行分段，得到m个曲线段；

根据插值精度K对每个曲线段进行分段，得到与每个曲线段对应的多个曲线分段；其中，每个曲线分段包括K个采样点数；

对每个曲线分段进行K阶牛顿插值表达式的计算，得到各曲线分段的插值表达式；

根据每个曲线分段的时间长度、目标采样点数以及所述插值表达式计算各个点位的插值，以实现对曲线分段的插值；

将已完成插值的曲线分段重新拼接，获得目标曲线。

优选地，根据曲线时间长度t及采样原始采样周期T，计算原始采样点数N＝t/T，得到原始曲线的所有点位M＝{Xn,Yn}。

优选地，根据原始输入曲线的曲线斜率对原始输入曲线进行分段，得到m个曲线段，具体包括：

计算原始输入曲线的曲线斜率；

当出现曲线斜率由正变为负，或由负变为正时，在曲线斜率为0处进行分段，以得到m个曲线段；其中，记m个曲线段为m＝{M1，M2，…，Mm}。

优选地，根据插值精度K对每个曲线段进行分段，得到与每个曲线段对应的多个曲线分段具体包括：

将每个曲线段Mm根据插值精度K分为Mmn＝Mm/K；其中，对于Mm的原始点数不能被K整除，则将最后一个曲线分段的插值精度设置为K’，K’为剩余点数。

优选地，各个曲线分段Mmn的插值表达式：

Nk(X)|Mmn＝

f(X0)+f(X0,X1)(X-X0)+…+f(X0,…,Xk)(X-X0)(X-X1)…(X-X(k-1))|k＝K’；其中，

f(X0,X1)＝[f(X1)-f(X0)]/(X1-X0)；

f(X0,X1,X2)＝[f(X0,X2)-f(X0,X1)]/(X2-X1)；

f[X0,X1,…Xk]＝f[X0,X1,…X(k-1)]-f[X1,…Xk]/[X0-Xk]；X_k为第K个采样点的横坐标值。

优选地，在Mmn中，设tmn为该段时间长度，则该曲线分段的目标点数n’＝tmn/T；则根据每个曲线分段的时间长度、目标采样点数以及所述插值表达式计算各个点位的插值具体为：

根据插值表达式以及第M段插值点位n’计算插值(Xn’,Yn’)|Mmn，Yn’|Mmn＝Nk(Xn’|M)|Mmn，表示插值后第M段的第n个点位的插值，插值后的曲线分段记为M’mn＝{Xn’,Yn’}。

优选地，将已完成插值的曲线分段重新拼接，获得目标曲线具体包括：

对各个曲线段Mm，将M’mn拼接，记为M’m；其中，设Xn’表示第M’mn段的第n’个点位,Xn表示第M’mn段最后一个点位；Xm表示拼接后的点位，则Xn’|M’mn等于第M’m(n-1)段最后一个点位值加n’,即Xm’|M’m＝Xn’|M’mn＝Xn|M’m(n-1)+n’，得到拼接后的M’m曲线点位M’m＝{Xm,Ym}

将M’m拼接为完整曲线M’；其中，设Xm’表示第Mm段的第m’个点位，Xm表示第Mm段的最后一个点位，XN表示拼接后的点位，则XN’|M’＝Xm’|M’m＝Xm|M’(m-1)+N’

得到拼接后的插值曲线M’＝{XN,YN}。

本发明实施例还提供了一种核电厂仪控系统测试曲线构建装置，其包括：

原始输入曲线获取单元，用于获取原始输入曲线，根据曲线时间长度t及原始采样周期T，计算原始采样点数N以及各个原始采样点的坐标；

目标设定单元，用于设定目标输入曲线的曲线参数；所述曲线参数包括目标采样周期T’、目标采样点数N’以及插值精度K；

第一分段单元，用于根据原始输入曲线的曲线斜率对原始输入曲线进行分段，得到m个曲线段；

第二分段单元，用于根据插值精度K对每个曲线段进行分段，得到与每个曲线段对应的多个曲线分段；其中，每个曲线分段包括K个采样点数；

插值表达式获取单元，用于对每个曲线分段进行K阶牛顿插值表达式的计算，得到各曲线分段的插值表达式；

插值单元，用于根据每个曲线分段的时间长度、目标采样点数以及所述插值表达式计算各个点位的插值，以实现对曲线分段的插值；

拼接单元，用于将已完成插值的曲线分段重新拼接，获得目标曲线。

本发明实施例还提供了一种核电厂仪控系统测试曲线构建设备，其包括存储器以及处理器，所述存储器内存储有计算机程序，所述计算机程序能够被所述处理器执行，以实现如上述的核电厂仪控系统测试曲线构建方法。

本发明实施例还提供了一种计算机可读存储介质，其存储有计算机程序，所述计算机程序能够被所述计算机可读存储介质所在设备的处理器执行，以实现如上述的核电厂仪控系统测试曲线构建方法。

本发明实施例通过对原始输入曲线进行分段及插值的方式来获得目标曲线，提高核电厂运行时的参数变化曲线的模拟精度，构建更接近真实运行环境的参数变化曲线，对反应堆保护系统进行统计测试，解决现有统计测试曲线构建中因数据采样周期大造成的失真度高、覆盖率低、趋势预判错误、峰值数据缺失等问题，充分模拟核电厂的运行参数曲线，提高统计测试的完备性、准确性、有效性，为核电厂仪控系统统计测试、可靠性评估等提供技术支撑。

附图说明

为了更清楚地说明本发明实施方式的技术方案，下面将对实施方式中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。

图1是现有的统计测试输入曲线示意图。

图2是现有的趋势预测曲线示意。

图3是统计测试曲线失真示意图。

图4是本发明第一实施例提供的核电厂仪控系统测试曲线构建方法的流程示意图。

图5是本发明第二实施例提供的核电厂仪控系统测试曲线构建装置的结构示意图。

具体实施方式

为使本发明实施方式的目的、技术方案和优点更加清楚，下面将结合本发明实施方式中的附图，对本发明实施方式中的技术方案进行清楚、完整地描述，显然，所描述的实施方式是本发明一部分实施方式，而不是全部的实施方式。基于本发明中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都属于本发明保护的范围。因此，以下对在附图中提供的本发明的实施方式的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施方式。基于本发明中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都属于本发明保护的范围。

请参阅图4，本发明第一实施例提供了一种核电厂仪控系统测试曲线构建方法，其可由核电厂仪控系统测试曲线构建设备来执行，特别的，由所述构建设备内的一个或者多个处理器来执行，以实现如下步骤：

S101，获取原始输入曲线，根据曲线时间长度t及原始采样周期T，计算原始采样点数N以及各个原始采样点的坐标。

在本实施例中，所述原始输入曲线即为原始的核电厂仪控系统的统计测试曲线。

在本实施例中，在获得所述原始输入曲线后，可对该曲线进行取点并获得各个原始采样点的坐标。

其中，原始采样点数N＝t/T，得到原始曲线的所有点位的坐标记为M＝{Xn,Yn}。

S102，设定目标输入曲线的曲线参数；所述曲线参数包括目标采样周期T’、目标采样点数N’以及插值精度K。

S103，根据原始输入曲线的曲线斜率对原始输入曲线进行分段，得到m个曲线段。

具体地，

首先，计算原始输入曲线的曲线斜率。

其中，曲线斜率的计算可采样求导方法或者其他方法，本发明不做具体限定。

然后，当出现曲线斜率由正变为负，或由负变为正时，在曲线斜率为0处进行分段，以得到m个曲线段；其中，记m个曲线段为m＝{M1，M2，…，Mm}。

S104，根据插值精度K对每个曲线段进行分段，得到与每个曲线段对应的多个曲线分段；其中，每个曲线分段包括K个采样点数。

具体地，将M中的每个段Mm根据插值精度K分为Mmn＝Mm/K，如M1＝{M11,M12,M13,…,M1n}；Mm＝{Mm1,Mm2,…,Mmn}。

考虑到由于Mm的原始点数存在不被K整除的可能，最后一段会出现点数小于K的情况，因此定义插值精度K’，即：

1)原段点数≥K时，K’＝K；

2)原段点数<K时，则K’＝剩余点数。

S105，对每个曲线分段进行K阶牛顿插值表达式的计算，得到各曲线分段的插值表达式。

其中，各个曲线分段Mmn的插值表达式：

Nk(X)|Mmn＝

f(X0)+f(X0,X1)(X-X0)+…+f(X0,…,Xk)(X-X0)(X-X1)…(X-X(k-1))|k＝K’；其中，

f(X0,X1)＝[f(X1)-f(X0)]/(X1-X0)；

f(X0,X1,X2)＝[f(X0,X2)-f(X0,X1)]/(X2-X1)；

f[X0,X1,…Xk]＝f[X0,X1,…X(k-1)]-f[X1,…Xk]/[X0-Xk]；X_k为第K个采样点的横坐标值。

S106，根据每个曲线分段的时间长度、目标采样点数以及所述插值表达式计算各个点位的插值，以实现对曲线分段的插值。

具体地，在Mmn中，设tmn为该段时间长度，则该曲线分段的目标点数n’＝tmn/T；则根据每个曲线分段的时间长度、目标采样点数以及所述插值表达式计算各个点位的插值具体为：

根据插值表达式以及第M段插值点位n’计算插值(Xn’,Yn’)|Mmn，Yn’|Mmn＝Nk(Xn’|M)|Mmn，表示插值后第M段的第n个点位的插值，插值后的曲线分段记为M’mn＝{Xn’,Yn’}。

S107，将已完成插值的曲线分段重新拼接，获得目标曲线。

具体地：

首先，对各个曲线段Mm，将M’mn拼接，记为M’m。

其中，设Xn’表示第M’mn段的第n’个点位,Xn表示第M’mn段最后一个点位；Xm表示拼接后的点位，则Xn’|M’mn等于第M’m(n-1)段最后一个点位值加n’,即Xm’|M’m＝Xn’|M’mn＝Xn|M’m(n-1)+n’，得到拼接后的M’m曲线点位M’m＝{Xm,Ym}

然后，将M’m拼接为完整曲线M’。

其中，设Xm’表示第Mm段的第m’个点位，Xm表示第Mm段的最后一个点位，XN表示拼接后的点位，则XN’|M’＝Xm’|M’m＝Xm|M’(m-1)+N’；

最后，得到拼接后的插值曲线M’＝{XN,YN}。

综上所述，本发明实施例通过对原始输入曲线进行分段及插值的方式来获得目标曲线，提高核电厂运行时的参数变化曲线的模拟精度，构建更接近真实运行环境的参数变化曲线，对反应堆保护系统进行统计测试，解决现有统计测试曲线构建中因数据采样周期大造成的失真度高、覆盖率低、趋势预判错误、峰值数据缺失等问题，充分模拟核电厂的运行参数曲线，提高统计测试的完备性、准确性、有效性，为核电厂仪控系统统计测试、可靠性评估等提供技术支撑。

请参阅图5，本发明第二实施例还提供了一种核电厂仪控系统测试曲线构建装置，其包括：

原始输入曲线获取单元210，用于获取原始输入曲线，根据曲线时间长度t及原始采样周期T，计算原始采样点数N以及各个原始采样点的坐标；

目标设定单元220，用于设定目标输入曲线的曲线参数；所述曲线参数包括目标采样周期T’、目标采样点数N’以及插值精度K；

第一分段单元230，用于根据原始输入曲线的曲线斜率对原始输入曲线进行分段，得到m个曲线段；

第二分段单元240，用于根据插值精度K对每个曲线段进行分段，得到与每个曲线段对应的多个曲线分段；其中，每个曲线分段包括K个采样点数；

插值表达式获取单元250，用于对每个曲线分段进行K阶牛顿插值表达式的计算，得到各曲线分段的插值表达式；

插值单元260，用于根据每个曲线分段的时间长度、目标采样点数以及所述插值表达式计算各个点位的插值，以实现对曲线分段的插值；

拼接单元270，用于将已完成插值的曲线分段重新拼接，获得目标曲线。

本发明第三实施例还提供了一种核电厂仪控系统测试曲线构建设备，其包括存储器以及处理器，所述存储器内存储有计算机程序，所述计算机程序能够被所述处理器执行，以实现如上述的核电厂仪控系统测试曲线构建方法。

本发明第四实施例还提供了一种计算机可读存储介质，其存储有计算机程序，所述计算机程序能够被所述计算机可读存储介质所在设备的处理器执行，以实现如上述的核电厂仪控系统测试曲线构建方法。

在本发明所提供的几个实施例中，应该理解到，所揭露的装置和方法，也可以通过其它的方式实现。以上所描述的装置和方法实施例仅仅是示意性的，例如，附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。

另外，在本发明各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。

所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，电子设备，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。

以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (5)

1.一种核电厂仪控系统测试曲线构建方法，其特征在于，

获取原始输入曲线，根据曲线时间长度t及原始采样周期T，计算原始采样点数N以及各个原始采样点的坐标，N=t/T，得到原始曲线的所有点位M={Xn,Yn}；

设定目标输入曲线的曲线参数；所述曲线参数包括目标采样周期T’、目标采样点数N’以及插值精度K；

根据原始输入曲线的曲线斜率对原始输入曲线进行分段，得到m个曲线段，具体包括：

计算原始输入曲线的曲线斜率；

当出现曲线斜率由正变为负，或由负变为正时，在曲线斜率为0处进行分段，以得到m个单调曲线段；其中，记m个曲线段为m={M1，M2，…，Mm}；

为每个曲线段设置插值精度K，根据插值精度K对每个曲线段进行分段，得到与每个曲线段对应的多个曲线分段； 具体包括：将每个曲线段Mm根据插值精度K分为Mmn=Mm/K；其中，对于Mm的原始点数不能被K整除，则将最后一个曲线分段的插值精度设置为K’，K’为剩余点数；

对每个曲线分段进行K阶牛顿插值表达式的计算，得到各曲线分段的插值表达式；

在Mmn中，设tmn为该段时间长度，则该曲线分段的目标点数n’=tmn/T；则根据每个曲线分段的时间长度、目标采样点数以及所述插值表达式计算各个点位的插值，具体为：

根据插值表达式以及第M段插值点位n’计算插值（Xn’,Yn’）|Mmn，Yn’|Mmn=NK(Xn’|M)|Mmn，表示插值后第M段的第n个点位的插值，插值后的曲线分段记为M’mn={Xn’, Yn’}；

将已完成插值的曲线分段重新拼接，获得目标曲线。

2.根据权利要求1所述的核电厂仪控系统测试曲线构建方法，其特征在于，将已完成插值的曲线分段重新拼接，获得目标曲线具体包括：

对各个曲线段Mm，将M’mn拼接，记为M’m；其中，设Xn’表示第M’mn段的第n’个点位, Xn表示第M’mn段最后一个点位；Xm表示拼接后的点位，则Xn’|M’mn等于第M’m(n-1)段最后一个点位值加n’, 即Xm’|M’m =Xn’|M’mn= Xn|M’m(n-1) + n’，得到拼接后的M’m曲线点位M’m={Xm,Ym}

将M’m拼接为完整曲线M’；其中，设Xm’表示第Mm段的第m’个点位，Xm表示第Mm段的最后一个点位，XN表示拼接后的点位，则XN’|M’= Xm’|M’m= Xm|M’(m-1) + N’

得到拼接后的插值曲线M’={XN, YN}。

3.一种核电厂仪控系统测试曲线构建装置，其特征在于，包括：

原始输入曲线获取单元，用于获取原始输入曲线，根据曲线时间长度t及原始采样周期T，计算原始采样点数N，N=t/T，得到原始曲线的所有点位M={Xn,Yn}；

目标设定单元，用于设定目标输入曲线的曲线参数；所述曲线参数包括目标采样周期T’、目标采样点数N’以及插值精度K；

第一分段单元，用于计算原始曲线上各点斜率并根据原始输入曲线的曲线斜率，对原始输入曲线进行分段，得到m个曲线段；当出现曲线斜率由正变为负，或由负变为正时，在曲线斜率为0处进行分段，以得到m个单调曲线段；其中，记m个曲线段为m={M1，M2，…，Mm}；

第二分段单元，用于根据插值精度K对每个曲线段进行分段，得到与每个曲线段对应的多个曲线分段；将每个曲线段Mm根据插值精度K分为Mmn=Mm/K；其中，对于Mm的原始点数不能被K整除，则将最后一个曲线分段的插值精度设置为K’，K’为剩余点数；

插值表达式获取单元，用于对每个曲线分段进行K阶牛顿插值表达式的计算，得到各曲线分段的插值表达式；

插值单元，用于根据每个曲线分段的时间长度、目标采样点数以及所述插值表达式计算各个点位的插值，以实现对曲线分段的插值；在Mmn中，设tmn为该段时间长度，则该曲线分段的目标点数n’=tmn/T，根据插值表达式以及第M段插值点位n’计算插值（Xn’,Yn’）|Mmn，Yn’|Mmn=NK(Xn’|M)|Mmn，表示插值后第M段的第n个点位的插值，插值后的曲线分段记为M’mn={Xn’, Yn’}；

拼接单元，用于将已完成插值的曲线分段重新拼接，获得目标曲线，对各个曲线段Mm，将M’mn拼接，记为M’m；其中，设Xn’表示第M’mn段的第n’个点位, Xn表示第M’mn段最后一个点位；Xm表示拼接后的点位，则Xn’|M’mn等于第M’m(n-1)段最后一个点位值加n’, 即Xm’|M’m =Xn’|M’mn= Xn|M’m(n-1) + n’，得到拼接后的M’m曲线点位M’m={Xm,Ym}，将M’m拼接为完整曲线M’；其中，设Xm’表示第Mm段的第m’个点位，Xm表示第Mm段的最后一个点位，XN表示拼接后的点位，则XN’|M’= Xm’|M’m= Xm|M’(m-1) + N’，得到拼接后的插值曲线M’={XN,YN}。

4.一种核电厂仪控系统测试曲线构建设备，其特征在于，包括存储器以及处理器，所述存储器内存储有计算机程序，所述计算机程序能够被所述处理器执行，以实现如权利要求1至2任意一项所述的核电厂仪控系统测试曲线构建方法。

5.一种计算机可读存储介质，其特征在于，存储有计算机程序，所述计算机程序能够被所述计算机可读存储介质所在设备的处理器执行，以实现如权利要求1至2任意一项所述的核电厂仪控系统测试曲线构建方法。

Images