CN106161298A - Non-linear ofdm system Iterative channel estimation based on compressed sensing - Google Patents
Non-linear ofdm system Iterative channel estimation based on compressed sensing Download PDFInfo
- Publication number
- CN106161298A CN106161298A CN201610512981.7A CN201610512981A CN106161298A CN 106161298 A CN106161298 A CN 106161298A CN 201610512981 A CN201610512981 A CN 201610512981A CN 106161298 A CN106161298 A CN 106161298A
- Authority
- CN
- China
- Prior art keywords
- channel
- linear
- channel estimation
- compressed sensing
- ofdm system
- 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.)
- Pending
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Noise Elimination (AREA)
Abstract
The present invention relates to a kind of non-linear ofdm system Iterative channel estimation based on compressed sensing, belong to wireless communication technology field, the purpose of this algorithm is for the OFDM condition of sparse channel estimation problem existed under the influence of non-linear distortion, system is carried out channel estimation, the precision of channel estimation of raising system, reduces the complexity of system.Being technically characterized in that of the present invention utilizes the dual openness of channel and nonlinear noise, first as observing matrix, pilot frequency information is compressed channel perception to estimate, gained channel information is regarded as observing matrix again and is compressed the estimation of perception non-linear distortion, and then signal is carried out nonlinear compensation, and loop iteration.The present invention considers the non-linear impact on system, estimates channel, be effectively increased the precision of channel estimation of non-linear ofdm system, improve operation efficiency under condition of sparse channel.
Description
Technical field
The invention belongs to wireless communication technology field, particularly to a kind of non-linear ofdm system based on compressed sensing repeatedly
For channel estimation method.
Background technology
OFDM (Orthogonal Frequency Division Multiplexing, OFDM) has frequency band
The advantages such as utilization rate symbol high, anti-crosstalk ability is strong, mitigating frequency-selective fading ability is strong, are the core in high-speed radiocommunication field
The heart.Channel estimation is one of key technology of ofdm system, and its precision directly affects the systematic function under fading channel.Channel is estimated
Meter and the two big key technologies that reduction peak-to-average force ratio is ofdm system.The precision that channel is estimated directly affects the system under fading channel
Performance, high peak-to-average power ratio cause signal through power amplifier produce non-linear distortion (nonlinear distortion,
NLD)。
At present, the existing researchs estimated about OFDM channel in a large number, wherein channel estimation studies based on auxiliary information is
For extensively, typical method has least square (least squares, LS), least mean-square error (minimum mean-squared
Error, MMSE), linear minimum mean-squared error (linear minimum mean-squared error, LMMSE) estimation etc..
In recent years, along with the proposition of compressive sensing theory, a large amount of scholars start to probe into and are applied to wireless communication field, and use it
OFDM channel in pilot aided is estimated to study.Compressed sensing is that a kind of Signal Compression for sparse signal is theoretical.Grind
Studying carefully and show, under condition of sparse channel, pilot tone based on equal number can reach preferably to estimate performance than traditional method.
But, existing channel estimation methods does not considers the nonlinear impact of ofdm signal mostly, it is believed that signal does not exists
Non-linear distortion or by reduce peak-to-average force ratio algorithm well solved.At present, the typical method reducing peak-to-average force ratio has pre-mistake
True method, partial transmission sequence method, Choose for user method etc., it is implemented in transmitters and majority has higher complexity.2003
Year, a kind of alternative manner based on signal reconstruction is suggested estimate in receivers and eliminate non-linear distortion.2005,
Tropp et al. proposes to be used for compressed sensing the impulsive noise of ofdm system and eliminates, and estimates impulsive noise based on pilot tone.2014
Year, A.Ghassemi et al. is it is further proposed that the non-linear distortion of ofdm signal is considered as sparse additive noise, available
It is estimated and compensates by compressed sensing in receivers.Said method can be avoided implementing to reduce peak-to-average force ratio in transmitter and calculate
Method, but all do not consider the impact that NLD is estimated by fading channel.
In non-linear ofdm system, cannot be carried out channel accurately based on the auxiliary information with non-linear distortion and estimate
Meter, and fading channel has a strong impact on NLD and estimates performance, channel is estimated and NLD estimates to influence each other, and mutually restriction becomes the two
The difficult point of Combined estimator.In recent years, the begin one's study channel of non-linear ofdm system of existing a small amount of scholar is asked with NLD Combined estimator
Topic.Such as based on iteration unified algorithm, carries out channel estimation based on LS and DFT interpolation, based on decision-feedback and signal reconstruction
Carry out NLD elimination.The estimated accuracy of the method is affected by pilot number and mistaken verdict, and needs known transmitter clipping operation
Prior information, receiver implement signal reconstruction add system complexity.The iterative joint of compressed sensing is combined based on LMMSE
Estimate also to be suggested, but LMMSE algorithm complex is higher, and need known channel prior information.Above-mentioned Combined estimator algorithm is the most not
Consider the sparse characteristic that channel is had under a lot of scenes.If channel is the most sparse, and pilot number is less than channel length, then
Pole is fallen sharply low by the channel estimating performance of LS, LMMSE.And under the influence of at home and abroad there is no at present about considering non-linear distortion
The research of Sparse Channel Estimation Algorithm.
Summary of the invention
It is an object of the invention to, for the OFDM condition of sparse channel estimation problem existed under the influence of non-linear distortion, utilize channel
Dual openness with NLD noise, proposes a kind of non-linear ofdm system Iterative channel estimation based on compressed sensing.Should
Algorithm has good estimation performance, it is achieved that the high-performance channel of non-linear ofdm system is estimated.
In order to reach above effect, the present invention non-linear ofdm system Iterative channel estimation based on compressed sensing,
Comprise the following steps:
Step 1: build the OFDM demodulation signal model existed when multipath channel affects with NLD;
Step 2: at receiving terminal, extracts pilot tone from the signal model of FFT and implements iterative estimate calculation based on pilot tone
Method.
Use non-linear ofdm system Iterative channel estimation based on compressed sensing, under condition of sparse channel, can be effective
Improve the precision of channel estimation of non-linear ofdm system, there is good estimation performance.
Accompanying drawing illustrates:
For clearer explanation inventive embodiments or technical scheme of the prior art, below will be to embodiment or existing
In technology description, the required accompanying drawing used is briefly described, and the accompanying drawing in describing below is only an enforcement of the present invention
Example, for those of ordinary skill in the art, do not pay creation laborious on the premise of, it is also possible to obtain it with reference to the accompanying drawings
His accompanying drawing.
Fig. 1 is the present invention non-linear ofdm system iterative channel estimation theory diagram;
Fig. 2 is that the channel of distinct methods under amplitude limit thresholding 4dB estimates mean square error;
Fig. 3 is the system bit error rate of distinct methods under amplitude limit thresholding 4dB;
Fig. 4 is that the channel under signal to noise ratio 30dB difference thresholding estimates mean square error;
Fig. 5 is the system bit error rate under signal to noise ratio 30dB difference thresholding.
Detailed description of the invention:
The purport of the present invention is to propose a kind of non-linear ofdm system Iterative channel estimation based on compressed sensing, profit
Estimate channel and non-linear distortion, and loop iteration by compressed sensing, can effectively reduce non-linear distortion and channel is estimated
Impact, reduce system complexity.
Below in conjunction with the accompanying drawings embodiment of the present invention is described in further detail.
One, the structure of ofdm system model
In an ofdm system, frequency domain signal X=[X0, X1..., XN-1]TTime-domain signal through IDFT modulation is represented by
Definition vector x=[x0, x1..., xN-1]T, then (1) formula is represented by matrix form
X=QX (2)
In formula, Q is N × N-dimensional inverse Fourier transform matrix,1≤m≤N, 1≤n≤N.
A length of N is added before N point datagCyclic Prefix form complete OFDM symbol through DAC export.
Signal produces NLD through power amplifier due to the restriction of its linearity, and the signal introducing NLD is represented by
x′n=β xn+dn (3)
In formula, β is for making xnWith dnIncoherent constant.
Definition vector d=[d0, d1..., dN-1]TFor the NLD introduced, it can regard sparse additive noise as.For analysis side
Just, making β in the present invention is 1, then the signal introducing NLD is represented by
X '=x+d=QX+d (4)
For making amplifier operation in the range of linearity, at transmitting terminal, signal is carried out amplitude limiting processing, then through the time domain of amplitude limit
Signal is represented by
In formula, A is amplitude limitation for normalization thresholding.
Affected by multipath channel, receive signal and be represented by
In formula,Represent cyclic convolution computing, z=[z0, z1..., zN-1]TRepresent additive white Gaussian noise, h=[h0,
h1..., hL-1, 01×(N-L+1)]T, wherein hl, l=0,1 ..., L-1 is channel impulse response (CIR).
It is represented by through DFT demodulation gained frequency-region signal
Definition vector Y=[Y0, Y1..., YN-1]T, then (7) formula is written as matrix form
Y=QHY=HX+HQHd+Z (8)
In formula, QHFor N × N-dimensional Fourier transform matrix, H=diag (H0, H1..., HN-1) it is channel frequency domain response vector
QHThe diagonal matrix of h deformation, Z=QHZ is the frequency domain form of white Gaussian noise.Formula (8) be exist simultaneously multipath channel with
OFDM demodulation signal model when NLD affects.
Two, receiving terminal extracts pilot tone and implements Iterative Method based on pilot tone
Specifically comprise the following steps that
Just estimate:
Nonlinear properties model (8) formula is deformed, vector X is rewritten as diagonal matrix form, and chooses pilot tone institute
Row, can obtain
In formula, i is iterations, Yp、ZpFor P × 1 dimensional vector, Xp'=diag (X1, X2..., XP) it is that P × P dimension is right
Angular moment battle array, P is pilot number.
Initialize i=0, and non-linear distortion now is regarded as noise, thus be integrated in noise item, then (9) formula can
It is written as
Based on compressed sensing principle, by known pilot information Xp' regard observing matrix as, willRegard orthogonal basis as, then can profit
With compressed sensing algorithm to the channel h in (10) formulaiSolve,
Algorithm completes the first estimation of channel.
Iterative estimate:
To the h estimatediCarry out Fourier transformation, obtain its frequency domain response Hi, by Hi(8) formula of substitution also chooses pilot tone row,
Can obtain
In formula, VpFor P × 1 dimensional vector.
Based on compressed sensing principle, willRegard observing matrix as, willRegard orthogonal basis as, utilize compressed sensing algorithm
To non-linear distortion d in (18) formulai+1Solve, can obtain
The d that will try to achievei+1(9) formula of substitution carries out nonlinear distortion compensation to signal, can obtain
In formula, UpFor p × 1 dimensional vector.
Recycling compressed sensing algorithm is to the h in (14) formulai+1Solve,
Algorithm completes interative computation for the first time.Circulation repetition (12), to the process of (15) formula, can carry out successive ignition fortune
Calculate, thus improve precision of channel estimation.The channel response utilizing iterative estimate to go out can be to carrying out nonlinear distortion compensation
Signal implements channel equalization.
Three, instance analysis
Build ofdm system for checking inventive algorithm performance, comprise 512 subcarriers, including 336 information
Carrier wave, 80 pilot sub-carriers, 96 null subcarrier, circulating prefix-length is 128, and 64QAM modulates, and 5/8LDPC encodes.Set up
Sparse multi-path channel model, channel length is 90, and wherein active path number is 12.
For embodying the performance of inventive algorithm, carrying out Performance comparision with existing algorithm in same system, wherein algorithm 1 is
A kind of Iterative Method, estimates channel based on LMMSE, estimates NLD noise based on compressed sensing;Algorithm 2 combines judgement for LS
The Joint iteration algorithm of feedback.
Fig. 2 is amplitude limit thresholding when being 4dB, and distinct methods channel under different signal to noise ratios estimates that mean square error (MSE) is bent
Line.It can be seen that along with the increase of signal to noise ratio, the estimated accuracy of each method is more and more higher.When signal to noise ratio is relatively low, the present invention calculates
The impact of method iterations is inconspicuous;And when signal to noise ratio is higher (>=17dB), along with the increase of iterations, precision of channel estimation
It is substantially better than the first estimation of not iteration, illustrates to consider necessity and the inventive algorithm that channel is estimated impact by non-linear distortion
For improving the superiority of nonlinear system channel estimating performance.Inventive algorithm and algorithm 1, algorithm 2 are in the feelings of equal iteration 3 times
Under condition, compared with assuming not consider the perfect channel estimation performance curve of NLD, inventive algorithm performance closest to.Due to letter
Road is condition of sparse channel, and channel length is more than pilot number, and the channel estimating performance of LS with LMMSE is poor.Simultaneously as son carries
Wave number is not the integral multiple of pilot number, and LS performance is inferior to LMMSE.
Fig. 3 is amplitude limit thresholding when being 4dB, distinct methods system bit error rate (BER) curve under different signal to noise ratios.Can
To find out, the channel estimating performance of nonlinear system all has certain gap with ideal linearity system, but inventive algorithm performance is
For close to ideal situation, and along with the increase of iterations, algorithm performance is significantly better than the first estimation of not iteration.In signal to noise ratio it is
During 30dB, the BER just estimated is 1 × 10-3, the BER that iteration is 3 times is 1 × 10-4。
Fig. 4 is signal to noise ratio when being 30dB, and distinct methods channel under different amplitude limit thresholdings estimates (MSE) curve.Permissible
Finding out, along with increasing of amplitude limit thresholding, the channel estimating performance of each method is become better and better.Under identical amplitude limit thresholding, this
At the beginning of bright algorithm, the performance of estimation, 2 iteration and 3 iteration gradually steps up.In three kinds of methods that 3 iteration are corresponding, the present invention calculates
Method best performance.
Fig. 5 is signal to noise ratio when being 30dB, distinct methods system (BER) curve under different amplitude limit thresholdings.Can from figure
To find out, along with the rising of threshold value, the systematic function performance that distinct methods is corresponding gradually steps up, and the property between various method
Can have same conclusion.When threshold value is 4dB, at the beginning of inventive algorithm, estimation, 2 iteration, the bit error rates of 3 iteration are divided
It is not 1 × 10-3、2×10-4、1×10-4, and it is significantly better than that algorithm 1, algorithm 2.
Sample result shows, this algorithm has good estimation performance, and under condition of sparse channel, the method can be effectively improved non-
The precision of channel estimation of linear ofdm system, good value estimated to ofdm system channel.
The foregoing is only presently preferred embodiments of the present invention, be not limiting as the present invention, all in the spirit and principles in the present invention
Within, any modification, equivalent substitution and improvement etc. made, should be included within the scope of the present invention.
Claims (3)
1. non-linear ofdm system Iterative channel estimation based on compressed sensing, comprises the following steps:
Step 1: build the OFDM demodulation signal model that there is multipath channel with non-linear distortion impact;
Step 2: at receiving terminal, extracts pilot tone from the signal model of FFT and implements Iterative Method based on pilot tone.
2. according to based on compressed sensing the non-linear ofdm system Iterative channel estimation described in right 1, it is characterised in that:
, there is OFDM demodulation signal model when multipath channel affects with non-linear distortion in step 1 simultaneously
Y=QHY=HX+HQHd+Z (1)
In formula, Y=[Y0, Y1..., YN-1]TFor the frequency domain response of receiving end signal, X=[X0, X1..., XN-1]TBelieve for transmitting terminal
Number frequency domain response, QHFor N × N-dimensional Fourier transform matrix, H=diag (H0, H1..., HN-1) it is channel frequency domain response vector
QHThe diagonal matrix of h deformation, h=[h0, h1..., hL-1, 01×(N-L+1)]T, wherein hl, l=0,1 ..., L-1 is channel
Shock response (CIR), Z=QHZ is the frequency domain form of white Gaussian noise.
3. according to based on compressed sensing the non-linear ofdm system Iterative channel estimation described in right 1, it is characterised in that:
Step 2, chooses the pilot tone row of transmitting terminal and receiving end signal,
In formula, i is iterations, Yp、ZpFor P × 1 dimensional vector, Xp'=diag (X1, X2..., XP) it is that P × P ties up angular moment
Battle array, P is pilot number.First by nonlinear noise merger noise item, compressed sensing is utilized to estimate the h of each iterationi, will estimate
HiBring formula (2) into, utilize compressed sensing to estimate di, by diBring formula (2) into and estimate the h of next iterationi+1, according to this
Plant the continuous iteration of mode to algorithmic statement.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610512981.7A CN106161298A (en) | 2016-06-30 | 2016-06-30 | Non-linear ofdm system Iterative channel estimation based on compressed sensing |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610512981.7A CN106161298A (en) | 2016-06-30 | 2016-06-30 | Non-linear ofdm system Iterative channel estimation based on compressed sensing |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106161298A true CN106161298A (en) | 2016-11-23 |
Family
ID=57349699
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610512981.7A Pending CN106161298A (en) | 2016-06-30 | 2016-06-30 | Non-linear ofdm system Iterative channel estimation based on compressed sensing |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106161298A (en) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106506430A (en) * | 2016-11-30 | 2017-03-15 | 黑龙江科技大学 | A kind of new algorithm of the compensation peak-to-average force ratio non-linear distortion based on compressed sensing technology |
CN113271269A (en) * | 2021-04-22 | 2021-08-17 | 重庆邮电大学 | Sparsity self-adaptive channel estimation method based on compressed sensing |
CN115086115A (en) * | 2022-06-10 | 2022-09-20 | 中国人民解放军空军工程大学 | Method for detecting and estimating nonlinear deformation signal caused by power amplifier |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102624658A (en) * | 2012-03-02 | 2012-08-01 | 清华大学 | Transmission method of time domain synchronous-orthogonal frequency division multiplexing (TDS-OFDM) based on theory of compressive sensing |
CN103152298A (en) * | 2013-03-01 | 2013-06-12 | 哈尔滨工业大学 | Blind signal reconstruction method based on distribution-type compressed sensing system |
CN103595414A (en) * | 2012-08-15 | 2014-02-19 | 王景芳 | Sparse sampling and signal compressive sensing reconstruction method |
CN104717162A (en) * | 2013-12-13 | 2015-06-17 | 天津工业大学 | OFDM ultra-wide band system nonlinear distortion restoring and channel estimation efficient uniting method |
-
2016
- 2016-06-30 CN CN201610512981.7A patent/CN106161298A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102624658A (en) * | 2012-03-02 | 2012-08-01 | 清华大学 | Transmission method of time domain synchronous-orthogonal frequency division multiplexing (TDS-OFDM) based on theory of compressive sensing |
CN103595414A (en) * | 2012-08-15 | 2014-02-19 | 王景芳 | Sparse sampling and signal compressive sensing reconstruction method |
CN103152298A (en) * | 2013-03-01 | 2013-06-12 | 哈尔滨工业大学 | Blind signal reconstruction method based on distribution-type compressed sensing system |
CN104717162A (en) * | 2013-12-13 | 2015-06-17 | 天津工业大学 | OFDM ultra-wide band system nonlinear distortion restoring and channel estimation efficient uniting method |
Non-Patent Citations (2)
Title |
---|
LI-JUN GE.ET: ""Joint Channel Estimation and Nonlinear Distortion Recovery Based on Compressed Sensing for OFDM Systems"", 《JOURNAL OF COMMUNICATIONS》 * |
MOHAMMAD MOHAMMADNIA-AVVAL.ET: ""Compressed Sensing Based Recovery of Nonlinearly Distorted OFDM Signals"", 《2011 IEEE INTERNATIONAL CONFERENCE ON COMMUNICATIONS》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106506430A (en) * | 2016-11-30 | 2017-03-15 | 黑龙江科技大学 | A kind of new algorithm of the compensation peak-to-average force ratio non-linear distortion based on compressed sensing technology |
CN106506430B (en) * | 2016-11-30 | 2019-10-11 | 黑龙江科技大学 | A kind of new algorithm of the compensation peak-to-average force ratio non-linear distortion based on compressed sensing technology |
CN113271269A (en) * | 2021-04-22 | 2021-08-17 | 重庆邮电大学 | Sparsity self-adaptive channel estimation method based on compressed sensing |
CN115086115A (en) * | 2022-06-10 | 2022-09-20 | 中国人民解放军空军工程大学 | Method for detecting and estimating nonlinear deformation signal caused by power amplifier |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101945066B (en) | Channel estimation method of OFDM/OQAM system | |
CN106789764A (en) | The transform domain quadratic estimate method of the denoising of joint Weighted Threshold and balanced judgement | |
CN101242388B (en) | Channel estimation method for high-speed single-carrier frequency domain balance ultra-wide broadband system | |
CN107181706B (en) | Frequency offset estimation and compensation method based on preamble symbol in constant envelope orthogonal frequency division multiplexing system | |
CN107360111A (en) | Impulsive Noise Mitigation Method in a kind of power line communication based on compressed sensing | |
CN107359906B (en) | Method for suppressing impulse noise in low-voltage power line communication system | |
CN106506430A (en) | A kind of new algorithm of the compensation peak-to-average force ratio non-linear distortion based on compressed sensing technology | |
CN104486267A (en) | Wavelet-denoising-based SC-FDE channel estimation method in short wave channel | |
CN110311872A (en) | A kind of underwater sound condition of sparse channel estimation method, system, equipment and storage medium | |
CN109729035A (en) | A kind of impulse noise suppression method in orthogonal frequency division multiplexi communication system | |
CN103873406A (en) | Inter-frame interference elimination method used for underwater sound orthogonal frequency-division multiplexing communication system | |
CN106161298A (en) | Non-linear ofdm system Iterative channel estimation based on compressed sensing | |
CN107682296A (en) | GFDM system high efficiency MMSE method of reseptances and device suitable for FSC | |
CN111355677A (en) | Multi-carrier underwater high-speed communication system based on filter bank | |
CN108737299A (en) | A kind of LMMSE channel estimation methods of low complex degree | |
CN100477651C (en) | High-performance OFDM channel estimation method based on combined pilot | |
CN108199988A (en) | Frequency domain GFDM low complex degree least mean-square error method of reseptances and receiver | |
CN102790746B (en) | Channel estimation method for OFDM (orthogonal frequency division multiplexing) system | |
Yang et al. | Delay-Doppler frequency domain-aided superimposing pilot OTFS channel estimation based on deep learning | |
CN101447969B (en) | Channel estimation method of multi-band orthogonal frequency division multiplexing ultra wide band system | |
CN103236993B (en) | A kind of channel estimation methods based on multipath delay profiles | |
CN102065035B (en) | Channel estimation method of multi-band orthogonal frequency-division multiplexing ultra-wideband system | |
CN105119856B (en) | Single carrier frequency domain equalization method based on sub-block analysis | |
CN115426224B (en) | Channel estimation method and system based on OTFS (optical transport plane) signal | |
CN108540187B (en) | Method and device for improving performance of nonlinear MU MIMO system based on MMSE (minimum mean square error) criterion |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20161123 |
|
WD01 | Invention patent application deemed withdrawn after publication |