CN116032310B - Signal self-adaptive detection reconstruction method based on channelized filtering - Google Patents

Signal self-adaptive detection reconstruction method based on channelized filtering Download PDF

Info

Publication number
CN116032310B
CN116032310B CN202310137163.3A CN202310137163A CN116032310B CN 116032310 B CN116032310 B CN 116032310B CN 202310137163 A CN202310137163 A CN 202310137163A CN 116032310 B CN116032310 B CN 116032310B
Authority
CN
China
Prior art keywords
signal
channel
filter
liquid crystal
crystal display
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.)
Active
Application number
CN202310137163.3A
Other languages
Chinese (zh)
Other versions
CN116032310A (en
Inventor
张羽
程岳云
余凯
李彦涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xi'an Hanbon Electronic Technology Co ltd
Original Assignee
Xi'an Hanbon Electronic Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xi'an Hanbon Electronic Technology Co ltd filed Critical Xi'an Hanbon Electronic Technology Co ltd
Priority to CN202310137163.3A priority Critical patent/CN116032310B/en
Publication of CN116032310A publication Critical patent/CN116032310A/en
Application granted granted Critical
Publication of CN116032310B publication Critical patent/CN116032310B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Abstract

The invention discloses a signal self-adaptive detection reconstruction method based on channelized filtering. And obtaining the amplitude of the down-conversion output signal of each channel by increasing or decreasing a stepping weighting method, obtaining the initial position of the channel with the signal according to the amplitude of the down-conversion output signal of each channel and a sliding window method, finding out that the same signal is decomposed in an adjacent channel by taking the initial position as a reference, occupying more channels with energy, and obtaining the channel number corresponding to each signal. And performing relevant modulation on the decomposed signals of the channels occupied by the same signal to obtain modulated decomposed signals. And obtaining a reconstruction signal according to a filtering rule of the comprehensive filtering by the modulated decomposition signal.

Description

Signal self-adaptive detection reconstruction method based on channelized filtering
Technical Field
The invention belongs to the field of signal processing, and particularly relates to a signal self-adaptive detection reconstruction method based on channelized filtering.
Background
The signal identification is widely applied to navigation communication, electronic countermeasure, investigation and detection and the like. The common digital signal identification is to perform time domain and frequency domain processing on the signal, extract the signal characteristics, then perform corresponding modulation on the signal, and then forward. On this basis, a polyphase filtering technique based on channelization is then created, which divides the digital receiver band into several equal parts, allowing the signal to be identified in a very wide band atmosphere. The common method is to calculate the amplitude of the signal of each channel, then identify the signal according to the amplitude, and modulate different signals. However, the conventional method has the following disadvantages: first, when the signal amplitude is calculated, the power operation is used, so that the consumption of resources is relatively large. Secondly, in the process of identifying the cross-channel signals, the signals cannot be identified effectively based on simple amplitude judgment, so that the efficiency of modulating and reconstructing the signals later is affected.
Disclosure of Invention
The invention aims to overcome the defects and provide a signal self-adaptive detection reconstruction method based on channelized filtering, which can calculate the amplitude of a down-conversion output signal according to a calculation method adopting an increasing and decreasing step weighting method. And according to the amplitude of the down-conversion output signal and by combining a channel merging method, the starting position and the ending position of the signal are effectively extracted, then a self-adaptive recognition method of the channel merging amplitude is adopted, after the signal occupying the most channels is determined, the signal is effectively recognized, then the decomposed sub-channel signals are modulated according to the starting position and the ending position of each signal, and the modulated sub-channel signals are reconstructed by using a comprehensive filter.
In order to achieve the above object, the method comprises the following steps:
s1, receiving a pulse signal;
s2, designing and analyzing a prototype filter and a comprehensive prototype filter;
s3, filtering the pulse signals by adopting an analysis prototype filter to obtain a plurality of filtered signals;
s4, calculating the amplitude of the filtered signal by adopting an increasing and decreasing step weighting method;
s5, judging whether a channel with corresponding amplitude has an initial signal or not by adopting a sliding window method, and obtaining a channel number set corresponding to an actual channel;
s6, modulating the filtered signals according to the channel number set to obtain modulated signals;
s7, adopting a comprehensive prototype filter to reconstruct and filter the modulated signal to obtain a reconstructed signal.
In S1, the expression of the pulse signal is:
Figure SMS_1
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure SMS_3
representing the pulse signal actually received,/->
Figure SMS_5
，/>
Figure SMS_7
，/>
Figure SMS_9
Is a set of natural numbers that are used to determine the number of the user,
Figure SMS_11
for slope, +>
Figure SMS_13
，/>
Figure SMS_14
Is bandwidth, when->
Figure SMS_2
Then point frequency is indicated, < >>
Figure SMS_4
For pulse width +.>
Figure SMS_6
For the sampling frequency +.>
Figure SMS_8
Is->
Figure SMS_10
Divided by->
Figure SMS_12
The remainder of (2).
S2, analyzing the prototype filter
Figure SMS_15
Is +.>
Figure SMS_16
The establishment method comprises the following steps:
establishing a low-pass analysis-based filter
Figure SMS_17
：
Figure SMS_18
Wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure SMS_20
，/>
Figure SMS_22
is a natural number set->
Figure SMS_23
，/>
Figure SMS_24
Is a positive integer greater than 2, +.>
Figure SMS_25
Indicating the length of the low-pass analysis-based filter, < >>
Figure SMS_26
Is an integer multiple of 8D or +.>
Figure SMS_27
，/>
Figure SMS_19
Is a positive integer greater than 4, +.>
Figure SMS_21
Is the sampling frequency;
order the
Figure SMS_28
，/>
Figure SMS_29
Is a low-pass analysis-based filter, +.>
Figure SMS_30
Windowing function for low-pass analysis-based filter,/->
Figure SMS_31
Normalization factor for low-pass analysis-based filter,/->
Figure SMS_32
。
In S3, the specific method for filtering the pulse signal by adopting the analysis prototype filter is as follows:
for analysis prototype filter
Figure SMS_33
Extracting and zero inserting to obtain analysis filter group +.>
Figure SMS_34
Let->
Figure SMS_35
The method comprises the following steps:
Figure SMS_36
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure SMS_37
is the slope;
at the same time to pulse signals
Figure SMS_38
Go->
Figure SMS_39
The method of extracting and inserting zero is as follows:
Figure SMS_40
then extracting and inserting zero pulse signal
Figure SMS_41
Analysis filter bank after extraction and zero insertion>
Figure SMS_42
Convolving to obtain a filtered signal +.>
Figure SMS_43
：
Figure SMS_44
Wherein the method comprises the steps of
Figure SMS_45
Representing convolution,/->
Figure SMS_46
Representing the decimation/interpolation multiple, +.>
Figure SMS_47
，/>
Figure SMS_48
，/>
Figure SMS_49
Is a positive integer greater than 2, +.>
Figure SMS_50
Is less than->
Figure SMS_51
Natural number of (3);
is provided with
Figure SMS_52
Is->
Figure SMS_53
The signal before the D-point inverse fast fourier transform is needed at the moment is expressed as follows:
Figure SMS_54
for a pair of
Figure SMS_55
Do->
Figure SMS_56
Inverse fast fourier transform of the point, get +.>
Figure SMS_57
Expressed as:
Figure SMS_58
order the
Figure SMS_59
Is->
Figure SMS_60
The signal after time down-conversion is represented as follows:
Figure SMS_61
obtaining
Figure SMS_62
。
S4, calculating the amplitude of the filtered signal by adopting an increasing and decreasing step weighting method comprises the following specific steps:
order the
Figure SMS_64
，/>
Figure SMS_65
Indicate->
Figure SMS_67
Down-converted output signal of individual channels, +.>
Figure SMS_69
Is->
Figure SMS_71
Real part of->
Figure SMS_73
Is->
Figure SMS_74
Imaginary part of->
Figure SMS_63
Is->
Figure SMS_66
At->
Figure SMS_68
Square of the amplitude of the moment, +.>
Figure SMS_70
Representing complex units, calculating the amplitude +.>
Figure SMS_72
The method of (2) is as follows:
in the first step, the first step is to provide,
Figure SMS_75
use->
Figure SMS_76
And->
Figure SMS_77
Comparing;
if it is
Figure SMS_78
Then->
Figure SMS_79
Doubling, add->
Figure SMS_80
Repeating the first step;
if it is
Figure SMS_81
，/>
Figure SMS_82
Halving, repeating the first step until +.>
Figure SMS_83
Entering a second step;
if it is
Figure SMS_84
Entering a second step;
second, calculating the amplitude
Figure SMS_85
The method comprises the following steps:
Figure SMS_86
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure SMS_87
is greater than->
Figure SMS_88
Step accumulated value, & gt>
Figure SMS_89
Is less than->
Figure SMS_90
Is initially of the step accumulated value of (a)
Figure SMS_91
，/>
Figure SMS_92
The accumulated step is represented, with an initial value of 1.
The specific method of S5 is as follows:
order the
Figure SMS_93
For channel marking, ++>
Figure SMS_94
Initial value->
Figure SMS_95
Judging the initial position of the channel; />
Figure SMS_96
Representing statistical parameters of signal channel, initially set +.>
Figure SMS_97
；/>
Figure SMS_98
Initially empty;
if there is a start signal, calculate the start position
Figure SMS_99
Corresponding end position->
Figure SMS_100
And channel number set +.>
Figure SMS_101
The method comprises the following steps:
the first step: judging channel mark
Figure SMS_102
If->
Figure SMS_103
Proceeding with the second step, otherwise->
Figure SMS_104
1, repeating the first step by self-adding;
and a second step of: calculate the first
Figure SMS_105
Start position of channel signal->
Figure SMS_106
Is provided with->
Figure SMS_108
For the current signal position, set +.>
Figure SMS_109
For decision threshold, value is taken according to variance of noise, let +.>
Figure SMS_110
For the amplitude +.>
Figure SMS_112
The judgment result of the threshold is carried out if the amplitude +.>
Figure SMS_113
Is greater than->
Figure SMS_107
Then->
Figure SMS_111
1, otherwise 0;
order the
Figure SMS_114
And->
Figure SMS_115
Respectively represent a left sliding window and a right sliding window, and the length is +.>
Figure SMS_116
The expression is as follows:
Figure SMS_117
Figure SMS_118
Figure SMS_119
Figure SMS_120
Figure SMS_121
Figure SMS_122
Figure SMS_123
Figure SMS_124
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure SMS_125
is greater than->
Figure SMS_126
；
Is provided with
Figure SMS_129
，/>
Figure SMS_130
Index of sliding window value for left sliding window or right sliding window, if +.>
Figure SMS_131
And->
Figure SMS_132
Figure SMS_133
，/>
Figure SMS_134
，
Figure SMS_135
Indicating the channel->
Figure SMS_127
The start position of a signal is found +.>
Figure SMS_128
；
For a pair of
Figure SMS_136
Amplitude +.>
Figure SMS_137
The starting position of the signal obtained by processing>
Figure SMS_138
For->
Figure SMS_139
Amplitude of (a) of (b)
Figure SMS_140
Processing to obtain the initial position of the signal +.>
Figure SMS_141
；
Start position of signal
Figure SMS_143
Start position of signal->
Figure SMS_144
Start position of sum signal->
Figure SMS_145
The minimum value of (2) is marked +.>
Figure SMS_146
First->
Figure SMS_147
Channel and->
Figure SMS_148
Channel up to->
Figure SMS_149
Before not finding, the respective initial positions are not judged, all are + ->
Figure SMS_142
；
And a third step of: calculation of
Figure SMS_150
：
Calculate the first
Figure SMS_151
Channel and->
Figure SMS_152
Sum of the modulus values of the corresponding time instants of the channels +.>
Figure SMS_153
First->
Figure SMS_154
Channel and->
Figure SMS_155
Sum of the modulus values of the corresponding time instants of the channels +.>
Figure SMS_156
The calculation method is as follows:
Figure SMS_157
Figure SMS_158
wherein the method comprises the steps of
Figure SMS_159
；
Figure SMS_160
To get in->
Figure SMS_161
Judging result of line threshold, if +.>
Figure SMS_162
Is greater than->
Figure SMS_163
Then->
Figure SMS_164
1, otherwise 0;
order the
Figure SMS_165
And->
Figure SMS_166
Respectively represent a left sliding window and a right sliding window, and the length is +.>
Figure SMS_167
The value is an integer multiple of 8, and is expressed as follows:
Figure SMS_168
Figure SMS_169
Figure SMS_170
Figure SMS_171
Figure SMS_172
Figure SMS_173
Figure SMS_174
Figure SMS_175
Figure SMS_176
is greater than->
Figure SMS_177
；
For the first
Figure SMS_178
The channel is judged as follows;
is provided with
Figure SMS_179
If->
Figure SMS_180
Channel presence->
Figure SMS_181
Then the end position of the signal is marked +.>
Figure SMS_182
；
For a pair of
Figure SMS_183
Processing to obtain the end position of the signal>
Figure SMS_184
；
End position of signal
Figure SMS_185
End position of sum signal->
Figure SMS_186
The maximum value is marked as +.>
Figure SMS_187
；
Fourth step: calculate the first
Figure SMS_188
Comparison factor of individual channels->
Figure SMS_189
And->
Figure SMS_190
Comparison factor of individual channels->
Figure SMS_191
And selecting the channel number corresponding to the largest value as the center channel number of the corresponding signal, the method is as follows:
Figure SMS_192
Figure SMS_193
comparison of
Figure SMS_194
And->
Figure SMS_195
；
If the maximum value is
Figure SMS_196
Then->
Figure SMS_197
，/>
Figure SMS_198
，/>
Figure SMS_199
At the same time make
Figure SMS_200
，/>
Figure SMS_201
；
If the maximum value is
Figure SMS_202
Then->
Figure SMS_203
，/>
Figure SMS_204
，/>
Figure SMS_205
At the same timeOrder the
Figure SMS_206
，/>
Figure SMS_207
，/>
Figure SMS_208
；
Recording device
Figure SMS_209
，/>
Figure SMS_210
，/>
Figure SMS_211
Add 1, if->
Figure SMS_212
And ending the judgment, otherwise, returning to the first step.
The specific method of S6 is as follows:
Figure SMS_213
and->
Figure SMS_214
Is a modulation function, and satisfies:
Figure SMS_215
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure SMS_216
index for signal start position->
Figure SMS_217
Indexing the signal end position;
then
Figure SMS_218
Wherein->
Figure SMS_219
Is less than->
Figure SMS_220
Natural number of (2)>
Figure SMS_221
Is a positive integer not greater than 3, +.>
Figure SMS_222
Is->
Figure SMS_223
The signal to be reconstructed which is subjected to amplitude modulation at the moment;
is provided with
Figure SMS_224
、/>
Figure SMS_225
And->
Figure SMS_226
Is 1 for each +.>
Figure SMS_227
The initial values are all 0;
the first step, the signal combination of the current channel is carried out, and the method is as follows:
if it is
Figure SMS_228
Then->
Figure SMS_229
Otherwise->
Figure SMS_230
Keeping the original value, and entering a second step;
in the second step, the second step is carried out,
Figure SMS_231
self-adding 1 if->
Figure SMS_232
Is greater than->
Figure SMS_233
Entering a third step, otherwise returning to the first step;
in the third step, the third step is that,
Figure SMS_234
self-adding 1 if->
Figure SMS_235
Is greater than->
Figure SMS_236
The modulation is completed.
In S7, the specific method of reconstruction filtering is as follows:
Figure SMS_237
for a pair of
Figure SMS_238
Time up-conversion processing result->
Figure SMS_239
Do->
Figure SMS_240
Fast fourier transform of the points, resulting in +.>
Figure SMS_241
Expressed as:
Figure SMS_242
obtaining
Figure SMS_243
；
Is provided with
Figure SMS_244
For integrated prototype filter->
Figure SMS_245
Extracting and zero inserting to obtain a reconstruction filter set +.>
Figure SMS_246
Order-making
Figure SMS_247
，/>
Figure SMS_248
The method comprises the following steps:
Figure SMS_249
by using
Figure SMS_250
And->
Figure SMS_251
Performing convolution
Figure SMS_252
For a pair of
Figure SMS_253
Go->
Figure SMS_254
Double interpolation to obtain +.>
Figure SMS_255
：
Figure SMS_256
Then reconstruct the signal
Figure SMS_257
The following are provided:
Figure SMS_258
wherein the method comprises the steps of
Figure SMS_259
Representing a natural number.
S2, synthesizing prototype filter
Figure SMS_260
Is +.>
Figure SMS_261
The establishment method comprises the following steps:
establishing a low-pass reconstruction basis filter
Figure SMS_262
：
Figure SMS_263
Wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure SMS_265
，/>
Figure SMS_266
is natural number (i.e.)>
Figure SMS_268
，/>
Figure SMS_269
Is a positive integer greater than 2, +.>
Figure SMS_270
Representing the length of the low-pass reconstructed basis filter, a->
Figure SMS_271
Is an integer multiple of 8D or +.>
Figure SMS_272
，/>
Figure SMS_264
Is a positive integer greater than 4, +.>
Figure SMS_267
Is the sampling frequency;
order the
Figure SMS_273
，/>
Figure SMS_274
Reconstructing the base filter for low pass, < >>
Figure SMS_275
A windowing function for a low-pass reconstructed basis filter,>
Figure SMS_276
reconstructing the normalization factor of the base filter for low pass,/->
Figure SMS_277
。
Compared with the prior art, the invention constructs the analysis prototype filter and the comprehensive prototype filter, extracts and interpolates the comprehensive prototype filter to obtain an analysis filter bank, and then respectively analyzes and filters the extracted input signals to obtain the filtering result of each channel. And obtaining the amplitude of the down-conversion output signal of each channel by increasing or decreasing a stepping weighting method, obtaining the initial position of the channel of the signal according to the amplitude of the down-conversion output signal of each channel and a sliding window method, finding out that the same signal is decomposed in an adjacent channel by taking the initial position as a reference, occupying more channels by energy, and obtaining the channel number corresponding to each signal. And performing relevant modulation on the decomposed signals of the channels occupied by the same signal to obtain modulated decomposed signals. And obtaining a reconstruction signal according to a filtering rule of the comprehensive filtering by the modulated decomposition signal. The invention can improve the efficiency of modulating and reconstructing the signals by taking the initial position of the signal channel as a reference, avoid the modulating process of different signals and greatly reduce the occupation of resources.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a signal time domain waveform at a signal-to-noise ratio of 5dB;
FIG. 3 is a signal spectrum at a signal-to-noise ratio of 5dB;
FIG. 4 is a graph of a spectrum corresponding to an analysis filter and a synthesis filter; wherein (a) is an analysis prototype filter and (b) is a synthesis prototype filter;
FIG. 5 is an amplitude plot of an input signal after amplitude modulation on each channel; wherein (a) is the amplitude decomposition of signal 1 at each channel; (b) decomposing the amplitude of the signal 2 in each channel;
FIG. 6 is a time domain diagram of the signal after modulation and reconstruction filtering;
fig. 7 is a spectrum of the signal after modulation and reconstruction filtering.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
Step one, generating pulse signals
Figure SMS_278
Let the bandwidth be->
Figure SMS_279
In Hz, sampling frequency is +.>
Figure SMS_280
(unit Hz):
Figure SMS_281
wherein the method comprises the steps of
Figure SMS_283
Representing the pulse signal actually received,/->
Figure SMS_284
，/>
Figure SMS_285
Indicating slope, & lt->
Figure SMS_286
，/>
Figure SMS_287
Indicates bandwidth, indicates pulse width, ">
Figure SMS_288
Is->
Figure SMS_289
Divided by->
Figure SMS_282
The remainder of (2).
Step two, designing the bandwidth as
Figure SMS_290
Is a prototype filter of (4)>
Figure SMS_291
，/>
Figure SMS_292
，/>
Figure SMS_293
Is->
Figure SMS_294
Positive integers greater than 2. Design Bandwidth is +.>
Figure SMS_295
Is>
Figure SMS_296
。
Design analysis prototype filter
Figure SMS_297
；
Establishing a low-pass analysis-based filter
Figure SMS_298
The expression is:
Figure SMS_299
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure SMS_301
，/>
Figure SMS_302
is a natural number set->
Figure SMS_303
，/>
Figure SMS_304
Is a positive integer greater than 2. />
Figure SMS_305
Indicating the length of the low-pass analysis-based filter, < >>
Figure SMS_306
Is an integer multiple of 8D or +.>
Figure SMS_307
，/>
Figure SMS_300
Typically a positive integer greater than 4.
Order the
Figure SMS_308
，/>
Figure SMS_309
Is a low-pass analysis-based filter, +.>
Figure SMS_310
Windowing function for low-pass analysis-based filter,/->
Figure SMS_311
For normalization factor->
Figure SMS_312
,/>
Figure SMS_313
Representing the filter length;
design of integrated prototype filter
Figure SMS_314
：
Establishing a low-pass reconstruction basis filter
Figure SMS_315
When the expression is:
Figure SMS_316
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure SMS_318
，/>
Figure SMS_319
is a natural number set->
Figure SMS_321
，/>
Figure SMS_322
Is a positive integer greater than 2, +.>
Figure SMS_323
Representing the length of the low-pass reconstructed basis filter, a->
Figure SMS_324
Is an integer multiple of 8D or +.>
Figure SMS_325
，/>
Figure SMS_317
Is a positive integer greater than 4, +.>
Figure SMS_320
For sampling frequency。
Order the
Figure SMS_326
，/>
Figure SMS_327
Reconstructing the base filter for low pass, < >>
Figure SMS_328
A windowing function for a low-pass reconstructed basis filter,>
Figure SMS_329
reconstructing the normalization factor of the base filter for low pass,/->
Figure SMS_330
；
Step three, adopting analysis prototype filter
Figure SMS_331
For pulse signal->
Figure SMS_332
Performing analysis filtering algorithm to obtain +.>
Figure SMS_333
The filtered signal->
Figure SMS_334
Wherein->
Figure SMS_335
Is less than->
Figure SMS_336
The natural number of (2) is expressed as follows:
is provided with
Figure SMS_337
Represents the decimation/interpolation multiple (wherein +.>
Figure SMS_338
) Analytical prototype filter->
Figure SMS_339
Extracting and zero inserting to obtain analysis filter group +.>
Figure SMS_340
Let->
Figure SMS_341
The method comprises the following steps:
Figure SMS_342
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure SMS_343
is the slope;
at the same time to pulse signals
Figure SMS_344
Go->
Figure SMS_345
The method of extracting and inserting zero is as follows:
Figure SMS_346
then extracting and inserting zero pulse signal
Figure SMS_347
Analysis filter bank after extraction and zero insertion>
Figure SMS_348
And (3) performing convolution:
Figure SMS_349
wherein the method comprises the steps of
Figure SMS_350
Representing convolution,/->
Figure SMS_351
Representing the decimation/interpolation multiple, +.>
Figure SMS_352
，/>
Figure SMS_353
，/>
Figure SMS_354
Is a positive integer greater than 2, +.>
Figure SMS_355
Is less than->
Figure SMS_356
Natural number of (a) is provided.
Is provided with
Figure SMS_357
Is->
Figure SMS_358
The signal before the D-point inverse fast fourier transform is needed at the moment is expressed as follows:
Figure SMS_359
for a pair of
Figure SMS_360
Do->
Figure SMS_361
Inverse fast fourier transform of the point, get +.>
Figure SMS_362
Expressed as:
Figure SMS_363
order the
Figure SMS_364
Is->
Figure SMS_365
The signal after time down-conversion is represented as follows:
Figure SMS_366
obtaining
Figure SMS_367
。
Step four, calculating by increasing and decreasing step weighting method
Figure SMS_368
Amplitude +.>
Figure SMS_369
The calculation method is as follows:
order the
Figure SMS_371
，/>
Figure SMS_372
Indicate->
Figure SMS_373
Down-converted output signal of individual channels, +.>
Figure SMS_375
Is->
Figure SMS_376
Real part of->
Figure SMS_377
Is->
Figure SMS_378
Imaginary part of->
Figure SMS_370
Representing plural units->
Figure SMS_374
Is->
Figure SMS_379
At->
Figure SMS_380
The square of the amplitude of the moment, the amplitude +.>
Figure SMS_381
The method of (2) is as follows:
is provided with
Figure SMS_382
Represents the accumulated value, initially +.>
Figure SMS_383
。/>
Figure SMS_384
Representing accumulation steps, wherein the initial value is 1;
a first step of,
Figure SMS_385
Use->
Figure SMS_386
And->
Figure SMS_387
A comparison is made.
Case one: if it is
Figure SMS_388
Then->
Figure SMS_389
Doubling, add->
Figure SMS_390
The first step is repeated. And a second case: if->
Figure SMS_391
，/>
Figure SMS_392
Halving, repeating the first step. Up to->
Figure SMS_393
And (5) entering a second step. Case three: if it is
Figure SMS_394
And (5) entering a second step.
Second, calculating the amplitude
Figure SMS_395
The method comprises the following steps:
Figure SMS_396
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure SMS_397
is greater than->
Figure SMS_398
Step accumulated value, & gt>
Figure SMS_399
Is less than->
Figure SMS_400
Is initially of the step accumulated value of (a)
Figure SMS_401
，/>
Figure SMS_402
The accumulated step is represented, with an initial value of 1.
Step five, to
Figure SMS_403
Judging to determine if the corresponding channel has initial signal to obtain the channel number set +.>
Figure SMS_404
The method comprises the following steps:
order the
Figure SMS_406
For channel marking, ++>
Figure SMS_407
Initial value->
Figure SMS_409
Judging the initial position of the channel; />
Figure SMS_410
Representing statistical parameters of signal channel, initially set +.>
Figure SMS_411
；/>
Figure SMS_412
Initially empty; if there is a start signal, calculate the start position +.>
Figure SMS_413
Corresponding end position->
Figure SMS_405
And channel number set +.>
Figure SMS_408
The method comprises the following steps:
the first step: judging channel mark
Figure SMS_414
If->
Figure SMS_415
Proceeding with the second step, otherwise->
Figure SMS_416
1, repeating the first step by self-adding;
and a second step of: calculate the first
Figure SMS_418
Start position of channel signal->
Figure SMS_420
Is provided with->
Figure SMS_421
For the current signal position, set +.>
Figure SMS_422
For decision threshold, value is taken according to variance of noise, let +.>
Figure SMS_423
For the amplitude +.>
Figure SMS_424
The judgment result of the threshold is carried out if the amplitude +.>
Figure SMS_425
Is greater than->
Figure SMS_417
Then->
Figure SMS_419
1, otherwise 0;
order the
Figure SMS_426
And->
Figure SMS_427
Respectively represent a left sliding window and a right sliding window, and the length is +.>
Figure SMS_428
The value is generally an integer multiple of 8, and is expressed as follows:
Figure SMS_429
Figure SMS_430
Figure SMS_431
Figure SMS_432
Figure SMS_433
Figure SMS_434
Figure SMS_435
Figure SMS_436
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure SMS_437
is greater than->
Figure SMS_438
；
Is provided with
Figure SMS_440
，/>
Figure SMS_442
Index of sliding window value for left sliding window or right sliding window, if +.>
Figure SMS_443
And->
Figure SMS_444
Figure SMS_445
，/>
Figure SMS_446
，
Figure SMS_447
Indicating the channel->
Figure SMS_439
Find a signalStart position of->
Figure SMS_441
；
Similarly, for
Figure SMS_449
The starting position of the signal can be obtained +.>
Figure SMS_450
For->
Figure SMS_452
The starting position of the signal can be obtained +.>
Figure SMS_454
. Start position of signal->
Figure SMS_456
Start position of signal->
Figure SMS_458
Start position of sum signal->
Figure SMS_459
The minimum value of (2) is marked +.>
Figure SMS_448
First->
Figure SMS_451
Channel and->
Figure SMS_453
Channel up to->
Figure SMS_455
Before not finding, the respective initial position determination is not performed, and the consistency is regarded as +.>
Figure SMS_457
。
And a third step of: calculation of
Figure SMS_460
：
Calculate the first
Figure SMS_461
Channel and->
Figure SMS_462
Sum of the modulus values of the corresponding time instants of the channels +.>
Figure SMS_463
First->
Figure SMS_464
Channel and->
Figure SMS_465
Sum of the modulus values of the corresponding time instants of the channels +.>
Figure SMS_466
The calculation method is as follows:
Figure SMS_467
Figure SMS_468
wherein the method comprises the steps of
Figure SMS_469
。
For the first
Figure SMS_471
Channel and->
Figure SMS_472
Sum of the modulus values of the corresponding time instants of the channels +.>
Figure SMS_475
，/>
Figure SMS_476
To get in->
Figure SMS_477
Judging result of line threshold, set +.>
Figure SMS_478
To get in->
Figure SMS_479
Judging result of line threshold, if +.>
Figure SMS_470
Is greater than->
Figure SMS_473
Then->
Figure SMS_474
1, otherwise 0.
Order the
Figure SMS_480
And->
Figure SMS_481
Respectively represent a left sliding window and a right sliding window, and the length is +.>
Figure SMS_482
The value is generally an integer multiple of 8, and is expressed as follows:
Figure SMS_483
Figure SMS_484
Figure SMS_485
Figure SMS_486
Figure SMS_487
Figure SMS_488
Figure SMS_489
Figure SMS_490
here, the
Figure SMS_491
Must be greater than +.>
Figure SMS_492
；
For the first
Figure SMS_493
The channel makes the following determination.
Is provided with
Figure SMS_494
If->
Figure SMS_495
Channel presence->
Figure SMS_496
Then the end position of the signal is marked +.>
Figure SMS_497
。
Similarly, pair
Figure SMS_498
Processing to obtain the end position of the signal>
Figure SMS_499
。
End position of signal
Figure SMS_500
End position of sum signal->
Figure SMS_501
The maximum value is marked as +.>
Figure SMS_502
。
Fourth step: calculate the first
Figure SMS_503
Comparison factor of individual channels->
Figure SMS_504
And->
Figure SMS_505
Comparison factor of individual channels->
Figure SMS_506
And selecting the channel number corresponding to the largest value as the center channel number of the corresponding signal, the method is as follows:
Figure SMS_507
Figure SMS_508
the following is the case:
case one: if the maximum value is
Figure SMS_509
，
Then
Figure SMS_510
，/>
Figure SMS_511
，/>
Figure SMS_512
At the same time let->
Figure SMS_513
，/>
Figure SMS_514
，/>
Figure SMS_515
Parameters are statistically calculated for the signal tracks.
And a second case: if the maximum value is
Figure SMS_516
。
Then
Figure SMS_517
，/>
Figure SMS_518
，/>
Figure SMS_519
At the same time let->
Figure SMS_520
，/>
Figure SMS_521
，/>
Figure SMS_522
。
Recording device
Figure SMS_523
，/>
Figure SMS_524
，/>
Figure SMS_525
Add 1, if->
Figure SMS_526
And ending the judgment, otherwise, returning to the first step.
Step six, according to the channel number set
Figure SMS_527
For->
Figure SMS_528
Performing corresponding modulation to obtain a modulated signal
Figure SMS_529
The treatment method comprises the following steps: />
Is provided with
Figure SMS_530
And->
Figure SMS_531
Is a modulation function and satisfies
Figure SMS_532
Wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure SMS_533
index for signal start position->
Figure SMS_534
Indexing the signal end position;
then
Figure SMS_535
Wherein->
Figure SMS_536
Is less than->
Figure SMS_537
Natural number of (2)>
Figure SMS_538
Is a positive integer not greater than 3, +.>
Figure SMS_539
Is->
Figure SMS_540
The signal to be reconstructed which is subjected to amplitude modulation at the moment comprises the following specific steps:
is provided with
Figure SMS_541
、/>
Figure SMS_542
And->
Figure SMS_543
Is 1 for each +.>
Figure SMS_544
The initial values are all 0;
the first step: the method for combining the signals of the current channel comprises the following steps:
if it is
Figure SMS_545
The->
Figure SMS_546
Otherwise->
Figure SMS_547
The original value is maintained. Then if->
Figure SMS_548
Less than 2->
Figure SMS_549
And 1 is added. And repeating the step, otherwise, entering the second step.
And a second step of:
Figure SMS_550
self-adding 1 if->
Figure SMS_551
Is greater than->
Figure SMS_552
And entering a third step, otherwise returning to the first step.
And a third step of:
Figure SMS_553
self-adding 1 if->
Figure SMS_554
Is greater than->
Figure SMS_555
The loop is jumped out.
Step seven, using integrated prototype filter
Figure SMS_556
For modulated signal->
Figure SMS_557
Performing reconstruction filtering to obtain a reconstruction signal +.>
Figure SMS_558
The method comprises the following steps:
using integrated prototype filters
Figure SMS_559
For modulated signal->
Figure SMS_560
The method for performing reconstruction filtering is as follows:
Figure SMS_561
for a pair of
Figure SMS_562
Time up-conversion processing result->
Figure SMS_563
Do->
Figure SMS_564
Fast fourier transform of the points, resulting in +.>
Figure SMS_565
Expressed as:
Figure SMS_566
obtaining
Figure SMS_567
；
Is provided with
Figure SMS_568
For integrated prototype filter->
Figure SMS_569
Extracting and zero inserting to obtain a reconstruction filter set +.>
Figure SMS_570
Order-making
Figure SMS_571
，/>
Figure SMS_572
The method comprises the following steps:
Figure SMS_573
reuse of
Figure SMS_574
And->
Figure SMS_575
Performing convolution
Figure SMS_576
/>
For a pair of
Figure SMS_577
Go->
Figure SMS_578
Double interpolation to obtain +.>
Figure SMS_579
：
Figure SMS_580
Then reconstruct the signal
Figure SMS_581
The following are provided:
Figure SMS_582
wherein the method comprises the steps of
Figure SMS_583
Representing a natural number.
Based on the above steps, a reconstructed signal after the modulation is obtained
Figure SMS_584
Examples:
referring to fig. 7, simulations were performed under Matlab 2007.
Input signal 1: pulse width 5us, pulse period 25us, bandwidth 50MHz, signal type is linear frequency modulation, digital intermediate frequency 93.75MHz, snr1=5 dB;
input signal 2: pulse width 5us, pulse period 25us, bandwidth 50MHz, signal type of chirped, digital intermediate frequency 203.125MHz, snr2=5 dB;
the synthesized input signal time domain diagram is shown in fig. 2, and the amplitude spectrum is shown in fig. 3:
reconstruction analysis prototype filter
Figure SMS_585
And integrated prototype filter->
Figure SMS_586
Is designed according to the design of (3).
If necessary, the length of the low-pass analysis base filter and the length of the low-pass reconstruction base filter are set to be the same, and the low-pass analysis base filter and the low-pass reconstruction base filter are both
Figure SMS_587
Windowing function of low-pass analysis-based filter +.>
Figure SMS_588
A hamming window is selected. Simultaneous make->
Figure SMS_589
Respectively obtain analysis filter group->
Figure SMS_590
（/>
Figure SMS_591
) And reconstruction Filter set->
Figure SMS_592
（/>
Figure SMS_593
）。
Analyzing prototype filters
Figure SMS_594
And integrated prototype filter->
Figure SMS_595
As shown in fig. 4a and 4 b:
amplitude of down-converted output signal after analysis of prototype filter
Figure SMS_596
(/>
Figure SMS_597
) As shown in fig. 4 a:
threshold parameters for judging the starting position and the ending position of each signal
Figure SMS_598
Length of sliding window->
Figure SMS_599
。
Referring to fig. 5a and 5b, let the modulation amplitude parameter of the input signal 1 be
Figure SMS_600
Amplitude modulation parameter of input signal 2 +.>
Figure SMS_601
. And modulates only the signal between the start position and the end position.
Modulating and comprehensively filtering to obtain a reconstructed signal
Figure SMS_602
The corresponding amplitude spectrum is shown in fig. 6. />

