專利名稱：多輸入多輸出系統(tǒng)中的串行干擾消除方法
技術(shù)領(lǐng)域：
本發(fā)明涉及信號檢測技術(shù)，特別涉及多輸入多輸出(MIMO)系統(tǒng)中多 串行干擾消除方法。
背景技術(shù)：
在多輸入多輸出系統(tǒng)中，用戶信號經(jīng)復(fù)用后從多根發(fā)射天線發(fā)射出去， 在接收端通過多根接收天線進行接收。基于最小均方誤差(MMSE)或迫零 (ZF)算法的垂直貝爾實驗室分層空時結(jié)構(gòu)(VBLAST)是一種低復(fù)雜度的 串行4企測方法。
其中，利用ZF算法進行串行干擾消除的基本思想是接收端在接收到 的信號中，按照一定的順序依次檢測各個發(fā)送信號，在檢測出每個發(fā)送信號 后，將該發(fā)送信號對其它發(fā)送信號的干擾消除掉，然后在消除干擾后的信號 中繼續(xù)檢測下一個發(fā)送信號，依此循環(huán)，直到檢測出所有的發(fā)送信號。利用 MMSE算法進行串行干擾抵消的信號檢測時，其基本思想與ZF算法相同， 同樣是按照一定順序檢測各個發(fā)送信號，與ZF算法的區(qū)別在于，在檢測出 每個發(fā)送信號并進行干擾消除時，除將檢測出的發(fā)送信號本身對其它發(fā)送信 號造成的干擾外，進一步消除各種噪聲對發(fā)送信號造成的干擾?？梢?，由于 考慮了其它噪聲干擾的影響，MMSE算法相比于ZF算法具有更好的檢測性 能。
在上述采用ZF或MMSE等串行干擾消除的接收方法時，發(fā)送信號的串 行檢測順序?qū)τ谙到y(tǒng)接收性能的影響很大。目前，采用VBLAST方案結(jié)合 ZF或MMSE進行串行干擾消除時，確定信號檢測順序的方式為在每才全測 出一個發(fā)送信號后，確定下一個要檢測的發(fā)送信號時，首先計算剩余的待檢測發(fā)送信號的信噪比，選擇信噪比最大的發(fā)送信號作為下一個檢測的發(fā)送信 號。在每次確定下一個檢測的發(fā)送信號過程中，均需要進行矩陣的求逆或偽 逆運算，而矩陣逆運算的復(fù)雜度很高，并且，隨著發(fā)射天線數(shù)的增加，需要 檢測的發(fā)送信號增加，檢測過程中對矩陣求逆或偽逆的次數(shù)也相應(yīng)增加，從
而導(dǎo)致整個VBLAST信號檢測的計算復(fù)雜度也成倍增加。

