CN103973627A - Full-rate distributed type multi-antenna bi-directional wireless cooperative relay transmission method - Google Patents

Abstract

A full-rate distributed type multi-antenna bi-directional wireless cooperative relay transmission method is suitable for a universal bi-directional relay communication network, wherein two communication source nodes and all relay nodes of the universal bi-directional relay communication network are respectively provided with a plurality of antennas. Information transmission between the nodes is carried out in an orthogonal frequency division multiplexing mode. A data communication process comprises two stages, wherein firstly, the communication source nodes conduct joint encoding on symbols to be transmitted in all sub-carriers and a plurality of time slots respectively and broadcast the encoded symbols to the relay nodes; secondly, each relay node conducts circular delay on received modulating signals in a time domain firstly, then the antennas on each relay node transmit processed signals to the communication source nodes sequentially in a time division working mode, and according to encoding and transmitting criterions, one source node obtains data transmitted by the other source node after eliminating signal interference of the former source node. The full-rate distributed type multi-antenna bi-directional wireless cooperative relay transmission method has the advantages that the expansibility is good, an encoding process is simple and easy to design, full-rate transmission can be achieved, the detection complexity is low, and very high application value is achieved.

Description

A kind of full rate distributed multi-antenna double-direction radio cooperating relay transmission method

Technical field

The present invention relates to a kind of full rate distributed multi-antenna double-direction radio cooperating relay transmission method, be applicable to the general scene that two communication source nodes and all via nodes all assemble many antennas, can realize full velocity transmission, cataloged procedure is simple and be easy to design, the detection model of receiving terminal and detection complexity all do not change with the increase of via node number, are conducive to reduce hardware implementation cost.

Background technology

Along with the fast development of radio communication service, high-quality mobile communication business requires to have caused the increase rapidly to capability of wireless communication system demand in world wide.But then, available wireless frequency spectrum is very limited, and how under limited frequency band condition, the spectrum efficiency that improves wireless communication system has become at present both at home and abroad the competitively heat subject of research.

In with Publication about Document (1) and document (2): (1) " Distributed MIMO Technologies in Cooperative Wireless Networks ", IEEE Communications Magazine, vol.49, no.5, pp.78-82, May2011; (2) " Spectral Efficiency of Distributed MIMO Systems ", IEEE Journal on Selected Areas in Communications, vol.31, no.10, pp.2112-2127, October2013, Ma and Wang etc. studies confirm that, distributed multiple-input and multiple-output (Distributed MIMO:Distributed Multiple Input Multiple Output) technology has become one of effective ways that can significantly improve capability of wireless communication system and the availability of frequency spectrum.Wherein, distributed wireless cooperating relay transmission method can also be utilized the radio relay communication feature of via node, and relaying and amplification, increase diversity performance, the reduction that realizes signal transmitted gross power, expanded the coverage of wireless network, has a extensive future.

In recent years, researcher mainly concentrates on two large fields to the research work of wireless distributed MIMO two-way communication theory both at home and abroad.First precoding technique, as with Publication about Document (3) to document (5): (3) " Joint Optimization for One and Two-Way MIMO AF Multiple-Relay Systems ", IEEE Transactions on Wireless Communications, vol.9, no.12, pp.3671-3681, December2010; (4) " Sum-Rate Maximization for Two-Way MIMO Amplify-and-Forward Relaying Systems ", IEEE69th Vehicular Technology Conference, VTC Spring, pp.1-5, April2009; (5) " A General Framework of Precoding Design for Multiple Two-way Relaying Communications ", IEEE Transactions on Signal Processing, vol.61, no.6, pp.1531-1535, March2013, can be by reaching the object that increases whole two-way wireless communication network capacity and transmittability to the precoding of communication source node modulation symbol to be transmitted and multiple via nodes place forward signal; But count and node antenna number all in larger MIMO Signal with Distributed Transmit Antennas at communication section, carry out the via node of precoding and not only will know the channel condition information relevant with this node (CSI:Channel State Information), sometimes also need to know the relevant CSI of other nodes, implementation complexity is high, is difficult to accomplish real-time processing.In addition, always there is certain error in the CSI valuation that each communication node obtains and actual value, and therefore, in real application systems, the cooperating relay scheme based on precoding technique is often difficult to reach theoretic heap(ed) capacity and optimal bit error performance.

On the other hand, distributed space time frequency coding can obtain multiple diversity gain potential in distributed multi-antenna system, improve error performance under transmitting terminal is not known the condition of any channel condition information, and its application in bidirectional communication network has obtained broad research equally.As in Publication about Document (6) and document (7): (6) " New Double Layer Space-Time Block Code for Distributed4 × 2MIMO Systems ", IEEE Wireless Communications and Networking Conference, pp.232-235, April2012; (7) " Distributed MIMO Coding Scheme with Low Decoding Complexity for Future Mobile TV Broadcasting ", Electronics Letters, vol.48, no.17, pp.1079-1081, August2012, can be by the stacking construction of the empty time-code of multiple small sizes be gone out to large-sized empty time-code, and then the diversity gain between antenna in acquisition MIMO Signal with Distributed Transmit Antennas, although its coding thinking is simple, be not easy to be generalized to general scene; In recent years, distributed Space-time-frequency Coding Technology based on different coding criterion and method for designing constantly comes forth, as with Publication about Document (8) to document (12): (8) " Adaptive Distributed Space-Time Coding Based on Adjustable Code Matrices for Cooperative MIMO Relaying Systems ", IEEE Transactions on Communications, vol.61, no.7, pp.2692-2703, July2013; (9) " Distributed Space Time Coding for Wireless Two-Way Relaying ", IEEE Transactions on Signal Processing, vol.61, no.4, pp.980-991, February2013; (10) " Distributed Concatenated Alamouti Codes for Two-Way Relaying Networks ", IEEE Wireless Communications Letters, vol.1, no.3, pp.197-200, June2012; (11) Novel Distributed Quasi-Orthogonal Space-Time Block Codes for Two-Way Two-Antenna Relay Networks ", IEEE Transactions on Wireless Communications, vol.12, no.9, pp.4338-4349, September 2013; (12) " Full-Rate Distributed Space-Time Codes for Cooperative Communications ", IEEE Transactions on Wireless Communications, vol.7, no.7, pp.2446-2451, July2008.Investigation discovery, the distributed space time frequency coding method majority of having announced is both at home and abroad all for a certain or a few application requirements design, and versatility is not strong, and also there is the problem such as coding, detection complexity height in some schemes.