Claims (8)

1. The signal self-adaptive detection reconstruction method based on the channelized filtering is characterized by comprising the following steps of:
s1, receiving a pulse signal;
s2, designing and analyzing a prototype filter and a comprehensive prototype filter;
s3, filtering the pulse signals by adopting an analysis prototype filter to obtain a plurality of filtered signals;
s4, calculating the amplitude of the filtered signal by adopting an increasing and decreasing step weighting method, wherein the specific method comprises the following steps of:
order the
Figure QLYQS_3
，/>
Figure QLYQS_5
Indicate->
Figure QLYQS_8
Down-converted output signal of individual channels, +.>
Figure QLYQS_2
Is->
Figure QLYQS_7
Real part of->
Figure QLYQS_10
Is->
Figure QLYQS_11
Is used to determine the imaginary part of (c),
Figure QLYQS_1
is->
Figure QLYQS_6
At->
Figure QLYQS_9
Square of the amplitude of the moment, +.>
Figure QLYQS_12
Representing complex units, calculating the amplitude +.>
Figure QLYQS_4
The method of (2) is as follows:
in the first step, the first step is to provide,
Figure QLYQS_13
use->
Figure QLYQS_14
And->
Figure QLYQS_15
Comparing;
if it is
Figure QLYQS_16
Then->
Figure QLYQS_17
Doubling, add->
Figure QLYQS_18
Repeating the first step;
if it is
Figure QLYQS_19
，/>
Figure QLYQS_20
Halving, repeating the first step until +.>
Figure QLYQS_21
Entering a second step;
if it is
Figure QLYQS_22
Entering a second step;
second, calculating the amplitude
Figure QLYQS_23
The method comprises the following steps:
Figure QLYQS_24
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure QLYQS_25
is greater than->
Figure QLYQS_26
Step accumulated value, & gt>
Figure QLYQS_27
Is less than->
Figure QLYQS_28
Is initially of the step accumulated value of (a)
Figure QLYQS_29
，/>
Figure QLYQS_30
Representing accumulation steps, wherein the initial value is 1;
s5, judging whether a channel with corresponding amplitude has an initial signal or not by adopting a sliding window method, and obtaining a channel number set corresponding to an actual channel;
s6, modulating the filtered signals according to the channel number set to obtain modulated signals;
s7, adopting a comprehensive prototype filter to reconstruct and filter the modulated signal to obtain a reconstructed signal.
2. The method for adaptive detection and reconstruction of signals based on channelization filtering as claimed in claim 1, wherein in S1, the expression of the pulse signal is:
Figure QLYQS_31
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure QLYQS_34
representing the pulse signal actually received,/->
Figure QLYQS_37
，/>
Figure QLYQS_41
，/>
Figure QLYQS_33
Is a natural number set->
Figure QLYQS_36
For slope, +>
Figure QLYQS_40
，/>
Figure QLYQS_43
Is bandwidth, when->
Figure QLYQS_32
Then point frequency is indicated, < >>
Figure QLYQS_39
For pulse width +.>
Figure QLYQS_42
For the sampling frequency +.>
Figure QLYQS_44
Is->
Figure QLYQS_35
Divided by
Figure QLYQS_38
The remainder of (2).
3. The method for adaptive detection reconstruction of a signal based on channelization filtering as claimed in claim 1, wherein in S2, the prototype filter is analyzed
Figure QLYQS_45
Is +.>
Figure QLYQS_46
The establishment method comprises the following steps:
establishing a low-pass analysis-based filter
Figure QLYQS_47
：
Figure QLYQS_48
Wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure QLYQS_50
，/>
Figure QLYQS_54
is a natural number set->
Figure QLYQS_56
，/>
Figure QLYQS_51
Is a positive integer greater than 2, +.>
Figure QLYQS_52
Indicating the length of the low-pass analysis-based filter, < >>
Figure QLYQS_55
Is an integer multiple of 8D or +.>
Figure QLYQS_57
，/>
Figure QLYQS_49
Is a positive integer greater than 4, +.>
Figure QLYQS_53
Is the sampling frequency;
order the
Figure QLYQS_58
，/>
Figure QLYQS_59
Is a low-pass analysis-based filter, +.>
Figure QLYQS_60
Windowing function for low-pass analysis-based filter,/->
Figure QLYQS_61
Normalization factor for low-pass analysis-based filter,/->
Figure QLYQS_62
。
4. A signal adaptive detection reconstruction method based on channelized filtering according to claim 3, wherein in S3, the specific method of filtering the pulse signal by using the analysis prototype filter is as follows:
for analysis prototype filter
Figure QLYQS_63
Extracting and zero inserting to obtain analysis filter group +.>
Figure QLYQS_64
Let->
Figure QLYQS_65
The method comprises the following steps:
Figure QLYQS_66
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure QLYQS_67
is the slope;
at the same time to pulse signals
Figure QLYQS_68
Go->
Figure QLYQS_69
The method of extracting and inserting zero is as follows:
Figure QLYQS_70
then extracting and inserting zero pulse signal
Figure QLYQS_71
Analysis filter bank after extraction and zero insertion>
Figure QLYQS_72
Convolving to obtain a filtered signal +.>
Figure QLYQS_73
：
Figure QLYQS_74
Wherein the method comprises the steps of
Figure QLYQS_75
Representing convolution,/->
Figure QLYQS_76
Representing the decimation/interpolation multiple, +.>
Figure QLYQS_77
，/>
Figure QLYQS_78
，/>
Figure QLYQS_79
Is a positive integer greater than 2, +.>
Figure QLYQS_80
Is less than->
Figure QLYQS_81
Natural number of (3);
is provided with
Figure QLYQS_82
Is->
Figure QLYQS_83
The signal before the D-point inverse fast fourier transform is needed at the moment is expressed as follows:
Figure QLYQS_84
for a pair of
Figure QLYQS_85
Do->
Figure QLYQS_86
Inverse fast fourier transform of the point, get +.>
Figure QLYQS_87
Watch (Table)The method is shown as follows:
Figure QLYQS_88
order the
Figure QLYQS_89
Is->
Figure QLYQS_90
The signal after time down-conversion is represented as follows:
Figure QLYQS_91
obtaining
Figure QLYQS_92
。
5. The method for adaptively detecting and reconstructing a signal based on channelized filtering as set forth in claim 4, wherein the specific method of S5 is as follows:
order the
Figure QLYQS_93
For channel marking, ++>
Figure QLYQS_94
Initial value->
Figure QLYQS_95
Judging the initial position of the channel; />
Figure QLYQS_96
Representing statistical parameters of signal channel, initially set +.>
Figure QLYQS_97
；/>
Figure QLYQS_98
Initially empty;
if there is a start signal, calculate the start position
Figure QLYQS_99
Corresponding end position->
Figure QLYQS_100
And channel number set +.>
Figure QLYQS_101
The method comprises the following steps:
the first step: judging channel mark
Figure QLYQS_102
If->
Figure QLYQS_103
Proceeding with the second step, otherwise->
Figure QLYQS_104
1, repeating the first step by self-adding;
and a second step of: calculate the first
Figure QLYQS_106
Start position of channel signal->
Figure QLYQS_108
Is provided with->
Figure QLYQS_112
For the current signal position, set +.>
Figure QLYQS_107
For decision threshold, value is taken according to variance of noise, let +.>
Figure QLYQS_110
For the amplitude +.>
Figure QLYQS_111
The judgment result of the threshold is carried out if the amplitude +.>
Figure QLYQS_113
Greater than
Figure QLYQS_105
Then->
Figure QLYQS_109
1, otherwise 0;
order the
Figure QLYQS_114
And->
Figure QLYQS_115
Respectively represent a left sliding window and a right sliding window, and the length is +.>
Figure QLYQS_116
The expression is as follows:
Figure QLYQS_117
Figure QLYQS_118
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure QLYQS_119
is greater than->
Figure QLYQS_120
；
Is provided with
Figure QLYQS_122
，/>
Figure QLYQS_126
Index of sliding window value for left sliding window or right sliding window, if
Figure QLYQS_128
And->
Figure QLYQS_123
Figure QLYQS_125
，/>
Figure QLYQS_127
，
Figure QLYQS_129
Indicating the channel->
Figure QLYQS_121
The start position of a signal is found +.>
Figure QLYQS_124
；
For a pair of
Figure QLYQS_130
Amplitude +.>
Figure QLYQS_131
The starting position of the signal obtained by processing>
Figure QLYQS_132
For->
Figure QLYQS_133
Amplitude of (a) of (b)
Figure QLYQS_134
Processing to obtain the initial position of the signal +.>
Figure QLYQS_135
；
Start position of signal
Figure QLYQS_138
Start position of signal->
Figure QLYQS_139
Start position of sum signal->
Figure QLYQS_141
The minimum value of (2) is recorded as
Figure QLYQS_137
First->
Figure QLYQS_140
Channel and->
Figure QLYQS_142
Channel up to->
Figure QLYQS_143
Before not finding, the respective initial positions are not judged, all are + ->
Figure QLYQS_136
；
And a third step of: calculation of
Figure QLYQS_144
：
Calculate the first
Figure QLYQS_145
Channel and->
Figure QLYQS_146
Sum of the modulus values of the corresponding time instants of the channels +.>
Figure QLYQS_147
First->
Figure QLYQS_148
Channel and->
Figure QLYQS_149
Sum of the modulus values of the corresponding time instants of the channels +.>
Figure QLYQS_150
The calculation method is as follows:
Figure QLYQS_151
wherein the method comprises the steps of
Figure QLYQS_152
；
Figure QLYQS_153
To get in->
Figure QLYQS_154
Judging result of line threshold, if +.>
Figure QLYQS_155
Is greater than->
Figure QLYQS_156
Then->
Figure QLYQS_157
1, otherwise 0;
order the
Figure QLYQS_158
And->
Figure QLYQS_159
Respectively represent a left sliding window and a right sliding window, and the length is +.>
Figure QLYQS_160
The value is an integer multiple of 8, and is expressed as follows:
Figure QLYQS_161
Figure QLYQS_162
Figure QLYQS_163
is greater than->
Figure QLYQS_164
；
For the first
Figure QLYQS_165
The channel is judged as follows;
is provided with
Figure QLYQS_166
If->
Figure QLYQS_167
Channel presence->
Figure QLYQS_168
Then the end position of the signal is marked +.>
Figure QLYQS_169
；
For a pair of
Figure QLYQS_170
Processing to obtain the end position of the signal>
Figure QLYQS_171
；
End position of signal
Figure QLYQS_172
End position of sum signal->
Figure QLYQS_173
The maximum value is marked as +.>
Figure QLYQS_174
；
Fourth step: calculate the first
Figure QLYQS_175
Comparison factor of individual channels->
Figure QLYQS_176
And->
Figure QLYQS_177
Comparison factor of individual channels->
Figure QLYQS_178
And selecting the channel number corresponding to the largest value as the center channel number of the corresponding signal, the method is as follows:
Figure QLYQS_179
comparison of
Figure QLYQS_180
And->
Figure QLYQS_181
；
If the maximum value is
Figure QLYQS_182
Then->
Figure QLYQS_183
，/>
Figure QLYQS_184
，/>
Figure QLYQS_185
At the same time let->
Figure QLYQS_186
，
Figure QLYQS_187
；
If the maximum value is
Figure QLYQS_188
Then->
Figure QLYQS_189
，/>
Figure QLYQS_190
，/>
Figure QLYQS_191
At the same time make
Figure QLYQS_192
，/>
Figure QLYQS_193
，/>
Figure QLYQS_194
；
Recording device
Figure QLYQS_195
，/>
Figure QLYQS_196
，/>
Figure QLYQS_197
Add 1, if->
Figure QLYQS_198
And ending the judgment, otherwise, returning to the first step.
6. The method for adaptively detecting and reconstructing a signal based on channelized filtering as set forth in claim 5, wherein the specific method of S6 is as follows:
Figure QLYQS_199
and->
Figure QLYQS_200
Is a modulation function, and satisfies:
Figure QLYQS_201
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure QLYQS_202
index for signal start position->
Figure QLYQS_203
Indexing the signal end position;
then
Figure QLYQS_204
Wherein->
Figure QLYQS_205
Is less than->
Figure QLYQS_206
Natural number of (2)>
Figure QLYQS_207
Is a positive integer not greater than 3,
Figure QLYQS_208
is->
Figure QLYQS_209
The signal to be reconstructed which is subjected to amplitude modulation at the moment;
is provided with
Figure QLYQS_210
、/>
Figure QLYQS_211
And->
Figure QLYQS_212
Is 1 for each +.>
Figure QLYQS_213
The initial values are all 0;
the first step, the signal combination of the current channel is carried out, and the method is as follows:
if it is
Figure QLYQS_214
Then->
Figure QLYQS_215
Otherwise->
Figure QLYQS_216
Keeping the original value, and entering a second step;
in the second step, the second step is carried out,
Figure QLYQS_217
self-adding 1 if->
Figure QLYQS_218
Is greater than->
Figure QLYQS_219
Entering a third step, otherwise returning to the first step;
in the third step, the third step is that,
Figure QLYQS_220
self-adding 1 if->
Figure QLYQS_221
Is greater than->
Figure QLYQS_222
The modulation is completed.
7. The method for adaptively detecting and reconstructing a signal based on channelized filtering as set forth in claim 6, wherein in S7, the specific method for reconstructing the filtering is as follows:
Figure QLYQS_223
for a pair of
Figure QLYQS_224
Time up-conversion processing result->
Figure QLYQS_225
Do->
Figure QLYQS_226
Fast fourier transform of the points, resulting in +.>
Figure QLYQS_227
Expressed as:
Figure QLYQS_228
obtaining
Figure QLYQS_229
；
Is provided with
Figure QLYQS_230
For integrated prototype filter->
Figure QLYQS_231
Extracting and zero inserting to obtain a reconstruction filter set +.>
Figure QLYQS_232
Let->
Figure QLYQS_233
，
Figure QLYQS_234
The method comprises the following steps:
Figure QLYQS_235
by using
Figure QLYQS_236
And->
Figure QLYQS_237
Performing convolution
Figure QLYQS_238
For a pair of
Figure QLYQS_239
Go->
Figure QLYQS_240
Double interpolation to obtain +.>
Figure QLYQS_241
：
Figure QLYQS_242
Then reconstruct the signal
Figure QLYQS_243
The following are provided:
Figure QLYQS_244
wherein the method comprises the steps of
Figure QLYQS_245
Representing a natural number.
8. The method for adaptive detection and reconstruction of signal based on channelization filtering as claimed in claim 1, wherein in S2, the prototype filter is synthesized
Figure QLYQS_246
Is +.>
Figure QLYQS_247
The establishment method comprises the following steps:
establishing a low-pass reconstruction basis filter
Figure QLYQS_248
：
Figure QLYQS_249
Wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure QLYQS_250
，/>
Figure QLYQS_254
is natural number (i.e.)>
Figure QLYQS_256
，/>
Figure QLYQS_251
Is a positive integer greater than 2, +.>
Figure QLYQS_253
Representing the length of the low-pass reconstructed basis filter, a->
Figure QLYQS_257
Is an integer multiple of 8D or +.>
Figure QLYQS_258
，/>
Figure QLYQS_252
Is a positive integer greater than 4, +.>
Figure QLYQS_255
Is the sampling frequency;
order the
Figure QLYQS_259
，/>
Figure QLYQS_260
Reconstructing the base filter for low pass, < >>
Figure QLYQS_261
A windowing function for a low-pass reconstructed basis filter,>
Figure QLYQS_262
reconstructing the normalization factor of the base filter for low pass,/->
Figure QLYQS_263
。
CN202310137163.3A 2023-02-20 2023-02-20 Signal self-adaptive detection reconstruction method based on channelized filtering Active CN116032310B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310137163.3A CN116032310B (en) 2023-02-20 2023-02-20 Signal self-adaptive detection reconstruction method based on channelized filtering

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310137163.3A CN116032310B (en) 2023-02-20 2023-02-20 Signal self-adaptive detection reconstruction method based on channelized filtering

