CN101909038A - Method and device for estimating blind signal bandwidth - Google Patents
Method and device for estimating blind signal bandwidth Download PDFInfo
- Publication number
- CN101909038A CN101909038A CN 201010274934 CN201010274934A CN101909038A CN 101909038 A CN101909038 A CN 101909038A CN 201010274934 CN201010274934 CN 201010274934 CN 201010274934 A CN201010274934 A CN 201010274934A CN 101909038 A CN101909038 A CN 101909038A
- Authority
- CN
- China
- Prior art keywords
- dft
- bandwidth
- long
- initial
- value
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
The invention discloses a method and a device for estimating a blind signal bandwidth. The method comprises the following steps of: transforming a received signal according to a preset discrete Fourier transformation (DFT) point number; acquiring a mid-value bandwidth DFT point number corresponding to a bandwidth mid-value in a signal bandwidth range required by a system and calculating to obtain an initial sliding window length according to the mid-value bandwidth DFT point number; smoothing a DFT transformation result by using the initial sliding window length to obtain an estimated value of an initial bandwidth DFT point number; comparing the estimated value of the initial bandwidth DFT point number with the mid-value bandwidth DFT point number and selecting a window length factor according to a comparison result; acquiring a variable sliding window length processed by the window length factor; processing the DFT transformation result according to the variable sliding window length; and acquiring an estimated value of a signal bandwidth. In the method for estimating the blind signal bandwidth disclosed by the embodiment of the invention, the variable window length is adjusted by increasing the size of an initial estimated bandwidth of an input signal and the input signal is further smoothed, so that an estimated result is more accurate and accuracy is improved.
Description
Technical field
The present invention relates to the cognitive radio technology field, relate in particular to a kind of blind signal bandwidth method of estimation and device.
Background technology
In cognitive radio system, the usually frequency band range of given system handles, the just bandwidth of system receiver when design system.Generally, have a plurality of signals and enter simultaneously in the system receiver, the modulation system difference of these signals, the amount of bandwidth of signal are also different.System only knows usually and received signal, and do not know the information such as amount of bandwidth, modulation type, carrier frequency of the signal that receives.In order better to carry out subsequent treatment such as parameter Estimation, Modulation Identification, demodulation to the received signal, estimate these signals bandwidth separately respectively after the signal that these need be received simultaneously separates usually.When prior art is carried out the signal bandwidth estimation, usually the signal after employing will separate carries out discrete Fourier transform (DFT), and utilize predefined smoothing windows length that signal is carried out smoothly at frequency domain, the method that a certain thresholding intercepting obtains the estimated bandwidth of signal is set then.
But, when utilizing default smoothing windows progress row level and smooth, if processed signal bandwidth is less, the bandwidth of handling the back signal so can be by excessive broadening, if and signal bandwidth is bigger, signal after the processing can leave more burr, estimates that therefore the bandwidth that obtains exists than mistake, and accuracy is low.
Summary of the invention
In view of this, the invention provides a kind of blind signal bandwidth method of estimation and device, utilize prior art to estimate that there is than mistake the problem that accuracy is low in the bandwidth that obtains with solution.Its concrete scheme is as follows:
A kind of blind signal bandwidth method of estimation comprises:
Received signal counted according to default discrete Fourier transform DFT carry out conversion;
The intermediate value bandwidth DFT that obtains the bandwidth intermediate value correspondence in the signal bandwidth scope of system requirements counts, and counts according to described intermediate value bandwidth DFT that to calculate initial smoothing windows long;
Utilize the described DFT transformation results of the long smoothing processing of described initial smoothing windows, obtain the initial bandwidth DFT estimated value of counting;
Count estimated value and described intermediate value bandwidth DFT of described initial bandwidth DFT counted and compare, select the long factor of window according to comparative result, the span of the long factor of described window is for the inverse of the estimated value of counting greater than described initial bandwidth DFT, less than 1;
It is long to obtain the count variable smoothing windows of estimated value correspondence of initial bandwidth DFT after the long factor of described window is handled;
According to the described DFT transformation results of described variable smoothing windows long process;
Obtain with described variable smoothing windows long process after the corresponding signal bandwidth estimated value of result.
Preferably, describedly count according to described intermediate value bandwidth DFT that to calculate initial smoothing windows long:
Determine that the sixth that described intermediate value bandwidth DFT counts is that initial smoothing windows is long.
Preferably, the described described DFT transformation results of the long smoothing processing of described initial smoothing windows of utilizing obtains the initial bandwidth DFT estimated value process of counting and comprises:
Result after the described smoothing processing is carried out normalization;
Choose the data after the thresholding intercepting normalization, obtain the initial bandwidth DFT estimated value of counting.
Preferably, count estimated value and described intermediate value bandwidth DFT of described initial bandwidth DFT counted and compares, select the process of the long factor of window to comprise according to comparative result:
The described initial bandwidth DFT estimated value and 1 of counting is deducted numerical value reciprocal multiplication after the described threshold value;
Whether more described multiplied result counts greater than described intermediate value bandwidth DFT, if, then choose value greater than median in the long factor span of described window, if not, then select value less than median in the long factor span of described window.
Preferably, the process of the described acquisition signal bandwidth estimated value corresponding with the result after the described variable smoothing windows long process comprises:
The DFT that obtains after the described variable smoothing windows long process counts;
Count according to described DFT and to calculate its corresponding signal bandwidth estimated value.
A kind of blind signal bandwidth estimation unit comprises:
The DFT converter unit is used for received signal counted according to default discrete Fourier transform DFT and carries out conversion;
The initial long acquiring unit of smoothing windows, the intermediate value bandwidth DFT that is used to obtain the interior bandwidth intermediate value correspondence of the signal bandwidth scope of system requirements counts, and counts according to described intermediate value bandwidth DFT that to calculate initial smoothing windows long;
The initial bandwidth DFT estimated value acquiring unit of counting is used to utilize the described DFT transformation results of the long smoothing processing of described initial smoothing windows, obtains the initial bandwidth DFT estimated value of counting;
The long factor selected cell of window, be used for count estimated value and described intermediate value bandwidth DFT of described initial bandwidth DFT counted and compare, select the long factor of window according to comparative result, the span of the long factor of described window is for the inverse of the estimated value of counting greater than described initial bandwidth DFT, less than 1;
The long acquiring unit of variable smoothing windows, it is long to be used to obtain the count variable smoothing windows of estimated value correspondence of initial bandwidth DFT after the long factor of described window is handled;
Variable smoothing windows long process unit is used for according to the described DFT transformation results of described variable smoothing windows long process;
Signal bandwidth estimated value acquiring unit, be used to obtain with described variable smoothing windows long process after the corresponding signal bandwidth estimated value of result.
Preferably, the described initial bandwidth DFT estimated value acquiring unit of counting comprises:
The normalization unit is used for the result after the described smoothing processing is carried out normalization;
The thresholding interception unit is used to choose the data after the thresholding intercepting normalization, obtains the initial bandwidth DFT estimated value of counting.
Preferably, the long factor selected cell of described window comprises:
The unit that multiplies each other is used for the described initial bandwidth DFT estimated value and 1 of counting is deducted numerical value reciprocal multiplication after the described threshold value;
The long factor determining unit of window, be used for more described multiplied result and whether count, if then choose value greater than median in the long factor span of described window greater than described intermediate value bandwidth DFT, if not, then select to grow the value of median in the factor span less than described window.
Preferably, described signal bandwidth estimated value acquiring unit comprises:
The DFT acquiring unit of counting, the DFT that is used to obtain after the described variable smoothing windows long process counts;
Signal bandwidth estimated value computing unit is used for counting according to described DFT and calculates its corresponding signal bandwidth estimated value.
From above-mentioned technical scheme as can be seen, the disclosed a kind of blind signal bandwidth method of estimation of the embodiment of the invention, long by increasing the variable window of regulating according to the initial estimation amount of bandwidth of input signal, input signal is carried out further smoothing processing, make estimated result more accurate, improved accuracy.
Description of drawings
In order to be illustrated more clearly in the embodiment of the invention or technical scheme of the prior art, to do to introduce simply to the accompanying drawing of required use in embodiment or the description of the Prior Art below, apparently, accompanying drawing in describing below only is some embodiments of the present invention, for those of ordinary skills, under the prerequisite of not paying creative work, can also obtain other accompanying drawing according to these accompanying drawings.
Fig. 1 is the flow chart of the disclosed blind signal bandwidth method of estimation of the embodiment of the invention;
Fig. 2 is the flow chart of the disclosed another blind signal bandwidth method of estimation of the embodiment of the invention;
Fig. 3 is the structural representation of the disclosed blind signal bandwidth estimation unit of the embodiment of the invention.
Embodiment
Below in conjunction with the accompanying drawing in the embodiment of the invention, the technical scheme in the embodiment of the invention is clearly and completely described, obviously, described embodiment only is the present invention's part embodiment, rather than whole embodiment.Based on the embodiment among the present invention, those of ordinary skills belong to the scope of protection of the invention not making the every other embodiment that is obtained under the creative work prerequisite.
Increased in a kind of blind signal bandwidth method of estimation disclosed by the invention and utilized the variable window length of regulating input signal to be carried out the process of further smoothing processing according to the initial estimation amount of bandwidth of input signal, make estimated result more accurate, improved accuracy.Its embodiment is as described below:
The flow process of the disclosed a kind of blind signal bandwidth method of estimation of present embodiment comprises as shown in Figure 1:
Step S11, received signal counted according to default discrete Fourier transform DFT carry out conversion;
Counting of DFT conversion can be chosen according to the principle of the following stated.Usually the signal bandwidth scope of cognitive radio system demand is [B
Min, B
Max], the system estimation accuracy requirement is given as d
n, sample frequency is f
s, sample frequency f
sGiven by the processing bandwidth that system is designed, satisfy nyquist sampling theorem, discrete Fourier transform (DFT) is counted and is N, and the discrete Fourier transform (DFT) frequency resolution is
For making the discrete Fourier transform (DFT) frequency resolution satisfy system requirements, need make the discrete Fourier transform (DFT) frequency resolution smaller or equal to system estimation accuracy requirement d
i≤ d
n, then N can be by sample frequency fs and discrete Fourier transform (DFT) frequency resolution d
iDetermine
Work as f
sSelect and system estimation accuracy requirement d
nGive regularly, can by
Promptly
Calculate minimum N value
When N greater than N
MinThe time, can select the size of N arbitrarily.
Step S12, the intermediate value bandwidth DFT that obtains the bandwidth intermediate value correspondence in the signal bandwidth scope of system requirements count, and count according to described intermediate value bandwidth DFT that to calculate initial smoothing windows long;
Carry out in the prior art smoothing windows long choose the time, observe little bandwidth signal usually and use short window long level and smooth, big bandwidth signal uses long window length principle of smoothing.Under the condition of known signal bandwidth, find that by a large amount of emulation the length of smoothing windows is 1/4 to 1/8 o'clock of signal bandwidth DFT points N, can obtain smooth effect preferably with this window length to signal smoothing on frequency domain.Because working as the signal of pre-treatment is Blind Source, signal bandwidth the unknown is so select the intermediate value B of system bandwidth scope for use
Mid, obtain its corresponding intermediate value bandwidth DFT and count, count according to this then and calculate the long M of initial smoothing windows
1
Step S13, utilize the described DFT transformation results of the long smoothing processing of described initial smoothing windows, obtain the initial bandwidth DFT estimated value of counting;
According to the method for calculating estimated bandwidth in the prior art,, obtain the initial bandwidth DFT estimated value of counting in this step according to the result behind the long smoothing processing DFT of initial smoothing windows.
Step S14, count estimated value and described intermediate value bandwidth DFT of described initial bandwidth DFT counted and compare, select the long factor of window according to comparative result, the span of the long factor of described window is for the inverse of the estimated value of counting greater than described initial bandwidth DFT, less than 1;
In this step, count estimated value and intermediate value bandwidth DFT of initial bandwidth DFT counted and compare, and then can judge, the bandwidth of current estimated signals is greater than bandwidth intermediate value in the signal bandwidth scope of system requirements or less than this intermediate value, choose the long factor of window according to the result, so that initial bandwidth DFT is counted according to being handled by the actual conditions of estimated signal.
Set the long M of variable window among the present invention
2By the long factor of a certain window and N
B1Multiplying each other obtains, and Zui Xiao window length is 1 in theory, the long factor 1/N of its corresponding window
B1Maximal window is long to be N
B1, the long factor of corresponding window is 1, the theoretical scope of the long factor of window is (1/N
B1, 1).In the above-mentioned steps, smooth effect preferably can be obtained, so can be to choose the long factor of window in (1/8,1/4) in scope also usually in 1/4 to 1/8 o'clock of counting for signal bandwidth DFT of the length of the smoothing windows of choosing.
Step S15, to obtain the count variable smoothing windows of estimated value correspondence of initial bandwidth DFT after the long factor of described window is handled long;
Can embody the bandwidth range characteristics of current processed signal through the initial bandwidth DFT point value after the long factor processing of window.
Step S16, the described DFT transformation results of the described variable smoothing windows long process of foundation;
Utilize variable smoothing windows long process DFT transformation results, its smooth effect is better.
The corresponding signal bandwidth estimated value of result after step S17, acquisition and the described variable smoothing windows long process.
More accurate according to the signal bandwidth estimated value that this step obtains.
Increased in the blind signal bandwidth method of estimation disclosed by the invention and utilized the variable window length of regulating input signal to be carried out the process of further smoothing processing according to the initial estimation amount of bandwidth of input signal, make estimated result more accurate, improved accuracy.
Another blind signal bandwidth method of estimation flow process disclosed by the invention comprises as shown in Figure 2:
Step S21, received signal counted according to default discrete Fourier transform DFT carry out conversion;
Step S22, the intermediate value bandwidth DFT that obtains the bandwidth intermediate value correspondence in the signal bandwidth scope of system requirements count, and count according to described intermediate value bandwidth DFT that to calculate initial smoothing windows long;
Calculate B
MidShared points N in N point DFT
Mid, operational formula is
It is the number of point shared among the result of signal bandwidth intermediate value after the system handles bandwidth is carried out the DFT conversion.Because the length of smoothing windows is 1/4 to 1/8 o'clock of signal bandwidth DFT points N, can obtain smooth effect preferably with this window length to signal smoothing on frequency domain, so choose 1/6 in the present embodiment.
Step S23, the described described DFT transformation results of the long smoothing processing of described initial smoothing windows of utilizing;
Step S24, the result after the described smoothing processing is carried out normalization;
Step S25, choose the data after the thresholding intercepting normalization, obtain the initial bandwidth DFT estimated value of counting;
Signal is after normalization, and the theoretical scope of its intercepting thresholding value is between 0 to 1, and threshold T hr and intercepting are worth relation in direct ratio.If counting, the shared DFT of signal true bandwidth is N
Br, by emulation, bandwidth intercepting value N
B1, thresholding Thr, signal true bandwidth points N
BrBut the following N of triadic relation's approximate expression
Br* (1-Thr) ≈ N
B1If it is higher that thresholding is chosen, intercepting back bandwidth value is less, and error is bigger; If thresholding hangs down the intercepting value near actual value, but affected by noise bigger.Thr can set according to actual conditions, and general span is between 0.4 to 0.7.
Step S26, the described initial bandwidth DFT estimated value and 1 of counting is deducted numerical value reciprocal multiplication after the described threshold value;
Initial bandwidth DFT in this step count estimated value be one by the DFT that obtains after the long level and smooth intercepting of the fixed window intercepting value N that counts
B1, the intermediate value bandwidth DFT that this intercepting value can not be directly corresponding with the bandwidth intermediate value counts and compares, and the value that compares should be able to be reflected in the signal bandwidth of fixing initial smoothing windows after level and smooth, by N
B1* 1/ (1-Thr) obtains.Bandwidth intercepting value N in step S25
B1With thresholding Thr, signal true bandwidth points N
BrThe following N of triadic relation's approximate expression
Br* (1-Thr) ≈ N
B1, obtain a shared DFT points N of signal bandwidth after front window is long level and smooth for making the shared DFT of the initialize signal estimated value of counting
BrApproximation, so described initial bandwidth estimated value and 1 is deducted numerical value reciprocal multiplication after the described threshold value.Its corresponding formula is N
B1* 1/ (1-Thr) is by expression formula N
Br* (1-Thr) ≈ N
B1Derivation obtains N
Br≈ N
B1* 1/ (1-Thr).
Whether step S27, more described multiplied result count greater than described intermediate value bandwidth DFT, if, execution in step S28a then, if not, execution in step S28b then;
Step S28a, choose value greater than median in the long factor span of described window;
The value of median in the factor span is grown in step S28b, selection less than described window;
Step S29, to obtain the count variable smoothing windows of estimated value correspondence of initial bandwidth DFT after the long factor of described window is handled long;
Respectively according to different situations, the result of value that utilization is chosen and step S26 multiplies each other, it is long to obtain variable smoothing windows, why variable, be because its signal bandwidth according to actual estimated has been carried out corresponding adjustment, window is long unfixing, to adapt to the needs of unlike signal bandwidth, guarantees result's accuracy.
Step S210, the described DFT transformation results of the described variable smoothing windows long process of foundation;
Utilize the long M of variable smoothing windows
2It is level and smooth that the DFT value of signal is carried out secondary, and level and smooth result is carried out normalization.Because variable smoothing windows length is according to the signal bandwidth of reality and difference, so the result after the processing also can be according to the difference of signal and difference.
DFT after step S211, the described variable smoothing windows long process of acquisition counts;
Choose Thr and the level and smooth result of normalization is intercepted, obtain the DFT points N of the bandwidth estimation value correspondence of signal as thresholding
B2
Step S212, count according to described DFT and to calculate its corresponding signal bandwidth estimated value.
According to the relation of the resolution of the counting of the shared DFT of bandwidth value=bandwidth * DFT, the DFT that obtains counted is converted to the bandwidth estimation value, is final signal bandwidth estimated value with the product of DFT resolution.
Blind signal bandwidth method of estimation disclosed by the invention, describe the concrete long factor of window of utilizing in detail and handle the initial bandwidth DFT estimated value process of counting, its concrete processing procedure is according to the situation of different processed signals and difference, its result also more can embody the bandwidth characteristics of different processed signals, initial bandwidth DFT after utilize the handling result after the variable smoothing windows progress row secondary smoothing processing that estimated value obtains that counts, more accurate, estimation precision is higher, and can adapt to multiple signal bandwidth scope, enlarged its range of application.
The embodiment of the invention does not limit the method that the DFT that utilizes the bandwidth correspondence in the intermediate treatment process counts and calculates, relatively waits, also can be according to the relation of the resolution of the counting of the shared DFT of bandwidth value=bandwidth * DFT, being converted into bandwidth value handles, just because the numerical value after the discrete Fourier transform (DFT) is discrete values, needing also in the subsequent treatment that bandwidth value is converted to its corresponding DFT counts, so compare with the foregoing description, operation is comparatively complicated.
Application implementation process disclosed by the invention comprises: at first according to the signal bandwidth scope [B of system requirements
Min, B
Max], the intermediate value of getting bandwidth range is designated as B
Mid, calculate B
MidShared points N among the result after the signal that will receive carries out N point DFT
Mid, operational formula is
Setting is with N
Mid1/6 long as initial smoothing windows, then the signal of handling is carried out N point DFT conversion, it is level and smooth that the data X after the conversion (n) does the long slip for M of window, X
1(k) be data after level and smooth, length is N-M+1, level and smooth back data X
1(k) divided by its maximum max (X
1(k)), obtain data X after the normalization
2(k), setting threshold value is 0.5, with X
2(k) value greater than 0.5 number as N
B1With with N
B1Corresponding bandwidth range is divided into two sections, greater than B
MidPart and less than B
MidPart; Calculate B
MidThe DFT points N
MidWith twice N
B1And N
MidRelatively, greater than N
MidThe long M of Shi Erci smoothing windows
2By N
B1With the factor 1/6 multiplication generation mutually, less than N
MidThe long M of Shi Erci smoothing windows
2By N
B1With the factor 1/8 multiplication generation mutually.It is long to set variable smoothing windows then, utilizes variable smoothing windows progress row smoothing processing once more, finally obtains comparatively accurate bandwidth estimation value.
The invention discloses a kind of blind signal bandwidth estimation unit, its structure as shown in Figure 3, comprise: count estimated value acquiring unit 33, window of DFT converter unit 31, the initial long acquiring unit 32 of smoothing windows, initial bandwidth DFT grown factor selected cell 34, the long acquiring unit 35 of variable smoothing windows, variable smoothing windows long process unit 36 and signal bandwidth estimated value acquiring unit 37, wherein:
DFT converter unit 31 is used for received signal counted according to default discrete Fourier transform DFT and carries out conversion; The intermediate value bandwidth DFT that the initial long acquiring unit 32 of smoothing windows is used to obtain the interior bandwidth intermediate value correspondence of the signal bandwidth scope of system requirements counts, and counts according to described intermediate value bandwidth DFT that to calculate initial smoothing windows long; The just initial bandwidth DFT estimated value acquiring unit 33 of counting is used to utilize the described DFT transformation results of the long smoothing processing of described initial smoothing windows, obtains the initial bandwidth DFT estimated value of counting; The long factor selected cell 34 of window, be used for count estimated value and described intermediate value bandwidth DFT of described initial bandwidth DFT counted and compare, select the long factor of window according to comparative result, the span of the long factor of described window is for the inverse of the estimated value of counting greater than described initial bandwidth DFT, less than 1; The long acquiring unit 35 of variable smoothing windows, the variable smoothing windows that is used to obtain the initial bandwidth estimated value correspondence after the long factor of described window is handled is long; Variable smoothing windows long process unit 36 is used for according to the described DFT transformation results of described variable smoothing windows long process; Signal bandwidth estimated value acquiring unit 37 be used to obtain with described variable smoothing windows long process after the corresponding signal bandwidth estimated value of result.
The initial bandwidth DFT estimated value acquiring unit 33 of counting comprises: normalization unit 331 is used for the result after the described smoothing processing is carried out normalization; Thresholding interception unit 332 is used to choose the data after the thresholding intercepting normalization, obtains the initial bandwidth DFT estimated value of counting.
The long factor selected cell 34 of window comprises: the unit 341 that multiplies each other is used for the described initial bandwidth DFT estimated value and 1 of counting is deducted numerical value reciprocal multiplication after the described threshold value; Whether the long factor determining unit 342 of window is used for more described multiplied result and counts greater than described intermediate value bandwidth DFT, if, then choose value, if not, then select value less than median in the long factor span of described window greater than median in the long factor span of described window.
The process of each unit execution concrete operations is as described below in this device:
The signal that DFT converter unit 31 receives system is counted according to default discrete Fourier transform DFT and is carried out conversion, the signal that obtains dispersing, the initial long acquiring unit 32 of smoothing windows obtains intermediate value bandwidth DFT according to the bandwidth intermediate value in the signal bandwidth scope of system requirements and counts, and count according to described intermediate value bandwidth DFT that to calculate initial smoothing windows long, usually choose that intermediate value bandwidth DFT counts 1/8th to 1/4th, the initial bandwidth DFT estimated value acquiring unit 33 of counting utilizes discrete signal after the described DFT conversion of the long smoothing processing of described initial smoothing windows, obtain the initial bandwidth DFT estimated value of counting, the long factor selected cell 34 of window is counted count estimated value and intermediate value bandwidth DFT of initial bandwidth DFT to compare, select the long factor of window according to comparative result, the span of the long factor of described window is the inverse of the estimated value of counting greater than described initial bandwidth DFT, less than 1, and utilize the long factor of the window that obtains to handle initial bandwidth DFT and count, the variable smoothing windows that the long acquiring unit 35 of variable smoothing windows obtains its correspondence according to the initial bandwidth estimated value after handling through the long factor of described window is long; Variable smoothing windows long process unit 36 is according to the described DFT transformation results of variable smoothing windows long process, and is last, signal bandwidth estimated value acquiring unit 37 obtain with described variable smoothing windows long process after the corresponding signal bandwidth estimated value of result.
Each embodiment adopts the mode of going forward one by one to describe in this specification, and what each embodiment stressed all is and the difference of other embodiment that identical similar part is mutually referring to getting final product between each embodiment.For the disclosed device of embodiment, because it is corresponding with the embodiment disclosed method, so description is fairly simple, relevant part partly illustrates referring to method and gets final product.
The professional can also further recognize, the unit and the algorithm steps of each example of describing in conjunction with embodiment disclosed herein, can realize with electronic hardware, computer software or the combination of the two, for the interchangeability of hardware and software clearly is described, the composition and the step of each example described prevailingly according to function in the above description.These functions still are that software mode is carried out with hardware actually, depend on the application-specific and the design constraint of technical scheme.The professional and technical personnel can use distinct methods to realize described function to each specific should being used for, but this realization should not thought and exceeds scope of the present invention.
The method of describing in conjunction with embodiment disclosed herein or the step of algorithm can directly use the software module of hardware, processor execution, and perhaps the combination of the two is implemented.Software module can place the storage medium of any other form known in random asccess memory (RAM), internal memory, read-only memory (ROM), electrically programmable ROM, electrically erasable ROM, register, hard disk, moveable magnetic disc, CD-ROM or the technical field.
To the above-mentioned explanation of the disclosed embodiments, make this area professional and technical personnel can realize or use the present invention.Multiple modification to these embodiment will be conspicuous concerning those skilled in the art, and defined herein General Principle can realize under the situation that does not break away from the spirit or scope of the present invention in other embodiments.Therefore, the present invention will can not be restricted to these embodiment shown in this article, but will meet and principle disclosed herein and features of novelty the wideest corresponding to scope.
Claims (9)
1. a blind signal bandwidth method of estimation is characterized in that, comprising:
Received signal counted according to default discrete Fourier transform DFT carry out conversion;
The intermediate value bandwidth DFT that obtains the bandwidth intermediate value correspondence in the signal bandwidth scope of system requirements counts, and counts according to described intermediate value bandwidth DFT that to calculate initial smoothing windows long;
Utilize the described DFT transformation results of the long smoothing processing of described initial smoothing windows, obtain the initial bandwidth DFT estimated value of counting;
Count estimated value and described intermediate value bandwidth DFT of described initial bandwidth DFT counted and compare, select the long factor of window according to comparative result, the span of the long factor of described window is for the inverse of the estimated value of counting greater than described initial bandwidth DFT, less than 1;
It is long to obtain the count variable smoothing windows of estimated value correspondence of initial bandwidth DFT after the long factor of described window is handled;
According to the described DFT transformation results of described variable smoothing windows long process;
Obtain with described variable smoothing windows long process after the corresponding signal bandwidth estimated value of result.
2. method according to claim 1 is characterized in that, describedly counts according to described intermediate value bandwidth DFT that to calculate initial smoothing windows long:
Determine that the sixth that described intermediate value bandwidth DFT counts is that initial smoothing windows is long.
3. method according to claim 2 is characterized in that, the described described DFT transformation results of the long smoothing processing of described initial smoothing windows of utilizing obtains the initial bandwidth DFT estimated value process of counting and comprises:
Result after the described smoothing processing is carried out normalization;
Choose the data after the thresholding intercepting normalization, obtain the initial bandwidth DFT estimated value of counting.
4. method according to claim 3 is characterized in that, count estimated value and described intermediate value bandwidth DFT of described initial bandwidth DFT counted and compare, and selects the process of the long factor of window to comprise according to comparative result:
The described initial bandwidth DFT estimated value and 1 of counting is deducted numerical value reciprocal multiplication after the described threshold value;
Whether more described multiplied result counts greater than described intermediate value bandwidth DFT, if, then choose value greater than median in the long factor span of described window, if not, then select value less than median in the long factor span of described window.
5. according to any described method among the claim 1-4, it is characterized in that the process of the signal bandwidth estimated value that described acquisition is corresponding with the result after the described variable smoothing windows long process comprises:
The DFT that obtains after the described variable smoothing windows long process counts;
Count according to described DFT and to calculate its corresponding signal bandwidth estimated value.
6. a blind signal bandwidth estimation unit is characterized in that, comprising:
The DFT converter unit is used for received signal counted according to default discrete Fourier transform DFT and carries out conversion;
The initial long acquiring unit of smoothing windows, the intermediate value bandwidth DFT that is used to obtain the interior bandwidth intermediate value correspondence of the signal bandwidth scope of system requirements counts, and counts according to described intermediate value bandwidth DFT that to calculate initial smoothing windows long;
The initial bandwidth DFT estimated value acquiring unit of counting is used to utilize the described DFT transformation results of the long smoothing processing of described initial smoothing windows, obtains the initial bandwidth DFT estimated value of counting;
The long factor selected cell of window, be used for count estimated value and described intermediate value bandwidth DFT of described initial bandwidth DFT counted and compare, select the long factor of window according to comparative result, the span of the long factor of described window is for the inverse of the estimated value of counting greater than described initial bandwidth DFT, less than 1;
The long acquiring unit of variable smoothing windows, it is long to be used to obtain the count variable smoothing windows of estimated value correspondence of initial bandwidth DFT after the long factor of described window is handled;
Variable smoothing windows long process unit is used for according to the described DFT transformation results of described variable smoothing windows long process;
Signal bandwidth estimated value acquiring unit, be used to obtain with described variable smoothing windows long process after the corresponding signal bandwidth estimated value of result.
7. device according to claim 6 is characterized in that, the described initial bandwidth DFT estimated value acquiring unit of counting comprises:
The normalization unit is used for the result after the described smoothing processing is carried out normalization;
The thresholding interception unit is used to choose the data after the thresholding intercepting normalization, obtains the initial bandwidth DFT estimated value of counting.
8. device according to claim 7 is characterized in that, the long factor selected cell of described window comprises:
The unit that multiplies each other is used for the described initial bandwidth DFT estimated value and 1 of counting is deducted numerical value reciprocal multiplication after the described threshold value;
The long factor determining unit of window, be used for more described multiplied result and whether count, if then choose value greater than median in the long factor span of described window greater than described intermediate value bandwidth DFT, if not, then select to grow the value of median in the factor span less than described window.
9. according to any described device among the claim 6-8, it is characterized in that described signal bandwidth estimated value acquiring unit comprises:
The DFT acquiring unit of counting, the DFT that is used to obtain after the described variable smoothing windows long process counts;
Signal bandwidth estimated value computing unit is used for counting according to described DFT and calculates its corresponding signal bandwidth estimated value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010274934 CN101909038B (en) | 2010-09-08 | 2010-09-08 | Method and device for estimating blind signal bandwidth |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 201010274934 CN101909038B (en) | 2010-09-08 | 2010-09-08 | Method and device for estimating blind signal bandwidth |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101909038A true CN101909038A (en) | 2010-12-08 |
| CN101909038B CN101909038B (en) | 2012-12-12 |
Family
ID=43264358
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 201010274934 Expired - Fee Related CN101909038B (en) | 2010-09-08 | 2010-09-08 | Method and device for estimating blind signal bandwidth |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101909038B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105916179A (en) * | 2016-05-30 | 2016-08-31 | 北京邮电大学 | Bandwidth estimation method and bandwidth estimation device |
| CN109507495A (en) * | 2018-10-17 | 2019-03-22 | 华北水利水电大学 | It is a kind of to become the long quasi- simultaneous interconnecting measurement method of parameters of window |
| CN112214895A (en) * | 2020-10-14 | 2021-01-12 | 中国船舶重工集团公司第七二四研究所 | Fourier multi-beam forming method with transform point number optimization design |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070070879A1 (en) * | 2005-09-27 | 2007-03-29 | Nec Corporation | Multi-user receiving apparatus converting SC-FDMA received signals of all users to signals in a frequency domain commonly |
| CN101485117A (en) * | 2006-07-05 | 2009-07-15 | 皇家飞利浦电子股份有限公司 | Bandwidth asymmetric communication system |
| CN101507219A (en) * | 2006-08-18 | 2009-08-12 | Nxp股份有限公司 | Time error estimation for data symbols |
| CN101682588A (en) * | 2007-03-27 | 2010-03-24 | 摩托罗拉公司 | Channel estimator for OFDM systems |
-
2010
- 2010-09-08 CN CN 201010274934 patent/CN101909038B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070070879A1 (en) * | 2005-09-27 | 2007-03-29 | Nec Corporation | Multi-user receiving apparatus converting SC-FDMA received signals of all users to signals in a frequency domain commonly |
| CN101485117A (en) * | 2006-07-05 | 2009-07-15 | 皇家飞利浦电子股份有限公司 | Bandwidth asymmetric communication system |
| CN101507219A (en) * | 2006-08-18 | 2009-08-12 | Nxp股份有限公司 | Time error estimation for data symbols |
| CN101682588A (en) * | 2007-03-27 | 2010-03-24 | 摩托罗拉公司 | Channel estimator for OFDM systems |
Non-Patent Citations (1)
| Title |
|---|
| 《清华大学学报(自然科学版)》 20040630 葛启宏等 《基于二维DFT 的多导频模式OFDM信道估计算法》 第44卷, 第6期 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105916179A (en) * | 2016-05-30 | 2016-08-31 | 北京邮电大学 | Bandwidth estimation method and bandwidth estimation device |
| CN105916179B (en) * | 2016-05-30 | 2019-07-26 | 北京邮电大学 | A kind of band width estimation method and device |
| CN109507495A (en) * | 2018-10-17 | 2019-03-22 | 华北水利水电大学 | It is a kind of to become the long quasi- simultaneous interconnecting measurement method of parameters of window |
| CN109507495B (en) * | 2018-10-17 | 2020-12-15 | 华北水利水电大学 | Variable-window-length quasi-synchronization grid-connected parameter measurement method |
| CN112214895A (en) * | 2020-10-14 | 2021-01-12 | 中国船舶重工集团公司第七二四研究所 | Fourier multi-beam forming method with transform point number optimization design |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101909038B (en) | 2012-12-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105743834B (en) | Reducing peak-to-average power ratio of OFDM signals | |
| CN114025379B (en) | Broadband multi-signal detection method, device and equipment | |
| CN110690931B (en) | Digital signal adaptive code rate estimation method and device based on multi-wavelet-base combination | |
| CN101909038B (en) | Method and device for estimating blind signal bandwidth | |
| US10715304B2 (en) | Modulation method and apparatus applicable to OvXDM system, and OvXDM system | |
| CN109490862A (en) | A kind of carrier frequency estimation method based on phase difference statistics spectrum | |
| CN106059973A (en) | Frequency offset estimation method and system | |
| US9197475B2 (en) | Base band processor with peak suppression function, transmitter with the same and method of transmitting signal | |
| CN111465051B (en) | SINR estimation method, system, electronic device and storage medium | |
| CN108900445B (en) | Method and device for estimating signal symbol rate | |
| CN111641420A (en) | Signal detection and acquisition method, device, receiver and storage medium | |
| US10805133B2 (en) | Method and apparatus for determining peak power, peak-to-average power ratio | |
| WO2021227217A1 (en) | Signal-to-noise ratio estimation method and apparatus, and electronic device and storage medium | |
| CN103297363A (en) | Symbol rate detecting device and symbol rate detecting method | |
| CN101917358B (en) | Method and device for estimating blind signal bandwidth | |
| CN110673139A (en) | Target detection radar system pulse compression coefficient generation method and device | |
| KR100457925B1 (en) | A frequency offset calculation method using log transform and linear approximation | |
| US6542101B1 (en) | Method and apparatus for performing analog-to-digital conversion using previous signal sample(s) | |
| CN107294625B (en) | Signal power estimation method and device and receiver | |
| CN112986922B (en) | Passive radar signal detection method based on sequential statistical filtering and binary detection | |
| CN108734188A (en) | A kind of clustering method, equipment and storage medium | |
| CN103873163A (en) | Sparseness self-adaptation compression spectrum sensing method based on asymptotic step length | |
| CN108169715B (en) | Method and system for determining in-phase channel phase imbalance degree and quadrature channel phase imbalance degree | |
| CN114124626B (en) | Signal noise reduction method and device, terminal equipment and storage medium | |
| KR101831198B1 (en) | Reduced complexity TWO-STEP TDOA/FDOA estimation method for communication signals |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121212 Termination date: 20200908 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |