Digitalization of inverting filter shaping circuit for nuclear pulse signals

Digitalization of inverting filter shaping circuit for nuclear pulse signals 1 2 corresp first-author Corresponding author, zhanghq821@163.com Corresponding author, zhanghq821@163.com zhanghq821@163.com https://orcid.org/0000-0001-9333-8585 2 2 https://orcid.org/0000-0002-5935-7747 2 State Key Laboratory of Nuclear Resources and Environment, East China University of Technology, Nanchang, 330013, China State Key Laboratory of Nuclear Resources and Environment, East China University of Technology, Nanchang, 330013, China School of Nuclear Science and Engineering, East China University of Technology, Nanchang 330013,China School of Nuclear Science and Engineering, East China University of Technology, Nanchang 330013,China In the design of filter shaping circuits for nuclear pulse signals, inverting filter shaping circuits perform better than non-inverting filter shaping circuits. Because these circuits facilitate changing the phase of a pulse signal, they are widely used in processing nuclear pulse signals. In this study, the transfer functions of four types of inverting filter shaping circuits, namely, the common inverting filter shaping, improved inverting filter shaping, multiple feedback low-pass filter shaping, and third-order multiple feedback low-pass filter shaping, in the Laplacian domain, are derived. We establish the numerical recursive function models and digitalize the four circuits, obtain the transfer functions in the Z domain, and analyze the filter performance and amplitude-frequency response characteristics in the frequency domain. Based on the actual nuclear pulse signal of the Si-PIN detector, we realize four types of inverting digital shaping. The results show that under the same shaping parameters, the common inverting digital shaping has better amplitude extraction characteristics, the third-order multiple feedback low-pass digital shaping has better noise suppression performance, and the multiple feedback digital shaping takes into account both pulse amplitude extraction and noise suppression performance. Nuclear pulse signal Inverting filter shaping circuit Digital shaping Amplitude-frequency response 1 Introduction 1 When nuclear pulse signals are processed, the output nuclear pulse signal of the pre-amplified needs to be filtered and shaped to improve the signal-to-noise ratio (SNR) of the system and change the pulse shape for subsequent processing. Several studies have focused on analog non-inverting filter shaping circuits for nuclear pulse signals. Evariste et al. designed and developed a set of amplifier filter shaping circuits for Si (Li), CdZnTe and CsI detectors, which included a fourth-order integration stage and a differentiation stage [ 1 1 ]. Gan et al. designed a low-noise nuclear pulse signal processing system, in which CR-RC was used to realize filter shaping [ 2 2 ]. Gao et al. designed CR-RC and Sallen-Key circuits based on CdZnTe and Si-PIN detectors, which were used to shaped into quasi-Gaussian [ 3 3 ]. Zhang et al. developed a low-noise integrated circuit based on CR-RC 2 , which exhibited excellent filtering effects [ 4 4 ]. The low-pass feedback terminal of the charge-sensitive preamplifier, based on RC low-pass feedback, was replaced by a junction field-effect transistor and an RC low-pass network for the high-resistance feedback resistor. The noise of the feedback termin al was considerably smaller than the thermal noise introduced by the high resistance feedback resistor [ 5 5 ]. In the work of Zeng et al. [ 6 6 ], a CR low-pass filter circuit was used to process the output signal of the high-purity germanium detector to filter out low-frequency interference. This circuit was equipped with a non-inverting follower and an adjustable gain amplifier circuit to amplify and filter the nuclear pulse signal. Wang et al. developed a compact 16-channel integrated charge-sensitive preamplifier to support the large-scale detector array used in modern nuclear physics experiments. A first-order RC filter circuit was employed to reduce the noise generated by the high voltage cable. Meanwhile, the resistor–capacitor feedback network formed a charge integrating and bleeding circuit that determined the energy sensitivity and pulse decay time [ 7 7 ]. The above circuits were designed and analyzed for their performance as analog non-inverting amplification filter shaping circuits for nuclear pulse signals. We summarize the major works reporting the digital non-inverting filter shaping method of nuclear pulse signals as follows. Zhang et al. used the digital trapezoidal shaping method to process the simulated nuclear signal and the actual sampled nuclear signal to obtain the best digital shaping parameters [ 8 8 ]. Zhang et al. [ 9 9 ] derived the transfer function of the filter shaping circuit based on Sallen-Key in the Laplacian and Z domains and improved its numerical recursive function model at various values of resistance and capacitance. A comparative analysis of performance indicators, such as amplitude extraction and noise suppression, led to the optimal design model and filter shaping parameters. The CR-RC m digital shaping method of nuclear pulse signal processing was implemented based on FPGA [ 10 10 ]. Zhou et al. derived the Sallen-Key numerical recursive function model using Kirchhoff’s current law. This was also the earliest digital shaping processing of nuclear pulse signals in the temporal domain [ 11 11 ]. Zhou et al. [ 12 12 ] also used the results of earlier studies to improve the digital Sallen-Key recursive function model and introduce the amplitude adjustable parameter. The transfer function of the digital Sallen-Key circuit was derived in the Z domain and the nuclear pulse signal was convolved with the impulse response function. However, the result, when applied to the processing of the double-exponential nuclear pulse signal, did not improve the real-time performance. This limited the further use of this method [ 13 13 ]. Zhang et al. constructed a set of digital processing platforms for nuclear pulse signals based on MATLAB, which realized the digital shaping processing of simulated nuclear signals and actual sampled nuclear signals under different shaping methods and shaping parameters [ 14 14 ]. In a previous study, the quality factor and cut-off frequency of the digital Sallen-Key shaping method were introduced and simultaneously compared with the trapezoidal shaping results [ 15 15 ]. It was observed that digital Gaussian shaping based on the Sallen-Key method exhibited better noise suppression performance, which was beneficial for pulse amplitude extraction. In contrast, digital Gaussian shaping based on CR-RC m had better counting rate characteristics, which aided pulse pile-up identification [ 16 16 ]. Zhang et al. studied the digitization of the pole-zero cancellation and Sallen-Key filter shaping circuits and obtained the numerical recursive function model and the amplitude-frequency response curve, respectively [ 17 17 ]. Liu et al. applied the CR-RC m digital shaping method to the γ -spectrum measurement system and obtained the relationship between the shaping parameters and energy resolution [ 18 18 ]. Regadio et al. designed and developed a set of adjustable digital shapers for particle physics, which could adjust the shaping algorithm and shaping parameters according to actual requirements, and implemented the algorithm in FPGA [ 19 19 ]. Valentin used the unfolding-synthesis technique in digital pulse processing systems by digitizing analog signals according to the unfolding system and converting the nuclear signal into the desired pulse shapes [ 20 20 ]. The adaptive digital pulse shaping technique was investigated for real-time signal processing, and a digital Sallen–Key low-pass filter was implemented to enhance the SNR and reduce baseline drifting in this system [ 21 21 ]. In comparison with the non-inverting digital filter shaping method, the inverting digital shaping method exhibits a better SNR, which aids in pulse amplitude extraction under the same shaping parameters. The digital model functions of the inverting filter circuits, which were obtained using the numerical solution of differential equations, confirmed that the inverting digital shaping method exhibited better noise suppression performance [ 22 22 ]. However, the amplitude-frequency response was not analyzed, and hence, the numerical recursive function model of the serialized inverting filter shaping circuits was not analyzed and studied. In this study, the digitalization of common inverting filter shaping circuits, improved inverting filter shaping circuits, multiple feedback low-pass filter shaping circuits and third-order multiple feedback low-pass filter shaping circuits is realized based on the preceding studies. The transfer function in the Laplacian domain, the amplitude-frequency response curve in the Z domain, and the numerical recursive function model are established. These are highly significant for the digital shaping processing of nuclear pulses and the optimal choice of parameters for nuclear pulse signal inverting filter shaping circuits. 2 Digitalization of inverting filter shaping circuits 1 Inverting filter shaping circuits for nuclear pulse signal processing include common inverting filter shaping circuits, improved inverting filter shaping circuits, multiple feedback low-pass filter shaping circuits, and third-order multiple feedback low-pass filter shaping circuits. 2.1 Digitalization of the common inverting filter shaping circuit 2 The common inverting filter shaping circuit is composed of an operational amplifier (OPAMP) A, resistors R 1 and R 2 and capacitor C. The circuit is shown in Fig. 1 Fig. 1 . As shown in Fig. 1 Fig. 1 , V i and V o are the input and output voltages, respectively, of the common inverting filter shaping circuit. The transfer function of the common inverting filter shaping circuit in the Laplacian domain is shown in Eq. (1). <math display="block" id="M1"><mrow><mi>H</mi><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msub><mi>V</mi><mtext>o</mtext></msub><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow><mrow><msub><mi>V</mi><mtext>i</mtext></msub><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow></mfrac><mo>=</mo><mfrac><mrow><mo>−</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>1</mn></msub><mi>C</mi></mrow></mfrac></mrow><mrow><mi>s</mi><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>2</mn></msub><mi>C</mi></mrow></mfrac></mrow></mfrac></mrow></math> When V i is converted into x i , V o is converted into y i , where R 1 = R 2 = R . According to differential numerical analysis, Eq. (1) leads to Eq. (2): <math display="block" id="M2"><mrow><mi>R</mi><mi>C</mi><mo>⋅</mo><mfrac><mrow><msub><mi>y</mi><mi>i</mi></msub><mo>−</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>1</mn></mrow></msub></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac><mo>+</mo><msub><mi>y</mi><mi>i</mi></msub><mo>=</mo><mo>−</mo><msub><mi>x</mi><mi>i</mi></msub><mn>.</mn></mrow></math> When x i is converted into a sequence of numbers, x [ n ], y i is converted into a sequence of numbers, y [ n ]. Thus, Eq. (3) can be obtained: <math display="block" id="M3"><mrow><mi>y</mi><mo stretchy="false">[</mo><mi>n</mi><mo stretchy="false">]</mo><mo>=</mo><mfrac><mrow><mi>k</mi><mo>⋅</mo><mi>y</mi><mo stretchy="false">[</mo><mi>n</mi><mo>−</mo><mn>1</mn><mo stretchy="false">]</mo><mo>−</mo><mi>x</mi><mo stretchy="false">[</mo><mi>n</mi><mo stretchy="false">]</mo></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mn>.</mn></mrow></math> When k = RC / Δ t (Δ t is sampling interval of ADC), the transfer function of the common inverting filter shaping circuit in the Z domain is as shown in Eq. (4): <math display="block" id="M4"><mrow><mi>H</mi><mo stretchy="false">(</mo><mi>z</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mo>−</mo><mn>1</mn></mrow><mrow><mi>k</mi><mo>+</mo><mn>1</mn><mo>−</mo><mi>k</mi><mo>⋅</mo><msup><mi>z</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow></mfrac><mn>.</mn></mrow></math> The frequency response function is: <math display="block" id="M5"><mrow><mi>H</mi><mo stretchy="false">(</mo><msup><mi>e</mi><mrow><mi>j</mi><mi>w</mi></mrow></msup><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mo>−</mo><mn>1</mn></mrow><mrow><mi>k</mi><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>j</mi><mi>w</mi></mrow></msup><mo stretchy="false">)</mo><mo>+</mo><mn>1</mn></mrow></mfrac><mn>.</mn></mrow></math> The modulus of Eq. (5) is shown in Eq. (6). <math display="block" id="M6"><mrow><mo>|</mo><mi>H</mi><mo stretchy="false">(</mo><msup><mi>e</mi><mrow><mi>j</mi><mi>w</mi></mrow></msup><mo stretchy="false">)</mo><mo>|</mo><mo>=</mo><mfrac><mn>1</mn><mrow><msqrt><mrow><mn>1</mn><mo>+</mo><mn>2</mn><mi>k</mi><mo>+</mo><mn>2</mn><msup><mi>k</mi><mn>2</mn></msup><mo>−</mo><mn>2</mn><mi>k</mi><mo stretchy="false">(</mo><mn>1</mn><mo>+</mo><mi>k</mi><mo stretchy="false">)</mo><mo>⋅</mo><mi>cos</mi><mo stretchy="false">(</mo><mi>w</mi><mo stretchy="false">)</mo></mrow></msqrt></mrow></mfrac><mn>.</mn></mrow></math> Eq. (6) may be used to obtain the amplitude- frequency response curves. 2.2 Digitalization of the improved inverting filter shaping circuit 2 The improved inverting filter shaping circuit is composed of an OPAMP A, resistors R 1 , R 2 and R 3 and capacitors C 1 and C 2 . The corresponding circuit is shown in Fig. 2 Fig. 2 . The transfer function in the Laplacian domain of the improved inverting filter shaping circuit shown in Fig. 2 Fig. 2 is given by Eq. (7). <math display="block" id="M7"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>H</mi><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow></mtd><mtd columnalign="left"><mrow><mtext>​</mtext><mo>=</mo><mtext>​</mtext><mfrac><mrow><msub><mi>V</mi><mtext>o</mtext></msub><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow><mrow><msub><mi>V</mi><mtext>i</mtext></msub><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow></mfrac></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mtext>​</mtext><mo>=</mo><mtext>​</mtext><mfrac><mrow><mo>−</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>1</mn></msub><msub><mi>C</mi><mn>1</mn></msub></mrow></mfrac><mi>s</mi><mo>−</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>1</mn></msub><msub><mi>R</mi><mn>2</mn></msub><msub><mi>C</mi><mn>1</mn></msub><msub><mi>C</mi><mn>2</mn></msub></mrow></mfrac></mrow><mrow><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>2</mn></msub><msub><mi>C</mi><mn>2</mn></msub></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>2</mn></msub><msub><mi>C</mi><mn>1</mn></msub></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>3</mn></msub><msub><mi>C</mi><mn>2</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mi>s</mi><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>2</mn></msub><msub><mi>R</mi><mn>3</mn></msub><msub><mi>C</mi><mn>1</mn></msub><msub><mi>C</mi><mn>2</mn></msub></mrow></mfrac></mrow></mfrac></mrow></mtd></mtr></mtable></mrow></math> When R 1 = R 2 = R 3 = R , C 1 = C 2 = C , according to the differential numerical analysis, Eq. (7) changes into Eq. (8). <math display="block" id="M8"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><msup><mrow><mo stretchy="false">(</mo><mi>R</mi><mi>C</mi><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>⋅</mo><mfrac><mrow><mfrac><mrow><msub><mi>y</mi><mi>i</mi></msub><mtext>​</mtext><mo>−</mo><mtext>​</mtext><msub><mi>y</mi><mrow><mi>i</mi><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>1</mn></mrow></msub></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mfrac><mrow><msub><mi>y</mi><mrow><mi>i</mi><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>1</mn></mrow></msub><mtext>​</mtext><mo>−</mo><mtext>​</mtext><msub><mi>y</mi><mrow><mi>i</mi><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>2</mn></mrow></msub></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>3</mn><mi>R</mi><mi>C</mi><mo>⋅</mo><mfrac><mrow><msub><mi>y</mi><mi>i</mi></msub><mtext>​</mtext><mo>−</mo><mtext>​</mtext><msub><mi>y</mi><mrow><mi>i</mi><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>1</mn></mrow></msub></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac><mtext>​</mtext><mo>+</mo><mtext>​</mtext><msub><mi>y</mi><mi>i</mi></msub></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>​</mtext><mo>=</mo><mtext>​</mtext></mrow></mtd><mtd columnalign="left"><mrow><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mi>R</mi><mi>C</mi><mo>⋅</mo><mfrac><mrow><msub><mi>x</mi><mi>i</mi></msub><mtext>​</mtext><mo>−</mo><mtext>​</mtext><msub><mi>x</mi><mrow><mi>i</mi><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>1</mn></mrow></msub></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac><mtext>​</mtext><mo>−</mo><mtext>​</mtext><msub><mi>x</mi><mi>i</mi></msub></mrow></mtd></mtr></mtable></mrow></math> Then Eq. (9) can be obtained as <math display="block" id="M9"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>y</mi><mo stretchy="false">[</mo><mi>n</mi><mo stretchy="false">]</mo><mo>=</mo></mrow></mtd><mtd columnalign="left"><mrow><mfrac><mrow><mo stretchy="false">(</mo><mn>2</mn><msup><mi>k</mi><mn>2</mn></msup><mo>+</mo><mn>3</mn><mi>k</mi><mo stretchy="false">)</mo><mo>⋅</mo><mi>y</mi><mo stretchy="false">[</mo><mi>n</mi><mo>−</mo><mn>1</mn><mo stretchy="false">]</mo><mo>−</mo><msup><mi>k</mi><mn>2</mn></msup><mo>⋅</mo><mi>y</mi><mo stretchy="false">[</mo><mi>n</mi><mo>−</mo><mn>2</mn><mo stretchy="false">]</mo></mrow><mrow><msup><mi>k</mi><mn>2</mn></msup><mo>+</mo><mn>3</mn><mi>k</mi><mo>+</mo><mn>1</mn></mrow></mfrac></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>−</mo><mfrac><mrow><mo>−</mo><mo stretchy="false">(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo><mo>⋅</mo><mi>x</mi><mo stretchy="false">[</mo><mi>n</mi><mo stretchy="false">]</mo><mo>−</mo><mi>k</mi><mo>⋅</mo><mi>x</mi><mo stretchy="false">[</mo><mi>n</mi><mo>−</mo><mn>1</mn><mo stretchy="false">]</mo></mrow><mrow><msup><mi>k</mi><mn>2</mn></msup><mo>+</mo><mn>3</mn><mi>k</mi><mo>+</mo><mn>1</mn></mrow></mfrac></mrow></mtd></mtr></mtable></mrow></math> When k = RC / Δ t , the transfer function of the improved inverting filter shaping circuit in the Z domain is shown in Eq. (10). <math display="block" id="M10"><mrow><mi>H</mi><mo stretchy="false">(</mo><mi>z</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mo>−</mo><mi>k</mi><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msup><mi>z</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo stretchy="false">)</mo><mo>−</mo><mn>1</mn></mrow><mrow><msup><mi>k</mi><mn>2</mn></msup><msup><mrow><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msup><mi>z</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>3</mn><mi>k</mi><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msup><mi>z</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo stretchy="false">)</mo><mo>+</mo><mn>1</mn></mrow></mfrac></mrow></math> The frequency response function is shown in Eq. (11): <math display="block" id="M11"><mrow><mi>H</mi><mo stretchy="false">(</mo><msup><mi>e</mi><mrow><mi>j</mi><mi>w</mi></mrow></msup><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mo>−</mo><mi>k</mi><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>j</mi><mi>w</mi></mrow></msup><mo stretchy="false">)</mo><mo>−</mo><mn>1</mn></mrow><mrow><msup><mi>k</mi><mn>2</mn></msup><msup><mrow><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>j</mi><mi>w</mi></mrow></msup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>3</mn><mi>k</mi><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>j</mi><mi>w</mi></mrow></msup><mo stretchy="false">)</mo><mo>+</mo><mn>1</mn></mrow></mfrac></mrow></math> The modulus of Eq. (11) is shown in Eq. (12). <math display="block" id="M12"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>|</mo><mi>H</mi><mo stretchy="false">(</mo><msup><mi>e</mi><mrow><mi>j</mi><mi>w</mi></mrow></msup><mo stretchy="false">)</mo><mo>|</mo></mrow></mtd><mtd columnalign="left"><mrow><mo>=</mo><msqrt><mrow><mfrac><mrow><mn>1</mn><mo>+</mo><mn>2</mn><mi>k</mi><mo>+</mo><mn>2</mn><msup><mi>k</mi><mn>2</mn></msup><mo>−</mo><mn>2</mn><mi>k</mi><mo stretchy="false">(</mo><mn>1</mn><mo>+</mo><mi>k</mi><mo stretchy="false">)</mo><mo>⋅</mo><mi>c</mi><mi>o</mi><mi>s</mi><mo stretchy="false">(</mo><mi>w</mi><mo stretchy="false">)</mo></mrow><mrow><mi>A</mi><mo>−</mo><mi>B</mi><mo>⋅</mo><mi>c</mi><mi>o</mi><mi>s</mi><mo stretchy="false">(</mo><mi>w</mi><mo stretchy="false">)</mo><mo>+</mo><mi>C</mi><mo>⋅</mo><mi>c</mi><mi>o</mi><mi>s</mi><mo stretchy="false">(</mo><mn>2</mn><mi>w</mi><mo stretchy="false">)</mo></mrow></mfrac></mrow></msqrt></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mi>A</mi><mo>=</mo><mn>1</mn><mo>+</mo><mn>6</mn><mi>k</mi><mo>+</mo><mn>20</mn><msup><mi>k</mi><mn>2</mn></msup><mo>+</mo><mn>18</mn><msup><mi>k</mi><mn>3</mn></msup><mo>+</mo><mn>6</mn><msup><mi>k</mi><mn>4</mn></msup></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mi>B</mi><mo>=</mo><mn>2</mn><mi>k</mi><mo stretchy="false">(</mo><mn>3</mn><mo>+</mo><mn>11</mn><mi>k</mi><mo>+</mo><mn>12</mn><msup><mi>k</mi><mn>2</mn></msup><mo>+</mo><mn>4</mn><msup><mi>k</mi><mn>3</mn></msup><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mi>C</mi><mo>=</mo><mn>2</mn><msup><mi>k</mi><mn>2</mn></msup><mo stretchy="false">(</mo><mn>1</mn><mo>+</mo><mn>3</mn><mi>k</mi><mo>+</mo><msup><mi>k</mi><mn>2</mn></msup><mo stretchy="false">)</mo></mrow></mtd></mtr></mtable></mrow></math> The amplitude-frequency response curves can be obtained using Eq. (12). 2.3 Digitalization of the multiple feedback low-pass filter shaping circuit 2 The multiple feedback low-pass filter shaping circuit consists of an OPAMP A, resistors R 1 , R 2 and R 3 and capacitors C 1 and C 2 . The corresponding circuit is shown in Fig. 3 Fig. 3 . In the Laplacian domain, the transfer function of the multiple feedback low-pass filter shaping circuit shown in Fig. 3 Fig. 3 is given by Eq. (13). <math display="block" id="M13"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>H</mi><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow></mtd><mtd columnalign="left"><mrow><mo>=</mo><mfrac><mrow><msub><mi>V</mi><mtext>o</mtext></msub><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow><mrow><msub><mi>V</mi><mtext>i</mtext></msub><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow></mfrac></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>=</mo><mfrac><mrow><mo>−</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>1</mn></msub><msub><mi>C</mi><mn>1</mn></msub><msub><mi>R</mi><mn>2</mn></msub><msub><mi>C</mi><mn>2</mn></msub></mrow></mfrac></mrow><mrow><msup><mi>s</mi><mn>2</mn></msup><mo>+</mo><mi>s</mi><mfrac><mn>1</mn><mrow><msub><mi>C</mi><mn>1</mn></msub></mrow></mfrac><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>1</mn></msub></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>2</mn></msub></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>3</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>2</mn></msub><msub><mi>C</mi><mn>1</mn></msub><msub><mi>R</mi><mn>3</mn></msub><msub><mi>C</mi><mn>2</mn></msub></mrow></mfrac></mrow></mfrac></mrow></mtd></mtr></mtable></mrow></math> When R 1 = R 2 = R 3 = R , C 1 = C 2 = C , according to differential numerical analysis, Eq. (13) can be changed into Eq. (14). <math display="block" id="M14"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><msup><mrow><mo stretchy="false">(</mo><mi>R</mi><mi>C</mi><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>⋅</mo><mfrac><mrow><mfrac><mrow><msub><mi>y</mi><mi>i</mi></msub><mo>−</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>1</mn></mrow></msub></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac><mo>−</mo><mfrac><mrow><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>1</mn></mrow></msub><mo>−</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>2</mn></mrow></msub></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mo>+</mo><mn>3</mn><mi>R</mi><mi>C</mi><mo>⋅</mo><mfrac><mrow><msub><mi>y</mi><mi>i</mi></msub><mo>−</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>1</mn></mrow></msub></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac><mo>+</mo><msub><mi>y</mi><mi>i</mi></msub><mo>=</mo><mo>−</mo><msub><mi>x</mi><mi>i</mi></msub></mrow></mtd></mtr></mtable></mrow></math> Further derivation yields Eq. (15) as <math display="block" id="M15"><mrow><mi>y</mi><mo stretchy="false">[</mo><mi>n</mi><mo stretchy="false">]</mo><mo>=</mo><mfrac><mrow><mo stretchy="false">(</mo><mn>2</mn><msup><mi>k</mi><mn>2</mn></msup><mo>+</mo><mn>3</mn><mi>k</mi><mo stretchy="false">)</mo><mo>⋅</mo><mi>y</mi><mo stretchy="false">[</mo><mi>n</mi><mo>−</mo><mn>1</mn><mo stretchy="false">]</mo><mo>−</mo><msup><mi>k</mi><mn>2</mn></msup><mo>⋅</mo><mi>y</mi><mo stretchy="false">[</mo><mi>n</mi><mo>−</mo><mn>2</mn><mo stretchy="false">]</mo><mo>−</mo><mi>x</mi><mo stretchy="false">[</mo><mi>n</mi><mo stretchy="false">]</mo></mrow><mrow><msup><mi>k</mi><mn>2</mn></msup><mo>+</mo><mn>3</mn><mi>k</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mn>.</mn></mrow></math> When k = RC / Δ t , the transfer function of the multiple feedback low-pass filter shaping circuit in the Z domain is shown in Eq. (16) <math display="block" id="M16"><mrow><mi>H</mi><mo stretchy="false">(</mo><msup><mi>e</mi><mrow><mi>j</mi><mi>w</mi></mrow></msup><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mo>−</mo><mn>1</mn></mrow><mrow><msup><mi>k</mi><mn>2</mn></msup><msup><mrow><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msup><mi>z</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>3</mn><mi>k</mi><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msup><mi>z</mi><mrow><mo>−</mo><mn>1</mn></mrow></msup><mo stretchy="false">)</mo><mo>+</mo><mn>1</mn></mrow></mfrac><mn>.</mn></mrow></math> The frequency response function is <math display="block" id="M17"><mrow><mi>H</mi><mo stretchy="false">(</mo><msup><mi>e</mi><mrow><mi>j</mi><mi>w</mi></mrow></msup><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mo>−</mo><mn>1</mn></mrow><mrow><msup><mi>k</mi><mn>2</mn></msup><msup><mrow><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>j</mi><mi>w</mi></mrow></msup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mo>+</mo><mn>3</mn><mi>k</mi><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>j</mi><mi>w</mi></mrow></msup><mo stretchy="false">)</mo><mo>+</mo><mn>1</mn></mrow></mfrac><mn>.</mn></mrow></math> The modulus of Eq. (17) is given by Eq. (18). <math display="block" id="M18"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mo>|</mo><mi>H</mi><mo stretchy="false">(</mo><msup><mi>e</mi><mrow><mi>j</mi><mi>w</mi></mrow></msup><mo stretchy="false">)</mo><mo>|</mo></mrow></mtd><mtd columnalign="left"><mrow><mo>=</mo><msqrt><mrow><mfrac><mn>1</mn><mrow><mi>A</mi><mo>−</mo><mi>B</mi><mo>⋅</mo><mi>cos</mi><mo stretchy="false">(</mo><mi>w</mi><mo stretchy="false">)</mo><mo>+</mo><mi>C</mi><mo>⋅</mo><mi>cos</mi><mo stretchy="false">(</mo><mn>2</mn><mi>w</mi><mo stretchy="false">)</mo></mrow></mfrac></mrow></msqrt></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mi>A</mi><mo>=</mo><mn>1</mn><mo>+</mo><mn>6</mn><mi>k</mi><mo>+</mo><mn>20</mn><msup><mi>k</mi><mn>2</mn></msup><mo>+</mo><mn>18</mn><msup><mi>k</mi><mn>3</mn></msup><mo>+</mo><mn>6</mn><msup><mi>k</mi><mn>4</mn></msup></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mi>B</mi><mo>=</mo><mn>2</mn><mi>k</mi><mo stretchy="false">(</mo><mn>3</mn><mo>+</mo><mn>11</mn><mi>k</mi><mo>+</mo><mn>12</mn><msup><mi>k</mi><mn>2</mn></msup><mo>+</mo><mn>4</mn><msup><mi>k</mi><mn>3</mn></msup><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mi>C</mi><mo>=</mo><mn>2</mn><msup><mi>k</mi><mn>2</mn></msup><mo stretchy="false">(</mo><mn>1</mn><mo>+</mo><mn>3</mn><mi>k</mi><mo>+</mo><msup><mi>k</mi><mn>2</mn></msup><mo stretchy="false">)</mo></mrow></mtd></mtr></mtable></mrow></math> The amplitude- frequency response curves can be obtained using Eq. (18). 2.4 Digitalization of the third-order multiple feedback low-pass filter shaping circuit 2 The third-order multiple feedback low-pass filter shaping circuit consists of an OPAMP A, resistors R 1 , R 2 , R 3 and R 4 and capacitors C 1 , C 2 and C 3 . The corresponding circuit is shown in Fig. 4 Fig. 4 . The transfer function in the Laplacian domain of the third-order multiple feedback low-pass filter shaping circuit shown in Fig. 4 Fig. 4 is given by Eq.(19). <math display="block" id="M19"><mtable columnalign="left"><mtr><mtd><mi>H</mi><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><msub><mi>V</mi><mtext>o</mtext></msub><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow><mrow><msub><mi>V</mi><mtext>i</mtext></msub><mo stretchy="false">(</mo><mi>s</mi><mo stretchy="false">)</mo></mrow></mfrac></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mrow><mo>−</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>1</mn></msub><msub><mi>R</mi><mn>2</mn></msub><msub><mi>R</mi><mn>3</mn></msub><msub><mi>C</mi><mn>1</mn></msub><msub><mi>C</mi><mn>2</mn></msub><msub><mi>C</mi><mn>3</mn></msub></mrow></mfrac></mrow><mrow><msup><mi>s</mi><mn>3</mn></msup><mo>+</mo><msup><mi>s</mi><mn>2</mn></msup><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>1</mn></msub><msub><mi>C</mi><mn>1</mn></msub></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>2</mn></msub><msub><mi>C</mi><mn>1</mn></msub></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>2</mn></msub><msub><mi>C</mi><mn>2</mn></msub></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>3</mn></msub><msub><mi>C</mi><mn>2</mn></msub></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>4</mn></msub><msub><mi>C</mi><mn>2</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mo>+</mo><mo>⋅</mo><mo>⋅</mo><mo>⋅</mo></mrow></mfrac></mtd></mtr><mtr><mtd><mi>s</mi><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>1</mn></msub><msub><mi>R</mi><mn>2</mn></msub><msub><mi>C</mi><mn>1</mn></msub><msub><mi>C</mi><mn>2</mn></msub></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>1</mn></msub><msub><mi>R</mi><mn>3</mn></msub><msub><mi>C</mi><mn>1</mn></msub><msub><mi>C</mi><mn>2</mn></msub></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>2</mn></msub><msub><mi>R</mi><mn>3</mn></msub><msub><mi>C</mi><mn>1</mn></msub><msub><mi>C</mi><mn>2</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mo>+</mo><mo>⋅</mo><mo>⋅</mo><mo>⋅</mo></mtd></mtr><mtr><mtd><mi>s</mi><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>2</mn></msub><msub><mi>R</mi><mn>4</mn></msub><msub><mi>C</mi><mn>1</mn></msub><msub><mi>C</mi><mn>2</mn></msub></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>1</mn></msub><msub><mi>R</mi><mn>4</mn></msub><msub><mi>C</mi><mn>1</mn></msub><msub><mi>C</mi><mn>2</mn></msub></mrow></mfrac><mo>+</mo><mfrac><mn>1</mn><mrow><msub><mi>R</mi><mn>3</mn></msub><msub><mi>R</mi><mn>4</mn></msub><msub><mi>C</mi><mn>2</mn></msub><msub><mi>C</mi><mn>3</mn></msub></mrow></mfrac></mrow><mo>)</mo></mrow><mo>+</mo><mo>⋅</mo><mo>⋅</mo><mo>⋅</mo></mtd></mtr><mtr><mtd><mfrac><mrow><msub><mi>R</mi><mn>1</mn></msub><mo>+</mo><msub><mi>R</mi><mn>2</mn></msub></mrow><mrow><msub><mi>R</mi><mn>1</mn></msub><msub><mi>R</mi><mn>2</mn></msub><msub><mi>R</mi><mn>3</mn></msub><msub><mi>R</mi><mn>4</mn></msub><msub><mi>C</mi><mn>1</mn></msub><msub><mi>C</mi><mn>2</mn></msub><msub><mi>C</mi><mn>3</mn></msub></mrow></mfrac></mtd></mtr></mtable></math> When R 1 = R 2 = R 3 = R 4 = R , C 1 = C 2 = C 3 = C , Eq. (19) can be changed into Eq. (20) using differential numerical analysis. <math display="block" id="M20"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><msup><mrow><mo stretchy="false">(</mo><mi>R</mi><mi>C</mi><mo stretchy="false">)</mo></mrow><mn>3</mn></msup><mtext>​</mtext><mo>⋅</mo><mtext>​</mtext><mfrac><mrow><mfrac><mrow><mfrac><mrow><msub><mi>y</mi><mi>i</mi></msub><mo>−</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>1</mn></mrow></msub></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mfrac><mrow><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>1</mn></mrow></msub><mo>−</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>2</mn></mrow></msub></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mfrac><mrow><mfrac><mrow><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>2</mn></mrow></msub><mo>−</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>3</mn></mrow></msub></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mfrac><mrow><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>3</mn></mrow></msub><mo>−</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>4</mn></mrow></msub></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac><mo>+</mo><mo>⋯</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mn>5</mn><mo stretchy="false">(</mo><mi>R</mi><mi>C</mi><msup><mo stretchy="false">)</mo><mn>2</mn></msup><mtext>​</mtext><mo>⋅</mo><mtext>​</mtext><mfrac><mrow><mfrac><mrow><msub><mi>y</mi><mi>i</mi></msub><mo>−</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>1</mn></mrow></msub></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac><mo>−</mo><mfrac><mrow><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>1</mn></mrow></msub><mo>−</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>2</mn></mrow></msub></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac><mo>+</mo><mo>⋯</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mn>6</mn><mi>R</mi><mi>C</mi><mtext>​</mtext><mo>⋅</mo><mtext>​</mtext><mfrac><mrow><msub><mi>y</mi><mi>i</mi></msub><mo>−</mo><msub><mi>y</mi><mrow><mi>i</mi><mo>−</mo><mn>1</mn></mrow></msub></mrow><mrow><mi>Δ</mi><mi>t</mi></mrow></mfrac><mo>+</mo><msub><mi>y</mi><mi>i</mi></msub><mo>=</mo><mo>−</mo><msub><mi>x</mi><mi>i</mi></msub></mrow></mtd></mtr></mtable></mrow></math> Upon further derivation, Eq. (21) can be obtained. <math display="block" id="M21"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>y</mi><mo stretchy="false">[</mo><mi>n</mi><mo stretchy="false">]</mo><mtext>​</mtext><mo>=</mo><mtext>​</mtext></mrow></mtd><mtd columnalign="left"><mrow><mfrac><mrow><mo stretchy="false">(</mo><mn>2</mn><msup><mi>k</mi><mn>3</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>10</mn><msup><mi>k</mi><mn>2</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>6</mn><mi>k</mi><mo stretchy="false">)</mo><mtext>​</mtext><mo>⋅</mo><mtext>​</mtext><mi>y</mi><mo stretchy="false">[</mo><mi>n</mi><mo>−</mo><mn>1</mn><mo stretchy="false">]</mo><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>5</mn><msup><mi>k</mi><mn>2</mn></msup><mtext>​</mtext><mo>⋅</mo><mtext>​</mtext><mi>y</mi><mo stretchy="false">[</mo><mi>n</mi><mo>−</mo><mn>2</mn><mo stretchy="false">]</mo></mrow><mrow><msup><mi>k</mi><mn>3</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>5</mn><msup><mi>k</mi><mn>2</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>6</mn><mi>k</mi><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>2</mn></mrow></mfrac><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mo>⋯</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mfrac><mrow><mo>−</mo><mn>2</mn><msup><mi>k</mi><mn>3</mn></msup><mo>⋅</mo><mi>y</mi><mo stretchy="false">[</mo><mi>n</mi><mo>−</mo><mn>3</mn><mo stretchy="false">]</mo><mo>−</mo><msup><mi>k</mi><mn>3</mn></msup><mo>⋅</mo><mi>y</mi><mo stretchy="false">[</mo><mi>n</mi><mo>−</mo><mn>4</mn><mo stretchy="false">]</mo><mo>−</mo><mi>x</mi><mo stretchy="false">[</mo><mi>n</mi><mo stretchy="false">]</mo></mrow><mrow><msup><mi>k</mi><mn>3</mn></msup><mo>+</mo><mn>5</mn><msup><mi>k</mi><mn>2</mn></msup><mo>+</mo><mn>6</mn><mi>k</mi><mo>+</mo><mn>2</mn></mrow></mfrac></mrow></mtd></mtr></mtable></mrow></math> When k = RC / Δ t , the transfer function of the third-order multiple feedback low-pass filter shaping circuit in the Z domain is shown in Eq. (22). <math display="block" id="M22"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>H</mi><mo stretchy="false">(</mo><mi>z</mi><mo stretchy="false">)</mo><mtext>​</mtext><mo>=</mo><mtext>​</mtext></mrow></mtd><mtd columnalign="left"><mrow><mfrac><mrow><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>1</mn></mrow><mrow><msup><mi>k</mi><mn>3</mn></msup><msup><mrow><mo stretchy="false">(</mo><mn>1</mn><mtext>​</mtext><mo>−</mo><mtext>​</mtext><msup><mi>z</mi><mrow><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>1</mn></mrow></msup><mo stretchy="false">)</mo></mrow><mn>3</mn></msup><mo stretchy="false">(</mo><mn>1</mn><mtext>​</mtext><mo>+</mo><mtext>​</mtext><msup><mi>z</mi><mrow><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>1</mn></mrow></msup><mo stretchy="false">)</mo><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>5</mn><msup><mi>k</mi><mn>2</mn></msup><msup><mrow><mo stretchy="false">(</mo><mn>1</mn><mtext>​</mtext><mo>−</mo><mtext>​</mtext><msup><mi>z</mi><mrow><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>1</mn></mrow></msup><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>6</mn><mi>k</mi><mo stretchy="false">(</mo><mn>1</mn><mtext>​</mtext><mo>−</mo><mtext>​</mtext><msup><mi>z</mi><mrow><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>1</mn></mrow></msup><mo stretchy="false">)</mo><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>2</mn></mrow></mfrac></mrow></mtd></mtr></mtable></mrow></math> The frequency response function is shown in Eq. (23), <math display="block" id="M23"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mi>H</mi><mtext>​</mtext><mo stretchy="false">(</mo><mi>z</mi><mo stretchy="false">)</mo><mtext>​</mtext><mtext>​</mtext><mo>=</mo><mtext>​</mtext><mfrac><mrow><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>1</mn></mrow><mrow><msup><mi>k</mi><mn>3</mn></msup><msup><mrow><mo stretchy="false">(</mo><mn>1</mn><mtext>​</mtext><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mtext>​</mtext><msup><mi>e</mi><mrow><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mi>j</mi><mi>w</mi></mrow></msup><mtext>​</mtext><mo stretchy="false">)</mo></mrow><mn>3</mn></msup><mo stretchy="false">(</mo><mn>1</mn><mtext>​</mtext><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mtext>​</mtext><msup><mi>e</mi><mrow><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mi>j</mi><mi>w</mi></mrow></msup><mtext>​</mtext><mo stretchy="false">)</mo><mtext>​</mtext><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mtext>​</mtext><mn>5</mn><msup><mi>k</mi><mn>2</mn></msup><mtext>​</mtext><msup><mrow><mo stretchy="false">(</mo><mtext>​</mtext><mn>1</mn><mtext>​</mtext><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mtext>​</mtext><msup><mi>e</mi><mrow><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mi>j</mi><mi>w</mi></mrow></msup><mtext>​</mtext><mo stretchy="false">)</mo></mrow><mn>2</mn></msup><mtext>​</mtext><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mtext>​</mtext><mn>6</mn><mi>k</mi><mo stretchy="false">(</mo><mn>1</mn><mtext>​</mtext><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mtext>​</mtext><msup><mi>e</mi><mrow><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mi>j</mi><mi>w</mi></mrow></msup><mtext>​</mtext><mo stretchy="false">)</mo><mtext>​</mtext><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>2</mn></mrow></mfrac></mrow></mtd></mtr></mtable></mrow></math> and its modulus is given by Eq. (24). <math display="block" id="M24"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><mtext>​</mtext><mo>|</mo><mtext>​</mtext><mi>H</mi><mtext>​</mtext><mo stretchy="false">(</mo><mtext>​</mtext><msup><mi>e</mi><mrow><mi>j</mi><mi>w</mi></mrow></msup><mtext>​</mtext><mo stretchy="false">)</mo><mtext>​</mtext><mo>|</mo><mtext>​</mtext></mrow></mtd><mtd columnalign="left"><mrow><mtext>​</mtext><mo>=</mo><mtext>​</mtext><mfrac><mn>1</mn><mrow><msqrt><mrow><mn>2</mn><mtext>​</mtext><mo>⋅</mo><mtext>​</mtext><mo stretchy="false">(</mo><mi>A</mi><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mtext>​</mtext><mi>B</mi><mtext>​</mtext><mo>⋅</mo><mtext>​</mtext><mi>cos</mi><mo stretchy="false">(</mo><mi>w</mi><mo stretchy="false">)</mo><mtext>​</mtext><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mi>C</mi><mtext>​</mtext><mtext>​</mtext><mo>⋅</mo><mtext>​</mtext><mi>cos</mi><mo stretchy="false">(</mo><mn>2</mn><mi>w</mi><mo stretchy="false">)</mo><mtext>​</mtext><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mtext>​</mtext><mi>D</mi><mtext>​</mtext><mtext>​</mtext><mo>⋅</mo><mtext>​</mtext><mi>cos</mi><mo stretchy="false">(</mo><mn>3</mn><mi>w</mi><mo stretchy="false">)</mo><mtext>​</mtext><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mtext>​</mtext><mi>E</mi><mtext>​</mtext><mtext>​</mtext><mo>⋅</mo><mtext>​</mtext><mi>cos</mi><mo stretchy="false">(</mo><mn>4</mn><mi>w</mi><mo stretchy="false">)</mo><mtext>​</mtext><mo stretchy="false">)</mo></mrow></msqrt></mrow></mfrac></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mi>A</mi><mtext>​</mtext><mo>=</mo><mtext>​</mtext><mn>5</mn><msup><mi>k</mi><mn>6</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>25</mn><msup><mi>k</mi><mn>5</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>57</mn><msup><mi>k</mi><mn>4</mn></msup><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>88</mn><msup><mi>k</mi><mn>3</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>46</mn><msup><mi>k</mi><mn>2</mn></msup><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>12</mn><mi>k</mi><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>2</mn></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mi>B</mi><mtext>​</mtext><mo>=</mo><mtext>​</mtext><mn>2</mn><mi>k</mi><mo stretchy="false">(</mo><mn>2</mn><msup><mi>k</mi><mn>5</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>10</mn><msup><mi>k</mi><mn>4</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>41</mn><msup><mi>k</mi><mn>3</mn></msup><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>58</mn><msup><mi>k</mi><mn>2</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>28</mn><mi>k</mi><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>6</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mi>C</mi><mtext>​</mtext><mo>=</mo><mtext>​</mtext><msup><mi>k</mi><mn>2</mn></msup><mo stretchy="false">(</mo><mn>4</mn><msup><mi>k</mi><mn>4</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>20</mn><msup><mi>k</mi><mn>3</mn></msup><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>37</mn><msup><mi>k</mi><mn>2</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>30</mn><mi>k</mi><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>10</mn><mo stretchy="false">)</mo></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mi>D</mi><mtext>​</mtext><mo>=</mo><mtext>​</mtext><mn>4</mn><msup><mi>k</mi><mn>6</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>20</mn><msup><mi>k</mi><mn>5</mn></msup><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>18</mn><msup><mi>k</mi><mn>4</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>4</mn><msup><mi>k</mi><mn>3</mn></msup></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow><mi>E</mi><mtext>​</mtext><mo>=</mo><mtext>​</mtext><msup><mi>k</mi><mn>6</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>5</mn><msup><mi>k</mi><mn>5</mn></msup><mtext>​</mtext><mo>−</mo><mtext>​</mtext><mn>6</mn><msup><mi>k</mi><mn>4</mn></msup><mtext>​</mtext><mo>+</mo><mtext>​</mtext><mn>2</mn><msup><mi>k</mi><mn>3</mn></msup></mrow></mtd></mtr><mtr columnalign="left"><mtd columnalign="left"><mrow/></mtd><mtd columnalign="left"><mrow/></mtd></mtr></mtable></mrow></math> Eq. (24) may be used to obtain the amplitude-frequency response curves. 2.5 Amplitude-frequency response curves of inverting filter shaping circuits 2 The amplitude-frequency response curves of common inverting filters at different shaping parameter values, according to Eq. (6), are shown in Fig. 5(a) Fig. 5(a) . Fig. 5(b) Fig. 5(b) shows the amplitude-frequency response curves of the improved inverting filters at different shaping parameter values according to Eq. (12). The amplitude-frequency response curves of the multiple feedback low-pass filters at different shaping parameter values, according to Eq. (18), are shown in Fig. 5(c) Fig. 5(c) . Similarly, Fig. 5(d) Fig. 5(d) shows the amplitude-frequency response curves of third-order multiple feedback low-pass filters at different shaping parameter values according to Eq. (24). Figure 5(a) Figure 5(a) shows that the common inverting filter shaping function model functions as a low-pass filter. When k is small, the amplitude-frequency response curve declines relatively slowly. As the value of k increases, the amplitude-frequency response curve declines faster. The closer the value of k is to zero, the smaller the cut-off frequency and the lower the noise suppression performance. Figure 5(b) Figure 5(b) shows that the cut-off frequency of the improved inverting filter is lower than that of the common inverting filter. Thus, it performs better as a low-pass filter. Higher values of the shaping parameter k indicate better filtering performance. From Fig. 5(c) Fig. 5(c) , we observe that the multiple feedback low-pass filter shaping function model behaves as a low-pass filter. A comparison with the common inverting filter reveals similar behavior. As shown in Fig. 5(d) Fig. 5(d) , the cut-off frequency of the third-order multiple feedback low-pass filter is lower than that of the multiple feedback low-pass filter, which indicates better low-pass filtering performance. The amplitude of this filter is the lowest among the four filters. According to Fig. 5 Fig. 5 , a large RC value, i.e., a large value of k leads to a small cut-off frequency of the amplitude-frequency response curve and better noise suppression performance by digital filter shaping. This is beneficial for amplitude extraction and energy resolution of the system, but it simultaneously widens the digital Gaussian output, produces smaller peaks, and causes a slower time response, which are not conducive for pulse pile-up identification. This affects pulse counting rate the system. On the contrary, a smaller value of k implies a larger cut-off frequency, which defeats the purpose of filter shaping and affects the effective extraction and analysis of pulse amplitudes. 3 Comparative analysis of digital shaping results 1 The time constant of the front-end circuit of the signal output of the Si-PIN detector system is 10 μs. Digital shaping is performed using different shaping methods and shaping parameters at a sampling rate of 40 MHz. Eqs. (3), (9), (15), and (21) are used to implement the four types of inverting digital filter shaping processes for nuclear pulse signals. 3.1 Shaping results with inverting digital filter 2 The original signal is shaped using the same parameter value ( k = 10). The numerical recursive modes of the common inverting filter shaping, improved inverting filter shaping, multiple feedback low-pass filter shaping, and third-order multiple feedback low-pass filter shaping are used to process nuclear pulse signals, and the results are shown in Fig. 6(a) Fig. 6(a) . Considering the noise suppression performance and pulse amplitude extraction, the digital shaping outputs of the multi-feedback low-pass filter shaping circuit and the third-order multi-feedback low-pass filter shaping circuit at different shaping parameter values are selected for comparison and analysis. For the same signal, the digital shaping results of the above two methods when k = 5, 10, 20, are shown in Fig. 6(b) Fig. 6(b) and Fig. 6(c) Fig. 6(c) , respectively. Figure 6(a) Figure 6(a) shows that at the same shaping parameter values, common inverting digital shaping output exhibits the largest amplitude, but the worst shaping symmetry and noise suppression performance. In comparison, improved inverting digital shaping performs better in terms of noise suppression and waveform symmetry. The third-order multi-feedback low-pass digital shaping output has the smallest pulse amplitude and the best noise suppression performance, whereas the multi-feedback low-pass digital shaping output exhibits the best noise suppression and pulse amplitude comprehensive performance. Figure 6(b) Figure 6(b) indicates that the pulse amplitude, noise suppression, and peak position information of the multiple feedback digital shaping output differ at different shaping parameter values. As k increases, the noise suppression performance improves, the pulse amplitude becomes smaller, and the peak position of the shaping output pulse shifts to the right. The above conclusion is consistent with the theoretical analysis and the results of the amplitude-frequency response curve. As shown in Fig. 6(c) Fig. 6(c) , the digital shaping results of the third-order multi-feedback low-pass filter resemble those in Fig. 6(b) Fig. 6(b) . However, in comparison with the multiple feedback digital shaping output, the third-order multiple feedback digital shaping output has an additional RC low-pass filter circuit at the front end, and therefore, under the same shaping parameter values, the third-order multiple feedback digital shaping output has better noise suppression performance. Additionally, the presence of the RC filter circuit lowers the amplitude of the output pulse amplitude slightly. Therefore, in actual inverse digital shaping, the multiple feedback low-pass filter shaping circuit is capable of appropriately increasing the shaping parameter values to obtain a better SNR, and the third-order multiple feedback filter shaping circuit can appropriately reduce the shaping parameter value. Thus, the premise of ensuring better noise suppression leads to a higher SNR. 3.2 Performance analysis of digital shaping methods 2 To compare and analyze the outputs of inverting digital shaping under the same shaping parameter values and different methods, or the same shaping methods and different shaping parameter values, the amplitude, symmetry, and SNR after inverting digital shaping are quantitatively evaluated. From Fig. 6 Fig. 6 , we observe that the output pulse amplitude of the inverting digital shaping has different degrees of attenuation for different shaping parameters. We quantitatively evaluate the amplitude attenuation index by using Eq. (25) to evaluate the amplitude attenuation. <math display="block" id="M25"><mrow><mi>d</mi><mo>=</mo><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><mfrac><mi>A</mi><mrow><msub><mi>A</mi><mn>0</mn></msub></mrow></mfrac><mo stretchy="false">)</mo><mo>×</mo><mn>100</mn><mi>%</mi></mrow></math> Here, A 0 and A are the pulse amplitudes before and after shaping, respectively, and d is the amplitude attenuation coefficient. A smaller d indicates lower pulse amplitude attenuation. Eq. (26) is used to quantitatively evaluate the symmetry of the digital shaping output. <math display="block" id="M26"><mrow><mi>δ</mi><mo>=</mo><mfrac><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover></mstyle><msup><mrow><mo stretchy="false">(</mo><msub><mi>x</mi><mrow><mn>0</mn><mo>−</mo><mi>i</mi></mrow></msub><mo>−</mo><msub><mi>x</mi><mrow><mn>0</mn><mo>+</mo><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></mrow><mn>2</mn></msup></mrow><mi>N</mi></mfrac></mrow></math> Here, i = 1, … N , is an integer, x 0- i is the i th point on the left side of the signal peak position, x 0+ i is the i th point on the right side of the signal peak position, and N is the range of symmetry, i. e., the number of points. δ is used to indicate the degree of symmetry and a smaller δ implies better symmetry. For example, the δ of a Gaussian pulse is 0. Eq. (27) is used to quantitatively evaluate the SNR of the digital shaping output. A higher SNR indicates better noise suppression performance. <math display="block" id="M27"><mrow><mi>S</mi><mi>N</mi><mi>R</mi><mo>=</mo><mfrac><mrow><msub><mi>V</mi><mrow><mtext>max</mtext></mrow></msub></mrow><mrow><msqrt><mrow><mrow><mo>(</mo><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover></mstyle><msup><mrow><mo stretchy="false">(</mo><msub><mi>V</mi><mi>i</mi></msub><mo>−</mo><msub><mi>V</mi><msup><mi>i</mi><mo>′</mo></msup></msub><mo stretchy="false">)</mo></mrow><mn>2</mn></msup></mrow><mo>)</mo></mrow><mo>/</mo><mi>N</mi></mrow></msqrt></mrow></mfrac></mrow></math> Here, V max is the maximum value of the output pulse amplitude, V i is the pulse amplitude value corresponding to the i th point, <math id="M28"><mrow><msub><mi>V</mi><msup><mi>i</mi><mo>′</mo></msup></msub></mrow></math> is the pulse amplitude value of the i th point obtained using polynomial fitting, and N is the number of calculation points. Table. 1 shows the amplitude, symmetry, and SNR corresponding to Fig. 6 Fig. 6 . In particular, Tables. 1(a), 1(b), and 1(c) correspond to Figs. 6(a), 6(b) Figs. 6(a), 6(b) , and 6(c), respectively. The amplitude attenuation is accurately estimated by using the third pulse in Fig. 5 Fig. 5 for calculating the amplitude attenuation parameters. Table. 1(a) Table. 1(a) shows that the common inverting digital shaping method has the smallest amplitude attenuation index (14.15%), but its SNR (38.69) and symmetry (120.35) are the worst. In comparison, the third-order multiple feedback low-pass digital shaping method has the best SNR (52.30) and symmetry (0.61), but its amplitude attenuation index is poor (37.29%), while the multiple feedback low-pass digital shaping method exhibits the best synthesis performance. 1 1 Shaping methods 1 1 d 1 1 δ 1 1 SNR 4 1 [Performance comparison of different shaping methods.] 1 1 Common inverting 1 1 14.15% 1 1 120.35 1 1 38.69 1 1 Improved inverting 1 1 24.68% 1 1 9.57 1 1 40.30 1 1 Multiple feedback 1 1 26.01% 1 1 3.45 1 1 46.23 1 1 Third-order multiple 1 1 37.29% 1 1 0.61 1 1 52.30 1 1 Shaping parameter 1 1 d 1 1 δ 1 1 SNR 4 1 [Performance comparison of the multiple feedback digital shaping method.] 1 1 k=5 1 1 17.44% 1 1 13.09 1 1 40.61 1 1 k=10 1 1 26.01% 1 1 3.45 1 1 46.23 1 1 k=20 1 1 38.18% 1 1 0.018 1 1 53.57 4 1 [Performance comparison of the third-order multiple feedback digital shaping method.] 1 1 k=5 1 1 25.45% 1 1 2.67 1 1 45.42 1 1 k=10 1 1 37.29% 1 1 0.61 1 1 52.30 1 1 k=20 1 1 48.95% 1 1 0.0062 1 1 63.02 As evident from Table. 1(b) Table. 1(b) , for the multiple feedback low-pass digital shaping output, the amplitude attenuation index increases from 17.44% to 38.18%, the symmetry index parameter decreases from 13.09 to 0.018, and the SNR index is increased from 40.61 to 53.57, as k increases. The value of k = 10 is recommended to obtain the best comprehensive performance. Table. 1(c) Table. 1(c) shows that for the third-order multiple feedback low-pass digital shaping output, the amplitude attenuation index increases from 25.45% to 48.95%, the symmetry index parameter decreases from 2.67 to 0.0062, and the SNR index increases from 45.42 to 63.02 as k increases. After comprehensive consideration, k = 5 is found to be the best choice. 4 Conclusion 1 The digitalization of the common inverting filter shaping circuit, improved inverting filter shaping circuit, multiple feedback low-pass filter shaping circuit, and third-order multiple feedback low-pass filtering shaping circuits were achieved. The transfer functions of the four types of circuits in the Laplacian and Z domains were obtained. An analysis of the filter performance showed that all the four types of filter shaping circuits exhibited low-pass filter performance. The digital shaping model of the third-order multiple feedback low-pass filter is preferred in terms of noise suppression and waveform symmetry performance. In contrast, in terms of comprehensive performance of amplitude extraction, noise suppression, and waveform symmetry, the digital shaping model of the multiple feedback low-pass filter is preferred. References W.T. Evariste , S.U. Hong , Q. Yi et al ., Design and simulation of Gaussian shaping amplifier made only with CMOS FET for FEE of particle detector . Nucl. Sci. Tech . 21 ( 5 ), 312 - 315 ( 2010 ). https://doi.org/10.13538/j.1001-8042/nst.21 https://doi.org/10.13538/j.1001-8042/nst.21 . Design and simulation of Gaussian shaping amplifier made only with CMOS FET for FEE of particle detector B. Gan , T.C. Wei , W. Gao , et al ., Design and performances of a low-noise and radiation-hardened readout ASIC for CdZnTe detectors . Journal of Semiconductors . 37 ( 6 ), 065007 ( 2016 ). https://doi.org/10.1088/1674-4926/37/6/065007 https://doi.org/10.1088/1674-4926/37/6/065007 Design and performances of a low-noise and radiation-hardened readout ASIC for CdZnTe detectors W. Gao , S. Li , Y. Duan et al ., Design and characterization of a low-noise front-end readout ASIC in 0.18 μm CMOS technology for CZT/Si-PIN detectors . IEEE Transact Nucl Sci . 65 ( 5 ), 1203 - 1211 ( 2018 ). https://doi.org/10.1109/TNS.2018.2826070 https://doi.org/10.1109/TNS.2018.2826070 Design and characterization of a low-noise front-end readout ASIC in 0.18 μm CMOS technology for CZT/Si-PIN detectors M.M. Zhang , Z.J. Chen , Y.C. Zhang et al ., Design and test results of a low-noise readout integrated circuit for high-energy particle detectors , Nucl. Sci. Tech . 21 ( 1 ), 44 - 48 ( 2010 ). https://doi.org/10.13538/j.1001-8042/nst.21.44-48 https://doi.org/10.13538/j.1001-8042/nst.21.44-48 Design and test results of a low-noise readout integrated circuit for high-energy particle detectors Y. Rong , G.Q. Zeng , L.Q. Ge et al ., Design of charge sensitive preamplifier based on RC low pass feedback . Nucl. Tech . 42 ( 9 ), 090402 ( 2019 ). https://doi.org/10.11889/j.0253-3219.2019.hjs.42.090402 https://doi.org/10.11889/j.0253-3219.2019.hjs.42.090402 . (in Chinese) Design of charge sensitive preamplifier based on RC low pass feedback G.Q. Zeng , S. Qing , X.P. OuYang et al ., Preliminary development of a new high efficiency and compact anti-Compton spectrometer . Nucl. Tech . 42 ( 2 ), 020401 ( 2019 ). https://doi.org/10.11889/j.0253-3219.2019.hjs.42.020401 https://doi.org/10.11889/j.0253-3219.2019.hjs.42.020401 . (in Chinese) Preliminary development of a new high efficiency and compact anti-Compton spectrometer D.X. Wang , C.J. Lin , L. Yang et al ., Compact 16-channel integrated charge-sensitive preamplifier module for silicon strip detectors . Nucl. Sci. Tech . 31 , 48 ( 2020 ). https://doi.org/10.1007/s41365-020-00755-0 https://doi.org/10.1007/s41365-020-00755-0 (in Chinese) Compact 16-channel integrated charge-sensitive preamplifier module for silicon strip detectors H.Q. Zhang , Z.D. Li , B. Tang , et al ., Optimal parameter choice of CR-RCm digital filter in nuclear pulse processing . Nucl. Sci. Tech . 30 , 108 ( 2019 ). https://doi.org/10.1007/s41365-019-0638-7 https://doi.org/10.1007/s41365-019-0638-7 . Optimal parameter choice of CR-RCm digital filter in nuclear pulse processing H.Q. Zhang , B. Tang , H.X. Wu , et al ., Study of Sallen-Key digital filters in nuclear pulse signal processing . Nucl. Sci. Tech . 30 , 145 ( 2019 ). https://doi.org/10.1007/s41365-019-0679-y https://doi.org/10.1007/s41365-019-0679-y Study of Sallen-Key digital filters in nuclear pulse signal processing Y.Y. Liu , J.L. Zhang , L.F. Liu , et al ., Implementation of real-time digital CR-RCm shaping filter on FPGA for gamma-ray spectroscopy . Nucl Instrum Methods Phys Res Sec. A 906 , 1 - 9 ( 2018 ). https://doi.org/10.1016/j.nima.2018.05.020 https://doi.org/10.1016/j.nima.2018.05.020 Implementation of real-time digital CR-RCm shaping filter on FPGA for gamma-ray spectroscopy J.B. Zhou , W. Zhou , J.R. Lei , et al ., Study of time-domain digital pulse shaping algorithms for nuclear signals . Nucl. Sci. Tech . 23 ( 3 ), 150 - 155 ( 2012 ). https://doi.org/10.13538/j.1001-8042/nst.23.150-155 https://doi.org/10.13538/j.1001-8042/nst.23.150-155 Study of time-domain digital pulse shaping algorithms for nuclear signals J.B. Zhou , W. Zhou , X. Hong , Improvement of digital S-K filter and its application in nuclear signal processing . Nucl. Sci. Tech . 24 ( 6 ), 060401 ( 2013 ). https://doi.org/10.13538/j.1001-8042/nst.2013.06.020 https://doi.org/10.13538/j.1001-8042/nst.2013.06.020 Improvement of digital S-K filter and its application in nuclear signal processing Q. Ge , L.Q. Ge , H.W. Yuan , et al ., A new digital Gaussian pulse shaping algorithm based on bilinear transformation . Nucl. Sci. Tech . 26 ( 1 ), 010402 ( 2015 ). https://doi.org/10.13538/j.1001-8042/nst.26.010402 https://doi.org/10.13538/j.1001-8042/nst.26.010402 . A new digital Gaussian pulse shaping algorithm based on bilinear transformation H.Q. Zhang , W.H. Lu , B. Tang . Digital shaping and performance analysis of nuclear pulse signal based on MATLAB . Nucl. Tech . 41 ( 10 ), 100401 ( 2018 ). https://doi.org/10.11889/j.0253-3219.2018.hjs.41.100401 https://doi.org/10.11889/j.0253-3219.2018.hjs.41.100401 . (in Chinese) Digital shaping and performance analysis of nuclear pulse signal based on MATLAB X. Hong , S.J. Ni , J.B. Zhou , et al ., Simulation study on Gaussian pulse shaping algorithm . Nucl. Tech . 39 ( 11 ), 110403 ( 2016 ). https://doi.org/10.11889/j.0253-3219.2016.hjs.39.110403 https://doi.org/10.11889/j.0253-3219.2016.hjs.39.110403 . (in Chinese) Simulation study on Gaussian pulse shaping algorithm Z.D. Li , H.Q. Zhang , J.Y. Liu , et al ., Implementation and analysis of Gaussian shaping method for digital nuclear pulse signal . Nucl. Tech . 42 ( 6 ), 060403 ( 2019 ). https://doi.org/10.11889/j.0253-3219.2019.hjs.42.060403 https://doi.org/10.11889/j.0253-3219.2019.hjs.42.060403 . (in Chinese) Implementation and analysis of Gaussian shaping method for digital nuclear pulse signal H.Q. Zhang , M.M. Yan , J.Y. Liu , Digitalization of filter shaping circuit for nuclear pulse signal . Nucl. Tech . 42 ( 7 ), 070402 ( 2019 ). https://doi.org/10.11889/j.0253-3219.2019.hjs.42.070402 https://doi.org/10.11889/j.0253-3219.2019.hjs.42.070402 . (in Chinese) Digitalization of filter shaping circuit for nuclear pulse signal Y.Y. Liu , J.L. Zhang , R. Zhou , et al ., Digitalization of CR-RCm filter . Nucl. Tech . 40 ( 6 ), 44 - 48 ( 2017 ). https://doi.org/10.11889/j.0253-3219.2017.hjs.40.060403 https://doi.org/10.11889/j.0253-3219.2017.hjs.40.060403 . (in Chinese) Digitalization of CR-RCm filter A. Regadio , S. Sanchez-Prieto , M. Prieto , et al ., Implementation of a real-time adaptive digital shaping for nuclear spectroscopy . Nucl Instrum Methods Phys. Res. Sec. A 735 , 297 - 303 ( 2014 ). https://doi.org/10.1016/j.nima.2013.09.063 https://doi.org/10.1016/j.nima.2013.09.063 Implementation of a real-time adaptive digital shaping for nuclear spectroscopy Valentin T. Jordanov . Unfolding-synthesis technique for digital pulse processing. Part 1: Unfolding. Nucl Instrum Methods Phys . Res. Sec. A 805 , 63 - 71 ( 2016 ). https://doi.org/10.1016/j.nima.2015.07.040 https://doi.org/10.1016/j.nima.2015.07.040 Unfolding-synthesis technique for digital pulse processing. Part 1: Unfolding. Nucl Instrum Methods Phys S. Saxena , A.I. Hawari . Investigation of FPGA-based real-time Adaptive digital pulse shaping for high-count-rate applications . IEEE Transact Nucl Sci . 64 ( 7 ), 1733 - 1738 ( 2017 ). https://doi.org/10.1109/TNS.2017.2692219 https://doi.org/10.1109/TNS.2017.2692219 Investigation of FPGA-based real-time Adaptive digital pulse shaping for high-count-rate applications J. Zhou , J.B. Zhou , J.R. Lei et al ., Simulation analysis of numerical model of inverting and noninverting amplifier filter circuit for nuclear signal in time domain . China Measurement & Test , 41 ( 9 ), 101 - 105 ( 2015 ). https://doi.org/10.11857/j.issn.1674-5124.2015.09.023 https://doi.org/10.11857/j.issn.1674-5124.2015.09.023 . (in Chinese) Simulation analysis of numerical model of inverting and noninverting amplifier filter circuit for nuclear signal in time domain 10.1007/s41365-020-00799-2 This work was supported by the National Key R & D Project (No. 2017YFF0106503), National Natural Science Foundation of China (Nos. 11665001 and 41864007). 2020-02-29 2020-07-22 2020-07-24 2020-08-28 2020-09-01 2021-04-10 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2020 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-020-00799-2 http://dx.doi.org/10.1007/s41365-020-00799-2 . 2020 Nuclear Science and Techniques 9 31