TW200814563A - Method for covariance matrix update - Google Patents

Method for covariance matrix update Download PDF

Info

Publication number
TW200814563A
TW200814563A TW096125879A TW96125879A TW200814563A TW 200814563 A TW200814563 A TW 200814563A TW 096125879 A TW096125879 A TW 096125879A TW 96125879 A TW96125879 A TW 96125879A TW 200814563 A TW200814563 A TW 200814563A
Authority
TW
Taiwan
Prior art keywords
matrix
square root
sample
sample set
covariance
Prior art date
Application number
TW096125879A
Other languages
English (en)
Inventor
Ali S Khayrallah
Original Assignee
Ericsson Telefon Ab L M
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 Ericsson Telefon Ab L M filed Critical Ericsson Telefon Ab L M
Publication of TW200814563A publication Critical patent/TW200814563A/zh

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B1/00Details of transmission systems, not covered by a single one of groups H04B3/00 - H04B13/00; Details of transmission systems not characterised by the medium used for transmission
    • H04B1/69Spread spectrum techniques
    • H04B1/707Spread spectrum techniques using direct sequence modulation
    • H04B1/7097Interference-related aspects
    • H04B1/711Interference-related aspects the interference being multi-path interference
    • H04B1/7115Constructive combining of multi-path signals, i.e. RAKE receivers
    • H04B1/712Weighting of fingers for combining, e.g. amplitude control or phase rotation using an inner loop
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B2201/00Indexing scheme relating to details of transmission systems not covered by a single group of H04B3/00 - H04B13/00
    • H04B2201/69Orthogonal indexing scheme relating to spread spectrum techniques in general
    • H04B2201/707Orthogonal indexing scheme relating to spread spectrum techniques in general relating to direct sequence modulation
    • H04B2201/7097Direct sequence modulation interference
    • H04B2201/709727GRAKE type RAKE receivers

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Noise Elimination (AREA)

Description