Publications (2)

Publication Number Publication Date
CN116032310A CN116032310A (en) 2023-04-28
CN116032310B true CN116032310B (en) 2023-07-04

Family

ID=86077609

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310137163.3A Active CN116032310B (en) 2023-02-20 2023-02-20 Signal self-adaptive detection reconstruction method based on channelized filtering

Country Status (1)

Country Link
CN (1) CN116032310B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108680786A (en) * 2018-04-11 2018-10-19 东南大学 A kind of pulse signal frequency domain adaptive filtering envelope extraction method
CN111835434A (en) * 2019-04-19 2020-10-27 深圳市鼎阳科技有限公司 Method and device for measuring broadband frequency response
US11184042B1 (en) * 2020-08-13 2021-11-23 Bae Systems Information And Electronic Systems Integration Inc. Partial band signal reconstruction using arbitrary numbers of synthesis channels

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6792057B2 (en) * 2002-08-29 2004-09-14 Bae Systems Information And Electronic Systems Integration Inc Partial band reconstruction of frequency channelized filters
EP2857860A4 (en) * 2012-06-05 2015-06-17 Panasonic Ip Man Co Ltd Signal processing device
CN109714146B (en) * 2018-01-16 2020-05-12 电子科技大学 Sliding window-based quadratic correlation frame synchronization method
CN109561035B (en) * 2018-12-21 2021-02-12 电子科技大学 Terahertz-based ultra-wideband single-channel signal blind detection method
CN110469326B (en) * 2019-08-02 2023-05-12 中石化石油工程技术服务有限公司 Self-adaptive correlation filtering method
CN112350691A (en) * 2020-11-02 2021-02-09 南京理工大学 Geomagnetic weak signal detection algorithm based on multi-stage filtering extraction and two-stage WOLA structure
CN113341378B (en) * 2021-04-13 2022-08-16 北京理工大学 Self-adaptive channelized receiving method based on frequency spectrum differential entropy detection
CN113627398B (en) * 2021-10-11 2022-01-04 西安瀚博电子科技有限公司 Signal characteristic detection method based on adaptive reconstruction filtering

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108680786A (en) * 2018-04-11 2018-10-19 东南大学 A kind of pulse signal frequency domain adaptive filtering envelope extraction method
CN111835434A (en) * 2019-04-19 2020-10-27 深圳市鼎阳科技有限公司 Method and device for measuring broadband frequency response
US11184042B1 (en) * 2020-08-13 2021-11-23 Bae Systems Information And Electronic Systems Integration Inc. Partial band signal reconstruction using arbitrary numbers of synthesis channels

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Unscented Kalman Filter for frequency and amplitude estimation;Novanda;Filter amplitude conver;全文 *

