9] 由初始电流测量值
，初始0时刻延迟消除后电流值为 i(〇) =c,i(〇)；
[0120] 对于任意k+1 化=0, 1，...）时亥Ij，+ 1)=W') + 公,.v''(A')，而k+1 时刻延迟消 除后的电流值为+1) =C./iX(A' +1)。
[0121] 本实施例在信号采样周期为0. 16秒的情况下，噪声为测量电流信号的1%。
[0122] 实施例2;
[0123] 本实施例在实施例1的基础上做出了如下进一步限定；本实施例在有换挡的情况 下，还包括按如下的信号处理方法对原始信号进行处理；在换档区域内，假设中子通量保持 不变，然后反推中子通量密度产生的电流信号，再与探测器实际输出电流相减，得到换挡突 变分量；在换档区域外，探测器输出电流减去换挡突变分量，得到中子通量密度产生的电流 信号，然后再对此电流信号进行延迟消除处理。
[0124] 本实施例换档区域设计结构如下；在换档区域内化1《k《k,)，假设中子通量密 度不变，则有：
[0125] n(k+1) =n(k) (巧）
[0126]
[012引可W反推出锭自给能探测器电流信号为：
[0129] I(k+1) =Ja化+l)+cn化+1) (18)
[0130] 将反推电流（18)当成探测器实际输出电流，然后进行延迟消除；
[0131] 其中，在换档区域时间边界k,处，换档引起的电流偏置量可W由下式进行估算：
[0132] D二 /似-.i>似 (19)
[0133]辛作:)表示在k2时刻的探测器实际输出电流。
[0134] 在换档区域外，需要对探测器实际输出电流进行偏置补偿W抵消换档所带来的影 响。因此本实施例在换档区域外，将探测器实际输出电流加上式（19)表示的换档引起的电 流偏置量，得到中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除。
[0135] 如上所述，可较好的实现本发明。
【主权项】
1.基于H2/H C?混合滤波的铑自给能探测器信号延迟消除方法，其特征在于：包括以下 步骤： 步骤1、建立铑与热中子的核反应模型： 在反应堆瞬态工况下，通量的变化引起铑自给能中子探测器电流的变化并不同步，后 者较前者有一定的滞后，描述上述反应的具体公式如下：(1) 彳Λ、 I (t) = cn (t) + λ Jin1 (t) (3) 其中，niiahn^t)分别表示ltl4Rh和ltl4mRh直接引起的电荷量，n(t)表示探测器处热中 子通量对应的探测器平衡状态下的探测器电流，λ i、λ 2分别表示ltl4Rh和ltl4niRh的衰变常 数，c表示探测器电流的瞬时响应份额， &1、&2分别表示ltl4Rh和ltl4mRh引起的电流份额，I⑴ 表示铭自给能电流； 步骤2、采用直接变换建立核反应模型对应的离散状态方程： 令上⑴=X1JIi1 (t)代入式（1)、式（2)及式（3)中，将连续时间的常微分方程直接进 行离散化，并添加噪声项，得到以下离散状态方程：(4) 其中：，W(k)为过程噪声项，V(k)为测量噪声项， 初始值为C7)； 少骤3、俩疋铱目葙酡休测益电孤tf」瞬H、」_胆併棚： 步骤4、利用H2/H c?混合滤波器对铑自给能探测器电流信号作延迟消除： 对于一个离散控制过程系统，该系统可用一个状态方程来描述： X (k+1) = Ax (k) +B1W (k) +B2V (k) y (k) = Cx (k) +D1W (k) +D2V (k) (8) z (k) = Lx (k) 其中，x(k)为第k次采样点的η维状态向量，y(k)为第k次采样点的测量值，z(k)为 1维待求向量，L为l*n维矩阵； 假定系统是渐近稳定的，则对给定的常数γ >〇,要求设计一个渐近稳定的满阶线性滤 波器 x(k + I) = A.x(k) + B, ν(Λ) j 卜 (9) z(k)=crm 使得 x(k +1) - Af(/:) + Bixv(k) + B2v(k) ( ) z(k) - Cx{k) 是渐近稳定的，则对应于通道vif的滤波误差方差有一个上界，即~， 对应于通道w - 5的滤波误差向量满足|fI2 <丫|卜|2，其中对于给定的常数γ >〇以及traceX)，系统存在一个Η2/Η c?混合滤波器，当且仅当以下 的线性矩阵不等式成立其中R、X为待求解的对称正定矩阵，而S、Z、T为待求解的一般矩阵； 得到上述矩阵后，H2/H c?混合滤波器的相关矩阵表示如下： Af= (R-X) ^1S, Bf= (R-X) ^1Z, Cf= T (14)； 对于铑自给能探测器，由其离散状态方程可知方程（8)中的对应矩阵为：通过求解线性矩阵不等式（11)、（12)，可得H2/H C?混合滤波器矩阵Af、Bf、Cf，从 而可以由如下步骤获取消除延迟后任意时刻的探测器电流值：由初始电流测量值夕(〇) 可獨，初始〇时刻延迟消除后电流值为5(0) = C, i(0);对于任意 k+1 (k = 0, 1，...）时刻，+ = + ，而k+Ι时刻延迟消除后的电流值为 i(A+1) = C,i(A+1)"2. 根据权利要求1所述的基于H2/H c?混合滤波的铑自给能探测器信号延迟消除方 法，其特征在于，所述步骤三中确定铑自给能探测器电流的瞬时响应份额具体包括以下步 骤：在反应堆启动物理实验阶段，通过升/降反应堆功率形成功率台阶，记录相应的堆外探 测器信号实测值与铑自给能探测器信号实测值，其中，堆外探测器能够瞬时响应中子通量 的变化，相应的实测值可认为是真实的中子通量；通过调整瞬时响应份额的理论值给定N 个不同的瞬时响应份额预测值，再将堆外探测器信号实测值代入离散状态方程，可以得到N 组铑自给能探测器信号理论值，将理论值与铑自给能探测器信号实测值进行比较，取其中 符合程度最好的一组理论值相应的瞬时响应份额预测值为后续延迟消除所采用的瞬时响 应份额。3. 根据权利要求1或2所述的基于H2/H c?混合滤波的铑自给能探测器信号延迟消 除方法，其特征在于，在有换挡的情况下，还包括按如下的信号处理方法对原始信号进行处 理：在换档区域内，假设中子通量保持不变，然后反推中子通量密度产生的电流信号，再与 探测器实际输出电流相减，得到换挡突变分量；在换档区域外，探测器输出电流减去换挡突 变分量，得到中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除处理。4. 根据权利要求3所述的基于H2/H c?混合滤波的铑自给能探测器信号延迟消除方法， 其特征在于，所述的换档区域设计结构如下： 在换档区域内Gc1S k2)，假设中子通量密度不变，则有： n (k+1) = n (k) (15)可以反推出铑自给能探测器电流信号为： I(k+1) = Ja(k+l)+cn(k+l) (18) 将反推电流（18)当成探测器实际输出电流，然后进行延迟消除； 其中，在换档区域时间边界1^2处，换档引起的电流偏置量可以由下式进行估算： D = !(Ii2Dik1) (19) K&)表示在1^2时刻的探测器实际输出电流； 在换档区域外，将探测器实际输出电流加上式（19)表示的换档引起的电流偏置量，得 到中子通量密度产生的电流信号，然后再对此电流信号进行延迟消除。
【专利摘要】本发明公开了基于H2/H∞混合滤波的铑自给能探测器信号延迟消除方法，包括依次进行的以下步骤：步骤1、建立铑与热中子的核反应模型；步骤2、采用直接变换建立核反应模型对应的离散状态方程；步骤3、确定铑自给能探测器电流的瞬时响应份额；步骤4、利用H2/H∞混合滤波器对铑自给能探测器电流信号作延迟消除。本发明应用时能对铑自给能中子探测器的电流信号进行延迟消除处理，并能有效抑制噪声，使得铑自给能中子探测器在反应堆瞬态工况时也能正常使用，且当输入信号是一个具有有限能量的不确定信号时，铑自给能中子探测器也能正常应用。
【IPC分类】G06F19/00, G01T3/00, G21C17/108
【公开号】CN104882179
【申请号】CN201510166193
【发明人】龚禾林, 陈长, 彭星杰, 赵文博, 刘启伟, 李向阳, 李庆, 于颖锐
【申请人】中国核动力研究设计院
【公开日】2015年9月2日
【申请日】2015年4月9日