For the two-way cooperative relay system of wireless distributed MIMO that is applicable to two communication source nodes and all via nodes and all assembles many antennas, also rarely have at present and can realize full velocity transmission, cataloged procedure simply and be easy to the general implementation method announcement designing, receiving terminal detection complexity is low, the key technology that this field relates to is to be worth further investigation.

Summary of the invention

Technical problem to be solved by this invention is to provide a kind of can simplify code Design process, realize the transmission of full rate bi-directional relaying and the lower universal distributed multi-antenna cooperative relay transmission method of receiving terminal detection complexity.

The present invention solves the problems of the technologies described above by the following technical solutions: a kind of full rate distributed multi-antenna double-direction radio cooperating relay transmission method, be applicable to the general bi-directional relaying communication network that two communication source nodes and all via nodes all assemble many antennas, define this network by K via node (R ₁, R ₂..., R _k) two communication source node (T of assistance ₁and T ₂) completing bidirectional data communication, each communication source node and via node assemble respectively N root antenna and M root antenna, note the system into (N, M, K) by abridging, define two communication source node T ₁and T ₂between lack tie link, each communication node is all operated under half duplex mode of communication, internodal communication adopts OFDM (OFDM:Orthogonal Frequency Division Multiplexing) mode to complete, each time slot, the every corresponding OFDM modulation transmissions of antenna.The double-direction radio cooperating relay transmission method principal character relating to is: once complete bidirectional data communication process mainly comprises two stages.

In the first stage, two communication source node T ₁and T ₂respectively the symbol to be transmitted in each subcarrier, multiple time slot is carried out to suitable combined coding, and signal after coding is sent to multiple via nodes with broadcast mode.

In second stage, first each via node carries out cyclic delay operation to each ofdm modulation signal receiving separately in time domain, the circulation delay amount that each via node adopts is therebetween all not identical, in signal amplification forwarding process subsequently, many antennas that assemble on each via node are sent to two communication source nodes by its reception signal after treatment successively with time-division working method, and source node can obtain the data that the other side's source node transmits according to coding and forwarding criterion after elimination self signal disturbs.

First stage, corresponding concrete communication process can be described below: based on OFDM communications mode, and two communication source node T ₁and T ₂original modulation symbol matrix to be sent on m time slot, a p subcarrier is used respectively s _mand c (p) _m(p) represent.

In formula, n _ffor can be in order to transmit the maximum sub-carrier number of data, b _mnand d (p) _mn(p) be respectively communication source node T ₁and T ₂on current time slots, current carrier wave its transmission symbol on root antenna, source block to be sent is crossed over M time slot.

Send diversity, two communication source node T for obtaining ₁and T ₂respectively the symbol to be transmitted on current carrier wave, all time slots, all antennas is carried out to combined coding, obtain final modulation symbol matrix S (p) and the C (p) sending, size is NM.

$S\left(p\right)=[{\tilde{s}}_1\left(p\right),{\tilde{s}}_2\left(p\right),...,{\tilde{s}}_M\left(p\right)],{\tilde{s}}_m\left(p\right)={\mathrm{\ψ}}_m{[s_1\left(p\right),s_2\left(p\right),...,s_M\left(p\right)]}^T---\left(3\right)$

$C\left(p\right)=[{\tilde{c}}_1\left(p\right),{\tilde{c}}_2\left(p\right),...,{\tilde{c}}_M\left(p\right)],{\tilde{c}}_m\left(p\right)={\mathrm{\ψ}}_m{[c_1\left(p\right),c_2\left(p\right),...,c_M\left(p\right)]}^T---\left(4\right)$

In formula ,] ^tthe transposition of representing matrix or vector, ψ _mfor the orthogonal combined coding matrix of N (NM), in order to realize the combined coding to the original modulation symbol of communication source node; In current stage of communication, the modulation symbol matrix S (p) after combined coding and C (p) are respectively by two communication source node T ₁and T ₂with single antenna average transmit power P ₁and P ₂be broadcast to each via node simultaneously.

$R_k\left(p\right)=\sqrt{2P_1}{H}_k\left(p\right)S\left(p\right)+\sqrt{2P_2}{G}_k\left(p\right)C\left(p\right)+W_k\left(p\right)---\left(5\right)$

$H_k\left(p\right)={[h_k^1\left(p\right),h_k^2\left(p\right),...,h_k^M\left(p\right)]}^T,h_k^m\left(p\right)={[h_k^{m1}\left(p\right),h_k^{m2}\left(p\right),...,h_k^{\mathrm{mN}}\left(p\right)]}^T---\left(6\right)$

$G_k\left(p\right)={[g_k^1\left(p\right),g_k^2\left(p\right),...,g_k^M\left(p\right)]}^T,g_k^m\left(p\right)={[g_k^{m1}\left(p\right),g_k^{m2}\left(p\right),...,g_k^{\mathrm{mN}}\left(p\right)]}^T---\left(7\right)$

In formula, R _k(p) be the reception signal matrix of k via node, size is MM, H _kand G (p) _k(p) be respectively two communication source node T ₁and T ₂to the channel matrix of k via node, size is MN, and wherein each element is the channel fading coefficient between respective antenna with available average is 0, variance N ₀=1 multiple gaussian variable characterizes, W _k(p) for being added in the white complex gaussian noise matrix at k via node place, size is MM, wherein each element can be 0 by average, variance N ₀=1 multiple gaussian variable characterizes.

Concrete communication process corresponding to second stage can be described below: first each via node carries out circulation delay to each ofdm modulation signal receiving separately in time domain, and this operation will change the phase place of corresponding frequency-region signal; In signal amplification forwarding process subsequently, many antennas that assemble on each via node are sent to two communication source nodes by its reception signal after treatment successively with time-division working method; Wherein communication source node T ₂reception signal can describe by following Mathematical Modeling.

$Y_2\left(p\right)=a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{{\tilde{G}}_k\left(p\right)\mathrm{diag}(r_k^p\left(\mathrm{1,1}\right),r_k^p\left(\mathrm{2,2}\right),...,r_k^p(M,M))e^{-j\frac{2p}{N_F}p{D}_k}\right\}+{\tilde{W}}_2\left(p\right)---\left(8\right)$

$r_k^p(m,m)=\sqrt{2P_1}{\left[h_k^m\left(p\right)\right]}^T{\tilde{s}}_m\left(p\right)+{\sqrt{2P_2}\left[g_k^m\left(p\right)\right]}^T{\tilde{c}}_m\left(p\right)+w_k^p(m,m)---\left(9\right)$

