基于LUENBERGER形式H∞滤波的银自给能探测器信号延迟消除方法.pdf

本发明公开了基于Luenberger形式H∞滤波的银自给能探测器信号延迟消除方法，包括依次进行的以下步骤：步骤1、建立银与热中子的核反应模型；步骤2、采用去耦变换建立核反应模型对应的离散状态方程；步骤3、确定银自给能探测器电流的瞬时响应份额；步骤4、利用Luenberger形式的H∞滤波器对银自给能探测器电流信号作延迟消除。本发明应用时能对银自给能中子探测器的电流信号进行延迟消除处理，并能有效抑制噪声，使得银自给能中子探测器在反应堆瞬态工况时也能正常使用，且由于本发明采用了Luenberger形式的H∞滤波器，作延迟消除时无需预先知道外部扰动输入信号的统计特性。

权利要求书
1.  基于Luenberger形式H∞滤波的银自给能探测器信号延迟消除方法，其特征在于：包括以下步骤：
步骤1、建立银与热中子的核反应模型：
在反应堆瞬态工况下，通量的变化引起银自给能中子探测器电流的变化并不同步，后者较前者有一定的滞后，描述上述反应的具体公式如下：
d108Ag(t)dt=σ107107Agφ(t)-λ108Ag108(t)---(1)]]>
d110Ag(t)dt=σ109109Agφ(t)-λ110Ag110(t)---(2)]]>
I(t)＝(109K109σ109Ag+107K107σ107Ag)φ(t)
              (3)
+110K110λ110Ag(t)+108K108λ108Ag(t)
其中，108Ag(t)表示108Ag的核密度，110Ag(t)表示110Ag的核密度，107Ag表示107Ag的核密度，109Ag表示109Ag的核密度，φ(t)表示探测器处中子通量，107σ表示107Ag的中子俘获截面，109σ表示109Ag的中子俘获截面，108λ表示108Ag的β衰变常数，110λ表示110Ag的β衰变常数，107K表示107Ag俘获中子后产生电流的概率，109K表示109Ag俘获中子后产生电流的概率，108K表示108Ag发生β衰变后产生电流的概率，110K表示110Ag发生β衰变后产生电流的概率，I(t)表示SPND电流；
步骤2、采用去耦变换获取核反应模型对应的离散状态方程：
使用Laplace变换将动态模型化为：
dx1(t)dt=-λ110x1(t)+λ110qpx3(t)---(4)]]>
dx2(t)dt=-λ108x2(t)+λ108rpx3(t)---(5)]]>
I(t)＝p(x1(t)+x2(t)+x3(t))  (6)
其中p为瞬时电流份额，q为110Ag(t)发生β衰变对应的延迟电流份额，r为108Ag(t)发生β衰变对应的延迟电流份额；
将式(4)、(5)和(6)进行时间离散化处理，并且添加过程噪声项和测量噪声项，可以得到离散状态方程如下：
xk+1=100qp·(1-e-λ108Ts)e-λ108Ts0rp·(1-e-λTs110)0e-λ110Ts·xk+100·wk---(7)]]>
Ik＝[p p p]xk+vk  (8)
φk＝[1 0 0]xk  (9)
其中xk=φkx1kx2kT,]]>wk为过程噪声，vk为测量噪声，TS为采样时间；
初始值为
x0=φ0x10x20=I0prI0p2rqI0---(10);]]>
步骤3、确定银自给能探测器电流的瞬时响应份额：
在反应堆启动物理实验阶段，通过升/降反应堆功率形成功率台阶，记录相应的堆外探测器信号实测值与银自给能探测器信号实测值；堆外探测器能够瞬时响应中子通量的变化，相应的实测值可认为是真实的中子通量；通过调整瞬时响应份额的理论值给定N个不同的瞬时响应份额预测值，再将堆外探测器信号实测值代入离散状态方程，可以得到N组银自给能探测器信号理论值，将 理论值与银自给能探测器信号实测值进行比较，取其中符合程度最好的某组理论值相应的瞬时响应份额预测值为后续延迟消除所采用的瞬时响应份额；
步骤4、利用Luenberger形式的H∞滤波器对银自给能探测器电流信号作延迟消除：
对于一个离散控制过程系统，该系统可用一个状态方程来描述：
x(k+1)＝Ax(k)+Bw(k)
y(k)＝Cx(k)+Dw(k)  (8)
z(k)＝Lx(k)
其中，x(k)为第k次采样点的n维状态向量，w(k)包含了系统过程噪声以及系统观测白噪声，y(k)为第k次采样点的测量值，z(k)为1维待求向量，L为l*n维矩阵；
针对离散系统(8)，设计如下渐近稳定的满阶线性Luenberger滤波器
x^k+1=Ax^k+K(yk-Cx^k)z^k=Lx^k---(9)]]>
式(9)为H∞滤波器，对于一个给定的γ，当且仅当如下的矩阵不等式有解：
Y0ATY-CTWTLT0γ2IBTY-DTWT0YA-WCYB-WDY0L001≥0---(10)]]>
其中Y＝YT∈Rn×n，W∈Rn×r，J＝JT∈Rm×m，H∞滤波器的增益K＝Y-1W；
对于银自给能探测器，由其离散状态方程可知方程(8)中的对应矩阵为：
A=100qp·(1-e-λ108Ts)e-λ108Ts0rp·(1-e-λ110Ts)0e-λ110Ts]]>
B=100000]]>
C＝[p p p]
D＝[0 1]
L＝[1 0 0]
通过求解线性矩阵不等式(10)，可得H∞滤波器矩阵K，从而可以由如下步骤获取消除延迟后任意时刻的探测器电流值：
由初始电流测量值可得x^(0)=y^(0)pr·y^(0)p2rq·y^(0),]]>初始0时刻延迟消除后电流值为z^(0)=Cfx^(0);]]>
对于任意k+1(k＝0,1,...)时刻，而k+1时刻延迟消除后的电流值为z^(k+1)=Cfx^(k+1).]]>

