專利名稱:多天線無線通信系統(tǒng)、發(fā)送/接收處理單元及其聯(lián)合處理方法
技術領域:
本發(fā)明涉及一種多天線無線通信系統(tǒng)(包括單/多載波的多天線系統(tǒng))、發(fā)送/接收處理單元及其聯(lián)合處理方法。
背景技術:
隨著無線網(wǎng)絡、多媒體技術和因特網(wǎng)的逐漸融合,人們對無線通信業(yè)務的類型和質量的要求越來越高。為滿足無線多媒體和高速率數(shù)據(jù)傳輸?shù)囊?,需要開發(fā)新一代無線通信系統(tǒng)。其中多天線輸入和輸出(MIMO)和正交頻分復用(OFDM)相結合的MIMO OFDM技術受到廣泛關注。
MIMO和OFDM相結合的MIMO OFDM系統(tǒng)具有兩者的優(yōu)點,它既通過OFDM調制把頻率選擇性MIMO衰落信道分解成一組并行平坦衰落信道,又利用MIMO提高了系統(tǒng)容量。
采用MIMO技術能獲得空間分集或復用增益,在無需耗費額外功率和帶寬的條件下大大增加了系統(tǒng)容量,同時能顯著提高傳輸鏈路質量,適用于傳輸高速率的音、視頻等多媒體業(yè)務。目前,一些高速無線接入技術標準,如IEEE 802.11n、802.16等均采用了MIMO和OFDM相結合的技術方式。一般的MIMO OFDM系統(tǒng)中,為獲得最優(yōu)系統(tǒng)容量,采用特征值分解方法把信道矩陣分解為多個并行獨立的子信道,然后利用注水算法對每個子信道進行功率分配。但功率分配要求射頻功放器件的動態(tài)范圍較寬,以便能連續(xù)地分配功率,故實現(xiàn)困難。另一方面,為提高系統(tǒng)的誤比特率性能,系統(tǒng)采用了自適應調制,且當對每個天線支路采用了相同的調制方式時,該方法存在的缺點是,由于各個子信道增益不一樣,低信道增益的子信道將會影響系統(tǒng)的誤比特率。因此,本發(fā)明提出一種設計方法,它能獲得高的系統(tǒng)容量,同時降低系統(tǒng)的誤比特率。
參考文獻列表[1]A.J.Paularj等人發(fā)表在Proceedings of IEEE、2004年2月、第92卷、第2號、第198~218頁上的題為“An Overview of MIMOCommunications-A Key to Gigabit Wireless”的文章。
Y.Jiang、W.W.Hager和J.Li于2003年12月公開發(fā)表的題為“The Geometric Mean Decomposition,Linear Algebra and ItsApplications”的文章。
B.Hassibi于2000年11月發(fā)表在第34屆Asilomar Conf.第1255~1259頁上的題為“A Fast Square-Root Implementation forBLAST”的文章。
發(fā)明內容
本發(fā)明目的在于提供一種用于多天線通信系統(tǒng)的方法和裝置,使用該方法及其裝置能降低系統(tǒng)的誤比特率,同時提高容量。
本發(fā)明的思想是通過對信道矩陣進行幾何平均分解,使得天線各個子信道的增益相同。這樣,當每個天線支路采用相同的調制方式時,系統(tǒng)的誤比特率性能不會因某一子信道存在深度衰落而降低。且系統(tǒng)采用閉環(huán)的發(fā)射、接收聯(lián)合處理,降低了系統(tǒng)容量損失。
為了實現(xiàn)上述目的和思想,根據(jù)本發(fā)明的第一方面,提出了一種用在多天線無線通信系統(tǒng)中的聯(lián)合處理方法,包括在接收端,對信道矩陣進行幾何平均分解,得到線性變換矩陣Q、上三角矩陣R和預濾波矩陣P,其維數(shù)均為K×K,將預濾波矩陣P和線性變換矩陣Q反饋到發(fā)送端,利用線性變換矩陣Q的Hermitian轉置矩陣QH對接收信號進行線性變換;在發(fā)送端,利用預濾波陣P對信號進行線性預編碼,利用線性變換矩陣Q的Hermitian轉置矩陣QH對信號進行線性變換;以及在接收端,對串行干擾刪除信號進行檢測,恢復與各個子信道對應的天線端的發(fā)射信號。
優(yōu)選地,所述串行干擾刪除信號檢測步驟包括從第k個子信道開始,用該子信道接收信號減去第(k+1)個子信道的接收信號;然后,再除以上三角矩陣R的對角元rkk,恢復與第k個子信道對應的天線端的發(fā)射信號;更新k=k-1;并判斷k是否小于1,如果k<1,則結束,若k≥1,則針對更新后的k繼續(xù)進行處理。
優(yōu)選地,所述聯(lián)合處理方法還可以包括在發(fā)送端,發(fā)射導頻信號到信道;以及在接收端,估計信道矩陣的增益系數(shù)。
根據(jù)本發(fā)明的第二方面,提出了一種用在多天線無線通信系統(tǒng)中的聯(lián)合處理方法,包括發(fā)送端發(fā)射導頻信號到信道,接收端估計信道矩陣的增益系數(shù);對信道矩陣進行幾何平均分解,得到線性變換矩陣Q、上三角矩陣R和預濾波矩陣P,其維數(shù)均為K×K;將預濾波陣P和線性變換矩陣Q反饋到發(fā)送端;發(fā)送端利用預濾波陣P對信號進行線性預編碼;發(fā)送端利用線性變換矩陣Q的Hermitian轉置矩陣QH對信號進行線性變換,同時,接收端也利用線性變換矩陣Q的Hermitian轉置矩陣QH對接收信號進行線性變換;對串行干擾刪除信號進行檢測從第k個子信道開始,用該子信道接收信號減去第(k+1)個子信道接收信號后,再除以R的對角元rkk,恢復與第k個子信道對應的天線端的發(fā)射信號;以及更新k=k-1,并判斷k是否小于1,如果k<1,則結束;若k≥1,則轉上述串行干擾刪除信號檢測步驟繼續(xù)處理。
根據(jù)本發(fā)明的第三方面,提出了一種在多天線無線通信系統(tǒng)中使用的發(fā)送處理單元,包括預編碼模塊,利用從接收端發(fā)送過來的預濾波陣P對信號進行線性預編碼;以及線性變換模塊,利用從接收端發(fā)送過來的線性變換矩陣Q的Hermitian轉置矩陣QH對信號進行線性變換。
根據(jù)本發(fā)明的第四方面,提出了一種在多天線無線通信系統(tǒng)中使用的接收處理單元,包括信道分解模塊,對信道矩陣進行幾何平均分解,得到線性變換矩陣Q、上三角矩陣R和預濾波矩陣P,其維數(shù)均為K×K;線性變換模塊,利用線性變換矩陣Q的Hermitian轉置矩陣QH對接收信號進行線性變換;以及迭代處理模塊,對串行干擾刪除信號進行檢測,恢復與各個子信道對應的天線端的發(fā)射信號。
優(yōu)選地,所述迭代處理模塊從第k個子信道開始,用該子信道接收信號減去第(k+1)個子信道的接收信號;然后,再除以上三角矩陣R的對角元rkk,恢復與第k個子信道對應的天線端的發(fā)射信號;更新k=k-1;并判斷k是否小于1,如果k<1,則結束,若k≥1,則針對更新后的k繼續(xù)進行處理。
根據(jù)本發(fā)明的第五方面,提出了一種多天線無線通信系統(tǒng)及其操作方法,其特征在于包括根據(jù)本發(fā)明的發(fā)送處理單元和接收處理單元。
同時,本發(fā)明提出的聯(lián)合處理方法可用于所有的多天線通信系統(tǒng),如單載波的MIMO CDMA和多載波的MIMO OFDM等。
下面,將參照附圖,對本發(fā)明的優(yōu)選實施例進行詳細的描述,其中圖1示出了根據(jù)本發(fā)明的MIMO系統(tǒng)發(fā)送端的示意圖;圖2示出了根據(jù)本發(fā)明的MIMO系統(tǒng)接收端的示意圖;圖3示出了根據(jù)本發(fā)明的聯(lián)合處理方法的流程圖;以及圖4示出了根據(jù)本發(fā)明的MIMO系統(tǒng)的操作性能的仿真結果。
具體實施例方式
下面結合附圖對本發(fā)明作具體說明。應該指出,所描述的實施例僅是為了說明的目的,而不是對本發(fā)明范圍的限制。
考慮一個單用戶MIMO OFDM通信鏈路,其發(fā)送和接收端分別如圖1和2所示,發(fā)送端和接收端分別安裝NT幅發(fā)送天線101~101’和NR幅接收天線106~106’。在發(fā)送端,信息符號流輸入到復用器103,復用器103以空分復用方式把符號流送入每個發(fā)送天線支路端;在每支路端,逆快速傅里葉變換(IFFT)模塊102~102’對符號流進行OFDM調制;調制后的符號從每幅天線101~101’發(fā)送到信道。在接收端,每幅天線106~106’接收到的信號送入快速傅里葉變換(FFT)模塊107~107’,在FFT模塊107~107’內完成OFDM解調;信道估計模塊108利用每幅天線接收到的導頻序列估計每條傳播路徑的信道增益;解復用器109把各天線符號轉換為串行輸出符號流。發(fā)送、接收端的聯(lián)合處理模塊(104~105、110~112),利用信道估計值,進行一系列發(fā)射、接收處理,其中在發(fā)送端完成預處理操作,在接收端完成后處理操作。
下面分析發(fā)送、接收聯(lián)合處理的原理。對于一般的MIMO系統(tǒng),接收信號為Y=HX+n(1)式中,Y為NR×1維接收信號向量,發(fā)射信號向量X∈CNT×1]]>的方差為σx2I,白高斯噪聲向量n的方差為σn2I,故信噪比SNR=σx2/σn2.]]>信道矩陣H∈CNR×NT]]>的秩為K,H中的元素hji為發(fā)射天線i到接收天線j的信道衰落系數(shù)。
對于空分復用MIMO系統(tǒng),當發(fā)送端未知信道狀態(tài)信息(CSI)時,一般的接收信號檢測方法是串行干擾刪除算法。例如,從第NR幅天線開始,利用線性均衡(如匹配濾波、迫零(ZF)均衡、最小均方誤差(MMSE)均衡)等方法恢復該幅天線接收信號;從接收信號中減去第NR幅天線接收信號后,繼續(xù)從第NR-1幅天線開始,恢復該幅天線信號;依次進行處理,直到恢復最后一幅天線的接收信號。該處理方法屬于開環(huán)工作模式,存在一定的系統(tǒng)容量損失。當接收端獲得信道增益信息并反饋到發(fā)送端后,即收、發(fā)端均獲得信道狀態(tài)信息,則系統(tǒng)工作在閉環(huán)模式下。這時,可采用如下等效處理方法。
對信道矩陣H進行QR分解得到,H=QNT×KRK×K---(2)]]>式中,QNT×K為列正交陣,RK×K為上三角陣。故式(1)重寫為Y=QRX+n(3)式(3)兩邊同乘以QH,(·)H表示Hermitian轉置(復共軛),得到Y=RX+n(4)展開即為
根據(jù)式(5)進行串行信號檢測,即對所有的K個子信道開始,計算x‾i←(y‾i-Σj=i+1Krijx‾j)/rii,]]>并把它映射為解調符號。
可以看到,通過對信道矩陣H的QR分解,MIMO信道被分解為K個并行子信道,即yi=riixi+ni,i=1,2,…,K。(6)與矩陣的SVD分解相同,信道矩陣QR分解后的子信道仍為呈現(xiàn)隨機衰落特征,如前所述,深衰落子信道將降低系統(tǒng)整體誤比特率性能。為此,在本發(fā)明中,采用參考文獻[2]所提出的幾何平均分解(GMD)方法,對信道傳輸矩陣進行幾何平均分解,H=QNR×KRK×KPNT×KH---(7)]]>這里,半酉矩陣Q、P分別為在發(fā)射、接收端進行預濾波和后濾波處理;上三角陣R具有相等的對角元rii=λ,λ定義為λ‾=(Πk=1Kλk)1/K---(8)]]>其中,λk為對信道矩陣H進行特征值分解(SVD)得到的對角陣diag[γ1…γK]的非零對角元,即H的非零特征值??梢?,通過GMD分解,得到的矩陣R的對角元為H的特征值的幾何平均值,使各子信道增益相同,故可降低由于各子信道衰落不同對系統(tǒng)整體性能的影響。GMD的實現(xiàn)可通過對稱置換和一對Givens旋轉操作來完成。
當引入線性預編碼矩陣P對發(fā)射信號預處理后,接收信號為Y=HPX+n (9)把式(7)代入式(9)得到Y=QRPHPX+n~=QRX+n~---(10)]]>兩邊同乘QH得到與式(4)相同的表示式,Y=RX+n,由此,可利用前述的干擾刪除方法來恢復信號。
根據(jù)以上分析,參考圖3,對根據(jù)本發(fā)明的系統(tǒng)聯(lián)合處理方法進行詳細的描述。
在步驟S001,發(fā)射導頻信號到信道,接收端估計信道矩陣的增益系數(shù);在步驟S002,對信道矩陣進行幾何平均分解,得到矩陣Q、R和P,其維數(shù)均為K×K;在步驟S003,把預濾波陣P和Q反饋到發(fā)送端;在步驟S004,發(fā)送端利用預濾波陣P對信號進行線性預編碼;在步驟S005,發(fā)送端利用QH對信號作線性變換,同時,接收端也利用QH對接收信號作處理;在步驟S006,串行干擾刪除信號檢測從第k個子信道開始,用該子信道接收信號減去第(k+1)個子信道接收信號后,再除以rkk(R的對角元),恢復與第k個子信道對應的天線端的發(fā)射信號;在步驟S007,更新k=k-1,并判斷k是否小于1,如果k<1,則結束;若k≥1,則轉到步驟S006進行循環(huán)處理。
具體地講,MIMO OFDM系統(tǒng)中的聯(lián)合處理裝置用以實現(xiàn)上述步驟,發(fā)射、接收端聯(lián)合處理裝置如圖1、2中虛線框內各個模塊所示,包括以下模塊信道分解模塊110,對信道矩陣進行幾何平均分解;預編碼模塊104,利用矩陣分解得到的預編碼器P對發(fā)射信號進行線性預編碼;線性變換模塊105和111,分別利用矩陣分解得到的矩陣QH在發(fā)送、接收端對信號進行線性變換;迭代處理模塊112,執(zhí)行信號處理恢復信號的算法,并判斷是否已經(jīng)完成對所有子信道的信號恢復。
結合圖1、圖2和圖3,對本發(fā)明所提出的MIMO OFDM系統(tǒng)的操作描述如下。
復用器103以空分復用方式把符號流送入每個發(fā)送天線支路端。在每個支路端,IFFT模塊102~102’對符號流進行OFDM調制。調制后的符號從每幅天線101~101’發(fā)送到信道。發(fā)送端發(fā)射導頻信號到信道,接收端估計信道矩陣的增益系數(shù)。每幅天線106~106’接收到的信號送入FFT模塊107~107’。在FFT模塊107~107’內完成OFDM解調。信道估計模塊108計算信道矩陣增益系數(shù)。解調后的符號送入解復用器109恢復為串行輸出符號流。信道分解模塊110對信道矩陣進行幾何平均分解得到的矩陣Q、R和P,其維數(shù)均為K×K。接收端把預濾波陣P和Q反饋到發(fā)送端。預編碼模塊104利用預濾波陣P對信號進行線性預編碼。線性變換模塊105利用QH對信號進行線性變換,同時,線性變換模塊111也利用QH對接收信號進行線性變換。迭代處理模塊112進行串行干擾刪除信號的檢測從第k個子信道開始,用該子信道接收信號減去第(k+1)個子信道接收信號后,再除以rkk(R的對角元),恢復與第k個子信道對應的天線端的發(fā)射信號,更新k=k-1,并判斷k是否小于1,如果k<1,則處理結束,若k≥1,則針對更新后的k繼續(xù)進行處理。
圖4示出了根據(jù)本發(fā)明的MIMO系統(tǒng)的操作性能的仿真結果。仿真參數(shù)為發(fā)送、接收天線數(shù)均為4;OFDM子載波數(shù)64,其循環(huán)前綴16;BPSK調制;未編碼;平坦瑞利衰落信道。系統(tǒng)的誤比特率(BER)性能如圖4所示。
本發(fā)明中,通過對信道矩陣進行幾何平均分解,使得天線的各個子信道增益相同。這樣,當每個天線支路采用相同的調制方式時,系統(tǒng)的誤比特率性能不會因某一子信道存在深度衰落而受影響。且系統(tǒng)采用閉環(huán)的發(fā)射、接收的聯(lián)合處理,相對開環(huán)多天線處理,降低了系統(tǒng)容量損失。因此,本發(fā)明提出的設計方法能同時獲得系統(tǒng)容量和誤比特率性能的增益。
盡管已經(jīng)針對典型實施例示出和描述了本發(fā)明,本領域的普通技術人員應該理解,在不脫離本發(fā)明的精神和范圍的情況下,可以進行各種其他的改變、替換和添加。因此,本發(fā)明不應該被理解為被局限于上述特定實例,而應當由所附權利要求所限定。
權利要求
1.一種用在多天線無線通信系統(tǒng)中的聯(lián)合處理方法,包括在接收端,對信道矩陣進行幾何平均分解,得到線性變換矩陣Q、上三角矩陣R和預濾波矩陣P,其維數(shù)均為K×K,將預濾波矩陣P和線性變換矩陣Q反饋到發(fā)送端,利用線性變換矩陣Q的Hermitian轉置矩陣QH對接收信號進行線性變換;在發(fā)送端,利用預濾波陣P對信號進行線性預編碼,利用線性變換矩陣Q的Hermitian轉置矩陣QH對信號進行線性變換;以及在接收端,對串行干擾刪除信號進行檢測,恢復與各個子信道對應的天線端的發(fā)射信號。
2.根據(jù)權利要求1所述的聯(lián)合處理方法,其特征在于所述串行干擾刪除信號檢測步驟包括從第k個子信道開始,用該子信道接收信號減去第(k+1)個子信道的接收信號;然后,再除以上三角矩陣R的對角元rkk,恢復與第k個子信道對應的天線端的發(fā)射信號;更新k=k-1;并判斷k是否小于1,如果k<1,則結束,若k≥1,則針對更新后的k繼續(xù)進行處理。
3.根據(jù)權利要求1或2所述的聯(lián)合處理方法,其特征在于還包括在發(fā)送端,發(fā)射導頻信號到信道;以及在接收端,估計信道矩陣的增益系數(shù)。
4.一種用在多天線無線通信系統(tǒng)中的聯(lián)合處理方法,包括發(fā)送端發(fā)射導頻信號到信道,接收端估計信道矩陣的增益系數(shù);對信道矩陣進行幾何平均分解,得到線性變換矩陣Q、上三角矩陣R和預濾波矩陣P,其維數(shù)均為K×K;將預濾波陣P和線性變換矩陣Q反饋到發(fā)送端;發(fā)送端利用預濾波陣P對信號進行線性預編碼;發(fā)送端利用線性變換矩陣Q的Hermitian轉置矩陣QH對信號進行線性變換,同時,接收端也利用線性變換矩陣Q的Hermitian轉置矩陣QH對接收信號進行線性變換;對串行干擾刪除信號進行檢測,從第k個子信道開始,用該子信道接收信號減去第(k+1)個子信道接收信號后,再除以R的對角元rkk,恢復與第k個子信道對應的天線端的發(fā)射信號;以及更新k=k-1,并判斷k是否小于1,如果k<1,則結束;若k≥1,則轉上述串行干擾刪除信號檢測步驟繼續(xù)處理。
5.一種在多天線無線通信系統(tǒng)中使用的發(fā)送處理單元,包括預編碼模塊(104),利用從接收端發(fā)送過來的預濾波陣P對信號進行線性預編碼;以及線性變換模塊(105),利用從接收端發(fā)送過來的線性變換矩陣Q的Hermitian轉置矩陣QH對信號進行線性變換。
6.根據(jù)權利要求5所述的發(fā)送處理單元,其特征在于所述發(fā)送處理單元應用于多天線輸入和輸出碼分多址系統(tǒng)。
7.根據(jù)權利要求5所述的發(fā)送處理單元,其特征在于所述發(fā)送處理單元應用于多天線輸入和輸出正交頻分復用系統(tǒng)。
8.根據(jù)權利要求5到7之一所述的發(fā)送處理單元,其特征在于所述發(fā)送處理單元位于基站或移動終端中。
9.一種在多天線無線通信系統(tǒng)中使用的接收處理單元,包括信道分解模塊(110),對信道矩陣進行幾何平均分解,得到線性變換矩陣Q、上三角矩陣R和預濾波矩陣P,其維數(shù)均為K×K;線性變換模塊(111),利用線性變換矩陣Q的Hermitian轉置矩陣QH對接收信號進行線性變換;以及迭代處理模塊(112),對串行干擾刪除信號進行檢測,恢復與各個子信道對應的天線端的發(fā)射信號。
10.根據(jù)權利要求9所述的接收處理單元,其特征在于所述迭代處理模塊(112)從第k個子信道開始,用該子信道接收信號減去第(k+1)個子信道的接收信號;然后,再除以上三角矩陣R的對角元rkk,恢復與第k個子信道對應的天線端的發(fā)射信號;更新k=k-1;并判斷k是否小于1,如果k<1,則結束,若k≥1,則針對更新后的k繼續(xù)進行處理。
11.根據(jù)權利要求9或10所述的接收處理單元,其特征在于所述接收處理單元應用于多天線輸入和輸出碼分多址系統(tǒng)。
12.根據(jù)權利要求9或10所述的接收處理單元,其特征在于所述接收處理單元應用于多天線輸入和輸出正交頻分復用系統(tǒng)。
13.根據(jù)權利要求9到12之一所述的接收處理單元,其特征在于所述接收處理單元位于移動終端或基站中。
14.一種多天線無線通信系統(tǒng),其特征在于包括根據(jù)權利要求5所述的發(fā)送處理單元和根據(jù)權利要求9所述的接收處理單元。
15.根據(jù)權利要求14所述的多天線無線通信系統(tǒng)的操作方法,其特征在于包括步驟信道分解模塊(110)對信道矩陣進行幾何平均分解,得到線性變換矩陣Q、上三角矩陣R和預濾波矩陣P,其維數(shù)均為K×K;預編碼模塊(104)利用從接收端發(fā)送過來的預濾波陣P對信號進行線性預編碼;線性變換模塊(105)利用從接收端發(fā)送過來的線性變換矩陣Q的Hermitian轉置矩陣QH對信號進行線性變換;線性變換模塊(111)利用線性變換矩陣Q的Hermitian轉置矩陣QH對接收信號進行線性變換;以及迭代處理模塊(112)對串行干擾刪除信號進行檢測,恢復與各個子信道對應的天線端的發(fā)射信號。
16.一種多天線無線通信系統(tǒng),其特征在于包括根據(jù)權利要求5所述的發(fā)送處理單元和根據(jù)權利要求10所述的接收處理單元。
17.根據(jù)權利要求16所述的多天線無線通信系統(tǒng)的操作方法,其特征在于包括步驟信道分解模塊(110)對信道矩陣進行幾何平均分解,得到線性變換矩陣Q、上三角矩陣R和預濾波矩陣P,其維數(shù)均為K×K;預編碼模塊(104)利用從接收端發(fā)送過來的預濾波陣P對信號進行線性預編碼;線性變換模塊(105)利用從接收端發(fā)送過來的線性變換矩陣Q的Hermitian轉置矩陣QH對信號進行線性變換;線性變換模塊(111)利用線性變換矩陣Q的Hermitian轉置矩陣QH對接收信號進行線性變換;以及所述迭代處理模塊(112)從第k個子信道開始,用該子信道接收信號減去第(k+1)個子信道的接收信號;然后,再除以上三角矩陣R的對角元rkk,恢復與第k個子信道對應的天線端的發(fā)射信號;更新k=k-1;并判斷k是否小于1,如果k<1,則結束,若k≥1,則針對更新后的k繼續(xù)進行處理。
全文摘要
一般的多天線系統(tǒng)中,由于各個子信道的增益不同,低信道增益的子信道將降低整個系統(tǒng)的誤比特率性能。因此,本發(fā)明提出了一種用在多天線無線通信系統(tǒng)中的聯(lián)合處理方法,包括在接收端,對信道矩陣進行幾何平均分解,得到線性變換矩陣Q、上三角矩陣R和預濾波矩陣P,其維數(shù)均為K×K,將預濾波矩陣P和線性變換矩陣Q反饋到發(fā)送端,利用Q
文檔編號H04L1/02GK1838578SQ200510055980
公開日2006年9月27日 申請日期2005年3月24日 優(yōu)先權日2005年3月24日
發(fā)明者黎海濤, 李繼峰 申請人:松下電器產(chǎn)業(yè)株式會社