${\tilde{G}}_k\left(p\right)=[{\tilde{g}}_k^1\left(p\right),{\tilde{g}}_k^2\left(p\right),...,{\tilde{g}}_k^M\left(p\right)],{\tilde{g}}_k^m\left(p\right)={[{\tilde{g}}_k^{1m}\left(p\right),{\tilde{g}}_k^{2m}\left(p\right),...,{\tilde{g}}_k^{\mathrm{Nm}}\left(p\right)]}^T---\left(10\right)$

$a=\sqrt{\frac{2P_{3}M}{2P_{1}N+2P_{2}N+N_0}}---\left(11\right)$

P＝NP ₁+NP ₂+KMP ₃(12)

$P_1=P_2=\frac{P}{4N},P_3=\frac{P}{2\mathrm{KM}}---\left(13\right)$

In formula, diag (X ₁, X ₂..., X _m) expression employing X ₁, X ₂... X _mthe block diagonal matrix of composition, be that k via node place receives signal matrix R _k(p) m element on leading diagonal, for noise matrix W _k(p) m element on leading diagonal, be that k via node is to communication source node T ₂channel fading coefficient while receiving signal after forward process between n root reception antenna and m transmit antennas, a is power amplification coefficient power amplification ratio when signal after via node forward process, D _kbe the circulation delay amount that k via node adopts, on each via node, all antennas adopt identical circulation delay amount, and the circulation delay amount adopting between via node is all not identical, for being added in communication source node T ₂place receives the white complex gaussian noise matrix in signal matrix, P ₃for the average transmit power that single antenna on each via node adopts, definition SNR=10log ₁₀(P/N ₀) (dB) as the signal to noise ratio standard of weighing whole two-way cooperation communication system performance quality, communication source node T ₂must eliminate self interference signal in each time slot carrying out before input.

$\begin{array}{c}{Y}_2^m\left(p\right)=a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{2P_1}{\tilde{g}}_k^m\left(p\right){\left[h_k^m\left(p\right)\right]}^T{\tilde{s}}_m\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\right\}\\ +a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{2P_2}{\tilde{g}}_k^m\left(p\right){\left[g_k^m\left(p\right)\right]}^T{\tilde{c}}_m\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\right\}\\ +a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{{\tilde{g}}_k^m\left(p\right)w_k^p(m,m)e^{-j\frac{2p}{N_F}p{D}_k}\right\}+{\tilde{w}}_2^m\left(p\right)\end{array}---\left(15\right)$

${\tilde{Y}}_2^m\left(p\right)=Y_2^m\left(p\right)-a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{2P_2}{\tilde{g}}_k^m\left(p\right){\left[g_k^m\left(p\right)\right]}^T{\tilde{c}}_m\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\right\}---\left(16\right)$

$E_2^m\left(p\right)=a\underset{k-1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{2P_1}{\tilde{g}}_k^m\left(p\right){\left[h_k^m\left(p\right)\right]}^T{e}^{-j\frac{2p}{N_F}p{D}_k}\right\}---\left(17\right)$

In formula, with be respectively communication source node T ₂on m time slot, a p subcarrier, eliminate the equivalent received matrix and the channel matrix that after self signal disturbs, obtain, for in corresponding noise on m column element position, according to coding and forward criterion and obtain following equivalent communication model.

$\left[\begin{array}{c}{\tilde{Y}}_2^1\left(p\right)\\ {\tilde{Y}}_2^2\left(p\right)\\ .\\ .\\ .\\ {\tilde{Y}}_2^M\left(p\right)\end{array}\right]=\mathrm{diag}(E_2^1\left(p\right),E_2^2\left(p\right),...,E_2^M\left(p\right))\left[\begin{array}{c}{\mathrm{\ψ}}_1\\ {\mathrm{\ψ}}_2\\ .\\ .\\ .\\ {\mathrm{\ψ}}_M\end{array}\right]\left[\begin{array}{c}{\left[s_1\left(p\right)\right]}^T\\ {\left[s_2\left(p\right)\right]}^T\\ .\\ .\\ .\\ {\left[s_M\left(p\right)\right]}^T\end{array}\right]+\left[\begin{array}{c}{Z}_2^1\\ Z_2^2\\ .\\ .\\ .\\ Z_2^M\end{array}\right]---\left(18\right)$

$Z_2^m=a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{{\tilde{g}}_k^m\left(p\right){w_k^p(m,m)e}^{-j\frac{2p}{N_F}p{D}_k}\right\}+{\tilde{W}}_2^m\left(p\right)---\left(19\right)$

Based on above-mentioned equivalent communication model, communication source node T ₂can adopt joint-detection mode to obtain communication source node T simultaneously ₁the all modulation symbols that send on current subcarrier, in a M time slot, and obtain by communication source node T ₁and the space diversity gain that brings of many antennas on single via node, the space diversity between via node is converted to frequency diversity by circulation delay transmission means, and this kind of frequency diversity can be obtained by error correction coding.

In addition, be it can also be seen that by above-mentioned equivalent communication model, the full rate distributed multi-antenna double-direction radio cooperating relay transmission method the present invention relates to is applicable to the general scene that two communication source nodes and all via nodes all assemble many antennas, can realize full velocity transmission, cataloged procedure is simple and be easy to design, the detection model of receiving terminal and detection complexity all do not change with the increase of via node number, are conducive to reduce hardware implementation cost.

By communication source node T ₂send to communication source node T ₁data can be by similar approach at T ₁place obtains.

Repeating process corresponding to second stage is not limited to a certain implementation, can be used in conjunction with existing other coding methods favorable expandability; In the time that the antenna number of each via node assembling is even number, the mode that can adopt timesharing to transmit two orthogonal Space Time Coding of antenna Alamouti completes data retransmission.Now, the detailed process of second stage can be described below:

$Y_2\left(p\right)=\sqrt{\frac{a}{2}}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{{\tilde{G}}_k\left(p\right)\mathrm{diag}(A_1,A_2,...,A_{M/2})e^{-j\frac{2p}{D_F}p{D}_k}\right\}+{\tilde{W}}_2\left(p\right)---\left(20\right)$

In formula ,] ^*represent the conjugation of scalar, A _mencoder matrix during for Alamouti orthogonal space, the reliability, the reduction receiving terminal that are conducive to improve forwarding data detect the error rate.Encoding characteristics during according to orthogonal space, to communication source node T ₂reception signal matrix modification after can obtain following mode.

