CN109818888B - Group sparse underwater acoustic channel estimation method in pulse interference environment - Google Patents
Group sparse underwater acoustic channel estimation method in pulse interference environment Download PDFInfo
- Publication number
- CN109818888B CN109818888B CN201910228247.1A CN201910228247A CN109818888B CN 109818888 B CN109818888 B CN 109818888B CN 201910228247 A CN201910228247 A CN 201910228247A CN 109818888 B CN109818888 B CN 109818888B
- Authority
- CN
- China
- Prior art keywords
- sparsity
- error
- group
- time
- underwater acoustic
- 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
Links
Images
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention provides a group sparse underwater acoustic channel estimation method in a pulse interference environment. (1) Defining a new cost function; (2) calculating prior error and gain; (3) updating estimator tap coefficients and regularization parameters representing sparsity; (4) and continuously iterating until the error converges. The invention has the advantages that: (1) the p-norm of the estimation error is used as a cost function, so that the influence of pulse interference is effectively inhibited, and the robustness is strong; (2) the mixed norm is used as sparse constraint, and the group sparse characteristic of the channel is fully utilized; (3) the regularization parameters are updated in a self-adaptive mode, the size of applied constraint is adjusted continuously according to the change of the sparsity degree of the underwater acoustic channel, and the estimation performance is better.
Description
Technical Field
The invention relates to an underwater acoustic signal processing method, in particular to a group sparse underwater acoustic channel estimation method in a pulse interference environment based on an RLP algorithm.
Background
The underwater acoustic signals are severely damaged in the transmission process, and accurate estimation of the channel is necessary in order to accurately recover the transmitted signals at the receiving end. Many underwater acoustic channels exhibit significant group sparsity, which means that the channels are distributed in groups by a small number of non-zero elements in the time domain. In such channels, a channel estimation algorithm based on general sparsity does not have good estimation performance because the group sparsity characteristics are not fully utilized.
Noise generated by organisms, ice layer cracking and the like in an underwater acoustic channel has a lot of impulses in a time domain, a probability density function has a heavy tail characteristic and does not obey Gaussian distribution any more, and the noise is generally described by using a symmetrical stable alpha distribution (S alpha S) model. In gaussian noise, the performance is optimal by using the square of the channel estimation error as a cost function, but the performance is seriously reduced under an impulse interference environment. In order to obtain a more robust channel estimation method, the method under gaussian noise must be improved.
Disclosure of Invention
The invention aims to provide a group sparse underwater acoustic channel estimation method in a pulse interference environment, which has strong robustness and better estimation performance.
The purpose of the invention is realized as follows:
(1) defining a new cost function;
(2) calculating prior error and gain;
(3) updating estimator tap coefficients and regularization parameters representing sparsity;
(4) and continuously iterating until the error converges.
The present invention may further comprise:
1. the new cost function is defined by combining a mixed norm characterizing the sparse characteristics of the channel and an RLP algorithm resisting the influence of impulse interference.
2. The defining of the new cost function specifically includes:
the cost function is defined as:
wherein: n is the estimation time, lambda is the forgetting factor,estimate error for time i, yiFor the signal to be expected at time i,for channel estimates at time n-1, superscriptHDenotes the Hermite transposition, uiInput vector, γ, for time inFor the regularization parameter, f (e)i) In order to be a function of the loss,a convex function for representing the sparsity of the channel;
loss function f (e) of RLP algorithmi) Is composed of
Wherein: upper label*The expression is the conjugate of complex number, p is more than or equal to 1, xi and delta are threshold values, and the value is related to the signal-to-noise ratio;
for channels with group sparsity, choose/1,0And l2,0Mixed norm sparsity constraints, i.e.OrWherein G is the number of groups, beta is a parameter characterizing the size and range of the sparse constraint,to compriseG thiA vector of group taps.
The invention provides a group sparse underwater acoustic channel estimation method combining RLP and group sparse norm, aiming at the underwater acoustic channel estimation problem under the pulse interference environment. RLP defines different costs in different error ranges, so that the RLP has a good inhibiting effect on large errors generated by pulse interference and strong robustness; and mixed norm constraint is added, so that the invention can fully utilize the sparse natural characteristic of the underwater acoustic channel group to obtain lower steady-state estimation error.
The invention has the advantages that: (1) the p-norm of the estimation error is used as a cost function, so that the influence of pulse interference is effectively inhibited, and the robustness is strong; (2) the mixed norm is used as sparse constraint, and the group sparse characteristic of the channel is fully utilized; (3) the regularization parameters are updated in a self-adaptive mode, the size of applied constraint is adjusted continuously according to the change of the sparsity degree of the underwater acoustic channel, and the estimation performance is better.
Drawings
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a graph comparing the RLP algorithm with the recursive least squares algorithm (RLS);
FIG. 3 is a graph of RLP mean square error under different sparsity constraints;
fig. 4 is a table showing the updating process of the tap coefficient at each iteration of the present invention.
Detailed Description
The invention mainly comprises the following steps with reference to fig. 1:
(1) defining a system cost function by combining the channel sparse characteristic and the influence of pulse interference;
(2) calculating prior error and gain;
(3) updating estimator tap coefficients and regularization parameters representing sparsity;
(4) and continuously iterating until the error converges.
The group sparse channel estimation method in the present invention is described in detail below by way of example.
1、l1,0-RLP、l2,0Implementation of RLP
In order to overcome the influence of impulse interference and utilize the sparse characteristic of the channel, a cost function is defined as follows:
wherein n is the estimation time, lambda is the forgetting factor,estimate error for time i, yiFor the signal to be expected at time i,for channel estimates at time n-1, superscriptHDenotes the Hermite transposition, uiInput vector, γ, for time inFor the regularization parameter, f (e)i) In order to be a function of the loss,is a convex function characterizing the sparsity of the channel.
For the RLP algorithm, its penalty function f (e)i) Is composed of
Wherein, the upper label*The expression is the conjugate of complex number, p is more than or equal to 1, xi and delta are threshold values, and the value is related to the signal-to-noise ratio.
The RLP cost is different in different ranges, and large errors caused by impulse interference can be effectively resisted. For channels with group sparsity, choose/1,0And l2,0Mixed norm sparsity constraints, i.e.OrWherein G is the number of groups, beta is a parameter characterizing the size and range of the sparse constraint,to compriseG thiA vector of group taps.
Fig. 4 is a table showing the updating process of the tap coefficient at each iteration of the present invention.
2. Simulation research:
simulation conditions are as follows: the underwater sound sparse channel has the channel length of 64, the number of non-zero groups is 2, and each group comprises 4 taps. The non-zero positions are randomly distributed, and the sum of the squares of the amplitudes is 1. Adding impulse interference complying with the S alpha S distribution, wherein alpha is 1.6, and gamma is 0.07. The signal-to-noise ratio is 10 dB. The mean square error between the estimated value and the true value is used as a measure.
FIG. 2 is a graph showing the comparison of the RLP algorithm and the RLS in the steady-state error under the impulse interference environment. It can be seen that the traditional RLS algorithm cannot suppress the influence of pulse interference, and the steady-state error curve is not converged; the RLP algorithm in the invention can effectively inhibit the influence of impulse interference.
FIG. 3 shows the formula I1,0-RLP、l2,0Steady state error comparison of RLP with other sparse constraint algorithms. It can be seen that, compared to other algorithms, l in the present invention1,0-RLP、l2,0The RLP algorithm obtains lower steady-state error by fully utilizing the group sparsity of the channel, thereby proving the effectiveness of the invention.
Claims (1)
1. A group sparse underwater acoustic channel estimation method under a pulse interference environment is characterized by comprising the following steps: in order to obtain a more stable channel estimation method under the environment of pulse interference generated by the cracking of organisms and ice layers, the method under Gaussian noise is improved, and the method comprises the following steps:
(1) defining a new cost function; the method specifically comprises the following steps:
the cost function is defined as:
wherein: n is the estimation time, lambda is the forgetting factor,estimate error for time i, yiFor the signal to be expected at time i,for channel estimates at time n-1, superscriptHDenotes the Hermite transposition, uiInput vector, γ, for time inFor the regularization parameter, f (e)i) In order to be a function of the loss,convex function for characterizing channel sparsity, wherein loss function f (e) of RLP algorithmi) Comprises the following steps:
wherein: upper label*The expression is the conjugate of complex number, p is more than or equal to 1, xi and delta are threshold values, and the value is related to the signal-to-noise ratio;
for channels with group sparsity, choose/1,0And l2,0Mixed norm sparsity constraints, i.e.OrWherein G is the number of groups, beta is a parameter characterizing the size and range of the sparse constraint,to compriseG thiA vector of group taps;
(3) Updating estimator tap coefficientsAnd regularization parameters characterizing sparsityMn=λ-1(Mn-1-KnuH nMn-1);
(4) And continuously iterating until the error converges.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910228247.1A CN109818888B (en) | 2019-03-25 | 2019-03-25 | Group sparse underwater acoustic channel estimation method in pulse interference environment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910228247.1A CN109818888B (en) | 2019-03-25 | 2019-03-25 | Group sparse underwater acoustic channel estimation method in pulse interference environment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109818888A CN109818888A (en) | 2019-05-28 |
CN109818888B true CN109818888B (en) | 2022-01-14 |
Family
ID=66610243
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910228247.1A Active CN109818888B (en) | 2019-03-25 | 2019-03-25 | Group sparse underwater acoustic channel estimation method in pulse interference environment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109818888B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112653640B (en) * | 2020-12-18 | 2022-05-17 | 哈尔滨工程大学青岛船舶科技有限公司 | Pulse noise suppression underwater acoustic channel estimation method |
CN114218777B (en) * | 2021-12-08 | 2024-10-01 | 哈尔滨工程大学 | Polar region underwater ice sound channel estimation method |
CN114584433B (en) * | 2022-02-24 | 2024-05-28 | 哈尔滨工程大学 | Synchronous signal detection method in multi-path channel under impulse noise environment |
CN117131757B (en) * | 2023-07-10 | 2024-05-03 | 中国人民解放军海军工程大学 | Construction method and device of low-frequency nonlinear channel equalizer |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108919240A (en) * | 2018-04-23 | 2018-11-30 | 东南大学 | A kind of underwater acoustic target radiated noise modulation spectrum reconstruction method based on group sparsity structure |
CN109462427A (en) * | 2018-10-12 | 2019-03-12 | 南京信息工程大学 | A kind of MIMO underwater acoustic channel estimation method optimizing smooth L0 norm based on improved ADAPTIVE MIXED |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9384447B2 (en) * | 2014-05-22 | 2016-07-05 | The United States Of America As Represented By The Secretary Of The Navy | Passive tracking of underwater acoustic sources with sparse innovations |
CN108616314B (en) * | 2018-05-10 | 2020-01-21 | 宁波大学 | Underwater acoustic communication system impulse noise suppression method based on OFDM |
CN109257030B (en) * | 2018-10-22 | 2020-10-20 | 中原工学院 | Variable step length lpSparse system identification method of norm LMS algorithm |
CN109405954B (en) * | 2018-10-23 | 2022-01-14 | 哈尔滨工程大学 | UUV mobile underwater acoustic communication technology based on vector hydrophone array |
CN109474352B (en) * | 2018-12-24 | 2021-07-06 | 三亚哈尔滨工程大学南海创新发展基地 | Underwater sound orthogonal frequency division multiplexing communication method based on deep learning |
-
2019
- 2019-03-25 CN CN201910228247.1A patent/CN109818888B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108919240A (en) * | 2018-04-23 | 2018-11-30 | 东南大学 | A kind of underwater acoustic target radiated noise modulation spectrum reconstruction method based on group sparsity structure |
CN109462427A (en) * | 2018-10-12 | 2019-03-12 | 南京信息工程大学 | A kind of MIMO underwater acoustic channel estimation method optimizing smooth L0 norm based on improved ADAPTIVE MIXED |
Also Published As
Publication number | Publication date |
---|---|
CN109818888A (en) | 2019-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109818888B (en) | Group sparse underwater acoustic channel estimation method in pulse interference environment | |
CN105227512B (en) | A kind of impulsive noise method of estimation in OFDM underwater sound communication systems | |
CN108900443B (en) | A kind of underwater acoustic channel interference elimination method in underwater sound communication | |
CN105553898A (en) | Equalizer and feedback equalization method | |
CN107124379B (en) | Orthogonal wavelet normal-modulus blind equalization method based on improved wolf pack optimization | |
CN106506430A (en) | A kind of new algorithm of the compensation peak-to-average force ratio non-linear distortion based on compressed sensing technology | |
CN114338308A (en) | Shallow sea underwater sound MPSK signal blind demodulation method and system | |
CN109799484B (en) | Multipath inhibition method, system and storage medium for external radiation source radar system | |
CN108306837B (en) | Proportional MSER adaptive decision feedback equalization system and implementation method thereof | |
CN113300988A (en) | Inter-modal interference suppression method for low-frequency underwater acoustic communication | |
Alam et al. | Performance comparison of STFT, WT, LMS and RLS adaptive algorithms in denoising of speech signal | |
CN112462352B (en) | Line spectrum enhancement method suitable for low signal-to-noise ratio condition | |
CN102355435B (en) | Wavelet weighted multi-modulus blind equalization algorithm based on fractional lower order statistics (WT-FLOSWMMA) | |
Salman et al. | A sparse variable step-size LMS algorithm for impulsive noise | |
Salman et al. | A zero-attracting sparse lncosh adaptive algorithm | |
CN110677362B (en) | Complex domain underwater acoustic channel self-adaptive equalization method | |
CN116341109A (en) | Adaptive filtering method and system suitable for modeling of aircraft system | |
CN108667538A (en) | The combined estimation method of time delay and amplitude fading under complicated noise | |
CN110034747B (en) | Robust complex scale symbol adaptive filter | |
Youwen et al. | A variable multiple step-size LMS algorithm with l 0-norm | |
Zhang et al. | l 0-norm penalized shrinkage LMS algorithm based DFE for underwater acoustic communication | |
TWI385941B (en) | Method and apparatus for canceling channel interference | |
Hasan et al. | A Performance Comparison of Adaptive LMS, NLMS, RLS, and AFP Algorithms for Wireless Blind Channel Identification | |
CN111899751B (en) | Generalized mixed norm self-adaptive echo cancellation method for resisting saturation distortion | |
CN115118560B (en) | Single carrier frequency domain equalization method based on time domain adaptive threshold noise reduction |
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 |