2.  根据权利要求1所述的基于Luenberger形式H∞滤波的银自给能探测器信号延迟消除方法，其特征在于，在有换挡的情况下，还包括按如下的信号处理方法对原始信号进行处理：在换档区域内，假设中子通量保持不变，然后反推中子通量密度产生的电流信号，再与探测器实际输出电流相减，得到换挡突变分量；在换档区域外，探测器输出电流减去换挡突变分量，得到中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除处理。

3.  根据权利要求1所述的基于Luenberger形式H∞滤波的银自给能探测器信号延迟消除方法，其特征在于，所述的换档区域设计结构如下：
在换档区域内(k1≤k≤k2)，假设中子通量密度不变，则有：
n(k+1)＝n(k)  (11)
x1(k+1)=e-λ108Tsx1(k)+qp·(1-e-λTs108)n(k)---(12)]]>
x2(k+1)=e-λ110Tsx2(k)+rp·(1-e-λTs110)n(k)---(13)]]>
可以反推出银自给能探测器电流信号为：
I(k+1)＝p(n(k+1)+x1(k+1)+x2(k+1))  (14)
将反推电流(14)当成探测器实际输出电流，通过权利要求1中所述步骤进行延迟消除；
在换档区域时间边界k2处，换档引起的电流偏置量可以由下式进行估算：
D=I(k2)-y^(k2)---(15)]]>
其中表示在k2时刻的探测器实际输出电流；
在换档区域外，需要对探测器实际输出电流进行偏置补偿以抵消换档所带来的影响，将探测器实际输出电流加上式(15)表示的换档引起的电流偏置量，得到中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除。