$\begin{array}{c}\left[\begin{array}{c}{Y}_2^{2m-1}\left(p\right)\\ {\left[Y_2^{2m}\left(p\right)\right]}^{*}\end{array}\right]=\sqrt{\frac{a}{2}}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{2P_1}{F}_k^m\mathrm{diag}({\left[h_k^{2m-1}\left(p\right)\right]}^T,{\left[h_k^{2m}\left(p\right)\right]}^T)\left[\begin{array}{c}{\tilde{s}}_{2m-1}\left(p\right)\\ {\tilde{s}}_{2m}\left(p\right)\end{array}\right]\right\}\\ +\sqrt{\frac{a}{2}}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{2P_2}{F}_k^m\mathrm{diag}({\left[g_k^{2m-1}\left(p\right)\right]}^T,{\left[g_k^{2m}\left(p\right)\right]}^T)\left[\begin{array}{c}{\tilde{c}}_{2m-1}\left(p\right)\\ {\tilde{c}}_{2m}\left(p\right)\end{array}\right]\right\}\\ +\sqrt{\frac{a}{2}}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{F_k^m\left[\begin{array}{c}{w}_k^p(2m-\mathrm{1,2}m-1)\\ w_k^p(2m,2m)\end{array}\right]\right\}+\left[\begin{array}{c}{\tilde{W}}_2^{2m-1}\left(p\right)\\ {\left[{\tilde{W}}_2^{2m}\left(p\right)\right]}^{*}\end{array}\right]\end{array}---\left(22\right)$

$F_k^m=\left[\begin{array}{cc}{\tilde{g}}_k^{2m-1}\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}& {\tilde{g}}_k^{2m}\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\\ {\left[{\tilde{g}}_k^{2m}\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\right]}^{*}& -{\left[{\tilde{g}}_k^{2m-1}\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\right]}^{*}\end{array}\right]---\left(23\right)$

For above-mentioned mode, can be at communication source node T ₂after eliminating self signal and disturbing, joint-detection is out from communication source node T ₁initial data.By communication source node T ₂send to communication source node T ₁data by similar approach at T ₁place obtains.

The invention has the advantages that: the radio data transmission method relating to is applicable to general bi-directional relaying communication network, favorable expandability, cataloged procedure is simple and be easy to design, can realize full velocity transmission and detection complexity is low, has very high practical value.

Brief description of the drawings

Fig. 1 is the theory diagram of full rate distributed multi-antenna double-direction radio cooperating relay transmission method.

Fig. 2 is that while adopting different via nodes to count, the present invention announces the performance of BER comparison diagram of coding transmission scheme.

Embodiment

Describe the present invention below in conjunction with accompanying drawing.

As shown in Figure 1, the present invention adopts a kind of general bi-directional relaying communication network that is applicable to two communication source nodes and all via nodes and all assembles many antennas, defines this network by K via node (R ₁, R ₂..., R _k) two communication source node (T of assistance ₁and T ₂) completing bidirectional data communication, each communication source node and via node assemble respectively N root antenna and M root antenna, note the system into (N, M, K) by abridging, define two communication source node T ₁and T ₂between lack tie link, each communication node is all operated under half duplex mode of communication, internodal communication adopts OFDM (OFDM:Orthogonal Frequency Division Multiplexing) mode to complete, each time slot, the every corresponding OFDM modulation transmissions of antenna.The double-direction radio cooperating relay transmission method principal character relating to is: once complete bidirectional data communication process mainly comprises two stages.

In the first stage, based on OFDM communications mode, two communication source node T ₁and T ₂original modulation symbol matrix to be sent on m time slot, a p subcarrier is used respectively s _mand c (p) _m(p) represent.

In formula, n _ffor can be in order to transmit the maximum sub-carrier number of data, b _mnand d (p) _mn(p) be respectively communication source node T ₁and T ₂on current time slots, current carrier wave its transmission symbol on root antenna, source block to be sent is crossed over M time slot.

Send diversity, two communication source node T for obtaining ₁and T ₂respectively the symbol to be transmitted on current carrier wave, all time slots, all antennas is carried out to combined coding, obtain final modulation symbol matrix S (p) and the C (p) sending, size is NM.

$S\left(p\right)=[{\tilde{s}}_1\left(p\right),{\tilde{s}}_2\left(p\right),...,{\tilde{s}}_M\left(p\right)],{\tilde{s}}_m\left(p\right)={\mathrm{\ψ}}_m{[s_1\left(p\right),s_2\left(p\right),...,s_M\left(p\right)]}^T---\left(3\right)$

$C\left(p\right)=[{\tilde{c}}_1\left(p\right),{\tilde{c}}_2\left(p\right),...,{\tilde{c}}_M\left(p\right)],{\tilde{c}}_m\left(p\right)={\mathrm{\ψ}}_m{[c_1\left(p\right),c_2\left(p\right),...,c_M\left(p\right)]}^T---\left(4\right)$

In formula ,] ^tthe transposition of representing matrix or vector, ψ _mfor the orthogonal combined coding matrix of N (NM), in order to realize the combined coding to the original modulation symbol of communication source node; In current stage of communication, the modulation symbol matrix S (p) after combined coding and C (p) are respectively by two communication source node T ₁and T ₂with single antenna average transmit power P ₁and P ₂be broadcast to each via node simultaneously.

$R_k\left(p\right)=\sqrt{2P_1}{H}_k\left(p\right)S\left(p\right)+\sqrt{2P_2}{G}_k\left(p\right)C\left(p\right)+W_k\left(p\right)---\left(5\right)$

$H_k\left(p\right)={[h_k^1\left(p\right),h_k^2\left(p\right),...,h_k^M\left(p\right)]}^T,h_k^m\left(p\right)={[h_k^{m1}\left(p\right),h_k^{m2}\left(p\right),...,h_k^{\mathrm{mN}}\left(p\right)]}^T---\left(6\right)$

$G_k\left(p\right)={[g_k^1\left(p\right),g_k^2\left(p\right),...,g_k^M\left(p\right)]}^T,g_k^m\left(p\right)={[g_k^{m1}\left(p\right),g_k^{m2}\left(p\right),...,g_k^{\mathrm{mN}}\left(p\right)]}^T$

$\left(7\right)$

In formula, R _k(p) be the reception signal matrix of k via node, size is MM, H _kand G (p) _k(p) be respectively two communication source node T ₁and T ₂to the channel matrix of k via node, size is MN, and wherein each element is the channel fading coefficient between respective antenna with available average is 0, variance N ₀=1 multiple gaussian variable characterizes, W _k(p) for being added in the white complex gaussian noise matrix at k via node place, size is MM, wherein each element can be 0 by average, variance N ₀=1 multiple gaussian variable characterizes.