Also Published As

Publication number Publication date
CN116032310A (en) 2023-04-28

Similar Documents

Publication Publication Date Title
CN101882964B (en) De-noising method of transient electromagnetic detecting echo signal
Herff et al. Extracting features from time series
CN109257068B (en) Adaptive filtering method for power line communication signals
CN108880621B (en) Adaptive filtering method for power line communication signals
CN103199945A (en) Recognition method for cognizing radio signal modulation modes under low-signal-to-noise-ratio condition
CN112702293A (en) Signal parameter analysis method and system based on cyclic spectrum
CN113341378B (en) Self-adaptive channelized receiving method based on frequency spectrum differential entropy detection
Phinyomark et al. EMG denoising estimation based on adaptive wavelet thresholding for multifunction myoelectric control
CN108737302A (en) The symbol rate estimation method and its device of accidental resonance joint wavelet transformation under Low SNR
CN114422311B (en) Signal modulation recognition method and system combining deep neural network and expert priori features
CN112380939A (en) Deep learning signal enhancement method based on generation countermeasure network
CN107123431A (en) A kind of underwater sound signal noise-reduction method
CN116032310B (en) Signal self-adaptive detection reconstruction method based on channelized filtering
CN109450829B (en) Method and apparatus for estimating code rate of digital modulation signal
CN102664843A (en) Method and device for single-frequency interference detection and elimination
Daldal A novel demodulation method for quadrate type modulations using a hybrid signal processing method
CN109039379B (en) Power line communication signal filtering method
Türk et al. Classification of chaos-based digital modulation techniques using wavelet neural networks and performance comparison of wavelet families
CN115378776A (en) MFSK modulation identification method based on cyclic spectrum parameters
CN109586763B (en) Denoising method and denoising system for spread spectrum signal in electronic communication system
CN111814703B (en) HB-based signal joint feature extraction method under non-reconstruction condition
Sharabaty et al. Alpha and theta wave localisation using hilbert-huang transform: empirical study of the accuracy
Berthelot et al. Filtering-based analysis comparing the DFA with the CDFA for wide sense stationary processes
Hadhami et al. Speech denoising based on empirical mode decomposition and improved thresholding
Zhang et al. Simultaneous denoising and preserving of seismic signals by multiscale time-frequency peak filtering

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: A Signal Adaptive Detection and Reconstruction Method Based on Channelized Filtering

Effective date of registration: 20230814

Granted publication date: 20230704

Pledgee: Xi'an innovation financing Company limited by guarantee

Pledgor: XI'AN HANBON ELECTRONIC TECHNOLOGY CO.,LTD.

Registration number: Y2023980052035

PE01 Entry into force of the registration of the contract for pledge of patent right