WO1998038772A1 - Estimation de la reponse impulsionnelle d'un canal au moyen de la decomposition d'une valeur singuliere - Google Patents
Estimation de la reponse impulsionnelle d'un canal au moyen de la decomposition d'une valeur singuliere Download PDFInfo
- Publication number
- WO1998038772A1 WO1998038772A1 PCT/FI1998/000159 FI9800159W WO9838772A1 WO 1998038772 A1 WO1998038772 A1 WO 1998038772A1 FI 9800159 W FI9800159 W FI 9800159W WO 9838772 A1 WO9838772 A1 WO 9838772A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- cir
- matrix
- eigenvalues
- present
- diagonal matrix
- Prior art date
Links
- 238000000354 decomposition reaction Methods 0.000 title claims description 4
- 239000011159 matrix material Substances 0.000 claims abstract description 58
- 238000000034 method Methods 0.000 claims abstract description 15
- 238000004891 communication Methods 0.000 claims abstract description 5
- 230000001131 transforming effect Effects 0.000 claims abstract description 3
- 238000012549 training Methods 0.000 abstract description 9
- 230000000694 effects Effects 0.000 abstract description 5
- 230000003247 decreasing effect Effects 0.000 abstract description 3
- 238000003012 network analysis Methods 0.000 abstract description 2
- 239000013598 vector Substances 0.000 description 7
- 238000013459 approach Methods 0.000 description 5
- 230000001143 conditioned effect Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000007792 addition Methods 0.000 description 1
- 238000010420 art technique Methods 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
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/0212—Channel estimation of impulse response
-
- 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
- H04L25/0242—Channel estimation channel estimation algorithms using matrix methods
- H04L25/0248—Eigen-space methods
Definitions
- the present invention relates to the estimating of CIR (channel impulse response) and SINR (signal-to-i ⁇ terference-plus-noise ratio) in CDMA receivers.
- the present invention is directed to a method of determining the channel impulse response (CIR) of a communication system, such as the CIR of radio channels of a digital mobile radio network (GSM network).
- the present invention relates to determining the CIR based on the reception of a known training sequence.
- the present invention also has application in interference cancellation and use in CDMA receivers.
- Synchronisation bursts (S3) are a useful portion of the signal.
- the SB are transmitted on at least one channel from every base station, and they are transmitted in a regular pattern. Decoding of the GSM protocols is not necessary. Both the data in the SB and the pattern with which they occur is fixed and substantially identical for all base stations.
- the advantage of using SB for determining the CIR is that they represent a relatively long, noise-like predetermined transmitted signal. Typically, 64 bits are transmitted over a period of 237 ⁇ sec. The determination of the CIR thus requires sufficient synchronisation to the bursts in the received signal to enable the extraction of the SB which is then processed to determine the CIR.
- the processing is done by using estimation techniques.
- a known training sequence is transmitted S (t), and this is corrupted by a communications channel producing the received signal S ra (t).
- the problem in estimating the CIR is to determine tap-weights ⁇ of a FIR filter (which is used to approximate the CIR), so that the known S tx (t) after passing through the filter is as close as possible to the received signal S, x (t).
- the tap weights ⁇ i are determined by simple correlation as
- This algorithm relies on the noise-like properties of the transmitted signal whereby its autocorrelation function should have low time sidelobes.
- the sidelobe performance of this algorithm has been found to be limited due to the correlation properties of the "clean" part of the training sequence and the fact that only partial correlations are performed for later weights.
- US 5,473,632 produces a matrix.
- the problem really relates to what is done to the matrix to calculate the CIR.
- a "factor ⁇ " is added to the matrix before the matrix is inverted, and then after the
- factor ⁇ is added, the matrix is inverted. In practice it has been found that the inversion of the matrix causes problems in calculating the CIR.
- the present invention seeks to alleviate the problems experienced in determining CIR and SINR with prior art techniques.
- the present invention stems from the realisation that the problems associated with the prior art can be alleviated by initially transforming the matrix obtained by virtue of network analysis into a 'diagonal matrix' and then inverting the resultant diagonal matrix. Thus there is no "factor ⁇ " added.
- An advantage of using a diagonal matrix is that it is possible to remove or ignore some eigenvalues (such as some" of the small problematic eigenvalues) and in this way avoid the magnified errors and noise of US 5,473,632 when the matrix is inverted.
- Another advantage of using a diagonal matrix is attributable to the nature of a diagonal matrix and that is, the matrix can be inverted by inverting only the diagonal elements. This exposes the difficulty caused by small eigenvalues - they lead to very large elements in the inverse.
- the present application discloses a still further invention which is directed to alleviating the problem of the production of large side lobes.
- This further invention is considered an improvement on the invention noted above and is based on, not discarding small eigenvalues (which is the approach noted above) but instead to multiply the reciprocal of the eigenvalue in the inverse matrix by a decreasing number until the sidelobes are removed or so small that they have little effect.
- S/N OdB
- Figure 2 illustrates CIR estimated using the algorithm disclosed in US
- Figure 3 illustrates CIR estimated using SVD and including 10 largest eigenvalues
- Figure 4 illustrates CIR estimated using SVD and including 16 largest eigenvalues
- Figure 5 illustrates CIR estimated using SVD and including 25 largest eigenvalues
- Figure 6 illustrates CIR estimated using SVD and including 36 largest eigenvalues
- n 2 25 S/N - OdB
- Figure 10 illustrates CIR estimated using SVD and using ramping function with n-
- One method of estimating the tap weights is to determine the weights which predict a received signal most closely matching the measured signal.
- the desired FIR tap weights (ot ⁇ ) are estimated by requiring that the cumulative square error between the measured received signal and the estimate of the received signal (i.e. the known transmitted signal passed through the FIR approximating the channel) is minimised, i.e. minimise
- T is the MxN matrix of conjugated delayed signals t jj -T j j O ⁇ i ⁇ M O ⁇ j ⁇ N (2 , 6)
- US 5,473,632 discloses an algorithm which stabilises the inversion of the matrix and produces results reasonably resilient to noise.
- US 5,473,632 estimates the CIR by c - ((A + e!)-rT)R (2.9) where I is the identity matrix and ⁇ is a constant (they call a noise term). This equation equates to the equation disclosed in US 5,473,632 at col. 6, line 44. The results of using this algorithm follow:
- the matrix A (which is positive definite) may be represented by where U and V are orthogonal and S is diagonal. In fact the diagonal elements of S are the eigenvalues (all positive) of A and these are arranged in decreasing order: s 1>1 > s 2 , 2 .>s n , ⁇ (2.1 1 ) If the matrix is singular, some of these are zero. The inverse of A is then
- A-1 VS-1U' (2.12) where S- 1 is formed by inverting the (diagonal) elements of S. The effects of a poorly conditioned matrix are clear here, the (almost zero) eigenvalues are inverted and produce an extremely large contribution to the inverse. These amplify the noise.
- Figures 16 and 17 illustrate the results for some more complex CIR's. It has also been found that the present embodiment has the ability to enable the algorithm to discriminate between peak, and Figure 18 shows close peaks split with time separations greater than about 5 ⁇ s.
- Matrix inversion as described above is also required when estimating the channel impulse response vector h.
- ⁇ h ⁇ (2.14)
- ⁇ is the crosscorrelation matrix between locally generated signals
- ⁇ is the crosscorrelation matrix between locally generated signals and the received signal.
- the principles and method as described above are equally applied in this situation.
- CDMA is applied usually using optimum combining of signals received from multiple sensors (antenna array). Optimum combining maximises the signal-to-interference-plus-noise ratio (SINR).
- w the weight vector for antenna elements
- ⁇ a constant
- R the received interference-plus-noise correlation matrix
- u d * is the conjugate of the desired signal vector.
- the preferred implementation of the inventions disclosed is to generate a CIR using the matrix multiplication as specified in equation (2.8) and to then use SVD to determine the inverse of the matrix A.
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Noise Elimination (AREA)
- Mobile Radio Communication Systems (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
AU62175/98A AU6217598A (en) | 1997-02-25 | 1998-02-23 | Channel impulse response estimation using singular value decomposition |
JP53734098A JP2001513954A (ja) | 1997-02-25 | 1998-02-23 | 特異値分解を用いたチャンネルインパルス応答の推定 |
EP98904204A EP0963642A1 (fr) | 1997-02-25 | 1998-02-23 | Estimation de la reponse impulsionnelle d'un canal au moyen de la decomposition d'une valeur singuliere |
NO994072A NO994072L (no) | 1997-02-25 | 1999-08-24 | Estimering av kanalimpulsrespons ved Õ benytte dekomposisjon av singulær verdi |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
AU14889/97A AU688228B1 (en) | 1997-02-25 | 1997-02-25 | A system and method of estimating CIR |
AU14889/97 | 1997-02-25 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO1998038772A1 true WO1998038772A1 (fr) | 1998-09-03 |
Family
ID=3704956
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/FI1998/000159 WO1998038772A1 (fr) | 1997-02-25 | 1998-02-23 | Estimation de la reponse impulsionnelle d'un canal au moyen de la decomposition d'une valeur singuliere |
Country Status (6)
Country | Link |
---|---|
EP (1) | EP0963642A1 (fr) |
JP (1) | JP2001513954A (fr) |
CN (1) | CN1248370A (fr) |
AU (1) | AU688228B1 (fr) |
NO (1) | NO994072L (fr) |
WO (1) | WO1998038772A1 (fr) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001050694A1 (fr) * | 1999-12-30 | 2001-07-12 | Telit Mobile Terminals S.P.A. | Procede et systeme recepteur wcdma-utra/fdd |
GB2365715A (en) * | 2000-05-15 | 2002-02-20 | Univ Bristol | A circuit for a communications receiver |
EP1185002A2 (fr) * | 2000-09-02 | 2002-03-06 | Lg Electronics Inc. | Procédé de traitement de signaux dans des systèmes de communication avec une pluralité d'antennes |
GB2370469A (en) * | 2000-09-07 | 2002-06-26 | Nec Corp | Improvements in CDMA receivers |
WO2003036891A2 (fr) * | 2001-10-24 | 2003-05-01 | Zenith Electronics Corporation | Procede d'estimation de canal combinant une correlation et la methode des moindres carres |
KR100413096B1 (ko) * | 2000-04-06 | 2003-12-31 | 가부시키가이샤 엔.티.티.도코모 | Cdma 시스템에서의 통신 품질 측정 장치 및 방법 |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040081131A1 (en) * | 2002-10-25 | 2004-04-29 | Walton Jay Rod | OFDM communication system with multiple OFDM symbol sizes |
CN1300962C (zh) * | 2002-12-31 | 2007-02-14 | 上海贝尔阿尔卡特股份有限公司 | 正交频分复用系统中均衡快衰落信道的方法及装置 |
CN100539557C (zh) * | 2004-10-18 | 2009-09-09 | 威盛电子股份有限公司 | 全球移动通信系统预测信道脉冲响应的装置及方法 |
MX2007016577A (es) * | 2005-07-13 | 2008-03-07 | Exxonmobil Upstream Res Co | Metodo para predecir lo mejor y lo peor en un conjunto de soluciones no unicas. |
CN104950326B (zh) * | 2015-03-04 | 2016-09-28 | 中石化石油工程技术服务有限公司 | 基于目的层频谱的可控震源非线性扫描信号的设计方法 |
CN114018250B (zh) * | 2021-10-18 | 2024-05-03 | 杭州鸿泉物联网技术股份有限公司 | 惯性导航方法、电子设备、存储介质和计算机程序产品 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4945502A (en) * | 1988-12-27 | 1990-07-31 | Eastman Kodak Company | Digital image sharpening method using SVD block transform |
US5209237A (en) * | 1990-04-12 | 1993-05-11 | Felix Rosenthal | Method and apparatus for detecting a signal from a noisy environment and fetal heartbeat obtaining method |
US5473632A (en) * | 1991-10-31 | 1995-12-05 | Rohde & Schwarz Gmbh & Co. | Method of determining the complex pulse response of a radio system |
-
1997
- 1997-02-25 AU AU14889/97A patent/AU688228B1/en not_active Ceased
-
1998
- 1998-02-23 CN CN 98802802 patent/CN1248370A/zh active Pending
- 1998-02-23 EP EP98904204A patent/EP0963642A1/fr not_active Withdrawn
- 1998-02-23 JP JP53734098A patent/JP2001513954A/ja active Pending
- 1998-02-23 WO PCT/FI1998/000159 patent/WO1998038772A1/fr not_active Application Discontinuation
-
1999
- 1999-08-24 NO NO994072A patent/NO994072L/no not_active Application Discontinuation
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4945502A (en) * | 1988-12-27 | 1990-07-31 | Eastman Kodak Company | Digital image sharpening method using SVD block transform |
US5209237A (en) * | 1990-04-12 | 1993-05-11 | Felix Rosenthal | Method and apparatus for detecting a signal from a noisy environment and fetal heartbeat obtaining method |
US5473632A (en) * | 1991-10-31 | 1995-12-05 | Rohde & Schwarz Gmbh & Co. | Method of determining the complex pulse response of a radio system |
Non-Patent Citations (6)
Title |
---|
BARTON M ET AL: "Reduced-rank least squares channel estimation", IEEE TRANSACTIONS ON ACOUSTICS, SPEECH AND SIGNAL PROCESSING, vol. 38, no. 8, August 1990 (1990-08-01), NEW YORK, US, pages 1403 - 1410, XP002063707 * |
BARTON M: "Approximate sample variance analysis of exponential channel derived from SVD-based algorithm", IEE PROCEEDINGS I (COMMUNICATIONS, SPEECH AND VISION), vol. 136, no. 6, December 1989 (1989-12-01), UK, pages 381 - 384, XP002063708 * |
EDFORS O ET AL: "OFDM CHANNEL ESTIMATION BY SINGULAR VALUE DECOMPOSITION", 1996 IEEE 46TH. VEHICULAR TECHNOLOGY CONFERENCE, MOBILE TECHNOLOGY FOR THE HUMAN RACE ATLANTA, APR. 28 - MAY 1, 1996, vol. VOL. 2, no. CONF. 46, 28 April 1996 (1996-04-28), NEW YORK, US, pages 923 - 927, XP000593108 * |
RAHMAN J ET AL: "Deconvolution and total least squares in finding the impulse response of an electromagnetic system from measured data", IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, APRIL 1995, USA, vol. 43, no. 4, NEW YORK, US, pages 416 - 421, XP002063705 * |
VAN DER VEEN A L ET AL: "Singular value analysis of space-time equalization in the GSM mobile system", 1996 IEEE INTERNATIONAL CONFERENCE ON ACOUSTICS, SPEECH, AND SIGNAL PROCESSING CONFERENCE PROCEEDINGS (CAT. NO.96CH35903), 7 May 1996 (1996-05-07) - 10 May 1996 (1996-05-10), NEW YORK, USA, pages 1073 - 1076 vol. 2, XP002063706 * |
ZWILLINGER D.: "Standard Mathematical Tables and Formulae, 30th edition", 1996, CRC PRESS, BOCA RATON, US, XP002063709 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001050694A1 (fr) * | 1999-12-30 | 2001-07-12 | Telit Mobile Terminals S.P.A. | Procede et systeme recepteur wcdma-utra/fdd |
KR100413096B1 (ko) * | 2000-04-06 | 2003-12-31 | 가부시키가이샤 엔.티.티.도코모 | Cdma 시스템에서의 통신 품질 측정 장치 및 방법 |
GB2365715A (en) * | 2000-05-15 | 2002-02-20 | Univ Bristol | A circuit for a communications receiver |
GB2365715B (en) * | 2000-05-15 | 2004-04-28 | Univ Bristol | Circuit |
EP1185002A2 (fr) * | 2000-09-02 | 2002-03-06 | Lg Electronics Inc. | Procédé de traitement de signaux dans des systèmes de communication avec une pluralité d'antennes |
EP1185002A3 (fr) * | 2000-09-02 | 2002-08-07 | Lg Electronics Inc. | Procédé de traitement de signaux dans des systèmes de communication avec une pluralité d'antennes |
US6876693B2 (en) | 2000-09-02 | 2005-04-05 | Lg Electronics Inc. | Method for processing signal in communications system having plurality antennas |
GB2370469A (en) * | 2000-09-07 | 2002-06-26 | Nec Corp | Improvements in CDMA receivers |
GB2370469B (en) * | 2000-09-07 | 2004-10-27 | Nec Corp | Improvements in CDMA receivers |
WO2003036891A2 (fr) * | 2001-10-24 | 2003-05-01 | Zenith Electronics Corporation | Procede d'estimation de canal combinant une correlation et la methode des moindres carres |
WO2003036891A3 (fr) * | 2001-10-24 | 2003-07-31 | Zenith Electronics Corp | Procede d'estimation de canal combinant une correlation et la methode des moindres carres |
US7035353B2 (en) | 2001-10-24 | 2006-04-25 | Zenith Electronics Corporation | Channel estimation method blending correlation and least-squares based approaches |
Also Published As
Publication number | Publication date |
---|---|
NO994072D0 (no) | 1999-08-24 |
NO994072L (no) | 1999-08-24 |
AU688228B1 (en) | 1998-03-05 |
CN1248370A (zh) | 2000-03-22 |
EP0963642A1 (fr) | 1999-12-15 |
JP2001513954A (ja) | 2001-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6765969B1 (en) | Method and device for multi-user channel estimation | |
KR100757696B1 (ko) | 개선된 신호 포착과 프로세싱을 지닌 코드 분할 다중 접속 시스템 및 오퍼레이션 방법 | |
CN102342075A (zh) | 从一个频率和时间到另一个的信道推断 | |
Dlugos et al. | Acquisition of spread spectrum signals by an adaptive array | |
CN101779426A (zh) | 用于使用近似逆矩阵来估计信号参数的信号处理器 | |
WO2002060138A2 (fr) | Procede et dispositif d'estimation de canaux de rang reduit dans un systeme de communication | |
EP0963642A1 (fr) | Estimation de la reponse impulsionnelle d'un canal au moyen de la decomposition d'une valeur singuliere | |
CN109450489B (zh) | 一种扩频水声通信的导频序列干扰抵消方法 | |
CN1883131B (zh) | 用于确定信道估计的设备与方法 | |
EP2009825A1 (fr) | Dispositif de reception radio, dispositif de transmission radio, station de base radio, procede de reception et procede de transmission | |
JP2008236065A (ja) | 受信装置、無線送受信システム及び無線受信方法 | |
US20140269880A1 (en) | Communication system and method using subspace interference cancellation | |
US8259854B2 (en) | Channel estimation using common and dedicated pilots | |
CN100553243C (zh) | 波束空间正交频分复用调制系统自适应波束形成方法 | |
CN1303545A (zh) | 使用网络用户标识符的具有初始化的迭代投影 | |
EP1716683B1 (fr) | Reception de communications a compensation de frequence | |
JP4448847B2 (ja) | 複雑さを低減させたスライディングウィンドウ方式による等化器 | |
JP4133587B2 (ja) | 受信装置 | |
US6904107B1 (en) | Space-weighted communication path estimation | |
EP2250776A2 (fr) | Procédé et appareil d'estimation de canal pour un traitement de signaux de communication | |
Yu et al. | Adaptive array beamforming for cyclostationary signals | |
JP4350020B2 (ja) | 受信装置 | |
WO2000077993A1 (fr) | Manipulation de matrice pour l'estimation de canaux | |
Zidane et al. | Bit error rate analysis of mc-cdma systems with channel identification using higher order cumulants | |
Elkamchouchi et al. | Robust doa estimation and data detection for mobile communication in multipath environment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 98802802.6 Country of ref document: CN |
|
AK | Designated states |
Kind code of ref document: A1 Designated state(s): AL AM AT AU AZ BA BB BG BR BY CA CH CN CU CZ DE DK EE ES FI GB GE GH GM GW HU ID IL IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MD MG MK MN MW MX NO NZ PL PT RO RU SD SE SG SI SK SL TJ TM TR TT UA UG US UZ VN YU ZW |
|
AL | Designated countries for regional patents |
Kind code of ref document: A1 Designated state(s): GH GM KE LS MW SD SZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH DE DK ES FI FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA GN ML MR NE SN TD TG |
|
DFPE | Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101) | ||
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
WWE | Wipo information: entry into national phase |
Ref document number: 09371602 Country of ref document: US |
|
ENP | Entry into the national phase |
Ref document number: 1998 537340 Country of ref document: JP Kind code of ref document: A |
|
WWE | Wipo information: entry into national phase |
Ref document number: 1998904204 Country of ref document: EP |
|
WWP | Wipo information: published in national office |
Ref document number: 1998904204 Country of ref document: EP |
|
REG | Reference to national code |
Ref country code: DE Ref legal event code: 8642 |
|
WWW | Wipo information: withdrawn in national office |
Ref document number: 1998904204 Country of ref document: EP |