In second stage, first each via node carries out circulation delay to each ofdm modulation signal receiving separately in time domain, and this operation will change the phase place of corresponding frequency-region signal; In signal amplification forwarding process subsequently, many antennas that assemble on each via node are sent to two communication source nodes by its reception signal after treatment successively with time-division working method; Wherein communication source node T ₂reception signal can describe by following Mathematical Modeling.

$Y_2\left(p\right)=a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{{\tilde{G}}_k\left(p\right)\mathrm{diag}(r_k^p\left(\mathrm{1,1}\right),r_k^p\left(\mathrm{2,2}\right),...,r_k^p(M,M))e^{-j\frac{2p}{N_F}p{D}_k}\right\}+{\tilde{W}}_2\left(p\right)---\left(8\right)$

$r_k^p(m,m)=\sqrt{2P_1}{\left[h_k^m\left(p\right)\right]}^T{\tilde{s}}_m\left(p\right)+{\sqrt{2P_2}\left[g_k^m\left(p\right)\right]}^T{\tilde{c}}_m\left(p\right)+w_k^p(m,m)---\left(9\right)$

${\tilde{G}}_k\left(p\right)=[{\tilde{g}}_k^1\left(p\right),{\tilde{g}}_k^2\left(p\right),...,{\tilde{g}}_k^M\left(p\right)],{\tilde{g}}_k^m\left(p\right)={[{\tilde{g}}_k^{1m}\left(p\right),{\tilde{g}}_k^{2m}\left(p\right),...,{\tilde{g}}_k^{\mathrm{Nm}}\left(p\right)]}^T---\left(10\right)$

$a=\sqrt{\frac{2P_{3}M}{2P_{1}N+2P_{2}N+N_0}}---\left(11\right)$

P＝NP ₁+NP ₂+KMP ₃(12)

$P_1=P_2=\frac{P}{4N},P_3=\frac{P}{2\mathrm{KM}}---\left(13\right)$

In formula, diag (X ₁, X ₂..., X _m) expression employing X ₁, X ₂... X _mthe block diagonal matrix of composition, be that k via node place receives signal matrix R _k(p) m element on leading diagonal, for noise matrix W _k(p) m element on leading diagonal, be that k via node is to communication source node T ₂channel fading coefficient while receiving signal after forward process between n root reception antenna and m transmit antennas, a is power amplification coefficient power amplification ratio when signal after via node forward process, D _kbe the circulation delay amount that k via node adopts, on each via node, all antennas adopt identical circulation delay amount, and the circulation delay amount adopting between via node is all not identical, for being added in communication source node T ₂place receives the white complex gaussian noise matrix in signal matrix, P ₃for the average transmit power that single antenna on each via node adopts, definition SNR=10log ₁₀(P/N ₀) (dB) as the signal to noise ratio standard of weighing whole two-way cooperation communication system performance quality, communication source node T ₂must eliminate self interference signal in each time slot carrying out before input.

$\begin{array}{c}{Y}_2^m\left(p\right)=a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{2P_1}{\tilde{g}}_k^m\left(p\right){\left[h_k^m\left(p\right)\right]}^T{\tilde{s}}_m\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\right\}\\ +a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{2P_2}{\tilde{g}}_k^m\left(p\right){\left[g_k^m\left(p\right)\right]}^T{\tilde{c}}_m\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\right\}\\ +a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{{\tilde{g}}_k^m\left(p\right)w_k^p(m,m)e^{-j\frac{2p}{N_F}p{D}_k}\right\}+{\tilde{w}}_2^m\left(p\right)\end{array}---\left(15\right)$

${\tilde{Y}}_2^m\left(p\right)=Y_2^m\left(p\right)-a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{2P_2}{\tilde{g}}_k^m\left(p\right){\left[g_k^m\left(p\right)\right]}^T{\tilde{c}}_m\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\right\}---\left(16\right)$

$E_2^m\left(p\right)=a\underset{k-1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{2P_1}{\tilde{g}}_k^m\left(p\right){\left[h_k^m\left(p\right)\right]}^T{e}^{-j\frac{2p}{N_F}p{D}_k}\right\}---\left(17\right)$

In formula, with be respectively communication source node T ₂on m time slot, a p subcarrier, eliminate the equivalent received matrix and the channel matrix that after self signal disturbs, obtain, for in corresponding noise on m column element position, according to coding and forward criterion and obtain following equivalent communication model.

$\left[\begin{array}{c}{\tilde{Y}}_2^1\left(p\right)\\ {\tilde{Y}}_2^2\left(p\right)\\ .\\ .\\ .\\ {\tilde{Y}}_2^M\left(p\right)\end{array}\right]=\mathrm{diag}(E_2^1\left(p\right),E_2^2\left(p\right),...,E_2^M\left(p\right))\left[\begin{array}{c}{\mathrm{\ψ}}_1\\ {\mathrm{\ψ}}_2\\ .\\ .\\ .\\ {\mathrm{\ψ}}_M\end{array}\right]\left[\begin{array}{c}{\left[s_1\left(p\right)\right]}^T\\ {\left[s_2\left(p\right)\right]}^T\\ .\\ .\\ .\\ {\left[s_M\left(p\right)\right]}^T\end{array}\right]+\left[\begin{array}{c}{Z}_2^1\\ Z_2^2\\ .\\ .\\ .\\ Z_2^M\end{array}\right]---\left(18\right)$

$Z_2^m=a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{{\tilde{g}}_k^m\left(p\right){w_k^p(m,m)e}^{-j\frac{2p}{N_F}p{D}_k}\right\}+{\tilde{W}}_2^m\left(p\right)---\left(19\right)$

Based on above-mentioned equivalent communication model, communication source node T ₂can adopt joint-detection mode to obtain communication source node T simultaneously ₁the all modulation symbols that send on current subcarrier, in a M time slot, and obtain by communication source node T ₁and the space diversity gain that brings of many antennas on single via node, the space diversity between via node is converted to frequency diversity by circulation delay transmission means, can be obtained by error correction coding and decoding process.

By communication source node T ₂send to communication source node T ₁data can be by similar approach at T ₁place obtains.

In addition, repeating process corresponding to second stage is not limited to a certain implementation, can be used in conjunction with existing other coding methods favorable expandability; In the time that the antenna number of each via node assembling is even number, the mode that can adopt timesharing to transmit two orthogonal Space Time Coding of antenna Alamouti completes data retransmission.Now, the detailed process of second stage can be described below.

