Search tools Text Classification Chemistry Measure Numbers Full documents Title Abstract Claims All Any Exact Not Add AND condition These CPCs and their children These exact CPCs Add AND condition
Exact Exact Batch Similar Substructure Substructure (SMARTS) Full documents Claims only Add AND condition
Add AND condition
Application Numbers Publication Numbers Either Add AND condition

Trapping wave forming method of originating broadband signal with variable parameters

Abstract
translated from Chinese

本发明涉及一种可变参数的发端宽带信号的陷波成形方法,属于通信信号处理技术领域。它包括:在信号发送端,实时检测出各有用的窄带信号的个数n,及每个有用的窄带信号的带宽和频带分布,根据发送端各有用窄带信号的带宽和频带分布构造与各有用窄带信号相应的n个子陷波器,并将各n个子陷波器串联成一个级联陷波器,将宽带发送信号通过该级联陷波器进行陷波,得到具有与所述各有用的窄带信号频率相应的陷波的宽带信号,继续检测各有用的窄带信号,如果发现参数发生变化,更新陷波器参数。本发明的特点在于:对窄带信号进行实时检测,并根据检测结果实时更新陷波器参数,从而有效抑制宽带发射信号对参数可变的窄带信号产生的干扰。

The invention relates to a variable-parameter notch forming method of a transmitting broadband signal, which belongs to the technical field of communication signal processing. It includes: at the signal sending end, real-time detection of the number n of each useful narrowband signal, and the bandwidth and frequency band distribution of each useful narrowband signal, according to the bandwidth and frequency band distribution of each useful narrowband signal at the sending end N sub-notch filters corresponding to the narrowband signal, and each n sub-notch filters are connected in series to form a cascaded notch filter, and the broadband transmission signal is notched through the cascaded notch filter to obtain useful The notched wideband signal corresponding to the frequency of the narrowband signal continues to detect useful narrowband signals, and if any parameter changes are found, the parameters of the notch filter are updated. The present invention is characterized in that the narrowband signal is detected in real time, and the notch filter parameters are updated in real time according to the detection result, thereby effectively suppressing the interference caused by the wideband transmission signal to the narrowband signal with variable parameters.

Landscapes

Show more

CN101795248B

China

Other languages
Chinese
Inventor
肖立民
关安福
孙引
周世东
冯伟
许希斌
李云洲
张秀军
何飞
Current Assignee
Tsinghua University

Worldwide applications
2010 CN

Application CN201010033671XA events
2012-11-14
Application granted
Expired - Fee Related
Anticipated expiration

Description