说明书基于Luenberger形式H∞滤波的银自给能探测器信号延迟消除方法
技术领域
本发明涉及核反应堆功率分布在线监测系统所用的堆内银自给能中子探测器信号的处理技术，具体是基于Luenberger形式H∞滤波的银自给能探测器信号延迟消除方法。
背景技术
用作先进堆芯测量系统堆内探测器的银自给能中子探测器，其敏感材料银与中子反应产生的次生核素发生β衰变产生电流，稳态情况下该电流大小与所在位置通量成正比，因此通过测量银自给能探测器能够推知其所在位置中子通量。由于该类探测器电流主要成分是由次生核素β衰变产生的，在反应堆瞬态情况(中子通量水平变化的情况)下，该类探测器电流不能实时反映通量水平的变化，而是有一定的延迟，延迟时间参数与次生核素的β衰变一致。因此，利用银自给能中子探测器作中子测量装置的先进堆芯测量系统，为了保证中子通量测量的准确性，需要对银自给能探器的电流信号作延迟消除处理。
由于实际的测量过程中总伴随有噪声(过程噪声和测量噪声)，利用直接的数学反演方法作延迟消除会将探测器电流信号噪声放大，最大可放大到20倍，影响测量的精度。因此，在延迟消除处理过程中，需要有效抑制噪声的放大。
目前应用于银自给能探测器信号延迟的消除主要基于Kalman滤波器实现，其应用时必须假定系统的外部扰动输入信号是一个具有已知统计特性的白噪声信号，当输入信号是一个具有有限能量的不确定信号，其统计特性难以得到，该方法就难以应用。
发明内容
本发明的目的在于克服现有技术的不足，提供了一种基于Luenberger形式H∞滤波的银自给能探测器信号延迟消除方法，其应用时能对银自给能中子探测器的电流信号进行延迟消除处理，并能有效抑制噪声，使得银自给能中子探测器在反应堆瞬态工况时也能正常使用，且由于本发明采用了Luenberger形式的H2滤波器，作延迟消除时无需预先知道外部扰动输入信号的统计特性。
本发明解决上述问题主要通过以下技术方案实现：基于Luenberger形式H∞滤波的银自给能探测器信号延迟消除方法，其特征在于，包括以下步骤：
步骤1、建立银与热中子的核反应模型：
在反应堆瞬态工况下，通量的变化引起银自给能中子探测器电流的变化并不同步，后者较前者有一定的滞后，描述上述反应的具体公式如下：
dAg108(t)dt=σ107Ag107φ(t)-λ108Ag108(t)---(1)]]>
dAg110(t)dt=σ109Ag109φ(t)-λ110Ag110(t)---(2)]]>
I(t)=(K109σ109Ag109+K107σ107Ag107)φ(t)+K110λ110Ag110(t)+K108λ108Ag108(t)---(3)]]>
其中，108Ag(t)表示108Ag的核密度，110Ag(t)表示110Ag的核密度，107Ag表示107Ag的核密度，109Ag表示109Ag的核密度，φ(t)表示探测器处中子通量，107σ表示107Ag的中子俘获截面，109σ表示109Ag的中子俘获截面，108λ表示108Ag的β衰变常数，110λ表示110Ag的β衰变常数，107K表示107Ag俘获中子后产生电流的概率，109K表示109Ag俘获中子后产生电流的概率，108K表示108Ag发生β衰变后产生电流的概率，110K表示110Ag发生β衰变后产生电流的概率，I(t)表示SPND电流；
步骤2、采用去耦变换获取核反应模型对应的离散状态方程：
使用Laplace变换将动态模型化为：
dx1(t)dt=-λ110x1(t)+λ110qpx3(t)---(4)]]>
dx2(t)dt=-λ108x2(t)+λ108qpx3(t)---(5)]]>
I(t)＝p(x1(t)+x2(t)+x3(t))(6)
其中p为瞬时电流份额，q为110Ag(t)发生β衰变对应的延迟电流份额，r为108Ag(t)发生β衰变对应的延迟电流份额；
将式(4)、(5)和(6)进行时间离散化处理，并且添加过程噪声项和测量噪声项，可以得到离散状态方程如下：
xk+1=100qp·(1-e-λ108Ts)e-λ108Ts0rp·(1-e-λ110Ts)0e-λ110Ts·xk+100·wk---(7)]]>
Ik＝[p p p]xk+vk                    (8)
φk＝[1 0 0]xk                           (9)
其中xk=φkx1kx2kT,]]>wk为过程噪声，vk为测量噪声，TS为采样时间；
初始值为
x0=φ0x10x20=I0prI0p2rqI0---(10);]]>
步骤3、确定银自给能探测器电流的瞬时响应份额：
在反应堆启动物理实验阶段，通过升/降反应堆功率形成功率台阶，记录相应的堆外探测器信号实测值与银自给能探测器信号实测值；堆外探测器能够 瞬时响应中子通量的变化，相应的实测值可认为是真实的中子通量；通过调整瞬时响应份额的理论值给定N个不同的瞬时响应份额预测值，再将堆外探测器信号实测值代入离散状态方程，可以得到N组银自给能探测器信号理论值，将理论值与银自给能探测器信号实测值进行比较，取其中符合程度最好的某组理论值相应的瞬时响应份额预测值为后续延迟消除所采用的瞬时响应份额；
步骤4、利用Luenberger形式的H∞滤波器对银自给能探测器电流信号作延迟消除：
对于一个离散控制过程系统，该系统可用一个状态方程来描述：
x(k+1)＝Ax(k)+Bw(k)
y(k)＝Cx(k)+Dw(k)                                (8)
z(k)＝Lx(k)
其中，x(k)为第k次采样点的n维状态向量，w(k)包含了系统过程噪声以及系统观测白噪声，y(k)为第k次采样点的测量值，z(k)为1维待求向量，L为l*n维矩阵；
针对离散系统(8)，设计如下渐近稳定的满阶线性Luenberger滤波器
x^k+1=Ax^k+K(yk-Cx^k)z^k=Lx^k---(9)]]>
式(9)为H∞滤波器，对于一个给定的γ，当且仅当如下的矩阵不等式有解：
Y0ATY-CTWTLT0γ2IBTY-DTWT0YA-WCYB-WDY0L00I≥0---(10)]]>
其中Y＝YT∈Rn×n，W∈Rn×r，J＝JT∈Rm×m，H∞滤波器的增益K＝Y-1W；
对于银自给能探测器，由其离散状态方程可知方程(8)中的对应矩阵为：
A=100qp·(1-e-λ108Ts)e-λ108Ts0rp·(1-e-λ110Ts)0e-λ110Ts]]>
B=100000]]>
C＝[p p p]
D＝[0 1]
L＝[1 0 0]
通过求解线性矩阵不等式(10)，可得H∞滤波器矩阵K，从而可以由如下步骤获取消除延迟后任意时刻的探测器电流值：
由初始电流测量值可得x^(0)=y^(0)pr·y^(0)p2rq·y^(0),]]>初始0时刻延迟消除后电流值为z^(0)=Cfx^(0);]]>
对于任意k+1(k＝0,1,...)时刻，而k+1时刻延迟消除后的电流值为z^(k+1)=Cfx^(k+1).]]>
本发明应用时利用Luenberger形式的H∞滤波器原理，在延迟消除过程中，可以有效地抑制噪声的放大，噪声抑制效果越好，延迟效果会逐渐变差，因此，本发明应用时需适当调节参数使延迟消除效果与噪声抑制达到最佳平衡。
当需要探测较大动态范围的中子通量密度，相应的也需要检测大动态范围的电流信号，而这一问题便集中在了模数转换器上。为了适应大动态范围的电流的量化，银自给能探测器的模数转换器采样分档电阻，当电流信号在大范围 变化时，模数转换器就会发生电阻档位转换。由于各档位没有完全匹配，各档位之间的切换会造成输出信号的近似于阶跃的突变。
换挡引起的突变分量进入延迟消除模块后，会被严重放大，使得时域上的阶跃突变被严重放大，影响最终信号延迟消除的质量(突变部分信号的严重失真)。在换挡时间段里，信号的变化主要由换挡突变贡献，相对而言，由中子通量密度变化引起的电流信号变化可以忽略。
为了处理换档导致的信号突变问题，在有换挡的情况下，还包括按如下的信号处理方法对原始信号进行处理：在换档区域内，假设中子通量保持不变，然后反推中子通量密度产生的电流信号，再与探测器实际输出电流相减，得到换挡突变分量；在换档区域外，探测器输出电流减去换挡突变分量，得到中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除处理。
所述的换档区域设计结构如下：
在换档区域内(k1≤k≤k2)，假设中子通量密度不变，则有：
n(k+1)＝n(k)                              (11)
x1(k+1)=e-λ108Tsx1(k)+qp·(1-e-λ108Ts)n(k)---(12)]]>
x2(k+1)=e-λ110Tsx2(k)+rp·(1-e-λ110Ts)n(k)---(13)]]>
可以反推出银自给能探测器电流信号为：
I(k+1)＝p(n(k+1)+x1(k+1)+x2(k+1))                 (14)
将反推电流(14)当成探测器实际输出电流，通过权利要求1中所述步骤进行延迟消除；
在换档区域时间边界k2处，换档引起的电流偏置量可以由下式进行估算：
D=I(k2)-y^(k2)---(15)]]>
其中表示在k2时刻的探测器实际输出电流；
在换档区域外，需要对探测器实际输出电流进行偏置补偿以抵消换档所带来的影响，将探测器实际输出电流加上式(15)表示的换档引起的电流偏置量，得到中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除。
在换档区域外，需要对探测器实际输出电流进行偏置补偿以抵消换档所带来的影响，将探测器实际输出电流加上式(15)表示的换档引起的电流偏置量，得到中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除。
综上所述，本发明具有以下有益效果：
1本发明整体工序简单，便于实现，能对银自给能中子探测器的电流信号进行延迟消除处理，并能有效抑制噪声，使得银自给能中子探测器在反应堆瞬态工况时也能正常使用；本发明基于Luenberger形式的H∞滤波器实现，在输入信号是一个具有有限能量的不确定信号时也能正常应用；本发明应用时将滤波器设计转化为相应线性矩阵不等式计算，方便计算，可以方便地使用Matlab的LMI Toolbox进行求解。
2本发明解决了先进堆芯测量系统(核反应堆功率分布在线监测系统)所用的堆内银自给能中子探测器信号的延迟消除问题。利用H∞滤波器对银自给能中子探测器信号进行延迟消除、平滑、降噪处理，通过适当选取Luenberger形式的H∞滤波器参数，能够很好的达到信号延迟消除效果和噪声抑制效果的最佳平衡。本发明能够保证银自给能探测器电流信号直接用于先进堆芯测量系统后续环节，而不丧失准确度。
3本发明对银自给能中子探测器的电流信号进行延迟消除处理，响应时间(阶跃通量变化时，信号恢复到稳态电流的90％所需的时间)在2～5秒内。
4本发明对银自给能中子探测器的电流信号延迟消除过程中，对测量电流 信号进行降噪处理，噪声放大倍数(延迟消除处理后的电流相对误差与噪声之比)抑制在1～4倍。
5本发明能有效处理因硬件换挡造成的阶跃对延迟消除效果的影响。
附图说明
图1为本发明的银自给能中子探测器结构图
图2为本发明一个具体实施例的处理流程图；
图3为银与热中子核反应图。
附图中标记及相应的零部件名称：
1—发射极，2—绝缘层，3—收集极，4-导线，5-保护壳，6-绝缘电缆，7-电流线，8-本底线，9-密封管，10-电流输出端。
具体实施方式
下面结合实施例及附图，对本发明做进一步地的详细说明，但本发明的实施方式不限于此。
实施例：
如图1所示的银自给能中子探测器结构图，其中各个序号的零部件名称对应为：1—发射极，2—绝缘层，3—收集极，4-导线，5-保护壳，6-绝缘电缆，7-电流线，8-本底线，9-密封管，10-电流输出端，该银自给能中子探测器，其特性参数为：λ108＝ln2/2.42/60s-1＝0.0048s-1,λ110＝ln2/24.4s-1＝0.0284s-1，p＝0.09，q＝0.66，r＝0.25；图3为银与中子核反应原理过程图，对于图3的反应过程中，采用图1的装置进行测量。如图2所示，基于Luenberger形式H∞滤波的银自给能探测器信号延迟消除方法，包括依次进行的以下步骤：步骤1、建立银与热中子的核反应模型；步骤2、采用去耦变换建立核反应模型对应的离散状态方程；步骤3、确定银自给能探测器电流的瞬时响应份额；步骤4、利用 Luenberger形式的H∞滤波器对银自给能探测器电流信号作延迟消除。
本实施例建立银与热中子的核反应模型的具体实施步骤如下：如图2所示，在反应堆瞬态工况下，通量的变化引起银自给能中子探测器电流的变化并不同步，后者较前者有一定的滞后，描述上述反应的具体公式如下：
dAg108(t)dt=σ107Ag107φ(t)-λ108Ag108(t)---(1)]]>
dAg110(t)dt=σ109Ag109φ(t)-λ110Ag110(t)---(2)]]>
I(t)=(K109σ109Ag109+K107σ107Ag107)φ(t)+K110λ110Ag110(t)+K108λ108Ag108(t)---(3)]]>
其中，108Ag(t)表示108Ag的核密度，110Ag(t)表示110Ag的核密度，107Ag表示107Ag的核密度，109Ag表示109Ag的核密度，φ(t)表示探测器处中子通量，107σ表示107Ag的中子俘获截面，109σ表示109Ag的中子俘获截面，108λ表示108Ag的β衰变常数，110λ表示110Ag的β衰变常数，107K表示107Ag俘获中子后产生电流的概率，109K表示109Ag俘获中子后产生电流的概率，108K表示108Ag发生β衰变后产生电流的概率，110K表示110Ag发生β衰变后产生电流的概率，I(t)表示SPND电流。
本实施例采用去耦变换建立核反应模型对应的离散状态方程的具体实施步骤如下：
使用Laplace变换将动态模型化为：
dx1(t)dt=-λ110x1(t)+λ110qpx3(t)---(4)]]>
dx2(t)dt=-λ108x2(t)+λ108qpx3(t)---(5)]]>
I(t)＝p(x1(t)+x2(t)+x3(t))                   (6)
其中p为瞬时电流份额，q为110Ag(t)发生β衰变对应的延迟电流份额，r为108Ag(t)发生β衰变对应的延迟电流份额。
将式(4)、(5)和(6)进行时间离散化处理，并且添加过程噪声项和测量噪声项，可以得到离散状态方程如下：
xk+1=100qp·(1-e-λ108Ts)e-λ108Ts0rp·(1-e-λ110Ts)0e-λ110Ts·xk+100·wk---(7)]]>
Ik＝[p p p]xk+vk              (8)
φk＝[1 0 0]xk                     (9)
其中xk=φkx1kx2kT,]]>wk为过程噪声，vk为测量噪声，TS为采样时间；
初始值为
x0=φ0x10x20=I0prI0p2rqI0---(10)]]>
本实施例确定银自给能探测器电流的瞬时响应份额的具体实施步骤如下：
在反应堆启动物理实验阶段，通过升/降反应堆功率形成功率台阶，记录相应的堆外探测器信号实测值与银自给能探测器信号实测值。堆外探测器能够瞬时响应中子通量的变化，相应的实测值可认为是真实的中子通量。通过调整瞬时响应份额的理论值给定N个不同的瞬时响应份额预测值，再将堆外探测器信号实测值代入离散状态方程，可以得到N组银自给能探测器信号理论值，将 理论值与银自给能探测器信号实测值进行比较，取其中符合程度最好的某组理论值相应的瞬时响应份额预测值为后续延迟消除所采用的瞬时响应份额。
本实施例利用H∞滤波器对银自给能探测器电流信号作延迟消除的具体实施步骤如下：
对于一个离散控制过程系统，该系统可用一个状态方程来描述：
x(k+1)＝Ax(k)+Bw(k)
y(k)＝Cx(k)+Dw(k)                         (11)
z(k)＝Lx(k)
其中，x(k)为第k次采样点的n维状态向量，w(k)包含了系统过程噪声以及系统观测白噪声，y(k)为第k次采样点的测量值，z(k)为1维待求向量，L为l*n维矩阵；
针对离散系统(11)，设计如下渐近稳定的满阶线性Luenberger滤波器
x^k+1=Ax^k+K(yk-Cx^k)z^k=Lx^k---(12)]]>
式(12)为H∞滤波器，对于一个给定的γ，当且仅当如下的矩阵不等式有解：
Y0ATY-CTWTLT0γ2IBTY-DTWT0YA-WCYB-WDY0L00I≥0---(13)]]>
其中Y＝YT∈Rn×n，W∈Rn×r，J＝JT∈Rm×m，H∞滤波器的增益K＝Y-1W；
对于银自给能探测器，由其离散状态方程可知方程(11)中的对应矩阵为：
A=100qp·(1-e-λ108Ts)e-λ108Ts0rp·(1-e-λ110Ts)0e-λ110Ts]]>
B=100000]]>
C＝[p p p]
D＝[0 1]
L＝[1 0 0]
通过求解线性矩阵不等式(13)，可得H∞滤波器矩阵K，从而可以由如下步骤获取消除延迟后任意时刻的探测器电流值：
由初始电流测量值可得x^(0)=y^(0)pr·y^(0)p2rq·y^(0),]]>初始0时刻延迟消除后电流值为z^(0)=Cfx^(0);]]>
对于任意k+1(k＝0,1,...)时刻，而k+1时刻延迟消除后的电流值为z^(k+1)=Cfx^(k+1).]]>
实施例2：
本实施例在实施例1的基础上做出了如下进一步限定：在有换挡的情况下，所述步骤4采用以下方式来进行延迟消除：
在换档区域内(k1≤k≤k2)，假设中子通量密度不变，则有：
n(k+1)＝n(k)(14)
x1(k+1)=e-λ108Tsx1(k)+qp·(1-e-λ108Ts)n(k)---(15)]]>
x2(k+1)=e-λ110Tsx2(k)+rp·(1-e-λ110Ts)n(k)---(16)]]>
可以反推出银自给能探测器电流信号为：
I(k+1)＝p(n(k+1)+x1(k+1)+x2(k+1))(17)
将反推电流(17)当成探测器实际输出电流，通过步骤4进行延迟消除。
在换档区域时间边界k2处，换档引起的电流偏置量可以由下式进行估算：
D=I(k2)-y^(k2)---(18)]]>
其中表示在k2时刻的探测器实际输出电流。
在换档区域外，需要对探测器实际输出电流进行偏置补偿以抵消换档所带来的影响，将探测器实际输出电流加上式(18)表示的换档引起的电流偏置量，得到中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除。
以上所述，仅是本发明的较佳实施例，并非对本发明做任何形式上的限制，凡是依据本发明的技术实质上对以上实施例所作的任何简单修改、等同变化，均落入本发明的保护范围之内。