$Y_2\left(p\right)=\sqrt{\frac{a}{2}}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{{\tilde{G}}_k\left(p\right)\mathrm{diag}(A_1,A_2,...,A_{M/2})e^{-j\frac{2p}{D_F}p{D}_k}\right\}+{\tilde{W}}_2\left(p\right)---\left(20\right)$

In formula ,] ^*represent the conjugation of scalar, A _mencoder matrix during for Alamouti orthogonal space, the reliability, the reduction receiving terminal that are conducive to improve forwarding data detect the error rate.Encoding characteristics during according to orthogonal space, to communication source node T ₂reception signal matrix modification after can obtain following mode.

$\begin{array}{c}\left[\begin{array}{c}{Y}_2^{2m-1}\left(p\right)\\ {\left[Y_2^{2m}\left(p\right)\right]}^{*}\end{array}\right]=\sqrt{\frac{a}{2}}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{2P_1}{F}_k^m\mathrm{diag}({\left[h_k^{2m-1}\left(p\right)\right]}^T,{\left[h_k^{2m}\left(p\right)\right]}^T)\left[\begin{array}{c}{\tilde{s}}_{2m-1}\left(p\right)\\ {\tilde{s}}_{2m}\left(p\right)\end{array}\right]\right\}\\ +\sqrt{\frac{a}{2}}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{2P_2}{F}_k^m\mathrm{diag}({\left[g_k^{2m-1}\left(p\right)\right]}^T,{\left[g_k^{2m}\left(p\right)\right]}^T)\left[\begin{array}{c}{\tilde{c}}_{2m-1}\left(p\right)\\ {\tilde{c}}_{2m}\left(p\right)\end{array}\right]\right\}\\ +\sqrt{\frac{a}{2}}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{F_k^m\left[\begin{array}{c}{w}_k^p(2m-\mathrm{1,2}m-1)\\ w_k^p(2m,2m)\end{array}\right]\right\}+\left[\begin{array}{c}{\tilde{W}}_2^{2m-1}\left(p\right)\\ {\left[{\tilde{W}}_2^{2m}\left(p\right)\right]}^{*}\end{array}\right]\end{array}---\left(22\right)$

$F_k^m=\left[\begin{array}{cc}{\tilde{g}}_k^{2m-1}\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}& {\tilde{g}}_k^{2m}\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\\ {\left[{\tilde{g}}_k^{2m}\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\right]}^{*}& -{\left[{\tilde{g}}_k^{2m-1}\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\right]}^{*}\end{array}\right]---\left(23\right)$

For above-mentioned mode, can be at communication source node T ₂after eliminating self signal and disturbing, joint-detection is out from communication source node T ₁initial data.

The full rate distributed multi-antenna double-direction radio cooperating relay transmission method that the applicant announces the present invention by Computer Simulation is carried out bit error rate (BER:Bit Error Rate) performance test, and provides simulation result corresponding under multiple systems parameter configuration.Main test environment parameter setting: adopt the multipath channel models with frequency selective fading characteristic to produce the internodal transmission channel of various types of communication, adopt that code rate is 1/2, constraint length is 7, generator polynomial is (171,133) convolution code is carried out chnnel coding to the original bit stream of two communication source nodes, adopt OFDM mode to transmit baseband modulation signal, in order to transmit the maximum sub-carrier number N of data _f=256, the additivity white complex gaussian noise that introduce average in data transmission procedures at different levels and be 0, variance is 1, high frequency modulated carrier frequency value is made as 3.5GHz.

(1), while adopting different via nodes to count, the present invention announces the performance of BER comparison of coding transmission scheme

The simulation result that Fig. 2 provides be adopt at communication source node that QPSK modulation system, frequency selective fading channels number of path are 1, communication source node and all via node antenna number be set to N=M=2, count K by changing via node, test to obtain.Herein will be in second stage repeating process, while adopting single antenna time division way to complete bi-directional relaying transmission, corresponding implementation of the present invention is called method I; Corresponding implementation of the present invention is called method II when adopting during based on two antenna Alamouti orthogonal space that encoding block time division way completes bi-directional relaying transmission.

As seen from Figure 2, along with the increase of via node quantity, two kinds of implementations corresponding to the present invention all can obtain larger diversity gain and lower bit error rate; Under identical via node number and signal to noise ratio condition, method II is owing to having introduced the orthogonal Space Time Coding of two antenna Alamouti in repeating process, thereby can obtain coding gain and the error performance that ratio method I is higher.As can be seen here, double-direction radio cooperating relay transmission method disclosed by the invention is not limited to a certain implementation, can forward coding method with existing other and be used in conjunction with, favorable expandability.

In addition, the detection model of receiving terminal and detection complexity all do not change with the increase of via node number, be not only applicable to two communication source nodes and all via nodes and all assemble the general scene of many antennas, and be conducive to reduce hardware design complexity and realize cost.

(2) the performance of BER comparison while adopting different coding method

Adopt at communication source node that BPSK modulation system, frequency selective fading channels number of path are 1, communication source node and all via node antenna number are set to N=M=2, via node is counted K=2, encoded and pass-through mode by change, test, test result shows: in contrast to the wireless two-way cooperative relay network that communication source node does not adopt combined coding technology or do not adopt time-division pass-through mode, transmission method disclosed by the invention can obtain lower bit error rate, has more obvious error performance advantage.

The foregoing is only the general embodiment of the invention; not in order to limit the invention; any amendment of doing within all spirit in the invention and principle, be equal to and replace and improvement etc., within all should being included in the protection range of the invention.

Claims (4)

1. a full rate distributed multi-antenna double-direction radio cooperating relay transmission method, is applicable to two communication source nodes and all via nodes and all assembles the general bi-directional relaying communication network of many antennas, defines this network by K via node: R ₁, R ₂..., R _kassist two communication source node: T ₁and T ₂complete bidirectional data communication, each communication source node and via node assemble respectively N root antenna and M root antenna, note the system into (N, M, K) by abridging, define two communication source node T ₁and T ₂between lack tie link, each communication node is all operated under half duplex mode of communication, internodal communication adopts OFDM (OFDM:Orthogonal Frequency Division Multiplexing) mode to complete, each time slot, the every corresponding OFDM modulation transmissions of antenna, is characterized in that: comprise two stages:

In the first stage, two communication source node T ₁and T ₂respectively the symbol to be transmitted in each subcarrier, multiple time slot is carried out to combined coding, and signal after coding is sent to multiple via nodes with broadcast mode;

In second stage, first each via node carries out cyclic delay operation to each ofdm modulation signal receiving separately in time domain, the circulation delay amount that each via node adopts is therebetween all not identical, in signal amplification forwarding process subsequently, many antennas that assemble on each via node are sent to two communication source nodes by its reception signal after treatment successively with time-division working method, and source node can obtain the data that the other side's source node transmits according to coding and forwarding criterion after elimination self signal disturbs.