200814563 九、發明說明: 【發明所屬之技術領域】 本發明一般係關於多信號接收器, V又付疋$之,係關 於更新用於多信號接收器之處理參數。 【先前技術】 傳統的無線接收器常常使用一協方差 /斗、次』丨μ 千^式的雜訊及 /或貝料協方差資訊以抑制多個所接收信號(即一單一的戶 發射信號之多個樣本集、多個不同信號或其任何組合)之 間的干擾。干擾抑制接收器之範例包括碼片等化界之
Rake接收器 ' —般化RAKE(GRAKE)接收器、單一輸二多 個輸出接收器、多個輸入多個輸出接收器等。 夕 此項技術中熟知的係,干擾抑制接收器需要對該協方差 矩陣中的協方差進行精確追蹤。由於所接收的樣丄集數目 眾多,因此追蹤該等協方差常常需要極複雜的計算。此等 複雜的計算常常限制一無線接收器精確追蹤並使用信號協 方差之能力。例如’-無線接收器可以實施數目減少的協 方差更新來減小接收器複雜性。但是,此減小可能限制用 於處理所接收信號的協方差之精確度。因此,無線產業不 斷尋找改良干擾抑制接收器的操作之方法。 【發明内容】 本發明提供一種用以在向具有兩個或更多樣本集之一所 接收信號添加或從該所接收信號移除一樣本集時更新一協 方差矩陣之一平方根矩陣之方法。一範例性方法包括決定 針對一所添加樣本集之減損關聯及一雜訊功率。該等減損 122662.doc 200814563 關聯表示所添加樣本集中的值與已存在樣本集中的值 的父叉關聯。該方法進一步包括依據該等減損關聯及該雜 说功率來補償該平方根矩陣。接著藉由添加具有一從該雜 訊功率導出的對角元素與從該等減損關聯導出的非對角元 素之仃與-列來增大已補償的平方根矩陣。還說明 應設備。 7 另耗例性方法包括決定針對一所移除樣本集之減損關 聯及-雜訊功率。該等減損關聯表示所移除樣本集中的值 與其餘樣本集中的值之間的交又關聯。該方法進一步包括 從一三角矩陣及一對角矩陣移除對應於所移除樣本集之一 ^與一行。該等三角及對角矩陣對應於該協方差矩陣。接 者依據減小的三角矩陣、減小的對角矩陣、減損關聯及雜 訊功率來計算一減小的平方根矩陣。還說明一對應設備。 【實施方式】 本發明提供一種用以在向一般化rake(GRAKe)接收器 添加或從該一般化rAKE(grake)接收器移除樣本集(即, 指峰(finger))時因所添加或所移除的樣本集而更新協方差 矩陣之方法及設備。圖丨解說適用於本發明之一範例性無 線接收器1〇〇之一方塊圖。無線接收器1〇〇包括一 grake 接收杰lio、一接收裔前端112及一天線114。接收器前端 112向GRAKE接收器11〇輸出所接收信號广⑺。該等所接收 信號r⑺包含從藉由一或多個天線丨14(其與一或多個接收 器前端112相關聯)接收的無線信號獲得之樣本值串流。一 範例性前端112可以包括放大器、濾波器、混合器、數位 122662.doc 200814563 在產生適合藉由該GRAKE接收 號時需要此等元件。在某些具 包括從不同天線(未顯示)發射之成 每一所接收信號咐-般包括自—從不同方向並以不同 時間延遲到達該接收器剛的信號產生之—或多個樣本 集。GRAKE接收器11 〇將盘所Γ啼 肝/、所接收佗號r⑺相關聯的樣本集
組合以產生針對第所符號週期之一符號估計攻岣。更特定 言之,GRAKE接收器、110對所接收信號吻進行延遲、解 展、加權及組合以估計z(w),下面將進一步說明。
化器及/或其他電子元件 器110來處理的已取樣信 體貫施例中,還可以 分0 GRAKE接收器11〇包含複數個RAKE指峰12〇、組合哭 130、通道估計器mo、處理器15〇及路徑搜尋器丨⑼。 RAKE指峰120藉由使用包括於每一 RAKE指峰12〇中的一延 遲元件122及一相關器124來處理雜訊及資料信號。延遲元 件122將所接收信號厂⑺延遲一由處理器15〇定義之數量。 相關器124將所延遲的信號與一擴展碼關聯以產生一解展 值y。向組合态130輸入來自相關器124之解展值y =(力 組合器130包括加權元件132與加法器134。加權元件132 接收來自RAKE指峰120之解展信號y並藉由處理器1 〇所決 定的加權係數…,wk)來對其加權,下面將進一 步說明。加法器134在每一符號週期期間將加權的解展信 號加總以形成一符號估計z(m)。 處理器150與通道估計器140及路徑搜尋器16〇協作來決 122662.doc 200814563 定延遲δ = (& A,…,心)及對應的加權係數w。熟習此項技 術者習知,路徑搜尋器160識別針對RAKE指峰120之候選 延遲。通道估計器140估計該通道並向處理器15〇提供估計 的通道係數h==(/z7, …,/^),在此亦稱為通道估計h。依 據該等通道估計h及所接收信號r⑺,處理器150將適當的 延遲δ及加權係數w分別指派給該等raKE指峰120及組人 器 130。 在一傳統的GRAKE中,依據橫跨該等尺八!^指峰12〇之 雜δίΐ關聯末计异加權係數w。該等加權係數係給定為· w = R'1h? (1) 其中R係表示RAKE指峰120之間的交叉關聯之一協方 差矩陣R之倒數。此項技術中習知,可以依據解展先導符 號來計算協方差矩陣R。接著將符號估計z(w)計算如下·· 2(m) = wHy? (2) 其中炉表示加權係數向量之厄米共軛(Hermhian)。可將 訊號對雜訊比估計為: SNR = wHh. ( 3 ) 為確保GRAKE接收器110精確地處理所接收信號^⑺,處 理器150週期性更新該協方差矩陣此週期性更新可以係 依據該等通道估計h所發生的變化、RAKE指峰12〇之數目 等。等式(4)顯不用以依據一當前協方差矩陣r、一指數衰 減因數《及一錯誤向量e來產生一更新的協方差矩陣r,之一 傳統方法。 (4) 122662.doc 200814563 錯誤向量e表示預期符號估計與當前接收器輸出符號估 計攻之間的差。該指數衰減因數具有一介於〇與1之間的 值。 每次在處理器150更新協方差矩陣尺時,便計算出新的組 合權重w及一新的。如等式(1)所示,產生新的組合權 重w百先需要計算更新的協方差矩陣w之倒數,其固有的 係複倒數。此外,儘管藉由頻繁更新方差矩陣而允許 GRAKE接收11 〇更精確地追蹤通道動態狀況,增加更新 頻率亦增加該接收器丨丨〇之複雜性。 a本發明藉由使用該協方差矩陣尺之一平方根“來計算該 等加權係數w以減少該GRAKE接收器11〇之複雜性。該協 方差矩陣Κ係一正定矩陣。因此,可以使用Cholesky分解 來計异滿足以下等式之矩陣Μ : 可以將該矩陣Μ解釋為該協方差矩陣尺之「平方根」。一 般也Μ係了二角矩陣,其正實數元素位於對角線上。 該協方差矩陣R亦可以係表示為相關矩陣之一函 數,如下式所示: R = LDLH,(6) 其中L係-下三角矩陣,其中】,在該對角線上㈣係具有 正兀素對角料1由將等式(5)與等式⑹相比較, 可以定義之間的關係。等式⑺顯示結果。 M = LD1/2 (7) 專式(5)、(6)及(7)顯示協方差矩料與—平方根矩陣m、 122662.doc 200814563 一二角矩陣L及一對角矩陣D之關係如何。 由於Μ係一下三角矩陣,因此藉由使用%而非該協方差 矩陣R來計算該等加權係數W不太複雜。可以在兩個步驟 中計算該等加權係數。首先,依據Μ及h來計算一中間向 量V。該中間向量v係給定為:
Mv = h · ( 8 ) 在獲得該_間向量¥後,依據下式計算該等加權係數·· (、 MHw = v. (9) 所接收#號之SNR係給定為: SNR = vHy. (10) 組合的RAKE輸出係給定為: 2 (m) = wHy = vHIW'1y . (11) 由於針對一三角矩陣(即Μ)之矩陣操控不太複雜,因此 使用Μ來替代該協方差矩陣r使得傳統上複雜的graKE操 作簡化。 (J 處理器1 50還可以藉由使用該平方根矩陣Μ來替代該協 方差矩陣R來簡化更新操作。如等式(4)所示,可以與r、α 及e成函數關係來計算一更新的協方差矩陣R,。藉由使用 等式(5)至(7)所示R、Μ、L及D之間的關係,可以開發出用 • 以依據L、D及e的當前版本來更新該平方根矩陣Μ之規 則’如以下等式(12)至(16)所示。 首先,假定可以依據下式計算針對三角矩陣L之更新 L,。 L=LK? (12) 122662.doc 10- 200814563 其中κ係具有下式所定 義結才冓 之一三角矩陣·· Κ: Ρ!β' 1 Ρ3^1 P^2 1 (13) -Ρηβ、Ρηβ2 …ρ β’ 向量心,心...,〜)中的參數_〗義為: P = L-V (14) 了以依據以下程库爽斗曾 ί' Ο ^ 汁开用於等式(13)之參數凡及所更 新的對角矩陣D,之對角元素々。 =1/α for/ = 1 :/7-1 r /,1=9/0(/ of;' = a / γ. (1 5 ) β丨=Pi I a end dn' = dn+\pn\2/yn 應明白,此程序提供用以計算!),的元素之一範例性方 法。其他替代性方法亦適用。一旦計算出一更新的三角矩 陣L及一更新的對角矩陣d’,便可以使用以下等式(16)來 更新該平方根矩陣。 M,:LO,1/2=LKD丨1/2 (16) 用以更新與等式(12)至(16)相關聯的平方根矩陣之程序 可以係縮寫為(ir,L,,D,) = /?t;(L,D,e,a),其中表示根更新。等 式(12)至(16)之i?y程序表示在以下說明中重複用作一用於 本發明的構建組塊之一範例性及"程序。但是,應明白更 新該平方根矩陣之其他及ί/程序亦可以應用於本發明。 122662.doc -11 - 200814563 上面說明用以更新針對一恆定數目的RAKE指峰120(其 對應於一怪定數目的樣本集)之平方根矩陣Μ之一程序。但 是,本發明還允許在rAKE指峰12〇之數目改變時更新該平 方根矩陣Μ。藉由改變raKE指峰120之數目而改變該協方 差矩陣R。給定一協方差矩陣R,則可以在添加或新配置 一 RAK指修120時計算一增大的更新之協方差矩陣R,。可 以將增大的協方差矩陣R,寫為··
-2 Hi (17) R =: σ Ρ
Ur」, 其中σ2表示與該樣本集相關聯之雜訊功率,而p表示針 …子處於母舊延遲的樣本集,處於新延遲的樣本集之雜訊 協方差之一向量。依據上述L、D、M&R之間的關係,可 以依據如下等式1 8來直接計算一增大的更新之平方根矩陣 Mf : 在等式(18)中,Μ”表示依據等式(12)至(16)計算出之一 補仏的平方根矩陣,其中錯誤向量e係替換為雜訊協方 j向量P,而其中指數衰減因數α係替換為卜用以補 仏與等式(12)至(16)相關聯的平方根矩陣之程序可以係簡 UVi'L'D’,尋 5D,pM/(j2) (19) 當移除或解配置-RAKE指峰12〇時,可以將當前協方差 ㈣R減少-列與_行,如下面等式(2g)中所*,以產生 一減小的更新之協方差矩陣R,。 122662.doc -12- 200814563 可以依據三角矩陣L及對角矩陣D之減小的版本來直接 計算減小的更新之平方根矩陣M,。可以藉由從該三角矩陣 L移除對應於所移除的RAKE指峰12〇之列與行來產生一減 小的三角矩陣L”。同樣,可以藉由從該對角矩陣d移除對 應於所移除的RAKE指峰120之列與行來產生一減小的對角 矩陣D,,。接著可以依據等式(12)至⑽來計算減小的更新 之平方根矩陣M,,#中L係替換為減小的三角矩陣l",d 係替換為減小的對角矩陣D”,以口替換錯誤向量e,而苴中 以(1〜2)來#換餘衰減因數α。貞增大的平方根矩料算 相同,ρ與σ2分別對應於一雜訊協方差向量及雜訊功率(對 應於移除的RAKE指峰120)。用以計算與等式(12)至(16)相 關聯之減小的更新之平方根矩陣M,之程序可以係簡寫為: (M|,L,,D,),(L”,D”,P,-l/a2). (21) 總而言之,處理器150可以依據[、D、〇1及6之當前版本 來計算一特殊的三角矩陣K及一更新的平方根矩陣河,、三 角矩陣L,及對角矩陣D,。另外,處理器15〇可以藉由替代 協方差矩陣R而使用平方根矩陣M來適應變化的樣本集數 目。M、L、D及K之固有的算術特性使得處理器bo能夠 藉由使用較一般與更新R相關聯的程序更加簡單的程序來 更新Μ。處理器150可以接著替代M而使用M,來處理所接收 信號。因此,藉由替代R而使用M,從而明顯降低藉由處 理器1 50實施的更新與處理操作之複雜性。 圖2 —般解說用以依據一變化數目的樣本集來更新該平 122662.doc •13- 200814563 方根矩陣Μ之一範例性程序200。如上所述,GRAKE接收 ^§110藉由使用Μ來處理所接收信號(步驟2丨〇)。當具有 一新延遲之一新樣本集出現在所接收信號中(步驟22〇)時, 處理器150增大Μ(步驟230)以應對新的RAKE指峰12〇。 GRAKE接收器110接著將先前的平方根矩陣M替換為增大 的更新之平方根矩陣M,,以便將來對所接收信號〃⑺進行
Ο 處理(步驟21〇)。當一樣本集不再存在於所接收信號中(步 驟240)時’處理器150減小M(步驟25〇)以應對移除的ram 指峰mjRAKE接收器110接著將先前的平方根矩陣崎 換為減小的更新之平方根矩陣以,,以便將來對所接收信號 進行處理(步驟21〇)。 圖3針對圖2解說一範例性增大程序(步驟㈣之特定細 ^為增大該平方根矩陣M,處理器⑼決定針對所添加 ,集之減損關聯P及一雜訊功率σ2(步驟232)。該等減損 關聯Ρ表示所添加樣本集中的值與現存樣本集中的值之間 的父又關聯。依據所決定的減損關聯ρ及雜 理器150補償當前平方 處 根矩陣Μ” (步驟234)。處理。。丫 產生一已翁的平方 已 益150接著藉由依據等式(19)向 已補 <員的平方根矩陳Μ ”a y 方麻… 行與—列來增大已補償的平 =矩陣M’,’以產生增大的更新之平方根矩陣‘ 元素所添加的订與列具有一從該雜訊功率σ2導出之對角 =、:該等減損關•導出之非對角的元素。 、十對圖2顯示一範例彳生% 節。為減小該平方根二;:序(步驟25〇)之特定細 早Μ處理器150決定針對所移除 122662.doc •14- 200814563 樣本集之減損關聯ρΑ_雜訊功率―(步驟252)。該等減^ 關聯p表示所移除樣本集中的值與其餘樣本集中心值之二 的交又關聯。處理^5〇接著從該三角矩陣l及該對角矩: D移除對應於所移除樣本集之列與行,以分別形成—減小 的三角矩陣L,及一減小的對角矩陣D"(步驟254)。依據減小 的三角矩陣L"、減小的對角 2 ^ 減知關聯P及雜訊功 率σ ,處理器15〇依據箄式 %儺寺式(21)汁异減小的平方根矩陣μ,
Ο (步驟256)。 上面說明依據變Α的通道條件及/或可用於該接收哭之 :變化數目的樣本集或指峰,針對一平方根矩陣m之更新 私序應明白’每次更新該平方根矩卩車Μ時,更新的版本 會替換先前的版本而用於將來的處理及更新操作。 、針對平方根矩陣Μ之上述更新程序假定該協方差矩陣r 或該平方根矩陣Μ已經存在及/或已初始化。但是,此假定 不-定始終正確。了面說明用以初始化該平方根矩陣社 多個程序。 在—具體實施例中’藉由使用—迭代程序來初始化Μ, 該程序在時^㈣―設定為恆等矩陣之—倍數的平方根 矩陣,並在時間卜m初始化的平方根矩陣财 束依據此十亙等開始」具體實施例,首先將該協方差矩 陣R設定為該恆等矩陣之某一倍數£,如等式(22)所示。 R(i = 〇) = d (22) 因此,該初始協方差矩陣並非奇異。依據等式(22)及等 式(5)至(7) ’在牯間t=〇的平方根矩陣M、三角矩陣[及對角 122662.doc -15- 200814563 矩陣D係給定為: M(f = 0) = ^1/2| L(f = o) = i . (23) D(f = 〇) = £*! 為在-隨後時間,=r決定初始平方根矩陣m,處理器 1 5 0實施以下迭代步驟: for ί = 1:Γ ^ [Μ(〇>L(f),D(f)] = (L(i -1),D(f-1)?e(f)}αJ (24) 其中妨(·)表不#式(12)至(16)所#更新程序。處理器 150在上述後續處理及更新操作中使用所產±的初始化平 方根矩陣Μ。 替代的初始化具體實施例係依據該協方差矩陣r首先係 設定為一 零矩陣(R(i = 〇H〇D之假定。由於該零協方差矩 陣R係尺寸《,因此要經過„次更新才產生一全秩協方差矩 陣R。即,為產生秩次„之一全秩協方差矩陣r,處理器 150依據等式⑷更新協方差矩陣„次。下面說明用以在藉由 使用此「零開始」程序來初始化R時將該平方根矩陣Μ初 始化之兩個替代具體實施例。 在一具體實施例中,處理器15〇將藉由第”次更新產生的 協方差矩陣定義為初始協方差矩陣R(i = w),本文中稱為混 合程序。隨後,處理器15〇藉由使用全秩的初始化協方差 矩陣R(~«)來依據等式(5)至(7)產生該初始平方根矩陣M。 處理器150在上述後續處理及更新操作中使用所產生的初 始化平方根矩陣]V[。依據該混合程序,處理器丨5 〇首先藉 由使用上述零開始程序來決定一秩次^之協方差矩陣R。隨 122662.doc -16- 200814563 後,處理器150使用Cholesky分解來決定對應的秩次π之 L、D及]V[矩陣。處理器150接著從時間? = η + !至某一所 需時間ί = Γ使用秩次π之L、D及Μ矩陣來實施與該恒等開 始程序相關聯之平方根更新程序,如下式所示: for f = π +1: Τ [M(〇^(〇>D(f)] = R(7(L(i-l),D(f.l),e(i),a) (25) end 處理器150在上述後續處理及更新操作中使用所產生的 已初始化平方根矩陣]VI。 Ο 在另一混合具體實施例中,處理器150針對首度”次迭代 而定義一新的平方根矩陣M與一新的更新之平方根矩陣 ’其中考量每一迭代中該協方差矩陣R之特定結構,在 此稱為經修改的混合程序。例如,對於第一迭代,該協方 差矩陣R係藉由R(1) = ae(1)e(i)H定義之一秩次i之矩陣。對於此 協方差矩陣R(1),可以將第一迭代之三角矩陣l(i)、對角 矩陣D(l)及對應的平方根矩陣Μ〇)定義為: Ο L(1)= e2/ei 1 (26) len'e, enA^-1 1 ♦I2 D(1)=. (27) 0— ei M(l) = L(l)D(l)1/2^i 62hi ; .βη 0 • · * · 0 ··· 0 122662.doc -17- (28) 200814563 從以上專式可以看出,L( 1)係一秩次π之矩陣,但是 M( 1)係一秩次1之矩陣。此使得處理器15 〇能夠在產生該平 方根矩陣M(2)的下一迭代之更新程序中使用。一般而 言,由於在該對角線中的2,因此每一三角矩陣L⑺將係一 秩次η之矩陣。此確保L之可倒性以適應計算該平方根更新 之需要,如等式(14)所示。另外,每一對角矩陣〇⑺在最 後π至/個對角元素中將具有〇,從而將秩次?賦予其及該平
方根矩陣Μ⑴。下面提供針對此範例性初始化程序產生針 對Μ的每一更新之特定細節。 對於第ί次迭代,處理器15〇依據下式來計算更新的對角 矩陣D'之最先的/個對角元素 7ι -Μα for / = 1:f a,+N2 ⑽ 丨 β; - Ρι / a end Ο 如上所述,其餘對角元素係構造為零。 處理斋1 5G使用該等結果來構造—新的特殊對角矩陣κ 藉由將在第'項目卩夕卜的右下側帛目設定為0而相對於等 (13)中所定義者來修改κ,如等式(3〇)所示。 式 Ριβ\ 1 Ρζβ\ p3y8; 1 κ = : : :
Pt+'β、ρΜβ2 ... ρ“'β; ' ^/+2Α ρΜβι · · · pt+2fi*t ο 1 A A pjI ... Ρηβ* 〇 ... 〇 122662.doc -18- (30) 200814563 最後,處理器150依據等式(12)計算更新的三角矩陣L,而 依據等式(16)計算更新的平方根矩陣%,。總而言之,處理 器150依據下式在時間,=Γ產生初始化的平方根矩陣μ ·· L(1)自等式25 D(1)自等式26 MC〇自等式27 假設f = 2: π [M(i),L(i),D(i) = Rt/S(L(i - 1)’D(i - 1),e(i),a,i)] (3 1 ) 結束
假設 f = r? +1: T
[M(0,L ⑴,D(i) = Ra(L(i 一 1),D(M),e(i),a)] 結束, 其中表示等式(12)至(16)所示之更新程序,而 及仍(·)表不藉由等式(12)、(14)、(16)、(29)及(30)針對奇 /、矩陣產生之更新程序。處理器15 〇在上述後續處理及更 新操作中使用所產生的初始化平方根矩陣Μ。 儘管上面根據一GRAKE接收器u〇來說明本發明,但熟 習此項技術者會明白本發明還應用於其他類型的接收器, 包括碼片等化器。圖5顯示一範例性碼片等化器17〇。除該 組合器130及通道估計器14〇外,碼片等化器17〇還包括一 荨化器處理1 5 2與複數個連續的延遲元件丨72。等化界處 理器152依據所接收的信號,⑺及通道估計h來計算一平方 根矩陣Μ、複數個延遲δ及複數個加權因數w。延遲元件 1 72依據計异出的延遲3來連續延遲所接收的信號γ⑺。組 合裔130中的加權元件132依據計算出的加權因數…來給連 續延遲的樣本加權。加法器134組合該等加權樣本以產生 122662.doc •19- 200814563 、Ί出樣本z(…。當向所接收信號添加或從該信號移 ^樣本π(即_延遲元件輸出)時,等化器處理器W分別 增::減小該平方根矩陣Μ,如上所述。 > ^未月確顯不,但本發明還應用於處理對應於一所發 γ守間對準的信號及/或樣本集之其他接收器。與 /木為之_型無關,本發明增大/減小一對應於該接收器 、有之協方差矩陣R的平方根矩陣μ,以應對添加/移 ’、、遲樣本集及/或#號。範例性接收器包括·· raKE接 收器’其針對每一 RAKE指峰具有-樣本集;彡個使用者 谓測接收ϋ,其針對每—使用者具有至少—樣本集;多個 輸多個輸出(ΜΙΜ0)接收器,其針對每一輸入/輸出具有 樣本集,等等。例如,考量具有·/個發射天線與尺個接收 天線之一 ΜΙΜ0系統。在接收端,該ΜΙΜΟ系統包含J個子 ^個赉射天線中的每一天線對應一子接收器。該 等/個子接收器之每一子接收器與上述接收器完全相同, 而每一子接收器係連接至接收藉由該等^個發射天線發射 勺L號之個接收天線。第j個子接收器係專用於處理藉由 第/個遠端天線發射而由所有足個接收天線接收之信號。此 第J個子接收器藉由使用與一單一天線發射的信號相關聯 之程序來計算路徑延遲、通道估計等。如等式(4)所示,該 協方差矩陣R經由該錯誤向量0捕獲其他,丨個所發射信號 對弟J個信號的接收產生之影響。 上述無線接收器100可以係具體化於任何無線裝置中, 例如一基地台或一行動終端。本文所使用的術語「行動終 122662.doc -20- 200814563 端」包括—蜂巢式無線電話與 的盔绫梦W —壬 1订包括一無線電話收發器 ::"置仃動終端亦可稱為「遍佈式計算」裝置。 本發明還應心❹—雜㈣ 裝置。此等裝置包括有線接收5 I陣之,、純號處理 廣泛的通道埒佑而一 中所接收的信號經歷 : 月一般需要某一形式的等化。此等裝置還 下,兮聲^ 麥克風_料之聲音接收器、在此情況 而受關於感測器之間的雜訊協方差之知識 又1 *此,以上範例旨在解說而非限制。 與當=:可採料同於本文明確提出的方式之其他方式來 ^丁7明而不脫離本發明之本質特徵。本具體實施例在 斤有方面欲考量為具說明性而非限制性,並 含^附以專利_之含義及等效範_的所有變化 【圖式簡單說明】 收器之 圖1顯示包含—GRAKE接收器之—範例性無線接 一方塊圖。 圖 圖2顯示依據本發明針對一 範例性更新程序之一流程 圖3顯示用以增大該協方差矩陣之一範例性流程圖。 圖4顯示用以減小該協方差矩陣之一範例性流程圖。 圖5顯示包含一碼片等化器之一範例性無線接收器 方塊圖。 【主要元件符號說明】 1〇〇 無線接收器 u〇 GRAKE接收器 122662.doc >21 - 200814563 112 接收器前端 114 天線 120 RAKE指峰 122 延遲元件 124 相關器 130 組合器 132 加權元件 134 加法器 140 通道估計器 150 處理器 152 等化器處理器 160 路徑搜尋器 170 碼片等化器 172 延遲元件
U 122662.doc -22-

Claims (1)

  1. 200814563 、申請專利範圍: 1· —種增大表示複數個樣 矩陣之—平方栌拓陆 ]的又又關聯之-協方差 又巨陣之方法,藉以反— 集,該方法包含: 所添加樣本 ^針對該所添加樣本集之減損關聯及—雜訊功率; 陣二:等減損關聯及該雜訊功率來補償該平方根矩 Ο (J :二添加具有—從該雜訊功率導出的對角元 專減知關聯導出的非對角元素: 補償的平方根矩陣。 丁…〗來增大該已 2·们之方法’其中依據該等減損關聯及該雜訊功 =補㈣平方根矩陣進—步包含依據對應㈣協方差 4之二角矩陣及一對角矩陣來補償該平方根矩陣。 :、員1之方法’其中藉由添加具有-從該雜訊功率 的對角元素與從該等減損關聯導出的非對角元素之 一行與一列來增大該已補償的平方根矩陣包含·依據該 雜訊功率之一標準差來導出該對角元素之一值。 人 、长項1之方去,其中藉由添加具有一從該雜訊功率 導出的對角元素與從該等減損關聯導出的非對角元素之 一行與一列來增大該已補償的平方根矩陣包含:依據與 该雜訊功率之一標準差成函數關係而縮放的該等減損關 聯來導出該等非對角元素。 5·如凊求項1之方法,其中該複數個樣本集包含具有來自 一 RAKE或GRAKE接收器的複數個指峰之值的複數個樣 122662.doc 200814563 本集’而其中該所添加的樣本集包含來自針對該r繞 或該GRAKE接收器之—所添加指峰的樣本值。 6 ·
    8. 如請求们之方法’其中該複數個樣本集包含具有來自 一瑪片等化器的複數個分接頭延遲輸出之值的複數個樣 本集,而其中該所添加的樣本集包含來自用於該碼片等 化器之一所添加分接頭輸出的樣本值。 如請求们之方法,其中該複數個樣本集包含且有來自 一多信號接收器的值之樣本集,而其中該所添加的樣本 集包含來自針對該多信號接收器之_所添加的輸入 之樣本值。 σ〜 包S依據該增大的平方根 ,其中包括該所添加的樣 如請求項1之方法,其進一步 矩陣來處理新的複數個樣本集 本集。 9·=求項1之方法’其進_步包含在補償該平方根矩陣 前初始化該平方根矩陣。 1〇·如請求項9之方法,其中初始化該平方根矩陣包含: 將一第一協方差矩陣設定為等於一零矩陣; 依據"個依序產生的通道估計集來對該第一協方差矩 陣作β次迭代更新,以產生一秩次π之協方差矩陣; 決定該秩次η之協方差矩陣之該平方根以產生一第一 平方根矩陣;以及 依據迭代更新的通道估計來迭代更新該第一平方根矩 陣,以產生一初始平方根矩陣。 11 ·如凊求項9之方法,其中初始化該平方根矩陣包含: 122662.doc 200814563 定義秩次一之一協方差矩陣; 一平方根矩 依據該秩次一之協方差矩陣來決定一第 陣; ^依據-奇異矩陣更新程序來對該第—協方差矩陣作, 次更新,以產生一秩次„之平方根矩陣;以及 依據-非奇異矩陣更新程序來迭代更新該秩次”之平 方根矩陣,以產生一初始平方根矩陣。 & -種減小表示複數個樣本集之間㈣又關聯之—協方差 矩陣之—平方根矩陣,以由此反映該等樣本集之一樣本 集之移除之方法,該方法包含: 決定針對-所移除樣本集之減損關聯及一雜訊功率· 從對應於該協方差矩陣之一三角矩陣及一對角矩陣兩 者,移除對應於該所移除樣本集之_列與—行,以由此 獲得減小的三角及對角矩陣;以及 依據該減小的三角矩陣、該減小的對角矩陣、該等減 損關聯及該雜訊功率來計算一減小的平方根矩陣。 13.如凊求項12之方法’其中該複數個樣本集包含且有來自 — RAKE或嶋㈣收器的複數個指峰之值的複數個樣 本集’而其中該所移除樣本集包含來自針對該rak^ 忒grake接收器之一已解配置的指 ⑷如請求項-方法,其中該複數個樣本==有來自 -碼片等化器的分接頭延遲輸出之值的複數個樣本集, 而其中該所移除樣本集包含來自針對該碼片等化器之一 已解配置的分接頭延遲輸出之樣本值。 122662.doc 200814563
    15.如請求項12之方法, 一多信號接收器的值. 包含來自對該多信號. 本值。 步包含依據該減小的平方根 集,該新的複數個樣本集不 16·如請求項12之方法,其進一步自 矩陣來處理新的複數個樣本集, 包括該所移除樣本集。 17· 士明求項12之方法,其進一步包含在減小該平方根矩陣 前初始化該平方根矩陣。 18. 如請求項17之方法,其中初始化該平方根矩陣包含: 將一第一協方差矩陣設定為等於一零矩陣; 依據"個依序產生的通道估計集來對該第一協方差矩 陣作”次迭代更新’以產生一秩次π之協方差矩陣; 決定該秩次”之協方差矩陣之該平方根,以產生一第 一平方根矩陣;以及 依據迭代更新的通道估計來迭代更新該第一平方根矩 陣,以產生一初始平方根矩陣。 19. 如請求項π之方法,其中初始化該平方根矩陣包含: 定義秩次一之一協方差矩陣; 依據該秩次一之協方差矩陣來決定一第一平方根矩 陣; 依據一奇異矩陣更新程序來對該第一協方差矩陣作η 次迭代更新以產生一秩次77之平方根矩陣;以及 依據一非奇異矩陣更新程序來迭代更新該秩次Μ之平 122662.doc 200814563 方根矩陣,以產生一初始平方根矩陣。 20· —種用以增大表示複數個樣本集之間的交叉關聯之一協 方差矩陣之一平方根矩陣之處理器,藉以反映一所添加 樣本集,該處理器經組態用以: • 》夬定針對-所添加樣本集之減損關聯及-雜訊功率; 依據該等減損關聯及該雜訊功率來補償該平方根矩 陣;以及 Γ #由添加具有—從該雜訊功率導出的對角線元素與從 該等減損關聯導出的非對角線元素之一行與一 該已補償的平方根矩陣。 θ 儿-種包括如請求項20之處理器的從磁接收器,其中該 稷數個樣本集包含來自複數個Rake指峰之值,而发中 該所添加樣本集包含來自一所添加rake指峰之樣本 U 22· 一種包括如請求項2G之處理器的碼片等化器,其中該複 數個樣本隼包合且右决ή 、 '3八有末自一碼片等化器的複數個分接頭 延遲輸出之值的複數個樣本华 、 隹—入+ 玄枣集,而其中該所添加的樣本 集匕3來自-所添加分接頭延遲之樣本值。 23. -種包括如請求項2〇之處理器的多信號接收器 複數個樣本集包含具有來白夕 ^ 八百;自一多信號接收器的值之樣本 杲’而其中該所添加的檨 ^ — 集對應於針對該多信號接收 杰之一所添加的輸入信號。 24· —種減小表示複數個 ,^ 汞之間的父叉關聯之一協方差 矩陣之一平方根矩陣之 — 裔,错以反映該等樣本集之 122662.doc 200814563 樣本集之移除,該處理器經組態用以: 決定針對-所移除樣本集之減損關聯及一雜訊 從對應於該協方差矩陣之—三角矩陣及—對角矩陣_ 除對應該所移除樣本集之—列與—行,以由此獲得一減夕 小的二角矩陣及一減小的對角矩陣;以及 依據該減小的三角矩陣、該減小的對角矩陣、 損關聯及該雜訊功率來計算一減小的平方根矩陣。/ 25·:種包括如請求項24之處理器的⑽㈣接收器, 複數個樣本集包含來自複數個rake指峰之 該所移除樣本集包含來自—已解配置的RAM指峰^ 本值。 26· -種包括如請求項24之處理器的碼片等化器,其 數個樣本集包含具有來自該碼片等化器的複數個分接頭 遲輸出之值的複數個樣本集,而其中該所移除樣本集 包含來自一已解配置的分接頭延遲之樣本值。 /、 27· -種包括如請求項24之處理器的多信號接收器,其中該 複數個樣本集包含具有來自該多信號接收器的值之樣^ 集’而其中該所移除樣本集對應於對該多信號接收器不 再關注之一輸入信號。 122662.doc
TW096125879A 2006-09-07 2007-07-16 Method for covariance matrix update TW200814563A (en)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
US11/470,676 US7769080B2 (en) 2006-09-07 2006-09-07 Method for covariance matrix update

Publications (1)

Publication Number Publication Date
TW200814563A true TW200814563A (en) 2008-03-16

Family

ID=37616913

Family Applications (1)

Application Number Title Priority Date Filing Date
TW096125879A TW200814563A (en) 2006-09-07 2007-07-16 Method for covariance matrix update

Country Status (5)

Country Link
US (1) US7769080B2 (zh)
EP (1) EP2067267A1 (zh)
CN (1) CN101512918B (zh)
TW (1) TW200814563A (zh)
WO (1) WO2008037537A1 (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7929593B2 (en) * 2008-04-15 2011-04-19 Telefonaktiebolaget Lm Ericsson (Publ) Method and apparatus for successive interference subtraction with covariance root processing
US20090323777A1 (en) * 2008-06-26 2009-12-31 Yi-Pin Eric Wang Methods and Apparatus for Sharing Signal Correlation Data in a Receiver
JP4809407B2 (ja) * 2008-09-12 2011-11-09 株式会社エヌ・ティ・ティ・ドコモ 移動通信システム、無線ネットワーク装置及び方法
US8270453B2 (en) * 2008-11-03 2012-09-18 Telefonaktiebolaget L M Ericsson (Publ) System and method for parallel interference cancellation with covariance root processing
US8483260B2 (en) * 2010-09-29 2013-07-09 Telefonaktiebolaget Lm Ericsson (Publ) Covariance factor matrix reuse
CN102811082B (zh) * 2011-06-02 2017-04-12 中兴通讯股份有限公司 一种信号处理方法及装置
US8787426B2 (en) * 2011-09-28 2014-07-22 Telefonaktiebolaget Lm Ericsson (Publ) Finger placement in multi-stage interference cancellation
US8761323B2 (en) * 2011-09-28 2014-06-24 Telefonaktiebolaget Lm Ericsson (Publ) Impairment covariance and combining weight updates during iterative turbo interference cancellation reception
CN103532888B (zh) * 2012-07-02 2017-11-21 中兴通讯股份有限公司 数据处理装置以及消除干扰的方法
EP2871789B1 (en) * 2013-11-12 2017-01-04 Huawei Technologies Co., Ltd. Method for estimating covariance matrices and use thereof
US20160014619A1 (en) * 2014-07-09 2016-01-14 Qualcomm Incorporated Multiple cell joint detection and interference cancellation

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6363104B1 (en) * 1998-10-02 2002-03-26 Ericsson Inc. Method and apparatus for interference cancellation in a rake receiver
US7058147B2 (en) * 2001-02-28 2006-06-06 At&T Corp. Efficient reduced complexity windowed optimal time domain equalizer for discrete multitone-based DSL modems
US6738682B1 (en) * 2001-09-13 2004-05-18 Advances Micro Devices, Inc. Method and apparatus for scheduling based on state estimation uncertainties
US6947715B2 (en) * 2002-03-30 2005-09-20 Broadcom Corporation VOFDM receiver correlation matrix processing using factorization
EP1357714A1 (en) 2002-04-23 2003-10-29 Siemens Information and Communication Networks S.p.A. Method and device for channel estimation in digital radio communication systems
US7200631B2 (en) * 2003-01-10 2007-04-03 Lucent Technologies Inc. Method and apparatus for determining an inverse square root of a given positive-definite hermitian matrix
US7746967B2 (en) * 2003-04-30 2010-06-29 Intel Corporation Beam-former and combiner for a multiple-antenna system
FI20040182A0 (fi) * 2004-02-06 2004-02-06 Nokia Corp Tietojenkäsittelymenetelmä, korjain ja vastaanotin
US8045638B2 (en) * 2004-03-05 2011-10-25 Telefonaktiebolaget Lm Ericsson (Publ) Method and apparatus for impairment correlation estimation in a wireless communication receiver
EP1723744B1 (en) 2004-03-12 2014-12-24 Telefonaktiebolaget L M Ericsson (publ) A method and apparatus for received signal quality estimation
US7539240B2 (en) * 2004-03-12 2009-05-26 Telefonaftiebolaget Lm Ericsson (Publ) Method and apparatus for parameter estimation in a generalized rake receiver
US7536158B2 (en) * 2004-03-29 2009-05-19 Telefonaktiebolaget Lm Ericsson (Publ) Impairment correlation estimation in a spread spectrum system

Also Published As

Publication number Publication date
US7769080B2 (en) 2010-08-03
WO2008037537A1 (en) 2008-04-03
CN101512918A (zh) 2009-08-19
US20080063033A1 (en) 2008-03-13
CN101512918B (zh) 2012-12-12
EP2067267A1 (en) 2009-06-10

Similar Documents

Publication Publication Date Title
TW200814563A (en) Method for covariance matrix update
JP4369294B2 (ja) 雑音電力推定装置、雑音電力推定方法及び信号検出装置
JP4829239B2 (ja) スペクトル拡散受信機における遅延選択のための装置、方法およびコンピュータ読取り可能な記録媒体
TWI247549B (en) Code tracking loop with automatic power normalization
JP4838250B2 (ja) 白色化フィルタを生成する方法、受信機回路、受信機およびコンピュータ可読記憶媒体
TWI302409B (en) Advanced receiver with sliding window block linear equalizer
TWI305450B (en) Method and system for receiving a wireless code division multiple access (cdma) signal
KR20070007029A (ko) 부호 특정 결합을 이용한 ds-cdma 간섭 억제를 위한방법 및 장치
JP2008544689A (ja) 1つ以上の受信信号処理の工程のなかでチップサンプル相関を用いるための方法と装置
JP2009206809A (ja) 無線通信装置、イコライザ、イコライザ重み係数演算プログラムおよびイコライザ重み係数演算方法
WO2012041675A1 (en) Covariance factor matrix reuse
US7724841B2 (en) Method and apparatus for iteratively calculating channel response estimates
WO2003021972A2 (en) Multi-user detection
WO2004088867A1 (ja) 受信装置
Lohan Multipath delay estimators for fading channels with applications in CDMA receivers and mobile positioning
US8213486B2 (en) Method for simplifying the calculations for pre-whitening in a G-rake receiver
CN103988444B (zh) 非冗余均衡
TW200522623A (en) Reduced complexity sliding window based equalizer
JP3747405B2 (ja) アンテナベリフィケーション方法及びアンテナベリフィケーション処理装置
JP4815638B2 (ja) 無線通信システムのためのマルチユーザ等化方法及び装置
US7630429B2 (en) Equalizer co-efficient generation apparatus and method therefor
CN112688889B (zh) 一种无人机测控系统的信道估计方法和装置
CN101345562B (zh) 耙式接收器及其解展频方法
CN110518921B (zh) 多径合并的方法、装置和存储介质以及电子设备
Fan et al. An adaptive linear prediction algorithm for joint blind equalization and blind multiuser detection in CDMA