A kind of trap manufacturing process of the broadband signal of making a start of variable element
Technical field
The invention belongs to the signal of communication processing technology field, particularly a kind of trap manufacturing process of the broadband signal of making a start of variable element.
Background technology
In present communication system, wide-band spread spectrum communication system and narrow-band communication system coexistence, because band resource is limited, sometimes at transmitting terminal, wideband spread-spectrum signal and the more shared frequency ranges of useful narrow band signal meeting, promptly the two frequency spectrum has taken place overlapping.If be sent out simultaneously at transmitting terminal wideband spread-spectrum signal and useful narrow band signal, also can mix on the frequency spectrum of narrow band signal has the part broadband signal, and narrow band signal has just received the interference of broadband signal like this.In actual scene, narrow-band communication system is time dependent in different operating frequency upper signal channel condition, and the bandwidth of the therefore number of useful narrow band signal, and each useful narrow band signal and these parameters of frequency band distribution are time dependent.Existing method can only be carried out trap and is shaped to the useful narrow band signal of preset parameter, can't effectively suppress the interference of broadband emission signal to the narrow band signal generation of changeable parameters.
At present both at home and abroad to communication system suppress The Study of Interference more be how to suppress the influence of narrow-band interference signal for wideband received signal at receiving terminal.In these research work, considered receiving terminal the narrow-band interference signal that will suppress be the signal of fixed frequency, like power frequency 50Hz interference signal, thus the parameter of the trapper of the receiving terminal of design fix, do not need real-time adjustment.And send in the scene of signal coexistence in transmitting terminal broadband signal and arrowband; Narrow-band communication system is time dependent in different operating frequency upper signal channel condition; So number of useful narrow band signal; Bandwidth and these parameters of frequency band distribution of reaching each useful narrow band signal are time dependent, need the design effective method suppress the interference of broadband emission signal to the narrow band signal generation of changeable parameters.
Summary of the invention
The objective of the invention is when sending signal in transmitting terminal wideband sending signal and arrowband and exist simultaneously, wideband sending signal sends the problem that signal produces interference to the arrowband, has proposed a kind of trap manufacturing process of the broadband signal of making a start of variable element.The inventive method is flexible, convenient, and complexity is lower, is easy to realize.To the number n of narrow band signal, sub-trapper-these parameters of centre frequency ω of 3dB attenuation bandwidth B W and sub-trapper detect in real time; And according to testing result real-time update trap parameter, thereby effectively suppress of the interference of broadband emission signal to the narrow band signal generation of changeable parameters.
A kind of trap manufacturing process of the broadband signal of making a start of variable element is characterized in that, may further comprise the steps:
1) at signal sending end, detect the number n of current each useful narrow band signal, n is a positive integer, and the bandwidth BW of each useful narrow band signal 1, BW 2..., BW nDistribute with frequency band;
2) be connected into a cascade trapper according to the bandwidth of each useful narrow band signal of transmitting terminal and frequency band distributed structure and the corresponding n of this useful narrow band signal sub-trapper, and with the individual sub-trapper of this n;
3) wideband sending signal is carried out trap through this cascade trapper, obtain having broadband signal with said each useful corresponding trap of narrow band signal frequency;
4) detect the number n of each useful narrow band signal in real time, and the bandwidth BW of each useful narrow band signal 1, BW 2..., BW nDistribute with frequency band, change, forward step 2 to if find these parameters and the detected parameter of step 1)), otherwise, repeating step 4).
Said step 2) constructs and the corresponding n of each useful narrow band signal sub-trapper, and each n sub-trapper is connected into a cascade trapper, specifically comprise;
21) according to the bandwidth BW of each useful narrow band signal 1, BW 2..., BW nConfirm each sub-trapper-3dB attenuation bandwidth [BW 1BW 2BW n] T, it be one by BW 1, BW 2..., BW nA vector that constitutes;
22) distribute according to the frequency band of each useful narrow band signal of transmitting terminal and confirm the centre frequency [ω of each sub-trapper 1ω 2ω n] T, the method for confirming is: if the initial frequency of useful narrow band signal i is f I1, the termination frequency is f I2, the centre frequency of i sub-trapper then ω i = f i 1 + f i 2 2 , I=1,2 ... N;
23) confirm the sample frequency f of trapper s, f sOperating frequency according to the transmitting terminal digital information processing system confirms that its value is not higher than the maximum operating frequency of transmitting terminal digital information processing system;
24) calculate the normalization attenuation bandwidth and the centre frequency of each sub-trapper, normalized bandwidth is:
[ BW 1 ′ BW 2 ′ · · · BW n ′ ] T = 2 π * [ BW 1 BW 2 · · · BW n ] T f s , Normalized centre frequency is:
[ ω 1 ′ ω 2 ′ · · · ω n ′ ] T = 2 π * [ ω 1 ω 2 · · · ω n ] T f s ;
25) calculate the frequency parameter and the bandwidth parameter of each sub-trapper, frequency parameter does Sin θ 1 i = Sin ( ω i ′ - π 2 ) , I=1 wherein, 2 ... N, bandwidth parameter does Sin θ 2 i = 1 - Tan ( BW i ′ / 2 ) 1 + Tan ( BW i ′ / 2 ) = Tan ( π 4 - BW i ′ 2 ) , I=1 wherein, 2 ... N;
26) according to step 25) frequency parameter and each sub-trapper of bandwidth parameter structure of the sub-trapper that obtains, obtain the transfer function of n sub-trapper, the transfer function of i sub-trapper is:
H i ( z ) = 1 + sin θ 2 i 2 1 + 2 sin θ 1 i z - 1 + z - 2 1 + sin θ 1 i ( 1 + sin θ 2 i ) z - 1 + sin θ 2 i z - 2
I=1 wherein, 2 ... N, z -1Expression is to sampling clock cycle T of digital signal sequences time-delay of input s, T s = 1 f s ;
27) with step 26) n of a gained trapper series connection obtains a cascade trapper, and the transfer function of cascade trapper is: H (z)=H 1(z) * H 2(z) * ... * H n(z), go out digital cascade trapper according to the transfer function configures of cascade trapper then.
Characteristics of the present invention and effect:
In the inventive method, for the number n of sub-trapper, sub-trapper-3dB attenuation bandwidth BW iAnd the centre frequency ω of sub-trapper iBut flexible configuration (i=1 wherein, 2 all ... N).Therefore, when useful narrow band signal changes, number n that can the sub-trapper of real-time update, sub-trapper-3dB attenuation bandwidth BW iAnd the centre frequency ω of sub-trapper i, just can obtain new cascade trapper structure, thereby effectively suppress of the interference of broadband emission signal the narrow band signal generation of changeable parameters.The inventive method is flexible, convenient, and complexity is lower, is easy to realize.
The present invention can be used as and suppresses the method for broadband signal to the interference of the narrow band signal of changeable parameters in the Shortwave Communication System, also can be used as other communication system when transmitting terminal wideband spread-spectrum signal and the coexistence of useful narrow band signal and suppresses the method for broadband signal to the interference of the useful narrow band signal of changeable parameters.
Description of drawings
Fig. 1 is the flow chart of the generation method of trapper among the present invention;
Fig. 2 is the structural representation of the sub-trapper of the present invention's employing;
Fig. 3 is the structural representation of the cascade trapper of the present invention's employing;
Fig. 4 is the amplitude-frequency response of the cascade trapper in the embodiment of the invention.
Embodiment
Further explain method of the present invention below in conjunction with accompanying drawing and embodiment.
The trap manufacturing process of the broadband signal of making a start of a kind of variable element of the present invention, its flow process is as shown in Figure 1, it is characterized in that, may further comprise the steps:
1) at signal sending end, detect the number n of current each useful narrow band signal, n is a positive integer, and the bandwidth BW of each useful narrow band signal 1, BW 2..., BW n(described detection method is the known conventional technology in this area) distributes with frequency band;
2) be connected into a cascade trapper according to the bandwidth of each useful narrow band signal of transmitting terminal and frequency band distributed structure and the corresponding n of each useful narrow band signal sub-trapper, and with each n sub-trapper;
3) wideband sending signal is carried out trap through this cascade trapper, obtain having broadband signal with said each useful corresponding trap of narrow band signal frequency;
4) detect the number n of each useful narrow band signal in real time, and the bandwidth BW of each useful narrow band signal 1, BW 2..., BW nDistribute with frequency band, change, forward step 2 to if find these parameters and the detected parameter of step 1)), otherwise, repeating step 4).
Said step 2) constructs and the corresponding n of each useful narrow band signal sub-trapper, and each n sub-trapper is connected into a cascade trapper, specifically comprise;
21) according to the bandwidth BW of each useful narrow band signal 1, BW 2..., BW nConfirm each sub-trapper-3dB attenuation bandwidth [BW 1BW 2BW n] T, it be one by BW 1, BW 2..., BW nA vector that constitutes;
22) distribute according to the frequency band of each useful narrow band signal of transmitting terminal and confirm the centre frequency [ω of each sub-trapper 1ω 2ω n] T, the method for confirming is: if the initial frequency of arrowband i is f I1, the termination frequency is f I2, the centre frequency of i sub-trapper then ω i = f i 1 + f i 2 2 , I=1,2 ... N;
23) confirm the sample frequency f of trapper s, f sOperating frequency according to the transmitting terminal digital information processing system confirms that its value is not higher than the maximum operating frequency of transmitting terminal digital information processing system;
24) calculate the normalization attenuation bandwidth and the centre frequency of each sub-trapper, normalized bandwidth does
[ BW 1 ′ BW 2 ′ · · · BW n ′ ] T = 2 π * [ BW 1 BW 2 · · · BW n ] T f s , Normalized centre frequency does
[ ω 1 ′ ω 2 ′ · · · ω n ′ ] T = 2 π * [ ω 1 ω 2 · · · ω n ] T f s ;
25) calculate the frequency parameter and the bandwidth parameter of each sub-trapper, frequency parameter does Sin θ 1 i = Sin ( ω i ′ - π 2 ) , I=1 wherein, 2 ... N, bandwidth parameter does Sin θ 2 i = 1 - Tan ( BW i ′ / 2 ) 1 + Tan ( BW i ′ / 2 ) = Tan ( π 4 - BW i ′ 2 ) , I=1 wherein, 2 ... N;
26) according to step 25) frequency parameter and each sub-trapper of bandwidth parameter structure of the sub-trapper that obtains, obtain the transfer function of n sub-trapper, the transfer function of i sub-trapper is:
H i ( z ) = 1 + sin θ 2 i 2 1 + 2 sin θ 1 i z - 1 + z - 2 1 + sin θ 1 i ( 1 + sin θ 2 i ) z - 1 + sin θ 2 i z - 2
I=1 wherein, 2 ... N, z -1Expression is to sampling clock cycle T of digital signal sequences time-delay of input s, T s = 1 f s ;
27) with step 26) n of a gained trapper series connection obtains a cascade trapper, and the transfer function of cascade trapper is: H (z)=H 1(z) * H 2(z) * ... * H n(z), then according to the digital cascade trapper of the transfer function configures of cascade trapper.
Embodiment 1:
The method of present embodiment may further comprise the steps:
1) at signal sending end, adopt spectrum analyzer to detect make a start number n=3 of useful narrow band signal of current time, the bandwidth of useful narrow band signal 1 is BW 1=2kHz, frequency band is distributed as 0.349-0.351MHz, and the bandwidth of useful narrow band signal 2 is BW 1=3kHz, frequency band is distributed as 0.3985-0.4015MHz, and the bandwidth of useful narrow band signal 3 is BW 1=3.5kHz, frequency band distribution 0.69825-0.70175MHz;
2) be connected into a cascade trapper according to the bandwidth of each useful narrow band signal of transmitting terminal and frequency band distributed structure and the corresponding n of each useful narrow band signal sub-trapper, and with each n sub-trapper;
Said step 2) constructs and the corresponding n of each useful narrow band signal sub-trapper, and each n sub-trapper is connected into a cascade trapper, specifically comprise;
21) according to the bandwidth of each useful narrow band signal confirm each sub-trapper-3dB attenuation bandwidth [BW 1BW 2BW n] T, it be one by BW 1, BW 2..., BW nA vector that constitutes, [BW 1BW 2BW 3] T=[200030003500] THz;
22) distribute according to the frequency band of each useful narrow band signal of transmitting terminal and confirm the centre frequency [ω of each sub-trapper 1ω 2ω n] T, the method for confirming is: if the initial frequency of arrowband i is f I1, the termination frequency is f I2, the centre frequency of i sub-trapper then ω i = f i 1 + f i 2 2 , I=1,2 ... N can obtain through calculating:
[ ω 1 ω 2 ω 3 ] T = [ 0.349 + 0.351 2 0.3985 + 0.4015 2 0.69825 + 0.70175 2 ] T * 10 6 ;
= [ 0.350.40.7 ] T * 10 6 Hz
23) confirm the sample frequency f of trapper s, f sOperating frequency according to the transmitting terminal digital information processing system confirms that its value is not higher than the maximum operating frequency of transmitting terminal digital information processing system, confirms f s=2*10 6Hz;
24) calculate the normalization attenuation bandwidth and the centre frequency of each sub-trapper, normalized bandwidth does
[ BW 1 ′ BW 2 ′ BW 3 ′ ] T = 2 π * [ BW 1 BW 2 BW 3 ] T f s
= 2 π * [ 200030003500 ] T 2 * 10 6 ,
= [ 0.00630.00940.0110 ] T
Normalized centre frequency does
[ ω 1 ′ ω 2 ′ ω 3 ′ ] T = 2 π * [ ω 1 ω 2 ω 3 ] T f s
= 2 π * [ 0.350.40.7 ] T * 10 6 2 * 10 6 ;
= [ 1.00961.25662.1991 ] T
25) calculate the frequency parameter and the bandwidth parameter of each sub-trapper, frequency parameter does Sin θ 1 i = Sin ( ω i ′ - π 2 ) , I=1 wherein, 2 ... N, bandwidth parameter does Sin θ 2 i = 1 - Tan ( BW i ′ / 2 ) 1 + Tan ( BW i ′ / 2 ) = Tan ( π 4 - BW i ′ 2 ) , I=1 wherein, 2 ... N,
The frequency parameter and the bandwidth parameter of sub-trapper 1 are respectively Sin θ 11 = Sin ( ω 1 ′ - π 2 ) = Sin ( 1.0096 - π 2 ) = - 0.4540 ,
sin θ 21 = tan ( π 4 - BW 1 ′ 2 ) = tan ( π 4 - 0.0063 2 ) = 0.9937 ;
The frequency parameter and the bandwidth parameter of sub-trapper 2 are respectively
sin θ 12 = sin ( ω 2 ′ - π 2 ) = sin ( 1.2566 - π 2 ) = - 0.3090 ,
sin θ 22 = tan ( π 4 - BW 2 ′ 2 ) = tan ( π 4 - 0.0094 2 ) = 0.9906 ;
The frequency parameter and the bandwidth parameter of sub-trapper 3 are respectively
sin θ 13 = sin ( ω 3 ′ - π 2 ) = sin ( 2.1991 - π 2 ) = 0.5878 ,
sin θ 23 = tan ( π 4 - BW 3 ′ 2 ) = tan ( π 4 - 0.0110 2 ) = 0.9891 ;
26) according to step 25) frequency parameter and each sub-trapper of bandwidth parameter structure of the sub-trapper that obtains, obtain the transfer function of 3 sub-trappers, the transfer function of each sub-trapper is respectively:
H 1 ( z ) = 1 + sin θ 21 2 1 + 2 sin θ 11 z - 1 + z - 2 1 + sin θ 11 ( 1 + sin θ 21 ) z - 1 + sin θ 21 z - 2
= 1 + 0.9937 2 1 - 2 * 0.4540 * z - 1 + z - 2 1 - 0.4540 * ( 1 + 0.9937 ) z - 1 + 0.9937 z - 2
= 0.9969 - 0.9051 z - 1 + 0.9969 z - 2 1 - 0.9051 z - 1 + 0.9937 z - 2
H 2 ( z ) = 1 + sin θ 22 2 1 + 2 sin θ 12 z - 1 + z - 2 1 + sin θ 12 ( 1 + sin θ 22 ) z - 1 + sin θ 22 z - 2
= 0.9953 - 0.6151 z - 1 + 0.9953 z - 2 1 - 0.6151 z - 1 + 0.9906 z - 2
H 3 ( z ) = 1 + sin θ 23 2 1 + 2 sin θ 13 z - 1 + z - 2 1 + sin θ 13 ( 1 + sin θ 23 ) z - 1 + sin θ 23 z - 2 ;
= 0.9945 + 1.1691 z - 1 + 0.9945 z - 2 1 + 1.1691 z - 1 + 0.9891 z - 2
27) with step 26) n of a gained trapper series connection obtains a cascade trapper, and the transfer function of cascade trapper is:
H ( z ) = H 1 ( z ) * H 2 ( z ) * H 3 ( z )
= 0.9969 - 0.9051 z - 1 + 0.9969 z - 2 1 - 0.9051 z - 1 + 0.9937 z - 2 * 0.9953 - 0 . 6151 z - 1 + 0.9953 z - 2 1 - 0.6151 z - 1 + 0.9906 z - 2 * 0.9945 + 1.1691 z - 1 + 0.9945 z - 2 1 + 1.691 z - 1 + 0.9891 z - 2
Then according to the digital cascade trapper of the transfer function configures of cascade trapper.
3) wideband sending signal is carried out trap through this cascade trapper, obtain having broadband signal with said each useful corresponding trap of narrow band signal frequency;
4) the employing spectrum analyzer detects the number n of each useful narrow band signal in real time, and the bandwidth BW of each useful narrow band signal 1, BW 2..., BW nDistribute with frequency band, change, forward step 2 to if find these parameters and the detected parameter of step 1)), otherwise, repeating step 4).
It is the digital notch method that has adopted in the digital filtering theory that the present invention adopts cascade trapper structure.Each sub-trapper can be restrained a narrow-band interference signal.The structure of sub-trapper is as shown in Figure 2, wherein uses the trend of the wire list registration word burst of band arrow, sin θ 1The frequency parameter that is called as sub-trapper, sin θ 2The bandwidth parameter that is called as sub-trapper,
Figure G201010033671XD000610
Be digital adder,
Figure G201010033671XD000611
Under to wear a minus sign be digital subtractor,
Figure G201010033671XD000612
The expression digital multiplier, the multiple of multiplier be marked on multiplier directly over, The expression d type flip flop, it can be to digital signal sequences sampling clock cycle of time-delay of input T s = 1 f s , F wherein sBe sample frequency.The input and output signal burst of digital adder, subtracter, multiplier and d type flip flop is represented by the direction of arrow.According to the relevant knowledge of the structure and the Digital Signal Processing of sub-trapper shown in Figure 2, the transfer function that can obtain sub-trapper is:
H ( z ) = 1 + sin θ 2 2 1 + 2 sin θ 1 z - 1 + z - 2 1 + sin θ 1 ( 1 + sin θ 2 ) z - 1 + sin θ 2 z - 2
Z wherein -1Expression is to sampling clock cycle T of digital signal sequences time-delay of input sConversely, according to the transfer function H (z) of sub-trapper, also can design the structure of sub-trapper according to structure chart shown in Figure 2.The cascade trapper is together in series a plurality of sub-trappers exactly, if the transfer function of i sub-trapper is H i(z), the transfer function of cascade trapper is exactly H (z)=H 1(z) * H 2(z) * ... * H n(z).The structure of cascade trapper is as shown in Figure 3.Make z=e among the H (z) Jw', w ' is the frequency independent variable of digital signal, H (e so Jw') be the frequency response of trapper, A (w ')=| H (e Jw') | be the amplitude-frequency response of trapper.
The amplitude-frequency response A (w ') of the cascade trapper of present embodiment=| H (e Jw') | as shown in Figure 4, this trapper sample frequency is f s=2MHz, trap number n=3, trap centre frequency and notch are respectively:
ω 1=0.35MHz,BW 1=2kHz;ω 2=0.4MHz,BW 1=3kHz;ω 3=0.7MHz,BW 3=3.5kHz。Can see that from Fig. 4 this cascade trapper is at digital signal frequency w 1' 0.35 π, w 2' 0.4 π, w 3Three frequencies of '=0.7 π have produced three traps, and these several frequency corresponding simulating signal frequencies do w 1 = w 1 ′ f s 2 π = 0.35 MHz , w 2 = w 2 ′ f s 2 π = 0.4 MHz , w 3 = w 3 ′ f s 2 π = 0.7 MHz , Just be the centre frequency of useful narrow band signal.Further analysis can be found the also lucky bandwidth corresponding to corresponding useful narrow band signal of each trap-3dB attenuation bandwidth.
Can find out the trap manufacturing process of the broadband signal of making a start of a kind of variable element proposed by the invention by the foregoing description, flexibly, convenient, and complexity is lower, is easy to realize.When spectrum analyzer detect the number n of narrow band signal, sub-trapper-when these parameters of centre frequency ω of 3dB attenuation bandwidth BW and sub-trapper change; Can construct the cascade trapper again according to new measurement result, thereby effectively suppress of the interference of broadband emission signal the narrow band signal generation of changeable parameters.
The foregoing description just is used to specify the trap shaping Algorithm of the configurable broadband signal of making a start of the present invention; Concrete data wherein just arbitrarily are provided with for convenient explanation; Can not be in order to limit protection scope of the present invention; Promptly as long as implement by the described step of this claim, wherein any variation of data all should belong to protection category of the present invention.

Claims (1)
Hide Dependent

1. the trap manufacturing process of the broadband signal of making a start of a variable element is characterized in that, may further comprise the steps:
1) at signal sending end, detect the number n of current each useful narrow band signal, n is a positive integer, and the bandwidth BW of each useful narrow band signal 1, BW 2..., BW nDistribute with frequency band;
2) be connected into a cascade trapper according to the bandwidth of each useful narrow band signal of transmitting terminal and frequency band distributed structure and the corresponding n of each useful narrow band signal sub-trapper, and with this n sub-trapper;
3) wideband sending signal is carried out trap through this cascade trapper, obtain having broadband signal with said each useful corresponding trap of narrow band signal frequency;
4) detect the number n of each useful narrow band signal in real time, and the bandwidth BW of each useful narrow band signal 1, BW 2..., BW nDistribute with frequency band, change, forward step 2 to if find these parameters and the detected parameter of step 1)), otherwise, repeating step 4);
Said step 2) constructs and the corresponding n of each useful narrow band signal sub-trapper, and this n sub-trapper is connected into a cascade trapper, specifically comprise;
21) according to the bandwidth BW of each useful narrow band signal 1, BW 2..., BW nConfirm each sub-trapper-3dB attenuation bandwidth [BW 1BW 2BW n] T, it be one by BW 1, BW 2..., BW nA vector that constitutes;
22) distribute according to the frequency band of each useful narrow band signal of transmitting terminal and confirm the centre frequency [ω of each sub-trapper 1ω 2ω n] T, the method for confirming is: if the initial frequency of useful narrow band signal i is f I1, the termination frequency is f I2, the centre frequency of i sub-trapper then ω i = f i 1 + f i 2 2 , I=1,2 ... N;
23) confirm the sample frequency f of trapper s, f sOperating frequency according to the transmitting terminal digital information processing system confirms that its value is not higher than the maximum operating frequency of transmitting terminal digital information processing system;
24) calculate the normalization attenuation bandwidth and the centre frequency of each sub-trapper, normalized bandwidth is:
[ B W 1 ′ B W 2 ′ · · · B W n ′ ] T = 2 π * [ B W 1 B W 2 · · · B W n ] T f s , Normalized centre frequency is:
[ ω 1 ′ ω 2 ′ · · · ω n ′ ] T = 2 π * [ ω 1 ω 2 · · · ω n ] T f s ;
25) calculate the frequency parameter and the bandwidth parameter of each sub-trapper, frequency parameter does
Figure FSB00000867213200014
I=1 wherein, 2 ... N, bandwidth parameter does Sin θ 2 i = 1 - Tan ( B W i ′ / 2 ) 1 + Tan ( B W i ′ / 2 ) = Tan ( π 4 - B W i ′ 2 ) , I=1 wherein, 2 ... N;
26) according to step 25) frequency parameter and each sub-trapper of bandwidth parameter structure of the sub-trapper that obtains, obtain the transfer function of n sub-trapper, the transfer function of i sub-trapper is:
H i ( z ) = 1 + sin θ 2 i 2 1 + 2 sin θ 1 i z - 1 + z - 2 1 + sin θ 1 i ( 1 + sin θ 2 i ) z - 1 + sin θ 2 i z - 2
I=1 wherein, 2 ... N, z -1Expression is to sampling clock cycle T of digital signal sequences time-delay of input s,
27) with step 26) n of a gained trapper series connection obtains a cascade trapper, and the transfer function of cascade trapper is: H (z)=H 1(z) * H 2(z) * ... * H n(z), go out digital cascade trapper according to the transfer function configures of cascade trapper then.