2. full rate distributed multi-antenna double-direction radio cooperating relay transmission method according to claim 1, is characterized in that: the first stage, corresponding concrete communication process can be described below.

In the first stage, based on OFDM communications mode, two communication source node T ₁and T ₂original modulation symbol matrix to be sent on m time slot, a p subcarrier is used respectively s _mand c (p) _m(p) represent:

In formula, n _ffor can be in order to transmit the maximum sub-carrier number of data, b _mnand d (p) _mn(p) be respectively communication source node T ₁and T ₂on current time slots, current carrier wave its transmission symbol on root antenna, source block to be sent is crossed over M time slot;

Send diversity, two communication source node T for obtaining ₁and T ₂respectively the symbol to be transmitted on current carrier wave, all time slots, all antennas is carried out to combined coding, obtain final modulation symbol matrix S (p) and the C (p) sending, size is NM:

$S\left(p\right)=[{\tilde{s}}_1\left(p\right),{\tilde{s}}_2\left(p\right),...,{\tilde{s}}_M\left(p\right)],{\tilde{s}}_m\left(p\right)={\mathrm{\ψ}}_m{[s_1\left(p\right),s_2\left(p\right),...,s_M\left(p\right)]}^T---\left(3\right)$

$C\left(p\right)=[{\tilde{c}}_1\left(p\right),{\tilde{c}}_2\left(p\right),...,{\tilde{c}}_M\left(p\right)],{\tilde{c}}_m\left(p\right)={\mathrm{\ψ}}_m{[c_1\left(p\right),c_2\left(p\right),...,c_M\left(p\right)]}^T---\left(4\right)$

In formula ,] ^tthe transposition of representing matrix or vector, ψ _mfor the orthogonal combined coding matrix of N (N M), in order to realize the combined coding to the original modulation symbol of communication source node; In current stage of communication, the modulation symbol matrix S (p) after combined coding and C (p) are respectively by two communication source node T ₁and T ₂with single antenna average transmit power P ₁and P ₂be broadcast to each via node simultaneously:

$R_k\left(p\right)=\sqrt{{2P}_1}{H}_k\left(p\right)S\left(p\right)+\sqrt{{2P}_2}{G}_k\left(p\right)C\left(p\right)+W_k\left(p\right)---\left(5\right)$

$H_k\left(p\right)={[h_k^1\left(p\right),h_k^2\left(p\right),...,h_k^M\left(p\right)]}^T,h_k^m\left(p\right)={[h_k^{m1}\left(p\right),h_k^{m2}\left(p\right),...,h_k^{\mathrm{mN}}\left(p\right)]}^T---\left(6\right)$

$G_k\left(p\right)={[g_k^1\left(p\right),g_k^2\left(p\right),...,g_k^M\left(p\right)]}^T,g_k^m\left(p\right)={[g_k^{m1}\left(p\right),g_k^{m2}\left(p\right),...,g_k^{\mathrm{mN}}\left(p\right)]}^T---\left(7\right)$

In formula, R _k(p) be the reception signal matrix of k via node, size is MM, H _kand G (p) _k(p) be respectively two communication source node T ₁and T ₂to the channel matrix of k via node, size is MN, and wherein each element is the channel fading coefficient between respective antenna with available average is 0, variance N ₀=1 multiple gaussian variable characterizes, W _k(p) for being added in the white complex gaussian noise matrix at k via node place, size is MM, wherein each element can be 0 by average, variance N ₀=1 multiple gaussian variable characterizes.

3. full rate distributed multi-antenna double-direction radio cooperating relay transmission method according to claim 2, is characterized in that: concrete communication process corresponding to second stage can be described below.

In second stage, many antennas that assemble on each via node are sent to two communication source nodes, wherein communication source node T by its reception signal after treatment successively with single antenna time-division working method ₂reception signal can describe by following Mathematical Modeling:

$Y_2\left(p\right)=a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{{\tilde{G}}_k\left(p\right)\mathrm{diag}(r_k^p\left(\mathrm{1,1}\right),r_k^p\left(\mathrm{2,2}\right),...,r_k^p(M,M))e^{-j\frac{2p}{N_F}p{D}_k}\right\}+{\tilde{W}}_2\left(p\right)---\left(8\right)$

$r_k^p(m,m)=\sqrt{{2P}_1}{\left[h_k^m\left(p\right)\right]}^T{\tilde{s}}_m\left(p\right)+\sqrt{{2P}_2}{\left[g_k^m\left(p\right)\right]}^T{\tilde{c}}_m\left(p\right)+w_k^p(m,m)---\left(9\right)$

${\tilde{G}}_k\left(p\right)=[{\tilde{g}}_k^1\left(p\right),{\tilde{g}}_k^2\left(p\right),...,{\tilde{g}}_k^M\left(p\right)],{\tilde{g}}_k^m\left(p\right)={[{\tilde{g}}_k^{1m}\left(p\right),{\tilde{g}}_k^{2m}\left(p\right),...,{\tilde{g}}_k^{\mathrm{Nm}}\left(p\right)]}^T---\left(10\right)$

$a=\sqrt{\frac{2P_{3}M}{2P_{1}N+2P_{2}N+N_0}}---\left(11\right)$

P＝NP ₁+NP ₂+KMP ₃(12)

$P_1=P_2=\frac{P}{4N},P_3=\frac{P}{2\mathrm{KM}}---\left(13\right)$

In formula, diag (X ₁, X ₂, L, X _m) expression employing X ₁, X ₂, LX _mthe block diagonal matrix of composition, be that k via node place receives signal matrix R _k(p) m element on leading diagonal, for noise matrix W _k(p) m element on leading diagonal, be that k via node is to communication source node T ₂channel fading coefficient while receiving signal after forward process between n root reception antenna and m transmit antennas, a is power amplification coefficient power amplification ratio when signal after via node forward process, D _kbe the circulation delay amount that k via node adopts, on each via node, all antennas adopt identical circulation delay amount, and the circulation delay amount adopting between via node is all not identical, for being added in communication source node T ₂place receives the white complex gaussian noise matrix in signal matrix, P ₃for the average transmit power that single antenna on each via node adopts, definition SNR=10log ₁₀(P/N ₀) (dB) as the signal to noise ratio standard of weighing whole two-way cooperation communication system performance quality, communication source node T ₂must eliminate self interference signal in each time slot carrying out before input:

