一种半导体探测器信号的分析方法和系统与流程

1.本发明属于核辐射探测技术领域，具体涉及一种半导体探测器信号的分析方法和系统。


背景技术：

2.在核物理和粒子物理领域中，需要用电子学方法收集核辐射粒子通过探测器产生的电信号，然后对其进行处理和分析。由于探测器输出信号一般都比较小，通常在探测器与放大器之间加一预放大器，也就是前置放大器。由于输出的幅度反映了输入电荷的大小，且与输入电容无关，故称之为电荷灵敏前置放大器。
3.现有技术中，将半导体检测器输出的电荷信号近似为阶跃信号或者负指数增长的信号，然后基于数字变换构建前置放大器的系统函数表达式，最后根据构建好的放大器对实际探测器信号进行处理。
4.现有技术存在的问题：
5.通过现有技术处理实际探测器信号，处理后的信号存在上下震荡，物理意义存在缺陷。


技术实现要素：

6.针对现有技术中存在的处理后的输出信号存在上下震荡，物理意义存在缺陷的问题，本发明提出了一种半导体探测器信号的分析方法和系统，其目的为：得到稳定的核信号冲激成型结果。
7.为实现上述目的本发明所采用的技术方案是：提供一种半导体探测器信号的分析方法，包括：
8.s1：获取半导体探测器的电荷信号，将所述电荷信号经过前置放大器，生成前置放大器信号，并得到前置放大器信号的数字表达式；
9.s2：将所述前置放大器信号经过c-r系统，得到输出信号和输出信号的数字表达式；
10.s3:构建三级级联变换系统，设输入信号为z[n]，第一级输出信号为y[n]，第二级输出信号为x[n]，第三级输出信号为p[n]，得出三级级联变换系统的信号变换表达式；
[0011]
s4：基于系统级联的特性，将s2中得到的输出信号经过三级级联变换系统进行逐级反卷积，得到各级级联系统具体的数字表达式。
[0012]
较优的，本发明s1中：
[0013]
所述电荷信号如公式(1)所示：
[0014][0015]
其中，a表示脉冲幅度，t表示时间，τ2和3表示时间衰减常数，τ2＞τ3；
[0016]
所述前置放大器的卷积描述如公式(2)所述：
[0017]
y2(t)＝y1(t)*ε(t)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0018]
其中，y2(t)为y1(t)的积分；
[0019]
将所述电荷信号经过所述前置放大器得到所述前置放大器信号如公式(3)所示：
[0020][0021]
其中，a表示脉冲幅度，t表示时间，τ2和3表示时间衰减常数，τ2＞τ3。
[0022]
较优的，本发明s2中，所述c-r系统包括c-r微分电路，对所述c-r微分电路进行分析，得到cr系统的传递函数如公式(4)所示：
[0023][0024]
其中，c为c-r微分电路中的电容值、r为c-r微分电路中的电阻值；
[0025]
整理可得公式(5)：
[0026][0027]
将公式(5)进行拉普拉斯变换得到c-r系统的响应函数的时间域表达式，如公式(6)所示：
[0028]
h(t)＝δ(t)-ae-at ε(t)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0029]
其中，a＝1/rc；
[0030]
将所述前置放大器信号经过c-r系统，得到的输出信号如公式(7)所示：
[0031][0032]
其中，a表示脉冲幅度，t表示时间，τ1＝1/a，τ2和3表示时间衰减常数，τ2＞τ3；
[0033]
整理可得公式(8)：
[0034][0035][0036]
较优的，本发明s3中，第一级为rc-1
系统，第二级为cr-1
系统，第三级为数字微分系统，则第一级输出信号y[n]、第二级输出信号x[n]和第三级输出信号p[n]分别如公式(9)、公式(10)和公式(11)所示：
[0037]
y[n]＝rc-1
(z[n],m)
ꢀꢀꢀꢀ
(9)
[0038]
x[n]＝cr-1
(y[n],m)
ꢀꢀꢀꢀ
(10)
[0039]
p[n]＝dif(x[n])
ꢀꢀꢀꢀꢀꢀꢀ
(11)
[0040]
其中，根据rc-1
算子的数字解推导得到公式(12)：
[0041]
y[n]＝z[n] m
·
(z[n]-z[n-1])
ꢀꢀꢀꢀꢀꢀꢀ
(12)
[0042]
其中，m＝(rc)/δt；
[0043]
根据cr-1
算子的数字解推导得公式(13)：
[0044]
[0045]
其中，
[0046]
将公式(12)代入公式(13)得到公式(14)：
[0047][0048]
对公式(14)两端同时进行数字积分，整理得到公式(15)：
[0049][0050]
由于冲击信号是阶跃信号的一阶导p(n)＝x(n)
′
，所以推导得到公式(16)：
[0051][0052]
较优的，本发明所述s4具体为：
[0053]
将rc-1
系统和cr-1系统的组合定义为inv-cr系统，基于系统级联的概念，将所述输出信号经过inv-cr系统，生成如公式(3)所描述的inv-cr系统信号，得到inv-cr系统表达式如公式(17)所示：
[0054]
x[n]＝k*∑y[n] y[n]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(17)
[0055]
将inv-cr系统得到的inv-cr系统信号经过数字微分系统，得到如公式(1)所描述的数字微分系统信号，得到数字微分系统表达式如公式(18)所示；
[0056]
x[n]＝y[n]-y[n-1]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(18)
[0057]
将数字微分系统信号经过构建的三级级联变换系统，生成冲激信号，得到三级级联变换系统表达式如公式(19)所示：
[0058][0059]
本发明还提出一种半导体探测器信号的分析系统，包括：
[0060]
半导体探测器：用于获取电荷信号；
[0061]
前置放大器：用于生成前置放大器信号，并得到前置放大器信号的数字表达式；
[0062]
c-r系统：将所述前置放大器信号经过c-r系统，得到输出信号和输出信号的数字表达式；
[0063][0064]
三级级联变换系统：包括rc-1
系统、cr-1
系统和数字微分系统，用于将所述输出信号变换为冲激信号。
[0065]
较优的，本发明所述电荷信号如公式(1)所示：
[0066][0067]
其中，a表示脉冲幅度，t表示时间，τ2和3表示时间衰减常数，τ2＞τ3；
[0068]
所述前置放大器的卷积描述如公式(2)所述：
[0069]
y2(t)＝y1(t)*ε(t)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0070]
其中，y2(t)为y1(t)的积分；
[0071]
将所述电荷信号经过所述前置放大器得到所述前置放大器信号如公式(3)所示：
[0072][0073]
其中，a表示脉冲幅度，t表示时间，τ2和3表示时间衰减常数，τ2＞τ3。
[0074]
较优的，本发明c-r系统中，所述c-r系统包括c-r微分电路，对所述c-r微分电路进行分析，得到cr系统的传递函数如公式(4)所示：
[0075][0076]
其中，c为c-r微分电路中的电容值、r为c-r微分电路中的电阻值；
[0077]
整理可得公式(5)：
[0078][0079]
将公式(5)进行拉普拉斯变换得到c-r系统的响应函数的时间域表达式，如公式(6)所示：
[0080]
h(t)＝δ(t)-ae-at ε(t)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0081]
其中，a＝1/rc；
[0082]
将所述前置放大器信号经过c-r系统，得到的输出信号如公式(7)所示：
[0083][0084]
其中，a表示脉冲幅度，t表示时间，τ1＝1/a，τ2和3表示时间衰减常数，τ2＞τ3；
[0085]
整理可得公式(8)：
[0086][0087]-1-1较优的，本发明三级级联变换系统中，第一级为rc系统，第二级为cr系统，第三级为数字微分系统，则第一级输出信号y[n]、第二级输出信号x[n]和第三级输出信号p[n]分别如公式(9)、公式(10)和公式(11)所示：
[0088]
y[n]＝rc-1
(z[n],m)
ꢀꢀꢀꢀꢀꢀ
(9)
[0089]
x[n]＝cr-1
(y[n],m)
ꢀꢀꢀꢀꢀꢀ
(10)
[0090]
p[n]＝dif(x[n])
ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(11)
[0091]
其中，根据rc-1
算子的数字解推导得到公式(12)：
[0092]
y[n]＝z[n] m
·
(z[n]-z[n-1])
ꢀꢀꢀꢀꢀꢀꢀ
(12)
[0093]
其中，m＝(rc)/δt；
[0094]
根据cr-1
算子的数字解推导得公式(13)：
[0095][0096]
其中，
[0097]
将公式(12)代入公式(13)得到公式(14)：
[0098]
[0099]
对公式(14)两端同时进行数字积分，整理得到公式(15)：
[0100][0101]
由于冲激信号是阶跃信号的一阶导p(n)＝x(n)
′
，所以推导得到公式(16)：
[0102][0103]
较优的，本发明将rc-1
系统和cr-1
系统的组合定义为inv-cr系统，基于系统级联的概念，将所述输出信号经过inv-cr系统，生成如公式(3)所描述的inv-cr系统信号，得到inv-cr系统表达式如公式(17)所示：
[0104]
x[n]＝k*∑y[n] y[n]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(17)
[0105]
将inv-cr系统得到的inv-cr系统信号经过数字微分系统，得到如公式(1)所描述的数字微分系统信号，得到数字微分系统表达式如公式(18)所示；
[0106]
x[n]＝y[n]-y[n-1]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(18)
[0107]
将数字微分系统信号经过构建的三级级联变换系统，生成冲激信号，得到三级级联变换系统表达式如公式(19)所示：
[0108][0109]
相比现有技术，本发明的技术方案具有如下优点/有益效果：
[0110]
1.本发明得到的输出信号的数字表达式更加简便，使信号分析更加便捷。
[0111]
2.本发明通过设置三级级联变换系统，能够得到更加稳定的冲激成型结果。
附图说明
[0112]
为了更清楚地说明本发明实施方式的技术方案，下面将对实施方式中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
[0113]
图1是本发明实施例1的积分型电荷灵敏前置放大器连接示意图。
[0114]
图2是本发明实施例1的双指数信号模拟示意图。
[0115]
图3是本发明实施例1的双指数信号积分模拟示意图。
[0116]
图4是本发明实施例1的c-r微分电路示意图。
[0117]
图5是本发明实施例1的前置放大器信号经过c-r系统的模拟示意图。
[0118]
图6是本发明实施例1的三级级联变换系统示意图。
[0119]
图7是本发明实施例1的将输出信号经过三级级联变换系统进行处理示意图。
[0120]
图8是本发明实施例1的探测器信号冲激成型示意图(τ1＝1.875μs,τ2＝0.100μs,τ3＝0.025μs)。
[0121]
图9是本发明实施例1的高计数率时sdd探测器信号冲激成型示意图(τ1＝1.875μs,τ2＝0.100μs,τ3＝0.025μs)。
具体实施方式
[0122]
为使本发明目的、技术方案和优点更加清楚，下面对本发明实施方式中的技术方
案进行清楚、完整地描述，显然，所描述的实施方式是本发明的一部分实施方式，而不是全部的实施方式。基于本发明中的实施方式，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式，都属于本发明保护的范围。因此，以下提供的本发明的实施方式的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施方式。
[0123]
实施例1：
[0124]
如图所示，本发明提供了一种半导体探测器信号的分析方法，包括：
[0125]
s1：获取半导体探测器的电荷信号，将所述电荷信号经过前置放大器，生成前置放大器信号，并得到前置放大器信号的数字表达式；如图1所示，为本实施例1的积分型电荷灵敏前置放大器连接图，s1中：
[0126]
所述电荷信号如公式(1)所示，图2为电荷信号(即双指数信号)的模拟示意图：
[0127][0128]
其中，a表示脉冲幅度，t表示时间，τ2和3表示时间衰减常数，τ2＞τ3；
[0129]
所述前置放大器的卷积描述如公式(2)所述：
[0130]
y2(t)＝y1(t)*ε(t)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0131]
其中，y2(t)为y1(t)的积分；
[0132]
将所述电荷信号经过所述前置放大器得到所述前置放大器信号如公式(3)所示，图3为信号模拟示意图：
[0133][0134]
其中，a表示脉冲幅度，t表示时间，τ2和3表示时间衰减常数，τ2＞τ3。
[0135]
s2：将所述前置放大器信号经过c-r系统，得到输出信号和输出信号的数字表达式；c-r系统是半导体探测器电荷灵敏前置放大器输出信号的后续处理单元，其电路图如图4所示。
[0136]
s2中，所述c-r系统包括c-r微分电路，对所述c-r微分电路进行分析，得到cr系统的传递函数如公式(4)所示：
[0137][0138]
其中，c为c-r微分电路中的电容值、r为c-r微分电路中的电阻值；
[0139]
整理可得公式(5)：
[0140][0141]
将公式(5)进行拉普拉斯变换得到c-r系统的响应函数的时间域表达式，如公式(6)所示：
[0142]
h(t)＝δ(t)-ae-at ε(t)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0143]
其中，a＝1/rc；
[0144]
该公式(6)表明c-r系统的卷积响应函数是负指数信号与冲击信号的线性组合。
[0145]
将所述前置放大器信号经过c-r系统，得到的输出信号如公式(7)所示，如图5所示为输出信号的模拟示意图：
[0146][0147]
其中，a表示脉冲幅度，t表示时间，τ1＝1/a，τ2和3表示时间衰减常数，τ2＞τ3；
[0148]
整理可得公式(8)：
[0149][0150]
将前置放大器信号与c-r系统直接拉普拉斯变换后，可以把得到的变换信号写成公式(20)的形式，公式(20)化简后得到公式(21)的形式，将公式(21)再拉普拉斯逆变换，得到与公式(8)相同的信号。
[0151][0152][0153]
s3:构建三级级联变换系统，设输入信号为z[n]，第一级输出信号为y[n]，第二级输出信号为x[n]，第三级输出信号为p[n]，得出三级级联变换系统的信号变换表达式；s3中，第一级为rc-1
系统，第二级为cr-1
系统，第三级为数字微分系统，则第一级输出信号y[n]、第二级输出信号x[n]和第三级输出信号p[n]分别如公式(9)、公式(10)和公式(11)所示：
[0154]
y[n]＝rc-1
(z[n],m)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(9)
[0155]
x[n]＝cr-1
(y[n],m)
ꢀꢀꢀꢀꢀ
(10)
[0156]
p[n]＝dif(x[n])
ꢀꢀꢀꢀꢀꢀꢀꢀ
(11)
[0157]
其中，根据rc-1
算子的数字解推导得到公式(12)：
[0158]
y[n]＝z[n] m
·
(z[n]-z[n-1])
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(12)
[0159]
其中，m＝(rc)/δt；
[0160]
根据cr-1
算子的数字解推导得公式(13)：
[0161][0162]
其中，
[0163]
将公式(12)代入公式(13)得到公式(14)：
[0164][0165]
对公式(14)两端同时进行数字积分，整理得到公式(15)：
[0166][0167]
公式(15)描述的是双指数信号到阶跃信号数字积分的数字变换，
[0168]
由于冲击信号是阶跃信号的一阶导p(n)＝x(n)
′
，所以推导得到公式(16)：
[0169][0170]
公式(16)即为三级级联系统变换的双指数信号冲击信号变换算法。
[0171]
s4：基于系统级联的特性和概念，将s2中得到的输出信号经过三级级联变换系统进行逐级反卷积，得到各级级联系统具体的数字表达式。s4具体为：
[0172]
将rc-1
系统和cr-1
系统的组合定义为inv-cr系统，基于系统级联的概念，将所述输出信号经过inv-cr系统，生成如公式(3)所描述的inv-cr系统信号，得到inv-cr系统表达式如公式(17)所示：
[0173]
x[n]＝k*∑y[n] y[n]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(17)
[0174]
将inv-cr系统得到的inv-cr系统信号经过数字微分系统，得到如公式(1)所描述的数字微分系统信号，得到数字微分系统表达式如公式(18)所示；
[0175]
x[n]＝y[n]-y[n-1]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(18)
[0176]
将数字微分系统信号经过构建的三级级联变换系统，生成冲激信号，得到三级级联变换系统表达式如公式(19)所示：
[0177][0178]
本发明还提出一种半导体探测器信号的分析系统，包括：
[0179]
半导体探测器：用于获取电荷信号；
[0180]
前置放大器：用于生成前置放大器信号，并得到前置放大器信号的数字表达式；
[0181]
所述电荷信号如公式(1)所示：
[0182][0183]
其中，a表示脉冲幅度，t表示时间，τ2和3表示时间衰减常数，τ2＞τ3；
[0184]
所述前置放大器的卷积描述如公式(2)所述：
[0185]
y2(t)＝y1(t)*ε(t)
ꢀꢀꢀꢀꢀꢀ
(2)
[0186]
其中，y2(t)为y1(t)的积分；
[0187]
将所述电荷信号经过所述前置放大器得到所述前置放大器信号如公式(3)所示：
[0188][0189]
其中，a表示脉冲幅度，t表示时间，τ2和3表示时间衰减常数，τ2＞τ3。
[0190]
c-r系统：将所述前置放大器信号经过c-r系统，得到输出信号和输出信号的数字表达式；本发明c-r系统中，所述c-r系统包括c-r微分电路，对所述c-r微分电路进行分析，得到cr系统的传递函数如公式(4)所示：
[0191][0192]
其中，c为c-r微分电路中的电容值、r为c-r微分电路中的电阻值；
[0193]
整理可得公式(5)：
[0194][0195]
将公式(5)进行拉普拉斯变换得到c-r系统的响应函数的时间域表达式，如公式(6)所示：
[0196]
h(t)＝δ(t)-ae-at ε(t)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0197]
其中，a＝1/rc；
[0198]
将所述前置放大器信号经过c-r系统，得到的输出信号如公式(7)所示：
[0199][0200]
其中，a表示脉冲幅度，t表示时间，τ1＝1/a，τ2和3表示时间衰减常数，τ2＞τ3；
[0201]
整理可得公式(8)：
[0202][0203]
三级级联变换系统：包括rc-1
系统、cr-1
系统和数字微分系统，用于将所述输出信号变换为冲激信号。
[0204]
本发明三级级联变换系统中，第一级为rc-1
系统，第二级为cr-1
系统，第三级为数字微分系统，则第一级输出信号y[n]、第二级输出信号x[n]和第三级输出信号p[n]分别如公式(9)、公式(10)和公式(11)所示：
[0205]
y[n]＝rc-1
(z[n],m)
ꢀꢀꢀꢀ
(9)
[0206]
x[n]＝cr-1
(y[n],m)
ꢀꢀꢀꢀ
(10)
[0207]
p[n]＝dif(x[n])
ꢀꢀꢀꢀꢀꢀꢀ
(11)
[0208]
其中，根据rc-1
算子的数字解推导得到公式(12)：
[0209]
y[n]＝z[n] m
·
(z[n]-z[n-1])
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(12)
[0210]
其中，m＝(rc)/δt；
[0211]
根据cr-1
算子的数字解推导得公式(13)：
[0212][0213]
其中，
[0214]
将公式(12)代入公式(13)得到公式(14)：
[0215][0216]
对公式(14)两端同时进行数字积分，整理得到公式(15)：
[0217][0218]
由于冲激信号是阶跃信号的一阶导p(n)＝x(n)
′
，所以推导得到公式(16)：
[0219][0220]
然后将rc-1
系统和cr-1
系统的组合定义为inv-cr系统，基于系统级联的概念，将所述输出信号经过inv-cr系统，生成如公式(3)所描述的inv-cr系统信号，得到inv-cr系统表达式如公式(17)所示：
[0221]
x[n]＝k*∑y[n] y[n]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(17)
[0222]
其中k为inv-cr系统相关的系数，将inv-cr系统得到的inv-cr系统信号经过数字微分系统，得到如公式(1)所描述的数字微分系统信号，得到数字微分系统表达式如公式(18)所示；
[0223]
x[n]＝y[n]-y[n-1]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(18)
[0224]
将数字微分系统信号经过构建的三级级联变换系统，生成冲激信号，得到三级级联变换系统表达式如公式(19)所示：
[0225][0226]
以上仅是本发明的优选实施方式，应当指出的是，上述优选实施方式不应视为对本发明的限制，本发明的保护范围应当以权利要求所限定的范围为准。对于本技术领域的普通技术人员来说，在不脱离本发明的精神和范围内，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。