Anti-aliasing response cancellation signal processing method and anti-aliasing response cancellation filter
By calculating and generating anti-aliasing response cancellation filter bank, using an optimization algorithm to design the filter tap coefficient, the problem of difficulty in performing system response cancellation and broadband multi-sampling rate signal processing in the prior art is solved, and the response cancellation and anti-aliasing filter with low operational amount is realized, which improves the efficiency and effect of signal processing.
Patent Information
- Application Number
- CN202411654268.7
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-11-19
- Publication Date
- 2025-03-04
AI Technical Summary
The prior art is difficult to simultaneously perform system response cancellation and broadband multi-sampling rate signal processing, resulting in problems such as beam width widening and signal detection probability decrease.
By measuring the system response, the anti-aliasing response cancellation filter group is calculated and the tap coefficient of the filter is designed using an optimization algorithm to achieve the fusion of response cancellation and anti-aliasing filtering.
Without increasing the calculation amount, system response cancellation and anti-aliasing filtering are implemented, which improves the efficiency and effect of signal processing and ensures the processing effect of subsequent steps.
Smart Images

Figure CN119561520A_ABST
Abstract
Description
Technical Field
[0001] The invention relates to an anti-aliasing response cancellation signal processing method and an anti-aliasing response cancellation filter, belonging to the field of broadband signal processing, and specifically to an anti-aliasing response cancellation signal processing method and an anti-aliasing response cancellation filter for broadband multi-sampling rate signal processing. Background Art
[0002] In actual situations, affected by the amplitude-frequency response of components such as transducers and preamplifiers, the system's transmit / receive signals will produce serious distortion. If the system response offset is not performed, the processing effects of subsequent beamforming, signal detection and other steps will be affected, resulting in problems such as beam width widening and decreased signal detection probability. There are two ways to achieve system response offset. One is to use the frequency domain method to first segment the input signal and perform FFT, then perform response compensation through the frequency domain weighting method, and finally perform IFFT and splice the segmented signals into the original signal. The problem with this method is that it increases the amount of calculation in the system and causes the problem of discontinuous output signals. The other is to design a filter through the time domain method and use the amplitude-frequency response of the filter to offset the system response. The advantage of this method is that it does not require signal segmentation processing, avoiding the problem of discontinuous output signals.
[0003] For broadband signals, the frequency division band processing method is generally used in processing. According to the Nyquist sampling constant, the system sampling rate can recover all the information of the signal when it exceeds twice the maximum frequency of the signal. For low-frequency signals, the multi-sampling rate signal processing technology can be used to reduce the system's computing burden and improve the system's computing efficiency and real-time performance. Multi-sampling rate signal processing technology requires multiple anti-aliasing filtering and downsampling processing of the signal to obtain frequency division band distortion-free signals with different sampling rates. Among them, anti-aliasing filtering is generally implemented using a FIR low-pass filter group. Summary of the invention
[0004] In view of this, an embodiment of the present invention provides an anti-aliasing response cancellation signal processing method and an anti-aliasing response cancellation filter to solve the problem that the prior art cannot simultaneously perform system response cancellation and broadband multi-sampling rate signal processing.
[0005] A first aspect of an embodiment of the present invention provides an anti-aliasing response cancellation signal processing method, comprising:
[0006] Measure the system response and calculate the cancellation response required to generate the filter bank;
[0007] According to the cancellation response, an anti-aliasing response cancellation filter group is obtained by an optimization method;
[0008] Get input signal;
[0009] The input signal is processed through the anti-aliasing response cancellation filter group to obtain a multi-sampling rate output signal.
[0010] In one embodiment, the cancellation response includes a response cancellation frequency point and a response amplitude corresponding to the frequency point.
[0011] In one embodiment, the measuring system response and calculating the cancellation response required by the filter group include: determining the response cancellation frequency point according to the signal processing frequency band, obtaining the amplitude response of the target system at the response cancellation frequency point using a standard measuring instrument, and calculating and generating the cancellation response.
[0012] In one embodiment, determining the response cancellation frequency point includes:
[0013] According to the signal processing frequency band requirements, according to f n =f l +n*100, n=0, 1, 2, ..., N-1 determines the response offset frequency point;
[0014] Among them, the f n is the nth response cancellation frequency point (in Hz), where n is the number of the response cancellation frequency point, and f u is the upper limit of the signal processing frequency band, the f l is the lower limit of the signal processing frequency band, In response to the number of offset frequency points, round is a rounding symbol.
[0015] In one embodiment, the calculating generates a cancellation response, comprising:
[0016] According to the system expected response (generally full-pass response), determining the offset response;
[0017] Among them, the H d (f n ) is at the frequency point f n The cancellation response at a (f n ) is at the frequency point f n The system expected response at t (f n ) is at the frequency point f n The system response at f n is the nth response cancellation frequency point, where n is the number of the response cancellation frequency point.
[0018] In one embodiment, after the calculation of the cancellation response required by the filter bank, the method further comprises:
[0019] The cancellation response is linearly interpolated at a resolution of 1 Hz to obtain the interpolated response cancellation frequency point and the corresponding response amplitude.
[0020] In one embodiment, the anti-aliasing response cancellation filter group is obtained by an optimization method according to the cancellation response, including: dividing the signal processing frequency band into octaves, determining the number of filters in the filter group and their working sampling frequencies, and using an optimization algorithm to design the tap coefficients of the filters.
[0021] In one embodiment, the anti-aliasing response cancellation filter bank is represented by h m =[h m (1), h m (2), ..., h m (km)], m=1, 2,..., M;
[0022] Among them, the h m is the mth filter in the filter bank, and the h m (k) is the kth tap coefficient of the mth filter in the filter bank, where k m is the tap length of the mth filter, and M is the number of filters in the anti-aliasing response cancellation filter bank.
[0023] In one embodiment, the operating sampling frequency of the mth filter in the anti-aliasing response cancellation filter is
[0024] Among them, F s is the sampling frequency of the original acquired data, and Q is the downsampling multiple, corresponding to the octave.
[0025] In one embodiment, designing the tap coefficients of the filter by an optimization method includes:
[0026] according to
[0027]
[0028] sth m d m (f j )<∈,f j ∈F sm
[0029] Optimize and solve to get the mth filter h m ;
[0030] Among them, H dm (f n ) is the mth filter at frequency point f n The expected response at h m is the tap coefficient of the mth filter in the filter bank, ∈ is the maximum response amplitude of the specified stop band, F pm is the passband of the mth filter, F sm is the stopband of the mth filter, and the transition band between the passband and the stopband is generally set to 6.25% of the passband width;
[0031] The first filter in the anti-aliasing response cancellation filter bank has an expected response within the passband according to The calculated cancellation response, for the 2nd to Mth filters, has an expected response within the passband as an all-pass response.
[0032] In one embodiment, the step of processing the input signal through the anti-aliasing response cancellation filter bank to obtain a multi-sampling rate output signal includes:
[0033] Perform step-by-step filtering and downsampling according to the octave division of the signal processing frequency band;
[0034] Let x im is the mth level signal input (the original input signal obtained is x i1 ), x om is the m-th level multi-sampling rate signal output, then in one embodiment, the specific steps of the step-by-step filtering and downsampling are: using the m-th filter of the anti-aliasing response offset filter group to filter the input signal x of this level im Filter and set the filtered output signal to x om As the m-th level multi-sampling rate signal output, x om Downsample by downsampling factor Q to get the input x for the m+1th level processing i(m+1) .
[0035] A second aspect of an embodiment of the present invention provides an anti-aliasing response cancellation filter, comprising:
[0036] A cancellation response calculation module, used for measuring the system response and calculating the cancellation response required for generating the filter bank;
[0037] A filter bank generation module is used to obtain an anti-aliasing response cancellation filter bank through an optimization method according to the cancellation response;
[0038] A data acquisition module, used for acquiring input signals;
[0039] The anti-aliasing multi-sampling rate processing module is used to process the input signal through the anti-aliasing response cancellation filter group to obtain a multi-sampling rate output signal.
[0040] The present invention also provides a computer-readable storage medium on which a computer program is stored. The program is executed by a processor to implement the steps of the anti-aliasing response cancellation signal processing method described in the present invention.
[0041] The beneficial effects of the present invention include:
[0042] (1) Low computational complexity: When performing broadband signal processing, if response cancellation and anti-aliasing filtering are to be achieved simultaneously, two steps are generally required, namely, response cancellation processing is performed first, and then anti-aliasing filtering processing is performed. The present invention can integrate the two processing steps of response cancellation and anti-aliasing through an optimization algorithm to form an anti-aliasing response cancellation filter, thereby achieving one-step processing, and can achieve anti-aliasing and response cancellation simultaneously without increasing the computational complexity;
[0043] (2) High flexibility: When generating a filter bank, parameters such as tap length, cutoff frequency, transition band width, and stopband attenuation strength can be flexibly set within a certain range as required;
[0044] (3) Good processing effect: The present invention solves the problem that the system response causes distortion of the system's transmit / receive signals to a certain extent, and can ensure the processing effect of subsequent steps such as beam forming and signal detection. BRIEF DESCRIPTION OF THE DRAWINGS
[0045] Figure 1 A schematic diagram of the implementation flow of the signal processing method provided in an embodiment of the present invention.
[0046] Figure 2 This is the system response diagram used in the simulation.
[0047] Figure 3 Graph showing the effect of filter bank response cancellation for anti-aliasing response cancellation.
[0048] Figure 4 Error plot of filter bank response cancellation for antialiasing response cancellation.
[0049] Figure 5 It is the response diagram of each filter in the anti-aliasing response cancellation filter bank (the vertical axis is the amplitude value).
[0050] Figure 6 Response diagram of each filter in the anti-aliasing response cancellation filter bank (the vertical axis is decibel value).
[0051] Figure 7 This is a comparison chart between the spectrum after adding the system response effect and the standard input signal spectrum.
[0052] Figure 8 The spectrum of the first-stage output signal of the anti-aliasing response cancellation filter bank is compared with the spectrum of the standard input signal.
[0053] Fig. 9 The spectrum of the second-stage output signal of the anti-aliasing response cancellation filter bank is compared with the spectrum of the standard input signal.
[0054] Fig.10 The spectrum of the third-stage output signal of the anti-aliasing response cancellation filter bank is compared with the spectrum of the standard input signal.
[0055] Fig.11 An exemplary diagram of a signal processing device provided in an embodiment of the present invention. DETAILED DESCRIPTION
[0056] In order to illustrate the technical solution of the present invention, a specific embodiment is provided below for illustration.
[0057] Example 1
[0058] See also Figure 1 As shown, a response cancellation signal processing method for anti-aliasing includes the following steps:
[0059] Step 1, measure the system response and calculate the cancellation response required to generate the filter bank, including:
[0060] Determine the response cancellation frequency point according to the signal processing frequency band, use standard measuring instruments to obtain the amplitude response of the target system at the response cancellation frequency point, and calculate and generate the cancellation response; then perform linear interpolation processing on the cancellation response at a resolution of 1 Hz to obtain the interpolated response cancellation frequency point and the corresponding response amplitude;
[0061] The cancellation response includes a response cancellation frequency point and a response amplitude corresponding to the frequency point; and determining the response cancellation frequency point includes:
[0062] According to the signal processing frequency band requirements, according to f n =f l +n*100, n=0, 1, 2, ..., N-1 determines the response offset frequency point;
[0063] Among them, the f n is the nth response cancellation frequency point (in Hz), where n is the number of the response cancellation frequency point, and f u is the upper limit of the signal processing frequency band, the f l is the lower limit of the signal processing frequency band, The number of response offset frequency points, round is the rounding symbol;
[0064] The calculation to generate the offset response comprises:
[0065] According to the system's expected response requirements, determining the offset response;
[0066] Among them, H d (f n ) is at the frequency point f n The cancellation response ata (f n ) is at the frequency point f n The expected response of the system at t (f n ) is at the frequency point f n The system response at f n is the nth response cancellation frequency point, and n is the number of the response cancellation frequency point.
[0067] Step 2, according to the cancellation response, obtaining an anti-aliasing response cancellation filter group by an optimization method, including:
[0068] The method of obtaining an anti-aliasing response cancellation filter group by an optimization method based on the cancellation response includes: dividing the signal processing frequency band into octaves, determining the number of filters in the filter group and their working sampling frequencies, and using an optimization algorithm to design the tap coefficients of the filters;
[0069] The anti-aliasing response cancellation filter bank is denoted as h m =[h m (1), h m (2), ..., h m (k m )], m=1, 2, ..., M;
[0070] Among them, h m is the mth filter in the filter bank, h m (k) is the kth tap coefficient of the mth filter in the filter bank, k m is the tap length of the mth filter, and M is the number of filters in the anti-aliasing response cancellation filter bank.
[0071] The operating sampling frequency of the mth filter in the anti-aliasing response cancellation filter is
[0072] Among them, F s is the sampling frequency of the original acquired data, and Q is the downsampling multiple, corresponding to the octave.
[0073] Step 3: Get the input signal.
[0074] Step 4, processing the input signal through the anti-aliasing response cancellation filter group to obtain a multi-sampling rate output signal, including: performing step-by-step filtering and downsampling according to the octave division of the signal processing frequency band;
[0075] Let x im The original input signal obtained is x i1 , x omThe m-th level multi-sampling rate signal output, the specific steps of the step-by-step filtering and downsampling are: using the anti-aliasing response to cancel the m-th filter of the filter group to the input signal x of this level im Filter and set the filtered output signal to x om As the m-th level multi-sampling rate signal output, x om Downsample by downsampling factor Q to get the input x for the m+1th level processing i(m+1) .
[0076] Example 2
[0077] An anti-aliasing response cancellation filter comprising:
[0078] A cancellation response calculation module is used to measure the system response and calculate the cancellation response required by the filter bank; the cancellation response includes the response amplitude corresponding to the cancellation frequency point and the spectrum point;
[0079] The measuring system response and calculating the cancellation response required by the filter bank include: determining the response cancellation frequency point according to the signal processing frequency band, obtaining the amplitude response of the target system at the response cancellation frequency point by using a standard measuring instrument, and calculating and generating the cancellation response;
[0080] The step of determining the response offset frequency point comprises:
[0081] According to the signal processing frequency band requirements, according to f n =f l +n*100, n=0, 1, 2, ..., N-1 determines the response offset frequency point;
[0082] Among them, f n is the nth response cancellation frequency point (in Hz), n is the number of the response cancellation frequency point, and f u is the upper limit of the signal processing frequency band, f l is the lower bound of the signal processing frequency band, Round is the rounding symbol;
[0083] The calculation generates a cancellation response, including:
[0084] According to the system expected response (generally full-pass response), determining the offset response;
[0085] Among them, H d (f n ) is at the frequency point f n The cancellation response at a (f n ) is at the frequency point f n The expected response of the system at t (fn ) is at the frequency point f n The system response at f n is the nth response cancellation frequency point, and n is the number of the response cancellation frequency point.
[0086] The filter group generation module is used to obtain an anti-aliasing response cancellation filter group through an optimization method according to the cancellation response; the anti-aliasing response cancellation filter group includes multiple anti-aliasing response cancellation filters; the anti-aliasing response cancellation filter is completed by designing the tap coefficients of the filter using an optimization algorithm.
[0087] The anti-aliasing response cancellation filter bank is denoted as h m =[h m (1), h m (2), ..., h m (km)], m=1, 2,..., M;
[0088] Among them, h m is the mth filter in the filter bank, h m (k) is the kth tap coefficient of the mth filter in the filter bank, k m is the tap length of the mth filter, and M is the number of filters in the anti-aliasing response cancellation filter bank.
[0089] The operating sampling frequency of the mth filter in the anti-aliasing response cancellation filter is
[0090] Among them, F s is the sampling frequency of the original acquired data, and Q is the downsampling multiple, corresponding to the octave.
[0091] The data acquisition module is used to obtain input signals.
[0092] The anti-aliasing multi-sampling rate processing module is used to process the input signal through the anti-aliasing response cancellation filter group to obtain a multi-sampling rate output signal.
[0093] Example 3
[0094] According to the signal processing method and filter provided in Example 1 and Example 2, simulation verification is performed, and the simulation parameters are set as follows: the original sampling rate is 50kHz, the system operating frequency band is 0Hz-4kHz, the collected data signal duration is 1 second, and the linear frequency modulation signals in four frequency bands are included: 0.2kHz-0.7kHz, 1.2kHz-1.7kHz, 3.2kHz-3.7kHz, and 5.2kHz-5.7kHz. The system response is as follows Figure 2 As shown, the expected response of the system is an all-pass filter response, the tap length of each filter in the filter bank is 377, and the filter stopband attenuation is -60dB.
[0095] In the simulation, the system response is set as Figure 2 As shown. Set the system expected response to the full-pass response, calculate and interpolate to obtain the cancellation response within the working frequency band as follows Figure 3 Shown by red line.
[0096] According to the octave rule (corresponding to a downsampling factor of 2), the entire operating frequency band can be divided into three sub-bands: 0Hz-4kHz, 0kHz-2kHz, and 0kHz-1kHz. It can be determined that there are three filters in the anti-aliasing response cancellation filter group, and their operating sampling rates are 50kHz, 25kHz, and 12.5kHz, respectively. The filter expected response, stopband constraints, transition band constraints, filter tap length, and operating sampling rate are input into the optimization algorithm to obtain the tap coefficients of each filter in the anti-aliasing response cancellation filter group. It can be calculated that the response of the first filter is as follows: Figure 3 As shown by the blue line, the error between the expected response and the operating frequency band is as follows: Figure 4 As shown; the responses of the 1st, 2nd and 3rd filters are as follows Figure 5 and Figure 6 shown.
[0097] The spectrum of the standard input signal generated by the simulation and the result after adding the system response are shown in the figure below. Figure 7 The signal after adding the system response is processed by the anti-aliasing response cancellation filter group, and the spectrum of the first, second, and third level output signals of the anti-aliasing response cancellation filter group is compared with the spectrum of the standard input signal as shown in the figure. Figure 8 , Fig. 9 , Fig.10 shown.
[0098] from Figure 4 It can be seen from the simulation that the maximum error of response cancellation is about 0.6%, which can meet the requirements of engineering applications. Figure 5 and Figure 6 It can be seen that in the simulation, the first filter in the anti-aliasing response cancellation filter group can complete the anti-aliasing filtering and response cancellation functions, the second and third filters can complete the anti-aliasing filtering function, and the stopband response of the three filters is lower than -60dB. Figure 7 It can be seen from the figure that the spectrum is distorted after the system response is added, and its amplitude changes with frequency. Figure 8 , Fig. 9 , Fig.10 It can be seen from the figure that after the signal with the system response effect is processed by the anti-aliasing response cancellation filter group, its spectrum is consistent with the standard input signal spectrum, which proves that the anti-aliasing response cancellation filter group has the function of response cancellation. Figure 8 , Fig. 9 , Fig.10 It can also be seen that the frequency spectrum ranges of the 1st, 2nd and 3rd level output signals of the anti-aliasing response cancellation filter group are 0kHz-4kHz, 0kHz-2kHz and 0Hz-1kHz respectively, which proves that the anti-aliasing response cancellation filter group has anti-aliasing filtering function.
[0099] In summary, the anti-aliasing response cancellation signal processing method and the anti-aliasing response cancellation filter can meet the requirements of response cancellation and multi-sampling rate signal processing. The simulation verifies the effectiveness of the anti-aliasing response cancellation signal processing method and the anti-aliasing response cancellation filter.
[0100] The signal processing method measures the system response and calculates the cancellation response required by the filter group; obtains the anti-aliasing response cancellation filter group through an optimization method based on the cancellation response; obtains the input signal; processes the input signal through the anti-aliasing response cancellation filter group to obtain a multi-sampling rate output signal, thereby realizing system response cancellation and anti-aliasing filtering without increasing the amount of system calculation, solving the problem that the prior art cannot simultaneously perform system response cancellation and broadband multi-sampling rate signal processing.
Claims
1. A response cancellation signal processing method for anti-aliasing, characterized in that: The following steps are involved: Step 1, measure the system response and calculate the cancellation response required to generate the filter bank, including: Determine the response cancellation frequency point according to the signal processing frequency band, obtain the amplitude response of the target system at the response cancellation frequency point using a standard measuring instrument, and calculate and generate the cancellation response; Step 2, according to the cancellation response, obtaining an anti-aliasing response cancellation filter group by an optimization method, including: Divide the signal processing frequency band into octaves, determine the number of filters in the filter bank and their working sampling frequency, and use optimization algorithms to design the filter tap coefficients; Step 3, obtaining an input signal; Step 4: Process the input signal through the anti-aliasing response cancellation filter group to obtain a multi-sampling rate output signal.
2. The signal processing method according to claim 1, characterized in that: The cancellation response includes a response cancellation frequency point and a response amplitude corresponding to the frequency point; The determining of the response offset frequency point comprises: According to the signal processing frequency band requirements, according to f n =f l +n*100, n=0, 1, 2, ..., N-1 determines the response offset frequency point; Among them, f n is the nth response cancellation frequency point, n is the number of the response cancellation frequency point, f u is the upper limit of the signal processing frequency band, f l is the lower bound of the signal processing frequency band, Round is the rounding symbol; After the calculation to generate the cancellation response required by the filter group, the method further includes performing linear interpolation processing on the cancellation response at a resolution of 1 Hz to obtain the interpolated response cancellation frequency point and the corresponding response amplitude.
3. The anti-aliasing response cancellation signal processing method according to claim 2, characterized in that: In step 1, the calculation to generate the offset response includes: According to the requirement that the system expected response is a full-pass response, determining the offset response; Among them, H d (f n ) is at the frequency point f n The cancellation response at a (f n ) is at the frequency point f n The expected response of the system at t (f n ) is at the frequency point f n The system response at f n is the nth response cancellation frequency point, and n is the number of the response cancellation frequency point.
4. The anti-aliasing response cancellation signal processing method according to claim 3, characterized in that: In step 2: The anti-aliasing response cancellation filter bank is denoted as h m =[h m (1), h m (2), ..., h m (k m )], m=1, 2, ..., M; Among them, h m is the mth filter in the filter bank, h m (k) is the kth tap coefficient of the mth filter in the filter bank, k m is the tap length of the mth filter, and M is the number of filters in the anti-aliasing response cancellation filter bank.
5. The anti-aliasing response cancellation signal processing method according to claim 4, characterized in that: In step 2: The operating sampling frequency of the mth filter in the anti-aliasing response cancellation filter is Among them, F s is the sampling frequency of the original acquired data, and Q is the downsampling multiple corresponding to the octave.
6. The anti-aliasing response cancellation signal processing method according to claim 5, characterized in that: In step 2, the designing of the filter tap coefficients by an optimization method includes: according to Optimize and solve to get the mth filter h m ; Among them, H dm (f n ) is the mth filter at frequency point f n The expected response at h m is the tap coefficient of the mth filter in the filter bank, ∈ is the maximum response amplitude of the specified stop band, F pm is the passband of the mth filter, F sm is the stopband of the mth filter, and the transition band between the passband and the stopband is generally set to 6.25% of the passband width; The anti-aliasing response cancels the first filter in the filter bank, and its expected response in the passband is based on The calculated cancellation response, for the 2nd to Mth filters, has an expected response within the passband as an all-pass response.
7. The anti-aliasing response cancellation signal processing method according to any one of claims 1 to 6, characterized in that: Step 4 also includes: Perform step-by-step filtering and downsampling according to the octave division of the signal processing frequency band; Let x im The original input signal obtained is x i1 , x om The m-th level multi-sampling rate signal output, the specific steps of the step-by-step filtering and downsampling are: using the anti-aliasing response to cancel the m-th filter of the filter group to the input signal x of this level im Filter and set the filtered output signal to x om As the m-th level multi-sampling rate signal output, and x om Downsample by downsampling factor Q to get the input x for the m+1th level processing i(m+1) .
8. A computer-readable storage medium having a computer program stored thereon, characterized in that: The program is executed by a processor to implement the steps of the anti-aliasing response cancellation signal processing method as described in any one of claims 1 to 7.
9. An anti-aliasing response cancellation filter, characterized in that: It includes a cancellation response calculation module, a filter bank generation module, a data acquisition module and an anti-aliasing multi-sampling rate processing module, wherein: The cancellation response calculation module is used to measure the system response and calculate the cancellation response required by the filter group; the cancellation response includes the response cancellation frequency point and the response amplitude corresponding to the frequency point; the measurement of the system response and the calculation of the cancellation response required by the filter group include: determining the response cancellation frequency point according to the signal processing frequency band, obtaining the amplitude response of the target system at the response cancellation frequency point by using a standard measuring instrument, and calculating the cancellation response; The filter group generation module is used to obtain an anti-aliasing response cancellation filter group through an optimization method according to the cancellation response; the anti-aliasing response cancellation filter group includes multiple anti-aliasing response cancellation filters; the anti-aliasing response cancellation filter is completed by designing the tap coefficients of the filter using an optimization algorithm; the anti-aliasing response cancellation filter group is represented by h m =[ h m(1),h m (2), ..., h m (k m )], m=1, 2, ..., M; Among them, h m is the mth filter in the filter bank, h m (k) is the kth tap coefficient of the mth filter in the filter bank, k m is the tap length of the mth filter, M is the number of filters in the anti-aliasing response cancellation filter bank; The data acquisition module is used to obtain input signals; The anti-aliasing multi-sampling rate processing module is used to process the input signal through the anti-aliasing response cancellation filter group to obtain a multi-sampling rate output signal.
10. The anti-aliasing response cancellation filter according to claim 9, characterized in that: The step of determining the response offset frequency point comprises: According to the signal processing frequency band requirements, according to f n =f l +n*100, n=0, 1, 2, ..., N-1 determines the response offset frequency point; Among them, f n is the nth response cancellation frequency point, n is the number of the response cancellation frequency point, f u is the upper limit of the signal processing frequency band, f l is the lower bound of the signal processing frequency band, Round is the rounding symbol; The calculation generates a cancellation response, including: According to the requirement that the system expected response is a full-pass response, determining the offset response; Among them, H d (f n ) is at the frequency point f n The cancellation response at a (f n ) is at the frequency point f n The expected response of the system at t (f n ) is at the frequency point f n The system response at f n is the nth response cancellation frequency point, and n is the number of the response cancellation frequency point.