$\begin{array}{c}{Y}_2^m\left(p\right)=a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{{2P}_1}{\tilde{g}}_k^m\left(p\right){\left[h_k^m\left(p\right)\right]}^T{\tilde{s}}_m\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\right\}\\ +a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{{2P}_2}{\tilde{g}}_k^m\left(p\right){\left[g_k^m\left(p\right)\right]}^T{\tilde{c}}_m\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\right\}\\ +a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{{\tilde{g}}_k^m\left(p\right)w_k^p(m,m)e^{-j\frac{2p}{N_F}p{D}_k}\right\}+{\tilde{W}}_2^m\left(p\right)\end{array}---\left(15\right)$

${\tilde{Y}}_2^m\left(p\right)=Y_2^m\left(p\right)-a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{{2P}_2}{\tilde{g}}_k^m\left(p\right){\left[g_k^m\left(p\right)\right]}^T{\tilde{c}}_m\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\right\}---\left(16\right)$

$E_2^m\left(p\right)=a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{{2P}_1}{\tilde{g}}_k^m\left(p\right){\left[h_k^m\left(p\right)\right]}^T{e}^{-j\frac{2p}{N_F}p{D}_k}\right\}---\left(17\right)$

In formula, with be respectively communication source node T ₂on m time slot, a p subcarrier, eliminate the equivalent received matrix and the channel matrix that after self signal disturbs, obtain, for in corresponding noise on m column element position, according to coding and forward criterion and obtain following equivalent communication model:

$\left[\begin{array}{c}{\tilde{Y}}_2^1\left(p\right)\\ {\tilde{Y}}_2^2\left(p\right)\\ .\\ .\\ .\\ {\tilde{Y}}_2^M\left(p\right)\end{array}\right]=\mathrm{diag}(E_2^1\left(p\right),E_2^2\left(p\right),...,E_2^M\left(p\right))\left[\begin{array}{c}{\mathrm{\ψ}}_1\\ {\mathrm{\ψ}}_2\\ .\\ .\\ .\\ {\mathrm{\ψ}}_M\end{array}\right]\left[\begin{array}{c}{\left[s_1\left(p\right)\right]}^T\\ {\left[s_2\left(p\right)\right]}^T\\ .\\ .\\ .\\ {\left[s_M\left(p\right)\right]}^T\end{array}\right]+\begin{array}{}\end{array}\left[\begin{array}{c}{Z}_2^1\\ Z_2^2\\ .\\ .\\ .\\ Z_2^M\end{array}\right]---\left(18\right)$

$Z_2^m=a\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{{\tilde{g}}_k^m\left(p\right)w_k^p(m,m)e^{-j\frac{2p}{N_F}p{D}_k}\right\}+{\tilde{W}}_2^m\left(p\right)---\left(19\right)$

Based on above-mentioned equivalent communication model, communication source node T ₂adopt joint-detection mode to obtain communication source node T simultaneously ₁the all modulation symbols that send on current subcarrier, in a M time slot, and obtain by communication source node T ₁and the space diversity gain that brings of many antennas on single via node, the space diversity between via node is converted to frequency diversity by circulation delay transmission means, and this kind of frequency diversity can be obtained by error correction coding;

Full rate distributed multi-antenna double-direction radio cooperating relay transmission method is by communication source node T ₂send to communication source node T ₁data by similar approach at T ₁place obtains.

4. full rate distributed multi-antenna double-direction radio cooperating relay transmission method according to claim 2, it is characterized in that: in the time that the antenna number of each via node assembling is even number, the mode that can adopt timesharing to transmit two orthogonal Space Time Coding of antenna Alamouti completes data retransmission, now, the detailed process of second stage can be described below:

$Y_2\left(p\right)=\sqrt{\frac{a}{2}}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{{\tilde{G}}_k\left(p\right)\mathrm{diag}(A_1,A_2,...,A_{M/2})e^{-j\frac{2p}{N_F}p{D}_k}\right\}+{\tilde{W}}_2\left(p\right)---(2$

$0)$

In formula ,] ^*represent the conjugation of scalar, A _mencoder matrix during for Alamouti orthogonal space, the reliability, the reduction receiving terminal that are conducive to improve forwarding data detect the error rate.Encoding characteristics during according to orthogonal space, to communication source node T ₂reception signal matrix modification after can obtain following mode:

$\begin{array}{c}\left[\begin{array}{c}{Y}_2^{2m-1}\left(p\right)\\ {\left[Y_2^{2m}\left(p\right)\right]}^{*}\end{array}\right]=\sqrt{\frac{a}{2}}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{{2P}_1}{F}_k^m\mathrm{diag}({\left[h_k^{2m-1}\left(p\right)\right]}^T,{\left[h_k^{2m}\left(p\right)\right]}^T)\left[\begin{array}{c}{\tilde{s}}_{2m-1}\left(p\right)\\ {\tilde{s}}_{2m}\left(p\right)\end{array}\right]\right\}\\ +\sqrt{\frac{a}{2}}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{\sqrt{{2P}_2}{F}_k^m\mathrm{diag}({\left[g_k^{2m-1}\left(p\right)\right]}^T,{\left[g_k^{2m}\left(p\right)\right]}^T)\left[\begin{array}{c}{\tilde{c}}_{2m-1}\left(p\right)\\ {\tilde{c}}_{2m}\left(p\right)\end{array}\right]\right\}\\ +\sqrt{\frac{a}{2}}\underset{k=1}{\overset{K}{\mathrm{\Σ}}}\left\{F_k^m\left[\begin{array}{c}{w}_k^p(2m-\mathrm{1,2}m-1)\\ w_k^p(2m,2m)\end{array}\right]\right\}+\begin{array}{}\end{array}\left[\begin{array}{c}{\tilde{W}}_2^{2m-1}\left(p\right)\\ {\left[{\tilde{W}}_2^{2m}\left(p\right)\right]}^{*}\end{array}\right]\end{array}---\left(22\right)$

$F_k^m=\left[\begin{array}{cc}{\tilde{g}}_k^{2m-1}\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}& {\tilde{g}}_k^{2m}\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\\ {\left[{\tilde{g}}_k^{2m}\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\right]}^{*}& -{\left[{\tilde{g}}_k^{2m-1}\left(p\right)e^{-j\frac{2p}{N_F}p{D}_k}\right]}^{*}\end{array}\right]---\left(23\right)$

For above-mentioned mode, can be at communication source node T ₂after eliminating self signal and disturbing, joint-detection is out from communication source node T ₁initial data;

By communication source node T ₂send to communication source node T ₁data by similar approach at T ₁place obtains.