發(fā)明內(nèi)容
有鑒于此，本發(fā)明提供一種在MIMO系統(tǒng)中進行串行干擾抵消時確定 檢測順序的方法，能夠有效降低VBLAST方案中信號檢測的計算復(fù)雜度和 提高信號檢測性能。
為實現(xiàn)上述目的，本發(fā)明采用如下的技術(shù)方案
一種多輸入多輸出系統(tǒng)中的串行干擾消除方法，包括
a、 獲取發(fā)送端到接收端的信道傳輸矩陣好；
b、 根據(jù)步驟a中所述信道傳輸矩陣^確定當(dāng)前信道傳輸矩陣，設(shè)置加 權(quán)矩陣0和第 一排序矩陣P分別為M x M維和x w維的單位矩陣，所述M 和iV分別為所述當(dāng)前信道傳輸矩陣的行數(shù)和列數(shù)，并將當(dāng)前信道傳輸矩陣的 第一列作為當(dāng)前列，即&=1, A表示當(dāng)前列的列號；
c、 從當(dāng)前信道傳輸矩陣的當(dāng)前列到最后一列，在每列中提取第&個元 素到最后一個元素構(gòu)成該列的重構(gòu)列向量，計算每列的重構(gòu)列向量的范數(shù)， 確定范數(shù)最小的重構(gòu)列向量所在列的列號《，并對該列的重構(gòu)列向量進行豪 思霍德householder變換，生成向量v和標(biāo)量-;根據(jù)向量v和標(biāo)量"構(gòu)造歸 一化正交矩陣込；
d、 將所述込與當(dāng)前信道傳輸矩陣的乘積作為更新后的信道傳輸矩陣， 將矩陣g,與當(dāng)前加權(quán)矩陣0的乘積作為更新后的加權(quán)矩陣；將更新后的信 道傳輸矩陣的第A列和第《列交換，將當(dāng)前第一排序矩陣的第^列和第《列 交換；
e、 將當(dāng)前信道傳輸矩陣的下一列作為當(dāng)前列，即A^A+1,并判斷A:是否等于當(dāng)前信道傳輸矩陣的行數(shù)和列數(shù)中的最小值，若是，則執(zhí)行步驟f，
否則返回步驟c;
f、根據(jù)當(dāng)前加權(quán)矩陣和步驟b中所述當(dāng)前信道傳輸矩陣確定上三角矩 陣及，根據(jù)矩陣及對各個信號進行串行干擾消除。
較佳地，所述根據(jù)當(dāng)前加權(quán)矩陣和步驟b中所述當(dāng)前信道傳輸矩陣確定 上三角矩陣及為將當(dāng)前加權(quán)矩陣與步驟b中獲取的所述信道傳輸矩陣以及 當(dāng)前第 一排序矩陣的共軛轉(zhuǎn)置矩陣之乘積作為矩陣及。
較佳地，步驟b中所述根據(jù)步驟a中的信道傳輸矩陣^確定當(dāng)前信道 傳輸矩陣為將步驟a中所述信道傳輸矩陣丑作為當(dāng)前信道傳輸矩陣；
所述根據(jù)矩陣W對各個信號進行串行干擾消除為
在接收端，計算當(dāng)前加權(quán)矩陣與接收信號的乘積，利用該乘積結(jié)果和矩 陣及，對接收信號進行串行干擾消除，檢測發(fā)射信號。
較佳地，步驟b中所述根據(jù)信道傳輸矩陣^確定當(dāng)前信道傳輸矩陣包 括根據(jù)信道傳輸矩陣丑構(gòu)造擴展信道矩陣及e-[及T，cj2/]1，其中，(72是系 統(tǒng)噪聲的均方差，J是A^x^的單位矩陣；將擴展信道矩陣作為當(dāng)前信道傳 輸矩陣；
所述根據(jù)矩陣及對各個發(fā)射信號進行串行干擾消除為在接收端，對 接收到的信號進行擴展，得到擴展后的接收信號為re=[/,0MxlT]T,建立擴展 后的接收信號、當(dāng)前加權(quán)矩陣0、矩陣及、第一排序矩陣P和發(fā)射信號rf的 函數(shù)關(guān)系為^e=i iW-其中，及'=[tfMxiVo/]T， /是A^x7V,的單 位矩陣；根據(jù)矩陣及和所述函數(shù)關(guān)系，檢測發(fā)射信號d,并進行串行干擾消 除。
較佳地，步驟b中所述根據(jù)信道傳輸矩陣丑確定當(dāng)前信道傳輸矩陣包 括計算步驟a中獲取的信道傳輸矩陣丑的共軛轉(zhuǎn)置矩陣^11,將計算結(jié)果 ^M乍為當(dāng)前信道傳輸矩陣；
所述根據(jù)矩陣及對各個信號進行串行預(yù)干擾消除為
在發(fā)送端，計算當(dāng)前第一排序矩陣的共軛轉(zhuǎn)置矩陣與矩陣及的共軛轉(zhuǎn)置矩陣的乘積i^及H,從選擇的第一發(fā)射信號開始，根據(jù)/^及H依次計算并 消除其它發(fā)射信號對該發(fā)射信號的干擾，并將消除干擾后的發(fā)射信號進行發(fā) 射。較佳地，所述根據(jù)P"及H依次消除其它發(fā)射信號對該發(fā)射信號的干擾為 '' j^^d"，j),其中，爿為取模除數(shù)，并由選擇的調(diào)制方式?jīng)Q定，A'為P^力H中的第k行第i個元素，《為第i發(fā)射信號，4為消除 干擾后的第i發(fā)射信號。較佳地，當(dāng)發(fā)射天線數(shù)小于等于接收天線數(shù)時，步驟b中所述根據(jù)信道 傳輸矩陣F確定當(dāng)前信道傳輸矩陣包括對步驟a中獲取的信道傳輸矩陣丑進行減格LR算法計算，得到矩陣及r和A,使及r-丑尸r,并將及r作為當(dāng)前信道傳輸矩陣；在步驟b中進一步包括將第二排序矩陣尸'初始化為矩陣A;在步驟d中執(zhí)行所述將當(dāng)前第 一排序矩陣的第A列和第《列交換的操作 前，進一步包括將當(dāng)前第一排序矩陣P設(shè)置為單位矩陣；在步驟d和e之間進一步包括對當(dāng)前信道傳輸矩陣的分矩陣進行LR 計算，生成矩陣及t和尸t，根據(jù)戶t構(gòu)造Mx7Vt維的幺模矩陣込'，將當(dāng)前第二 排序矩陣、當(dāng)前第一排序矩陣和矩陣g/的乘積作為更新后的第二排序矩陣， 將當(dāng)前信道傳輸矩陣與矩陣g/的乘積作為更新后的信道傳輸矩陣；步驟f中所述根據(jù)當(dāng)前加權(quán)矩陣和步驟b中所述當(dāng)前信道傳輸矩陣確定 上三角矩陣及為將當(dāng)前加權(quán)矩陣、步驟a中獲取的所述信道傳輸矩陣丑、 當(dāng)前第二排序矩陣/>'的乘積作為上三角矩陣及。較佳地，所述當(dāng)前信道傳輸矩陣的分矩陣為當(dāng)前信道傳輸矩陣的第 A:+l行到第M行中的第A+l個元素到第7V個元素構(gòu)成的矩陣。較佳地，所述根據(jù)A構(gòu)造Mxi^維的幺模矩陣g/為將當(dāng)前信道傳輸 矩陣的第壯l行到第M行中的第A+l個元素到第7V個元素應(yīng)用LR減格算法構(gòu)成矩陣込'。較佳地，所述根據(jù)矩陣及對各個發(fā)射信號進行串行干擾消除為 在接收端，建立接收信號r、當(dāng)前第二排序矩陣戶'、當(dāng)前加權(quán)矩陣g、上三角矩陣及和發(fā)射信號d的函數(shù)關(guān)系為gr-及(尸')-V+01,根據(jù)所述函數(shù)關(guān)系，利用迫零ZF或MMSE算法檢測發(fā)射信號rf，并進行串行干擾消除。 較佳地，所述根據(jù)向量v和標(biāo)量"構(gòu)造歸一化正交矩陣込為在向量v的上面添加U-1)個零，構(gòu)成向量v',即v'=[^-1)xl,vT]T，再奸算/m^VV 將計算結(jié)果作為矩陣込。由上述技術(shù)方案可見，本發(fā)明中，獲取發(fā)送端到接收端的信道傳輸矩陣 丑，并基于househoulder變換進行排序QR分解獲取上三角矩陣及，根據(jù)矩 陣及對各個信號進行串行干擾消除。通過上述過程對信道傳輸矩陣進行排序 QR分解后，得到矩陣R,由于該矩陣R是一個上三角矩陣，因此可以利用 上三角矩陣特性依次對各個信號進行串行干擾消除，從而避免利用矩陣求逆 或偽逆的方式確定串行;f企測順序，進而簡化串行干擾消除的計算復(fù)雜度。


圖1為本發(fā)明實施例一中串行干擾消除的具體流程圖。 圖2為本發(fā)明實施例二中串行干擾消除的具體流程圖。 圖3為應(yīng)用實施例一和二的方法，與無排序的ZF和MMSE算法相比，進行信號檢測的性能比較示意圖。圖4為本發(fā)明實施例三中串行干擾消除的具體流程圖。 圖5為本發(fā)明實施例四中串行干擾消除的具體流程圖。 圖6為應(yīng)用實施例四的方法、背景技術(shù)中介紹的ZF算法、ZF算法與LR算法結(jié)合時，進行信號檢測的性能比較示意圖。
具體實施方式
為使本發(fā)明的目的、技術(shù)手段和優(yōu)點更加清楚明白，以下結(jié)合附圖，對本發(fā)明做進一步詳細說明。本發(fā)明的基本思想是對信道傳輸矩陣進行排序的QR分解，利用QR 分解后得到的上三角矩陣R進行串行干擾消除。在本發(fā)明中，進行排序的QR分解時，采用基于豪思霍德(householder) 變換的方式進行。在具體實施時，可以進行單純的排序QR分解獲取上三角 矩陣R，再利用該上三角矩陣進行發(fā)射信號檢測和串行干擾消除；另外，為 進一步改善系統(tǒng)的信號檢測性能，也可以通過將排序QR分解與減格(LR， Lattice Reduction)算法相結(jié)合，獲取上三角矩陣R,再利用該上三角矩陣進 行發(fā)射信號檢測和串行干擾消除。下面將對這兩種方式下的具體實施作詳細 描述。實施例一在本實施例中，通過單純的排序QR分解過程獲取上三角矩陣。在M 根發(fā)射天線、M根接收天線的MIMO系統(tǒng)中，圖1為本發(fā)明實施例一中串 行干擾消除的具體流程圖。如圖l所示，該方法包括步驟IOI,獲取發(fā)送端到接收端的信道傳輸矩陣丑。本實施例中，在接收端進行信道估計，從而獲取發(fā)送端到接收端的信道 傳輸矩陣H。接下來，通過步驟102~ 106對信道傳輸矩陣^進行QR分解，其中， 設(shè)該QR分解的矩陣的維數(shù)為Mx見本實施例中，M=7Vr, A^Vt。步驟102，將步驟101獲取的信道傳輸矩陣丑作為當(dāng)前信道傳輸矩陣， 對其進行QR分解；根據(jù)信道傳輸矩陣丑將加權(quán)矩陣g和第一排序矩陣P 分別初始化為單位矩陣，并將當(dāng)前信道傳輸矩陣的第 一列作為當(dāng)前列。具體地，將加權(quán)矩陣g初始化為單位矩陣，并且該單位矩陣的維數(shù)與 當(dāng)前信道傳輸矩陣的行數(shù)相等，即^ = /^將第一排序矩陣戶初始化為單位 矩陣，并且該單位矩陣的維數(shù)與當(dāng)前信道傳輸矩陣的列數(shù)相等，即尹=&。 利用變量A:代表當(dāng)前信道傳輸矩陣的當(dāng)前列列號，因此，本步驟中A:-1。步驟103,從當(dāng)前信道傳輸矩陣的當(dāng)前列到最后一列，在每列中提取第A:個元素到最后一個元素構(gòu)成該列的重構(gòu)列向量，計算每列的重構(gòu)列向量的范數(shù)；并確定計算結(jié)果中范數(shù)最小的重構(gòu)列向量，確定該重構(gòu)列向量所在列的列號《。信道傳輸矩陣為MxN維的矩陣，它表示發(fā)射端到接收端的信道傳輸特 性，其中的任意元素hij表示第y根發(fā)射天線到第/根接收天線的信道沖激響 應(yīng)，矩陣中當(dāng)前列的列向量表示第A根發(fā)射天線到接收端的信道傳輸特性。本步驟計算重構(gòu)列向量范數(shù)的方式為計算該向量中各個元素的模的平 方和， 再對該平方和取平方才艮， 即 Normall = (IVJ2+IU2+…十IV,12)1 其中，片A:,…,W。在上述計算的各個重構(gòu) 列向量范數(shù)中，選擇范數(shù)最小的重構(gòu)列向量所在列的列號，即《= min(NormALL(/::琳這里，重構(gòu)列向量的范數(shù)表征了該列對應(yīng)的發(fā)射天線到第k根至最后一 根接收天線的信道增益狀況，由于該增益狀況近似反映了信號的信噪比狀 況，因此，選擇范數(shù)最小的重構(gòu)列向量所在的列，也就是近似選擇了這些列 對應(yīng)的發(fā)射信號中信噪比最小的 一根發(fā)射天線。步驟104,對第《列的重構(gòu)列向量進行householder變換，生成向量v 和標(biāo)量A;根據(jù)向量v和標(biāo)量"構(gòu)造MxM維的歸一化正交矩陣込。本步驟對第q列的重構(gòu)列向量進行豪思霍德(householder)變換。householder變換中，輸入一個w維的列向量經(jīng)Householder(x)的變 換操作后生成w維的向量v和標(biāo)量"，并且這兩個量滿足5 = / -々 "為 一wx"維正交歸一矩陣，5* jc:norm(jc"d,其中，q為wxl的列向量，該 A的第一個元素為l,其他元素為0。根據(jù)householder變換的上述特性，可以利用向量v和標(biāo)量 〃 構(gòu)造Mx M維的歸一化正交矩陣込。具體地，首先在向量v的上面添加(k-l)個零， 構(gòu)成M維的向量v,即v'=
T，然后，再計算J,"v' '^將計算結(jié) 果作為矩陣込，根據(jù)householder變換的特性可得，矩陣込是個MxM維的 歸一化正交矩陣。步驟105,將矩陣込與當(dāng)前信道傳輸矩陣的乘積作為更新后的信道傳輸 矩陣，將矩陣込與當(dāng)前加權(quán)矩陣^的乘積作為更新后的加權(quán)矩陣；將更新 后的信道傳輸矩陣的第A:列和第《列交換，將當(dāng)前第一排序矩陣的第A:列和 第《列交換。本步驟利用構(gòu)造的矩陣込更新當(dāng)前信道傳輸矩陣和加權(quán)矩陣，具體地， 好=込*丑，g=g*2,。將更新后的信道傳輸矩陣的第A:列和第《列交換，將當(dāng) 前第一排序矩陣的第A:列和第《列交換，也就是將當(dāng)前列和步驟103中選擇 的范數(shù)最小的重構(gòu)列向量所在的列進行交換。如前所述，范數(shù)最小的重構(gòu)列 向量所在的列被近似認(rèn)為與具有最小信噪比的發(fā)射天線對應(yīng)，將該列與當(dāng)前 列交換后，即將該列置前。具體地，在第一次循環(huán)時，則在所有列中確定范數(shù)最小的重構(gòu)列 向量所在列的列號《，并將第一列與第《列交換，使第《列成為第一列；在 第二次循環(huán)時，在剩余的第二列到最后一列中確定范數(shù)最小的重構(gòu)列向量所 在列的列號《，并將第二列與第《列交換；如此循環(huán)下去，未被選擇的列中 范數(shù)最小的列置前，直到將所有的列選擇完畢。形成近似按照信噪比遞增的 順序排列的加權(quán)矩陣、信道傳輸矩陣和第 一 排序矩陣。步驟106，將當(dāng)前信道傳輸矩陣的下一列作為當(dāng)前列，即A:-A+1,并判 斷k是否等于當(dāng)前信道傳輸矩陣的行數(shù)和列數(shù)中的最小值，若是，則執(zhí)行步 驟107,否則返回步驟103。到本步驟前，已完成一次排序過程。本步驟中，將當(dāng)前信道傳輸矩陣的 下一列作為當(dāng)前列，即令k-k+l，進行下一次循環(huán)。在進行下一次循環(huán)前， 首先判斷k是否等于當(dāng)前信道傳輸矩陣的行數(shù)和列數(shù)中的最小值，若是，則 表明QR分解已完成，可以結(jié)束循環(huán)執(zhí)行步驟107，利用加權(quán)矩陣和信道傳 輸矩陣計算上三角矩陣；否則，表明QR分解仍未完成，返回步驟103,繼 續(xù)進行下一次循環(huán)。步驟107,將當(dāng)前加權(quán)矩陣的轉(zhuǎn)置矩陣、步驟101中獲得的所述當(dāng)前信 道傳輸矩陣以及當(dāng)前第 一排序矩陣的乘積作為矩陣及。經(jīng)過上述步驟101 ~ 106的操作對信道傳輸矩陣/T進行排序QR分解后， 得到加權(quán)矩陣g和第一排序矩陣尸，根據(jù)上述信道傳輸矩陣和QR分解后的 結(jié)果計算上三角矩陣及的方式為及=0及,(1)。由于矩陣R是個上三角矩陣，因此可以利用該矩陣對檢測信號進行串 行干擾消除和信號檢測。在本實施例中，采用串行迭代方法，通過上三角矩 陣進行串行干擾消除及信號檢測。具體地，由于^及d+w，由公式(l)可得 Z^gH^CP")-1 = gH及凡將其代入上式可得，r=gH^iW+/i,將該式左右同乘0 可得，gr=l iV+gw。將iW看作發(fā)射信號排序過程，由于通過上述QR分解已經(jīng)確定了加權(quán) 矩陣g,并且接收信號f已知，于是可以求得gr;又矩陣及為上三角矩陣， 因此，可以/人/W向量中的最后一個元素開始^r測，并在4企測完該元素后， 將該元素對其它元素的干擾消除掉，再對倒數(shù)第二個元素進行檢測，直到檢 測到第一個元素，從而能夠獲得/W的取值。然后，再將該向量左乘第一排 序矩陣的逆矩陣P"1即得到發(fā)射信號J。具體利用上述關(guān)系式進行信號檢測和串行干擾消除的步驟包括步驟108，在接收端，對利用濾波特性為g的濾波器對接收信號進行濾波，得到向量0"。步驟109,依據(jù)0*和及計算7W。步驟110，將向量i^左乘第一排序矩陣的逆矩陣FM尋到發(fā)射信號J。 經(jīng)過上述步驟108 ~ 110即完成了接收端的信號檢測。 至此，本實施例中進行串行干擾消除的方法流程結(jié)束。由上述流程可見， 本發(fā)明中，通過排序QR分解確定上三角矩陣及，并且根據(jù)上三角矩陣與發(fā) 射信號和接收信號的關(guān)系，進行發(fā)射信號的檢測和串行干擾消除，該信號檢 測過程中近似采用的信噪比遞減的順序進行，能夠獲得較佳的檢測性能。同 時，由于該檢測過程中避免利用矩陣求逆或偽逆過程，因此降低了信號檢測 的復(fù)雜度。在本實施例中，采用了單純的排序QR分解獲取及矩陣，并且利用ZF檢測器進行信號檢測。事實上，通過MMSE檢測器與上述排序QR分解過 程相結(jié)合，由于MMSE過程考慮了信號噪聲，因此能夠獲得更優(yōu)的檢測性 能。下面通過實施例二對這種實施方式進行描述。 實施例二本實施例中，仍然通過單純的排序QR分解獲取R矩陣。具體地，圖2 為本發(fā)明實施例二中串行干擾消除的具體流程圖。如圖2所示，該方法包括步驟201，獲取發(fā)送端到接收端的信道傳輸矩陣^，并構(gòu)造擴展信道矩 陣W^[及T,cj2/]1,將其作為當(dāng)前信道傳輸矩陣。本步驟中，接收端通過信道估計獲取發(fā)送端到接收端的信道傳輸矩陣 再根據(jù)信道傳輸矩陣/r確定當(dāng)前信道傳輸矩陣。具體確定當(dāng)前信道傳輸 矩陣的方式為根據(jù)信道傳輸矩陣丑構(gòu)造擴展信道矩陣」^=[〃7，(12/]7，其 中，(12是系統(tǒng)噪聲的均方差，/是7V^V,的單位矩陣。步驟202,對擴展信道矩陣，利用householder變換進行排序QR分解。本實施例中，將擴展信道矩陣作為當(dāng)前的信道傳輸矩陣，因此對擴展信 道矩陣進行排序QR分解，具體分解方式與實施例一中步驟102- 106的過 程相同。在步驟106中，當(dāng)判斷k與當(dāng)前信道傳輸矩陣的行數(shù)和列數(shù)的最小 值相等時，執(zhí)行步驟203。步驟203，將當(dāng)前加權(quán)矩陣、步驟201中確定的擴展信道矩陣以及當(dāng)前 第 一排序矩陣的共軛轉(zhuǎn)置矩陣之乘積作為矩陣及e。經(jīng)過上述步驟201 ~ 203的操作對擴展信道傳輸矩陣好e進行排序QR分 解后，得到加權(quán)矩陣g和第一排序矩陣尸，根據(jù)上述信道傳輸矩陣和QR分 解后的結(jié)果計算上三角矩陣及的方式為2 e=0e^, (2)。在本實施例中，采用MMSE檢測器，通過上三角矩陣進行串行干擾消除及信號檢測。具體地，由于F^W+"，由公式(2)可得仏-ge^ eCP")-^ geH|U> (3)。其中，在本實施例中，,-F1。步驟204,在接收端，擴展接收信號為re=[/,0MxlT]T,并根據(jù)gere和本步驟中，首先擴展接收信號為i e=[/,^xlT]T;根據(jù)接收信號、發(fā)射信號和信道傳輸矩陣的關(guān)系r=/w+",以及擴展后的接收信號和擴展信道傳輸 矩陣的表達式，可以得到擴展后的接收信號、擴展信道傳輸矩陣和發(fā)射信號間的關(guān)系，具體為re= ^rf+w-(ji V/ (4),其中，及'=[tf瓶&， / ]t, J是^xT^ 的單位矩陣。將公式(3)代入公式(4)可得re= geH/U/+ii-oi '</,將該式的等號左右兩邊同乘ge得到gje= J^/W+仏n- ge(J及W,由于/^i^/，因此，可以利用函數(shù)關(guān)系2jf及^^ - 0e(W,^"計算向量尸J的取值。步驟205 ,將向量iV/左乘第 一排序矩陣的逆矩陣戶"得到發(fā)射信號</。經(jīng)過上述步驟204 ~ 205即完成了接收端的信號檢測。至此，本實施例中進行串行干擾消除的方法流程結(jié)束。本實施例中，對 考慮了信道噪聲的擴展信道矩陣進行排序QR分解，從而確定上三角矩陣R。 并進一步根據(jù)擴展信道矩陣和擴展后接收信號的關(guān)系進行發(fā)射信號的檢測， 與實施例一中類似。同時，與實施例一相比，由于考慮到了信道噪聲，因此 具有更優(yōu)的檢測效果。為說明本發(fā)明的優(yōu)點，對應(yīng)用上述實施例一和實施例二的系統(tǒng)進行了仿 真實^r。圖3為應(yīng)用實施例一和二的方法，與不進行排序、單純采用ZF和 MMSE算法相比，進行信號檢測的性能比較示意圖。該系統(tǒng)性能仿真的仿真 參數(shù)為4根發(fā)射天線和4根接收天線的MIMO系統(tǒng)，系統(tǒng)信道為最大時延 為16碼片的瑞利衰落信道，系統(tǒng)采用256個子載波的OFDM技術(shù)消除頻率 選擇性信道產(chǎn)生的符號間干擾，在一個OFDM符號的傳輸時間內(nèi)信道是平 穩(wěn)的，在一幀中的OFDM符號間是相關(guān)的，不同幀間的OFDM符號是不相 關(guān)的，調(diào)制方式為4-QAM,接收端能夠進行準(zhǔn)確的信道估計。如圖3所示， 橫軸表示接收端檢測信號的信噪比，縱軸表示接收端檢測信號的誤碼率。其中，曲線301是無排序的ZF算法進行信號檢測時信噪比與誤碼率的 關(guān)系曲線；曲線302是應(yīng)用實施例一的方式進行信號檢測時信噪比與誤碼率 的關(guān)系曲線；曲線303是無排序的MMSE算法進行信號檢測時信噪比與誤 碼率的關(guān)系曲線；曲線304是應(yīng)用實施例二的方式進行信號檢測時信噪比與誤碼率的關(guān)系曲線。由曲線301和302的比較可見，本發(fā)明實施例一的方式與無排序的ZF 算法進行信號檢測相比，能夠帶來更好的檢測性能，在高SNR下，能夠帶 來約2dB的增益。又由曲線303和304的比較可見，本發(fā)明實施例二的方式 與無排序的MMSE算法進行信號檢測相比，也能夠帶來更好的檢測性能， 在誤碼率為10—2時，系統(tǒng)增益約為4dB。另外，由曲線302和304的比較可 見，本發(fā)明實施例二的方式與實施例一的方式相比，也能夠帶來一定的系統(tǒng) 檢測性能的增益，說明排序QR分解與MMSE算法結(jié)合能夠具有更優(yōu)的檢 測性能。在上述實施例一和實施例二中，均是通過單純的排序QR分解獲取上三 角矩陣，然后結(jié)合接收端進行的ZF或MMSE算法進行信號檢測。下面的實 施例三也將通過單純的排序QR分解獲取上三角矩陣R,但是利用該上三角 矩陣R的方式與實施例一和實施例二不同，將在發(fā)送端利用該矩陣R進行 串行干擾消除，具體實施方式
如下所述。實施例三在Nt根發(fā)射天線、Nr根接收天線的MIMO系統(tǒng)中，圖4為本發(fā)明實施 例三中串行干擾消除的具體流程圖。如圖4所示，該方法包括步驟401，獲取發(fā)送端到接收端的信道傳輸矩陣^。本實施例中，發(fā)送端獲取發(fā)送端到接收端的信道傳輸矩陣^，具體方式 可以為在接收端進行信道估計獲取發(fā)送端到接收端的信道傳輸矩陣丑，然 后接收端將估計結(jié)果反饋給發(fā)送端；或者在上下行信道特性相同的系統(tǒng)環(huán)境 中，發(fā)送端進行信道估計獲取接收端到發(fā)送端的信道傳輸矩陣，并將信道估 計結(jié)果作為發(fā)送端到接收端的信道傳輸矩陣及。步驟402，計算信道傳輸矩陣的共軛轉(zhuǎn)置矩陣^，基于householder變 換對信道傳輸矩陣f進行排序QR分解，生成加權(quán)矩陣g、第一排序矩陣 P和上三角矩陣及。本實施例中，將信道傳輸矩陣的共軛轉(zhuǎn)置矩陣^M乍為當(dāng)前信道傳輸矩陣。本步驟中對信道傳輸矩陣^進行排序QR分解的過程與實施例一中步 驟102 107相同，這里就不再贅述。進行QR分解后，加權(quán)矩陣0、第一 排序矩陣P、上三角矩陣及和信道傳輸矩陣^的關(guān)系為及-g^1, (5)步驟403，在發(fā)送端，根據(jù)發(fā)射信號確定湯姆林森-哈利西瑪 (Tomlinson-Harishima )預(yù)編碼后的發(fā)射信號J 。本步驟中，由上述公式(5)可得信道傳輸矩陣丑-i^及Hg,將其代入 r=^/+"可得1^P"及Hgi/ + n。選擇t^gH作為發(fā)送端的預(yù)編碼矩陣，則有r -/^^Hgf/rf + n = i^及Hj + n。其中，,及H是對矩陣及H的列進行重新排列后得到的矩陣，由于及H是個下三角矩陣，i^及H重新排列后的矩陣中，各個數(shù) 據(jù)流間的相互干擾仍和及H —樣，因此在發(fā)送端仍可以得到各個發(fā)射信號(《}間的干擾關(guān)系，其中，《為向量rf中的第i個元素，具體地，其它信號對第i個發(fā)射信號4的干擾為E二^A , A,i是矩陣^力H中的元素，Z' =1,...， M， A=1,...,A^， ^為Tomlinson-Harishima預(yù)編碼后的發(fā)射信號^的第A個元素。根據(jù)上述發(fā)射信號間的干擾，在發(fā)送端通過Tomlinson-Harishima預(yù)編 碼進行串行干擾消除，具體為j =j — y"i義，i廣mod",力，mod為取A,*模操作，取模除數(shù)大小為爿，由所采用的調(diào)制方式?jīng)Q定，從而得到 Tomlinson-Harishima預(yù)編碼后的發(fā)射信號^ 。至此，本實施例中進行預(yù)編碼的方法流程結(jié)束。在接收端，首先獲取矩 陣P"及h，具體獲取方式可以是接收端按照與發(fā)送端相同的方式完成QR 分解獲取P矩陣和R矩陣，然后計算i^i H,或者，也可以由發(fā)送端將矩陣 i^及H中的元素^w直接發(fā)送到接收端。然后，接收端對接收到的上述發(fā)射信 號利用相應(yīng)的矩陣^及h中的元素7 w進行均衡，其中，A:=7，...i:,〖為最大 發(fā)送流的數(shù)目，再對均衡結(jié)果進行與發(fā)送端相同的取模操作。經(jīng)過上述處理 后的接收信號即為檢測后的發(fā)射信號rf。在本實施例中，通過在發(fā)送端利用 上三角矩陣R進行Tomlinson-Harishima預(yù)編碼實現(xiàn)了發(fā)射信號間的串行干擾消除，該處理流程在發(fā)送端進行，相應(yīng)地在接收端的處理過程簡化，降低 了信號檢測時的計算復(fù)雜度，并且此方法沒有在接收端的串行干擾消除的錯 誤傳播帶來的性能損失。在上述三個實施例中，獲取上三角矩陣及的方式均是通過單純進行排序QR分解獲得的，在下面的實施例中，通過排序QR分解結(jié)合LR算法的 方式獲取上三角矩陣R。 實施例四在Nt根發(fā)射天線、Nr根接收天線的MIMO系統(tǒng)中，圖5為本發(fā)明實施 例四中串行干擾消除的總體流程圖。如圖5所示，該方法包括步驟501,獲取發(fā)送端到接收端的信道傳輸矩陣丑。本實施例中，在接收端進行信道估計，從而獲取發(fā)送端到接收端的信道 傳輸矩陣丑。接下來，通過步驟502 506對信道傳輸矩陣^進行QR分解，其中， 設(shè)該QR分解的矩陣的維數(shù)為M x 7V。本實施例中，M=iVr， 7V=A/"t。步驟502,對步驟501中獲取的信道傳輸矩陣丑進行LR計算，得到矩 陣及r和A,使」^ =丑^，并將好r作為當(dāng)前信道傳輸矩陣。步驟503,將A^xA/"r維的加權(quán)矩陣g和iVtx ATt維的第一排序矩陣p初 始化為單位矩陣，并將當(dāng)前信道傳輸矩陣及r的第一列作為當(dāng)前列，即A:-1， A:表示當(dāng)前列的列號，將第二排序矩陣尸初始化為矩陣i>r。步驟504,從當(dāng)前信道傳輸矩陣的當(dāng)前列到最后一列，在每列中提取第 A個元素到最后一個元素構(gòu)成該列的重構(gòu)列向量，計算每列的重構(gòu)列向量的 范數(shù)；并確定計算結(jié)果中范數(shù)最小的重構(gòu)列向量，確定該重構(gòu)列向量所在列 的列號《。步驟505，對第《列的重構(gòu)列向量進行householder變換，向量v和標(biāo)量 ";才艮據(jù)向量v和標(biāo)量yS構(gòu)造MxM維的歸一化正交矩陣込。上述步驟504- 505的操作與實施例一中步驟103 ~ 104的操作相同，這 里就不再贅述。步驟506，將矩陣込與當(dāng)前信道傳輸矩陣的乘積作為更新后的信道傳輸 矩陣，將矩陣込與當(dāng)前加權(quán)矩陣g的乘積作為更新后的加權(quán)矩陣；將更新 后的信道傳輸矩陣的第A列和第《列交換，將當(dāng)前第 一 排序矩陣/ 設(shè)置為單 位矩陣，并交換該單位矩陣的第A:列和第《列。步驟507,對當(dāng)前信道傳輸矩陣的分矩陣進行LR計算，生成矩陣" 和i>t,根據(jù)i>t構(gòu)造7Vt x 7Vt維的幺模矩陣込'。本步驟中，當(dāng)前信道傳輸矩陣的分矩陣為當(dāng)前信道傳輸矩陣的第A+l 行到第M行中的第A+l個元素到第W個元素構(gòu)成的矩陣，即1 :M，A + 1 :iV)。對上述形成的分矩陣進行LR計算生成矩陣及t和Pt。根據(jù)尸t構(gòu)造A/"rX7Vr維的歸一化正交矩陣g/為將kxk維的單位矩陣 作為左上角分矩陣，將A作為右下角分矩陣，其余位置填0,構(gòu)成A^xATt 維的幺模矩陣込'，即込=[/*,^(糾);0 ， A]。步驟508，將當(dāng)前第二排序矩陣、當(dāng)前第一排序矩陣和矩陣g/的乘積作為更新后的第二排序矩陣，將當(dāng)前信道傳輸矩陣與矩陣2/的乘積作為更新后的信道傳輸矩陣。本步驟中，更新第二排序矩陣，具體將當(dāng)前第二排序矩陣、當(dāng)前第一排 序矩陣和矩陣g/的乘積作為更新后的第二排序矩陣，即,=/>'*/>*込；更新信道傳輸矩陣，具體將當(dāng)前信道傳輸矩陣與矩陣込'的乘積作為更新后的 信道傳輸矩陣，即丑=^ *込。步驟509，將當(dāng)前信道傳輸矩陣的下一列作為當(dāng)前列，即&=&+1，并判 斷k是否等于當(dāng)前信道傳輸矩陣的列數(shù)，若是，則執(zhí)行步驟510,否則返回 步驟504。步驟510,將當(dāng)前加權(quán)矩陣、步驟501中獲取的信道傳輸矩陣丑、當(dāng)前 第二排序矩陣尸'的乘積作為上三角矩陣及。經(jīng)過上述步驟501 ~ 510的操作對信道傳輸矩陣H進行排序QR分解后， 得到加權(quán)矩陣g和第二排序矩陣尸'，根據(jù)上述信道傳輸矩陣和QR分解后的 結(jié)果計算上三角矩陣及的方式為及=0及尸 (6)。由于矩陣及是個上三角矩陣，因此可以利用該矩陣對檢測信號進行串行干擾消除和信號檢測。在本實施例中，采用ZF檢測器，通過上三角矩陣進行串行干擾消除及信號檢測。具體執(zhí)行如下步驟步驟511,在接收端，對接收信號進行濾波得到向量0r，并根據(jù)0r和及計算向量CP')"rf。由公式(6 )可得及二gH及CP')-1，將其代入r=^/+"可得，r^gH及CP')"iZ+",將該式左右同時左乘g可得，0^及CP')-V+gn。在該等式關(guān)系中，矩陣及為上三角矩陣，在ZF算法中噪聲項gw可以忽略，因此可以利用該等式計算向量CP')'V的取值。步驟512,將向量(尸')-V左乘第二排序矩陣尸得到發(fā)射信號</。經(jīng)過上述步驟511 ~ 512即完成了接收端的信號檢測。至此，本實施例中進行串行干擾消除的方法流程結(jié)束。本實施例中，通過排序QR分解與LR算法相結(jié)合的方式確定上三角矩陣R,并且根據(jù)上三角矩陣與發(fā)射信號和接收信號的關(guān)系，進行發(fā)射信號的檢測和串行干擾消除。為說明本實施例的優(yōu)點，對應(yīng)用上述實施例四的系統(tǒng)進行了仿真實驗。 圖6為應(yīng)用實施例四的方法、背景技術(shù)中介紹的ZF算法、ZF算法與LR算 法結(jié)合時，進行信號檢測的性能比較示意圖。該系統(tǒng)仿真的仿真參數(shù)與圖3 所示仿真結(jié)果的仿真參數(shù)相同。如圖6所示，橫軸表示接收端檢測信號的信 噪比，縱軸表示接收端檢測信號的誤碼率。其中，曲線601是采用背景技術(shù)中的無排序的ZF算法進行信號檢測時 信噪比與誤碼率的關(guān)系曲線；曲線602是采用本發(fā)明排序的QR分解的ZF 算法進行信號檢測時信噪比與誤碼率的關(guān)系曲線；曲線603是應(yīng)用實施例四 QR與LR算法結(jié)合的方式進行信號檢測時信噪比與誤碼率的關(guān)系曲線。由曲線601、 602和603的比較可見，應(yīng)用LR算法后，能夠帶來更優(yōu) 的信號檢測性能，而應(yīng)用本發(fā)明實施例四的方式進行信號檢測，能夠獲得最 佳的信號檢測性能。以上僅為本發(fā)明的較佳實施例而已，并非用于限定本發(fā)明的保護范圍。 凡在本發(fā)明的精神和原則之內(nèi)，所作的任何修改、等同替換、改進等，均應(yīng) 包含在本發(fā)明的保護范圍之內(nèi)。
權(quán)利要求
1、一種多輸入多輸出系統(tǒng)中的串行干擾消除方法，其特征在于，該方法包括a、獲取發(fā)送端到接收端的信道傳輸矩陣H；b、根據(jù)步驟a中所述信道傳輸矩陣H確定當(dāng)前信道傳輸矩陣，設(shè)置加權(quán)矩陣Q和第一排序矩陣P分別為M×M維和N×N維的單位矩陣，所述M和N分別為所述當(dāng)前信道傳輸矩陣的行數(shù)和列數(shù)，并將當(dāng)前信道傳輸矩陣的第一列作為當(dāng)前列，即k＝1，k表示當(dāng)前列的列號；c、從當(dāng)前信道傳輸矩陣的當(dāng)前列到最后一列，在每列中提取第k個元素到最后一個元素構(gòu)成該列的重構(gòu)列向量，計算每列的重構(gòu)列向量的范數(shù)，確定范數(shù)最小的重構(gòu)列向量所在列的列號q，并對該列的重構(gòu)列向量進行豪思霍德householder變換，生成向量v和標(biāo)量β；根據(jù)向量v和標(biāo)量β構(gòu)造歸一化正交矩陣Qt；d、將所述Qt與當(dāng)前信道傳輸矩陣的乘積作為更新后的信道傳輸矩陣，將矩陣Qt與當(dāng)前加權(quán)矩陣Q的乘積作為更新后的加權(quán)矩陣；將更新后的信道傳輸矩陣的第k列和第q列交換，將當(dāng)前第一排序矩陣的第k列和第q列交換；e、將當(dāng)前信道傳輸矩陣的下一列作為當(dāng)前列，即k＝k+1，并判斷k是否等于當(dāng)前信道傳輸矩陣的行數(shù)和列數(shù)中的最小值，若是，則執(zhí)行步驟f，否則返回步驟c；f、根據(jù)當(dāng)前加權(quán)矩陣和步驟b中所述當(dāng)前信道傳輸矩陣確定上三角矩陣R，根據(jù)矩陣R對各個信號進行串行干擾消除。
2、 根據(jù)權(quán)利要求1所述的方法，其特征在于，所述根據(jù)當(dāng)前加權(quán)矩陣和步驟b中所述當(dāng)前信道傳輸矩陣確定上三角矩陣 及為將當(dāng)前加權(quán)矩陣與步驟b中獲取的所述信道傳輸矩陣以及當(dāng)前第一排序 矩陣的共軛轉(zhuǎn)置矩陣之乘積作為矩陣i 。
3、 根據(jù)權(quán)利要求2所述的方法，其特征在于，步驟b中所述才艮據(jù)步驟a中的信道傳輸矩陣^確定當(dāng)前信道傳輸矩陣為 將步驟a中所述信道傳輸矩陣丑作為當(dāng)前信道傳輸矩陣； 所述根據(jù)矩陣及對各個信號進行串行干擾消除為在接收端，計算當(dāng)前加權(quán)矩陣與接收信號的乘積，利用該乘積結(jié)果和矩陣 及，對接收信號進行串行干擾消除，檢測發(fā)射信號。
4、 根據(jù)權(quán)利要求2所述的方法，其特征在于，步驟b中所述才艮據(jù)信道傳輸矩陣W確定當(dāng)前信道傳輸矩陣包括才艮據(jù)信道 傳輸矩陣丑構(gòu)造擴展信道矩陣及e-[及T，a2/]1，其中，02是系統(tǒng)噪聲的均方差， /是A^XA/",的單位矩陣；將擴展信道矩陣作為當(dāng)前信道傳輸矩陣；所述根據(jù)矩陣W對各個發(fā)射信號進行串行干擾消除為在接收端，對接收 到的信號進行擴展，得到擴展后的接收信號為re=[/,diT]T,建立擴展后的接 收信號、當(dāng)前加權(quán)矩陣0、矩陣及、第一排序矩陣P和發(fā)射信號rf的函數(shù)關(guān)系 為0^J iW-ff0i 7^iV,其中，及'=[~威。/]7， /是A^X7V",的單位矩陣；根據(jù) 矩陣及和所述函數(shù)關(guān)系，檢測發(fā)射信號d,并進行串行干擾消除。
5、 根據(jù)權(quán)利要求2所述的方法，其特征在于，步驟b中所述根據(jù)信道傳輸矩陣丑確定當(dāng)前信道傳輸矩陣包括將步驟a 中獲取的信道傳輸矩陣丑的共軛轉(zhuǎn)置矩陣^作為當(dāng)前信道傳輸矩陣；所述根據(jù)矩陣及對各個信號進行串行預(yù)干擾消除為在發(fā)送端，計算當(dāng)前第一排序矩陣戶的共軛轉(zhuǎn)置矩陣與矩陣及的共軛轉(zhuǎn)置 矩陣的乘積i^及h,從選擇的第一發(fā)射信號開始，根據(jù)/^及h依次計算并消除其 它發(fā)射信號對該發(fā)射信號的干擾，并將消除干擾后的發(fā)射信號進行發(fā)射。
6、 根據(jù)權(quán)利要求5所述的方法，其特征在于，所述根據(jù)P"及h依次消除其它發(fā)射信號對該發(fā)射信號的干擾為J=d-5T&A, ijmod",力，其中，爿為取模除數(shù)，并由選擇的調(diào)制方式?jīng)Q定，&為/^及h中的第k行第i個元素，《為第i發(fā)射信號，4為消除干擾后的第i發(fā)射信號。
7、 根據(jù)權(quán)利要求1所述的方法，其特征在于，當(dāng)發(fā)射天線數(shù)小于等于接收天線數(shù)時，步驟b中所述根據(jù)信道傳輸矩陣丑 確定當(dāng)前信道傳輸矩陣包括對步驟a中獲取的信道傳輸矩陣^進行減格LR 算法計算，得到矩陣fir和A,使好^ =丑^,并將Hr作為當(dāng)前信道傳輸矩陣；在步驟b中進一步包括將第二排序矩陣戶'初始化為矩陣iV，在步驟d中執(zhí)行所述將當(dāng)前第一排序矩陣的第A列和第《列交換的才喿作前， 進一步包括將當(dāng)前第一排序矩陣P設(shè)置為單位矩陣；在步驟d和e之間進一步包括對當(dāng)前信道傳輸矩陣的分矩陣進行LR計 算，生成矩陣試和尸t，根據(jù)A構(gòu)造A^X7Vt維的幺模矩陣込'，將當(dāng)前第二排序 矩陣、當(dāng)前第一排序矩陣和矩陣込'的乘積作為更新后的第二排序矩陣，將當(dāng)前 信道傳輸矩陣與矩陣込'的乘積作為更新后的信道傳輸矩陣；步驟f中所述根據(jù)當(dāng)前加權(quán)矩陣和步驟b中所述當(dāng)前信道傳輸矩陣確定上 三角矩陣及為將當(dāng)前加權(quán)矩陣、步驟a中獲取的所述信道傳輸矩陣及、當(dāng)前 第二排序矩陣/>'的乘積作為上三角矩陣及。
8、 根據(jù)權(quán)利要求7所述的方法，其特征在于，所述當(dāng)前信道傳輸矩陣的分 矩陣為當(dāng)前信道傳輸矩陣的第壯l行到第M行中的第A:+l個元素到第7V個元 素構(gòu)成的矩陣。
9、 根據(jù)權(quán)利要求7所述的方法，其特征在于，所述根據(jù)戶t構(gòu)造7VtXA/t維 的幺模矩陣込'為將當(dāng)前信道傳輸矩陣的第^+l行到第M行中的第fc+l個元 素到第W個元素應(yīng)用LR減格算法構(gòu)成矩陣込'。
10、 根據(jù)權(quán)利要求7所述的方法，其特征在于，所述根據(jù)矩陣及對各個發(fā) 射信號進行串行干擾消除為在接收端，建立接收信號r、當(dāng)前第二排序矩陣F、當(dāng)前加權(quán)矩陣g、上 三角矩陣及和發(fā)射信號rf的函數(shù)關(guān)系為^ =及(尸)-^+01，根據(jù)所述函數(shù)關(guān)系， 利用迫零ZF或MMSE算法檢測發(fā)射信號</,并進行串行干擾消除。
11、 根據(jù)權(quán)利要求l到IO中任一所述的方法，其特征在于，所述根據(jù)向量 v和標(biāo)量"構(gòu)造歸一化正交矩陣込為在向量v的上面添加(^-l )個零，構(gòu)成 向量v'，即v'=[^_1)xl,vT]T，再計算/3^V w,將計算結(jié)果作為矩陣込。
全文摘要
本發(fā)明公開了一種多輸入多輸出系統(tǒng)中的串行干擾消除方法，包括基于householder變換對信道傳輸矩陣進行QR排序分解獲得上三角矩陣R，根據(jù)矩陣R對各個發(fā)射信號進行預(yù)編碼的串行干擾消除或?qū)邮招盘栠M行串行干擾消除。由于該矩陣R是一個上三角矩陣，因此可以利用上三角矩陣特性依次對各個發(fā)射信號進行串行干擾消除，從而避免利用VBLAST對矩陣求逆或偽逆的方式確定串行檢測順序，進而簡化串行干擾消除的計算復(fù)雜度，提高信號檢測的性能。
文檔編號H04L1/06GK101409604SQ20071017564
公開日2009年4月15日 申請日期2007年10月9日 優(yōu)先權(quán)日2007年10月9日
發(fā)明者拉蓋施, 王大飛, 強 薛 申請人:鼎橋通信技術(shù)有限公司