CN100553187C - 时域同步正交频分复用系统的发射分集方法 - Google Patents

Abstract

本发明属于数字信息传输领域。步骤为：1)将频域输入序列按照子载波序号分为奇数和偶数子序列，长度均为N/2；2)将连续两帧的上述子序列分别作反离散傅里叶变换，将结果缓存；3)对缓存数据进行计算得到四个发射链路的时域信号；4)在四个发射链路的TDS-OFDM保护间隔内分别插入不同的PN序列作为帧头，将帧头和上一步得到的帧体分别组成各自发射链路的完整信号帧；5)将信号帧进行成形滤波、数模变换和前端处理，分别通过四个天线在预定的频道带宽中发射出去。本发明简单、快速、准确，保持了系统的传输效率，更适用于时间、频率双选择性信道，且支持“软失败”，增加了系统可靠性，同时也易于移植到其他系统中。

Description

时域同步正交频分复用系统的发射分集方法

技术领域

本发明属于数字信息传输技术领域，更具体地涉及一种时域同步正交频分复用(TimeDomain Synchronous OFDM，TDS-OFDM)系统中基于空时分组编码(Space-Time Block Code)的时频域联合发射分集方法。

背景技术

在复杂的无线环境中，周围的物体(如房屋、建筑物或树木等)对无线电波会起到反射的作用。这些障碍物会产生幅度衰减和相位延迟不同的反射波。如果发射一个调制信号，那么该发射信号的多个反射波就会从不同方向经过不同传播延迟到达接收天线。这些反射信号经过位于各处的接收机天线接收后，根据其随机相位的不同，对接收信号会起到加强或减弱的作用。由此会造成接收端信号的幅度变化，形成衰落。统计表明，在障碍物均匀的城市街道或森林环境中，信号包络起伏近似于满足Rayleigh分布，故多径快衰落又称为Rayleigh衰落。短期快衰落是由于收发信号双方的相对运动而产生：多径信号的存在造成时间扩散，从而引起传输信号的符号间干扰；而相对运动造成的多普勒效应会引起传输信号的相位迅速变化，在不同的测试环境下有不同的快衰落特性。

除了接收信号的瞬时值会出现快速Rayleigh衰落之外，场强中值也会出现缓慢变化。变化的原因主要有两个方面：一是由移动环境中的固定障碍物(如建筑物、山丘、森林等)的阴影效应引起的；二是由于气象条件的变化，导致大气相对介电常数的垂直梯度发生缓变，即电波折射系数随时间变化，从而多径传播到达固定接收点的信号时延也随之变化。这种由阴影效应和气象原因引起的信号变化，称为慢衰落。特别地，由气象原因引起的变化较小，通常忽略。慢衰落环境下的接收信号近似服从对数正态分布，变化幅度取决于障碍物状况、工作频率、变化速率、障碍物和接收机移动速度等。

在多径移动接收中，多径效应引起的时延扩展和多普勒效应引起的多普勒频率展宽同时存在，称为频率和时间选择性衰落信道。在这种信道下，接收信号的信噪比很不稳定，当信道处于深度衰落中时接收信噪比低，判决错误的概率就大，严重降低信号传输的可靠性。为了提高系统的抗衰落性能，可以采用各种信道均衡技术、正交频分复用(OFDM)多载波调制技术等。而分集技术是克服频率和时间选择性衰落的有效技术，它将相同信息经过几个不相关的衰落信道，然后对接收信号进行合成。因为几个信道同时处于深衰落的概率较低，因此可以达到平滑信道衰落，增加信噪比，改善接收机误码特性的目的。在数字电视地面广播网络中，由于分集技术降低了接收机的信噪比(SNR)门限要求，因此在同样的发射功率下，还可以扩展电视信号的覆盖范围。

传统的信号分集技术如时间分集、频率分集等，它们的作用实际上相当于信道编码中的重复编码加交织技术，虽然能够改善系统的误码性能，但由于相同信息的重复传送，往往要牺牲较大的的传输效率。另一个常用分集方式是采用多重天线进行空间分集，这种技术在发射端或接收端都可以实现，分别称为发射分集和接收分集。

其中，使用多个接收天线的接收分集是一种传统而有效的分集技术，它不需要牺牲传输效率，在接收端可以采用最大值切换、最大比例合并等简单方式完成多个接收信号的选择或合并，然后再按照常规方法进行译码和判决。2001年，法国Harris公司的研究人员进行了欧洲地面数字视频广播DVB-T的接收天线分集实验(Faria G.Mobile DVB-T using antennadiversity receivers.2001.Available：Http://www.broadcastpapers.com)，在各种复杂多径环境下，效果很好，测试样机的平均SNR门限下降约6dB，抗多普勒能力增加100％。该接收分集方案如图1所示，接收端采用两组独立的射频前端以及OFDM解调和信道估计模块，经过OFDM解调后，两路接收信号在第k个子载波上的样值分别为：

$\left\{\begin{array}{c}{R}_{\mathrm{Rx}1}\left(k\right)=H_{\mathrm{Rx}1}\left(k\right)X\left(k\right)+N_{\mathrm{Rx}1}\left(k\right)\\ R_{\mathrm{Rx}2}\left(k\right)=H_{\mathrm{Rx}2}\left(k\right)X\left(k\right)+N_{\mathrm{Rx}2}\left(k\right)\end{array}\right.$ (0≤k≤N-1)

其中，H_Rx1(k)和H_Rx2(k)分别为两条接收路径第k个子载波上的频率响应值，N_Rx1(k)和N_Rx2(k)则分别表示相应路径上的噪声。

假设两路接收信号经历了互不相关的信道衰落且信道估计结果正确，那么最大比率合并是最好的信号合并方式。将两组接收信号分别乘以其子载波频率响应值的共轭再相加，得到：

$R\left(k\right)=H_{\mathrm{Rx}1}^{*}\left(k\right)R_{\mathrm{Rx}1}\left(k\right)+H_{\mathrm{Rx}2}^{*}\left(k\right)R_{\mathrm{Rx}2}\left(k\right)$

$=({\left|H_{\mathrm{Rx}1}\left(k\right)\right|}^2+{\left|H_{\mathrm{Rx}2}\left(k\right)\right|}^2)X\left(k\right)+H_{\mathrm{Rx}1}^{*}\left(k\right)N_{\mathrm{Rx}1}\left(k\right)+H_{\mathrm{Rx}2}^{*}\left(k\right)N_{\mathrm{Rx}2}\left(k\right)$

其中，上式的合成信号与非分集系统得到的信号形式相同，可以直接用于译码和判决。显然，合并后的信号的信噪比高于两个支路信号的信噪比最大值，因此获得了分集增益。如果采用更多接收天线并按上述方式进行最大比率合并，还可以获得更大增益。

接收分集能够获得良好的效果，分集方法也很简单，但在应用于数字电视地面广播(DTTB)领域中时会受到一些限制。一是接收分集需要接收机有多套并行的射频前端处理，增加了接收机的成本和复杂度，这对于广播系统不合算；二是要使各路接收信号不相关，每两个接收天线的距离要为载波波长的10倍量级，在DTTB所处的VHF/UHF频段，此距离约为4～7m，这对于很多移动和便携式接收终端很难实现。相反地，对于发射机来说，上述限制都不成问题，因此，发射分集技术日益成为研究的热点。由于采用发射分集时，在多个发射天线间往往需要进行信号矢量的合理配置以尽量增加分集阶数(Diversity Order)，所以这项技术也被看成是放置在传统信道编码后的又一层“内码”，称为分集编码。其接收端的信号采样是多个发射信号的叠加，需要通过适当的处理方式进行分离和译码。

近年来已有很多发射分集的研究成果。在文献“Wittneben A.A new bandwidth efficienttransmit antenna modulation diversity scheme for linear digital modulation.in Proc.of IEEE ICC’93.Geneva，Switzerland：IEEE，1993.1630-1634”和“Winters J.The diversity gain of transmit diversityin wireless system with Rayleigh fading.in Proc.of IEEE ICC’94.New Orleans，LA：IEEE，1994.1121-1125”中采用的分集方式是将相同信号延时发射，形成一个“人为多径”并用类似Rake接收的方式加以合并。文献“Foschini G and Gans M.On limits of wireless communications in afading environment when using multiple antenna.Wireless Personal Communications，1998，6(3)：311-335”介绍了Blast系统中的空时分层码，它是将发射信号分成多路并分别进行传统的信道编码和交织。1998年，Tarokh在“Tarokh V，SeshadriN，and Calderbank A.Space-time codes forhigh data rate wireless communications：performance criterion and code construction.IEEE Trans.onInformation Theory，1998，44(2)：744-765”中介绍了空时格型编码(STTC)的概念，它将信道编码和天线分集联合进行设计，该方法在理论上可以获得最大限度的增益，但是需要改变整个发射系统的设计，而且即使在较少的发射天线数和低阶星座图的情况下，译码复杂度仍然很大。为了解决这个问题，Alamouti于1998年在其经典论文“Alamouti S.A simple transmitdiversity technique for wireless communications.IEEE Trans.on Select Areas in Communications，1998，16(8)：1451-1458”中提出了一种空时分组编码(STBC)方案，在两发射天线系统中应用，其编码构造和译码算法非常简单，同样可以获得分集增益。Tarokh等人在1999年将STBC推广至任意发射天线数的情况，对这种方案给出了理论分析和构造准则。上述这些STBC编码是基于正交(orthogonal)的结构提出的。在复数域上，假设具有k个符号(x₁，x₂，...，x_k)的p×n_T分组编码矩阵为G(x₁，x₂，...，x_k)，G中的元素满足：(1)每一项都是x_i、x_i ^*或其线性组合之一；(2)G^HG＝(|x₁|²+|x₂|²+...+|x_k|²)I_n，I_n为单位矩阵。此码字用于发射天线分集系统时，将连续k个输入符号按照G所示进行编码，矩阵的每一列都表示供一个发射天线发出的符号序列。由于编码矩阵G具有(2)所示的正交性质，使得接收端译码时可以将每个符号分开，对x_i分别进行译码，这样就大大降低了接收端的译码复杂度。

定义STBC码字的编码效率为R＝k/p，其中，k为输入符号数，p为编码延时。对于Alamouti提出的2天线发射分集系统，R＝1。但是当发射天线数目(n_T)多于2时，在Tarokh等人的文献中已经证明，基于正交性质设计的STBC码的编码效率小于1，即p＞k，因而损失了有效传输码率。这就意味着要使采用发射分集后的系统仍保持原有的单发射机系统的传输码率，必须增大原有系统所占用的带宽，这对带宽固定的系统(如DTTB)来说是一个很大的缺点，因此引入了基于准正交(quasi-orthogonal)性质设计的STBC码(“Jafarkhani H，A quasi-orthogonalspace-time block code.IEEE Trans.on Communications，2001，49(1)：1-4”)。准正交STBC码放松了正交性条件的约束，虽然降低了一些分集增益，但可以使得编码效率达到1，即使在带宽固定的系统中也能有效应用。在复数域上，仍然假设分组编码矩阵为G，此时p＝k。以Jafarkhani在文献中提出的4发射天线情况为例，对于准正交设计G中的元素，上述性质(1)保留，(2)改为

$G^{H}G=\left[\begin{array}{cccc}a& 0& 0& b\\ 0& a& -b& 0\\ 0& -b& a& 0\\ b& 0& 0& a\end{array}\right]$

其中， $a=\underset{i=1}{\overset{4}{\mathrm{\Σ}}}{\left|x_i\right|}^2,$ b＝2Re(x₁x₄ ^*-x₂x₃ ^*)。因此，对于矩阵中的各列向量(V_i，i＝1，2，3，4)来说，可以分为2组：(V₁，V₄)和(V₂，V₃)。每组内的向量之间不正交，但不同组中的向量之间是正交的。在接收端进行最大似然(ML)译码时，可以将判决公式按照组分为两部分，这样运算量仍然很小。

由于STBC(正交设计或准正交设计)编码方式具有快速译码的优点，所以很快得到了广泛研究，并迅速由平衰落单载波信道扩展到频率选择性衰落的OFDM信道中，形成了基于OFDM系统的空时分组码。由于在OFDM调制技术中涉及到时域和频域两种信号，因此STBC编码可以分别在时域和频域中进行。

如果在时域输入符号中进行STBC编码(STC-OFDM)，即将连续的k个OFDM符号按照一定的编码格式G所示进行编码(“Lee K and Williams D.A space-time coded transmitter diversitytechnique for frequency selective fading channels.in Proc.IEEE Sensor Array and MultichannelSignal Processing Workshop.Cambridge，MA：IEEE，2000.149-152”)。此时，为了保持传输矩阵G正交(或准正交)的性质，必须假设信道是准静态的，即在相邻的连续k个OFDM符号时间内信道保持不变。这个假设在快衰落信道下会产生很大误差，因而像在DTTB的移动接收等情况下并不适用。

STBC编码也可以在频域中进行(SFC-OFDM)，将相邻的k个子载波上的数据编码(“LeeK and Williams D.A space-frequency transmitter diversity technique for OFDM systems.in Proc.IEEE GLOBECOM’00.San Francisco，CA：IEEE，2000，1473-1477”)。这样，为了保持传输矩阵G正交(或准正交)的性质，需要假设相邻的k个子载波上的频率响应相同。虽然SFC-OFDM可以适用于快衰落信道，但同样由于信道假设带来的误差，在频率选择性衰落信道中并不适用。而实际存在的信道环境大多是双选择性(时间和频率)的，因此STC-OFDM和SFC-OFDM系统在实际应用中都会带来较大误差。

根据数字电视地面广播的工程特点，TDS-OFDM系统采用的发射分集方法应遵循以下设计原则：(1)不牺牲传输效率。分集方案应保证系统保持原有的信息吞吐能力，不能引入冗余，也就是说每个射频信道不会因为增加了分集编码而降低信息传输率，因而在设计多于2发射天线的分集系统时，要采用上述准正交的设计结构；(2)尽量不改动原有发射系统。使分集方案成为发射系统的一个可选“配件”，这样运营者可以根据本地信道特点和覆盖情况等客观实际决定是否需要采用分集，也最大限度降低了发射机成本；(3)接收机算法改动小，复杂度增加少。增加发射天线不可避免地会增加接收机的功能模块，如信道估计和信号合并器等；(4)支持“软失败(Soft Failure)”。所谓支持“软失败”，是指发射天线间的信号矢量配置应当保证，当一个接收路径因某种原因失效时，另一个接收路径仍能使系统正常接收，仅是牺牲了1/2的平均接收功率。因此，这种发射分集技术实际上还增加了系统的可靠性。

目前世界上数字电视地面广播传输标准主要有三种：美国的ATSC(高级电视系统委员会Advanced Television Systems Committee)、欧洲的DVB-T(地面数字视频地面广播Digital VideoTerrestrial Broadcasting-Terrestrial)和日本的ISDB-T(地面综合业务数字广播Integrated ServiceDigital Broadcasting-Terrestrial，ISDB-T)。我国自1994年起，也开始了高清晰度电视的研究工作。在此背景下，清华大学提出了地面数字多媒体广播(Digital Multimedia Broadcasting forTerrestrial，DMB-T)传输协议。

清华DMB-T中采用的时域同步正交频分复用(TDS-OFDM)调制属于多载波调制技术，但与欧洲DVB-T采用的编码正交频分复用(COFDM)技术不同，在TDS-OFDM系统中没有插入频域导频信号，而是在OFDM的保护间隔中以时域的方式插入了伪随机(PN)序列，用于帧同步、频率同步、定时同步、信道传输特性估计和跟踪相位噪声等。

为了实现快速和稳定的同步，清华大学提出的TDS-OFDM传输系统采用了分级帧结构。帧结构的基本单元称为信号帧，如图2所示。200/225个信号帧定义为一个帧群，512个帧群定义为一个超帧。帧结构的顶层称为日帧，由超帧组成。帧群中的每一个信号帧有唯一的帧号，它被编码在帧头的PN序列中。

TDS-OFDM传输系统的信号帧使用时域同步的正交频分复用调制，或者称为以PN序列为保护间隔的正交频分复用调制。一个信号帧由帧同步和帧体两部分组成，它们具有相同的基带符号率7.56MS/s(1/T)。一个信号帧可以作为一个正交频分复用(OFDM)块。一个OFDM块进一步分成一个保护间隔和一个离散傅里叶逆变换(IDFT)块。对于TDS-OFDM来说，帧同步PN序列作为OFDM的保护间隔，而帧体作为IDFT块，如图3所示。

关于DMB-T、TDS-OFDM的相关情况详见授权号为00123597.4名为“地面数字多媒体电视广播系统”、授权号为01115520.5名为“时域同步正交频分复用调制方法”、授权号为ZL01130659.9名为“地面数字多媒体电视广播系统中的帧同步产生方法”，以及授权号为01124144.6名为“正交频分复用调制系统中保护间隔的填充方法”等清华大学申请的中国发明专利。

为了在TDS-OFDM系统中实现发射天线分集，必须满足下列条件：

(1)发射天线之间保持足够的距离，以使到达接收机的各条传输信道统计独立；

(2)接收端在进行信道估计时能够准确估计出当前时刻每个信道的信道信息；

(3)通过合适的方法把接收到的多路信号分离出来，使其互不相关，然后将分离出的多路信号合并，获得最大的信噪比。

针对上述背景，本发明提出了一种针对TDS-OFDM系统的基于空时分组编码的时频域联合发射分集方法。

发明内容

本发明的目的在于提供一种时域同步正交频分复用(Time Domain Synchronous OFDM，TDS-OFDM)系统中基于空时分组编码(Space Time Block Code，STBC)的一种时频域联合发射分集方法。

本发明针对数字电视地面广播系统中的发射分集问题，提出了一种基于准正交结构STBC编码的时频域联合发射分集方法。本发明提出的发射分集方法框图如图4所示。

本发明所述的时域同步正交频分复用，即TDS-OFDM，时频域联合发射分集方法，其特征在于，所述方法是一种基于空时分组编码的时频域联合发射分集方法，它在专用数字集成电路中是按照以下步骤依次实现的：

步骤1.记频域输入序列为X(k，l)，其中k表示子载波序号，0≤k≤N-1，N为OFDM系统中的子载波数，l表示信号帧序号，将X(k，l)按照子载波序号分为奇数子序列X_o(k，l)和偶数子序列X_e(k，l)，它们的长度均为N/2；

步骤2.将X_o(k，l)和X_e(k，l)分别作N/2点反离散傅里叶变换，得到的时域序列为x₁ ^to(n，l)和x₁ ^te(n，l)；

步骤3.将连续两帧输入数据进行反离散傅里叶变换后的结果x₁ ^to(n，l)、x₁ ^te(n，l)、x₁ ^to(n，l+1)和x₁ ^te(n，l+1)存入到缓存中；

步骤4.然后将缓存中的数据按照下面所述四种情况进行不同运算，分别得到用于四个天线发射所需的时域信号：

(a)对于第一个发射天线Tx1，时域信号x_Tx1(n，l)、x_Tx1(n，l+1)为

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}1}(n,l)=[x_1^{\mathrm{te}}(n,l)+x_1^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}1}(n+N/2,l)=[x_1^{\mathrm{te}}(n,l)-x_1^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}1}(n,l+1)=[x_1^{\mathrm{te}}(n,l+1)+x_1^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\\ x_{\mathrm{Tx}1}(n+N/2,l+1)=[x_1^{\mathrm{te}}(n,l+1)-x_1^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\end{array}\right.$ 0≤n ≤N/2-1；

其中， $W_N^{-n}=e^{j\frac{2\mathrm{\π}}{N}n},$ N为OFDM系统中的子载波数；

(b)对于第二个发射天线Tx2，先将缓存中的数据经过空频编码后得到

$\left\{\begin{array}{c}{x}_2^{\mathrm{te}}(n,l)={x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l),\\ x_2^{\mathrm{to}}(n,l)={-x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l),\\ x_2^{\mathrm{te}}(n,l+1)={x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l+1),\\ x_2^{\mathrm{to}}(n,l+1)={-x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l+1),\end{array}\right.$ 0≤n ≤N/2-1；

其中，*表示复数共轭运算，(n)_N/2表示对n取模N/2运算，然后得到时域信号x_Tx2(n，l)、x_Tx2(n，l+1)为

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}2}(n,l)=[x_2^{\mathrm{te}}(n,l)+x_2^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}2}(n+N/2,l)=[x_2^{\mathrm{te}}(n,l)-x_2^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}2}(n,l+1)=[x_2^{\mathrm{te}}(n,l+1)+x_2^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\\ x_{\mathrm{Tx}2}(n+N/2,l+1)=[x_2^{\mathrm{te}}(n,l+1)-x_2^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\end{array}\right.$ 0≤n≤N/2-1；

(c)对于第三个发射天线Tx3，先将缓存中的数据经过空时编码后得到

$\left\{\begin{array}{c}{x}_3^{\mathrm{te}}(n,l)={x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l+1),\\ x_3^{\mathrm{to}}(n,l)={x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l+1),\\ x_3^{\mathrm{te}}(n,l+1)=-{x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l),\\ x_3^{\mathrm{to}}(n,l+1)={-x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l),\end{array}\right.$ 0≤n≤N/2-1；

然后得到时域信号x_Tx3(n，l)、x_Tx3(n，l+1)为

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}3}(n,l)=[x_3^{\mathrm{te}}(n,l)+x_3^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}3}(n+N/2,l)=[x_3^{\mathrm{te}}(n,l)-x_3^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}3}(n,l+1)=[x_3^{\mathrm{te}}(n,l+1)+x_3^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\\ x_{\mathrm{Tx}3}(n+N/2,l+1)=[x_3^{\mathrm{te}}(n,l+1)-x_3^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\end{array}\right.$ 0≤n≤N/2-1；

(d)对于第四个发射天线Tx4，缓存中的数据先经过空时编码，得到如(c)中所示的结果x₃ ^to(n，l)、x₃ ^te(n，l)、x₃ ^to(n，l+1)和x₃ ^te(n，l+1)，然后再经过空频编码得到

$\left\{\begin{array}{c}{x}_4^{\mathrm{te}}(n,l)={x_3^{\mathrm{to}}}^{*}({(-n)}_{N/2},l),\\ x_4^{\mathrm{to}}(n,l)={-x_3^{\mathrm{te}}}^{*}({(-n)}_{N/2},l),\\ x_4^{\mathrm{te}}(n,l+1)={x_3^{\mathrm{to}}}^{*}({(-n)}_{N/2},l+1),\\ x_4^{\mathrm{to}}(n,l+1)={-x_3^{\mathrm{te}}}^{*}({(-n)}_{N/2},l+1),\end{array}\right.$ 0≤n≤N/2-1；

最后得到时域信号x_Tx4(n，l)、x_Tx4(n，l+1)为

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}4}(n,l)=[x_4^{\mathrm{te}}(n,l)+x_4^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}4}(n+N/2,l)=[x_4^{\mathrm{te}}(n,l)-x_4^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}4}(n,l+1)=[x_4^{\mathrm{te}}(n,l+1)+x_4^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\\ x_{\mathrm{Tx}4}(n+N/2,l+1)=[x_4^{\mathrm{te}}(n,l+1)-x_4^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\end{array}\right.$ 0≤n≤N/2-1；

步骤5.根据时域同步正交频分复用系统信号帧帧头的长度420或945，生成四个不同的相应长度的PN序列；

步骤6.按照时域同步正交频分复用系统的信道帧结构，在发射天线Tx1、Tx2、Tx3和Tx4四个链路的TDS-OFDM保护间隔内分别插入上述不同的PN序列作为帧头，将帧头PN序列和步骤(4)得到的帧体x_Tx1(n，l)、x_Tx2(n，l)、x_Tx3(n，l)、x_Tx4(n，l)分别组成四个发射链路各自完整的信号帧；

步骤7.将上述完整的TDS-OFDM信号进行成形滤波和数模变换处理，然后经过包含频率上变换和功放在内的前端处理，最后分别通过天线Tx1、Tx2、Tx3和Tx4在预定的频道带宽中发射出去，完成发射天线分集。

本发明所述的时域同步正交频分复用，即TDS-OFDM，时频域联合发射分集方法，其特征在于，所述的时频域联合编码的等效空时分组编码结构用矩阵G表示为：

$G=\left[\begin{array}{cccc}X(2k,l)& X^{*}(2k+1,l)& X^{*}(2k,l+1)& X(2k+1,l+1)\\ X(2k+1,l)& {-X}^{*}(2k,l)& X^{*}(2k+1,l+1)& -X(2k,l+1)\\ X(2k,l+1)& X^{*}(2k+1,l+1)& {-X}^{*}(2k,l)& -X(2k+1,l)\\ X(2k+1,l+1)& -X^{*}(2k,l+1)& {-X}^{*}(2k+1,l)& X(2k,l)\end{array}\right],$ (0≤k≤N/2-1)；

其中，X(k，l)为频域输入序列，k表示子载波序号，l表示OFDM帧序号，这是一个准正交(quosi-orthogonal)设计的空时分组编码(STBC)结构。

本发明所述的时域同步正交频分复用，即TDS-OFDM，时频域联合发射分集方法，其特征在于，可以保持原有的单发射机链路(Tx1)基本不变，而只是在其中加入了缓存器，使得网络构造灵活。其他三个发射链路(Tx2、Tx3、Tx4)的时域信号只须对缓存中的数据进行简单处理即可得到，因此，在传输一帧OFDM信号的时间内只须做一次IDFT运算，运算复杂度很小。

本发明所述的时域同步正交频分复用，即TDS-OFDM，时频域联合发射分集方法，其特征在于，分集系统工作的鲁棒性强，即如果部分传输链路工作不正常，不用对原有系统进行任何修改，接收端仍然可以正常译码，并且误码性能至少不低于单发射机系统下的情况。在本发明中，如果任一个发射链路出现故障，那么其余三个发射链路仍可以组成准正交STBC结构，从而获得分集增益。当只有两个发射链路工作正常时，可以分为三种情况：1)Tx1和Tx3(或Tx2和Tx4)，可以组成一个2天线的正交STC-OFDM结构；2)Tx1和Tx2(或Tx3和Tx4)，可以组成一个2天线的正交SFC-OFDM结构；3)Tx1和Tx4(或Tx2和Tx3)。在1)和2)情况下，系统仍能获得分集增益。

同时，本发明所提出的时频域联合发射分集方法不失一般性，可以很方便地移植到其他多载波DTTB系统。本发明所述的发射分集方案并不排斥接收分集，在本发明中可以引入多个接收天线进行接收分集。

下面我们对本发明中提出的时频域联合发射分集方法的原理和性能进行分析。发射分集方法的应用系统结构框图如图5所示。

假设频域输入信号序列为X(k，l)，其中k表示子载波序号(0≤k≤N-1，N为OFDM系统中的子载波数)，l表示信号帧序号。在OFDM系统中，经过反离散Fourier变换(IDFT)，得到的时域信号x_Tx1(n，l)为

$x_{\mathrm{Tx}1}(n,l)=\frac{1}{N}\underset{k=0}{\overset{N}{\mathrm{\Σ}}}X(k,l)W_N^{-\mathrm{nk}}$ (0≤n≤N-1)

其中， $W_N^k=e^{-j\frac{2\mathrm{\π}}{N}k}.$

将X(k，l)按照子载波序号分为奇数子序列X_o(k，l)和偶数子序列X_e(k，l)，它们的长度均为N/2，若记X_o(k，l)和X_e(k，l)做N/2点IDFT变换的结果为x₁ ^to(n，l)和x₁ ^te(n，l)，则时域信号x_Tx1(n，l)可以改变形式表示为：

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}1}(n,l)=[x_1^{\mathrm{te}}(n,l)+x_1^{\mathrm{to}}(n,l)W_N^{-n}]/2\\ x_{\mathrm{Tx}1}(n+N/2,l)=[x_1^{\mathrm{te}}(n,l)-x_1^{\mathrm{to}}(n,l)W_N^{-n}]/2\end{array}\right.$ (0≤n≤N/2-1)

这个结果x_Tx1(n，l)就是用于第一个发射天线Tx1的时域信号。同时将x_Tx1(n，l)、x_Tx2(n，l)、x_Tx3(n，l)、x_Tx4(n，l)保存到缓存中。

同样地，用于第二个发射天线Tx2的时域信号x_Tx2(n，l)和x_Tx2(n，l+1)可以表示为：

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}2}(n,l)=[x_2^{\mathrm{te}}(n,l)+x_2^{\mathrm{to}}(n,l)W_N^{-n}]/2\\ x_{\mathrm{Tx}2}(n+N/2,l)=[x_2^{\mathrm{te}}(n,l)-x_2^{\mathrm{to}}(n,l)W_N^{-n}]/2\\ x_{\mathrm{Tx}2}(n,l+1)=[x_2^{\mathrm{te}}(n,l+1)+x_2^{\mathrm{to}}(n,l+1)W_N^{-n}]/2\\ x_{\mathrm{Tx}2}(n+N/2,l+1)=[x_2^{\mathrm{te}}(n,l+1)-x_2^{\mathrm{to}}(n,l+1)W_N^{-n}]/2\end{array}\right.$ (0≤n≤N/2-1)

公式中的x₂ ^te(n，l)，x₂ ^to(n，l)，x₂ ^te(n，l+1)和x₂ ^to(n，l+1)是由缓存中的信号经过空频编码(SFC)后得到的：

$\left\{\begin{array}{c}{x}_2^{\mathrm{te}}(n,l)={x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l)\\ x_2^{\mathrm{to}}(n,l)={-x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l)\\ x_2^{\mathrm{te}}(n,l+1)={x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l+1)\\ x_2^{\mathrm{to}}(n,l+1)={-x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l+1)\end{array}\right.$ (0≤n ≤N/2-1)

其中，^*表示复数共轭运算，(n)_N/2表示对n取模N/2运算。根据离散Fourier变换(DFT)的运算性质，可得x₂ ^te(n，l)，x₂ ^to(n，l)，x₂ ^te(n，l+1)和x₂ ^to(n，l+1)的N/2点DFF变换结果可以写为：

$\left\{\begin{array}{c}{x}_2^{\mathrm{te}}(n,l)\stackrel{\mathrm{DFT}(N/2)}{\→}{X}_o^{*}(k,l)\\ x_2^{\mathrm{to}}(n,l)\stackrel{\mathrm{DFT}(N/2)}{\→}-X_e^{*}(k,l)\\ x_2^{\mathrm{te}}(n,l+1)\stackrel{\mathrm{DFT}(N/2)}{\→}{X}_o^{*}(k,l+1)\\ x_2^{\mathrm{to}}(n,l+1)\stackrel{\mathrm{DFT}(N/2)}{\→}-X_e^{*}(k,l+1)\end{array}\right.$ (0≤n，k≤N/2-1)

因此，对于第二个发射链路，其等效频域输入信号为

X_Tx2＝[X^*(1)，-X^*(0)...X^*(2k+1)，-X^*(2k)...X^*(N-1)，-X^*(N-2)](0≤k≤N/2-1)

上式所示的性质与OFDM帧号无关，所以在括号中省略了第二项l。

对于第三个发射链路Tx3，要先将缓存中的信号经过空时编码(STC)，可得

$\left\{\begin{array}{c}{x}_3^{\mathrm{te}}(n,l)={x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l+1)\\ x_3^{\mathrm{to}}(n,l)={x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l+1)\\ x_3^{\mathrm{te}}(n,l+1)=-{x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l)\\ x_3^{\mathrm{to}}(n,l+1)={-x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l)\end{array}\right.$ (0≤n≤N/2-1)

同理，应用DFT变换的运算性质，x₃ ^te(n，l)、x₃ ^to(n，l)、x₃ ^te(n，l+1)和x₃ ^to(n，l+1)的N/2点DFT变换的结果可以表示为

$\left\{\begin{array}{c}{x}_3^{\mathrm{te}}(n,l)\stackrel{\mathrm{DFT}(N/2)}{\→}{X_e}^{*}(k,l+1)\\ x_3^{\mathrm{to}}(n,l)\stackrel{\mathrm{DFT}(N/2)}{\→}{X_o}^{*}(k,l+1)\\ x_3^{\mathrm{te}}(n,l+1)\stackrel{\mathrm{DFT}(N/2)}{\→}-{X_e}^{*}(k,l)\\ x_3^{\mathrm{to}}(n,l+1)\stackrel{\mathrm{DFT}(N/2)}{\→}-{X_o}^{*}(k,l)\end{array}\right.$ (0≤n，k≤N/2-1)

因此，第三个发射链路的等效频域输入信号为

$\left\{\begin{array}{c}{X}_{\mathrm{Tx}3}(k,l)=X^{*}(k,l+1)\\ X_{\mathrm{Tx}3}(k,l+1)={-X}^{*}(k,l)\end{array}\right.$ (0≤k≤N-1)

进一步，可得时域信号x_Tx3(n，l)、x_Tx3(n，l+1)为

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}3}(n,l)=[x_3^{\mathrm{te}}(n,l)+x_3^{\mathrm{to}}(n,l)W_N^{-n}]/2\\ x_{\mathrm{Tx}3}(n+N/2,l)=[x_3^{\mathrm{te}}(n,l)-x_3^{\mathrm{to}}(n,l)W_N^{-n}]/2\\ x_{\mathrm{Tx}3}(n,l+1)=[x_3^{\mathrm{te}}(n,l+1)+x_3^{\mathrm{to}}(n,l+1)W_N^{-n}]/2\\ x_{\mathrm{Tx}3}(n+N/2,l+1)=[x_3^{\mathrm{te}}(n,l+1)-x_3^{\mathrm{to}}(n,l+1)W_N^{-n}]/2\end{array}\right.$ (0≤n≤N/2-1)

然后，考虑第四个发射链路，将经过空时编码后得到的数据x₃ ^te(n，l)、x₃ ^to(n，l)、x₃ ^te(n，l+1)和x₃ ^to(n，l+1)再经过空频编码可得

$\left\{\begin{array}{c}{x}_4^{\mathrm{te}}(n,l)={x_3^{\mathrm{to}}}^{*}({(-n)}_{N/2},l)\\ x_4^{\mathrm{to}}(n,l)={-x_3^{\mathrm{te}}}^{*}({(-n)}_{N/2},l)\\ x_4^{\mathrm{te}}(n,l+1)={x_3^{\mathrm{to}}}^{*}({(-n)}_{N/2},l+1)\\ x_4^{\mathrm{to}}(n,l+1)={-x_3^{\mathrm{te}}}^{*}({(-n)}_{N/2},l+1)\end{array}\right.$ (0≤n≤N/2-1)

同样根据DFT变换运算的性质，x₄ ^te(n，l)、x₄ ^to(n，l)、x₄ ^te(n，l+1)和x₄ ^to(n，l+1)的N/2点DFT变换的结果可以表示为

$\left\{\begin{array}{c}{x}_4^{\mathrm{te}}(n,l)\stackrel{\mathrm{DFT}(N/2)}{\→}{X_{\mathrm{Tx}3}^o}^{*}(k,l)=X_o(k,l+1)\\ x_2^{\mathrm{to}}(n,l)\stackrel{\mathrm{DFT}(N/2)}{\→}{{-X}_{\mathrm{Tx}3}^e}^{*}(k,l)={-X}_e(k,l+1)\\ x_2^{\mathrm{te}}(n,l+1)\stackrel{\mathrm{DFT}(N/2)}{\→}{X_{\mathrm{Tx}3}^o}^{*}(k,l+1)={-X}_o(k,l)\\ x_2^{\mathrm{to}}(n,l+1)\stackrel{\mathrm{DFT}(N/2)}{\→}{{-X}_{\mathrm{Tx}3}^e}^{*}(k,l+1)=X_e(k,l)\end{array}\right.$ (0≤n，k≤N/2-1)

式中，X_Tx3 ^o(k，l)、X_Tx3 ^e(k，l)分别为x₃ ^to(n，l)、x₃ ^te(n，l)的N/2点DFT变换的结果。上式即为第四个发射链路的等效频域输入信号。

最后，得到时域信号x_Tx4(n，l)、x_Tx4(n，l+1)为

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}4}(n,l)=[x_4^{\mathrm{te}}(n,l)+x_4^{\mathrm{to}}(n,l)W_N^{-n}]/2\\ x_{\mathrm{Tx}4}(n+N/2,l)=[x_4^{\mathrm{te}}(n,l)-x_4^{\mathrm{to}}(n,l)W_N^{-n}]/2\\ x_{\mathrm{Tx}4}(n,l+1)=[x_4^{\mathrm{te}}(n,l+1)+x_4^{\mathrm{to}}(n,l+1)W_N^{-n}]/2\\ x_{\mathrm{Tx}4}(n+N/2,l+1)=[x_4^{\mathrm{te}}(n,l+1)-x_4^{\mathrm{to}}(n,l+1)W_N^{-n}]/2\end{array}\right.$ (0≤n≤N/2-1)

从上面的分析中我们可以得到四个发射链路的等效频域输入编码矩阵为

$G=\left[\begin{array}{cccc}X(2k,l)& X^{*}(2k+1,l)& X^{*}(2k,l+1)& X(2k+1,l+1)\\ X(2k+1,l)& {-X}^{*}(2k,l)& X^{*}(2k+1,l+1)& -X(2k,l+1)\\ X(2k,l+1)& X^{*}(2k+1,l+1)& {-X}^{*}(2k,l)& -X(2k+1,l)\\ X(2k+1,l+1)& -X^{*}(2k,l+1)& {-X}^{*}(2k+1,l)& X(2k,l)\end{array}\right]$ (0≤k≤N/2-1)

这是一个准正交(quosi-orthogonal)设计的STBC结构。

不失一般性，假设在接收端只有一个接收机，以使得分析简单。为了在接收端分别得到四个发射链路的信道估计结果，在TDS-OFDM系统的保护间隔内要分别插入不同的PN序列作为帧头。在下面的分析中，假设接收机可以得到每条发射链路的准确信道信息，各个发射天线的信号经过的信道是不相关的，并且是加性信道，用一个接收天线收到的信号是各个发射信号的叠加。经过OFDM解调后，接收信号在第2k、2k+1个子载波上的样值可以表示为：

$\left\{\begin{array}{c}{R}_l\left(2k\right)=H_{1,l}\left(2k\right)X(2k,l)+H_{2,l}\left(2k\right)X^{*}(2k+1,l)+\\ H_{3,l}\left(2k\right)X^{*}(2k,l+1)+H_{4,l}\left(2k\right)X(2k+1,l+1)+{\mathrm{\η}}_l\left(2k\right)\\ R_l(2k+1)=H_{1,l}(2k+1)X(2k+1,l)-H_{2,l}(2k+1)X^{*}(2k,l)+\\ H_{3,l}(2k+1)X^{*}(2k+1,l+1)-H_{4,l}(2k+1)X(2k,l+1)+{\mathrm{\η}}_l(2k+1)\\ R_{l+1}\left(2k\right)=H_{1,l+1}\left(2k\right)X(2k,l+1)+H_{2,l+1}\left(2k\right)X^{*}(2k+1,l+1)-\\ H_{3,l+1}\left(2k\right)X^{*}(2k,l)-H_{4,l+1}\left(2k\right)X(2k+1,l)+{\mathrm{\η}}_{l+1}\left(2k\right)\\ R_{l+1}(2k+1)=H_{1,l+1}(2k+1)X(2k+1,l+1)-H_{2,l+1}(2k+1)X^{*}(2k,l+1)-\\ H_{3,l+1}(2k+1)X^{*}(2k+1,l)+H_{4,l+1}(2k+1)X(2k,l)+{\mathrm{\η}}_{l+1}(2k+1)\end{array}\right.$

(0≤k≤N/2-1)

式中，H_i，l表示在第l个时间段内(一个时间段表示传输一个完整的OFDM帧所用的时间，下同)第i个传输链路的复值信号响应向量，η_l表示在第l个时间段内复值加性白高斯噪声(AWGN)向量。

假设在相邻的两个时间段以及相邻的两个子载波之间的信道响应近似相同，即

H_i(k)＝H_i，l(2k)≈H_i，l(2k+1)≈H_i，l+1(2k)≈H_i，l+1(2k+1)

(i＝1，2，3，40≤k≤N/2-1)

则接收信号可以简化表示为

$\left[\begin{array}{c}{R}_l\left(2k\right)\\ {R_l}^{*}(2k+1)\\ {R_{l+1}}^{*}\left(2k\right)\\ R_{l+1}(2k+1)\end{array}\right]=\left[\begin{array}{cccc}{H}_1\left(k\right)& H_2\left(k\right)& H_3\left(k\right)& H_4\left(k\right)\\ {{-H}_2}^{*}\left(k\right)& {H_1}^{*}\left(k\right)& {{-H}_4}^{*}\left(k\right)& {H_3}^{*}\left(k\right)\\ {{-H}_3}^{*}\left(k\right)& {{-H}_4}^{*}\left(k\right)& {H_1}^{*}\left(k\right)& {H_2}^{*}\left(k\right)\\ H_4\left(k\right)& {-H}_3\left(k\right)& {-H}_2\left(k\right)& H_1\left(k\right)\end{array}\right]\left[\begin{array}{c}X(2k,l)\\ X^{*}(2k+1,l)\\ X^{*}(2k,l+1)\\ X(2k+1,l+1)\end{array}\right]+\left[\begin{array}{c}{\mathrm{\η}}_l\left(2k\right)\\ {{\mathrm{\η}}_l}^{*}(2k+1)\\ {{\mathrm{\η}}_{l+1}}^{*}\left(2k\right)\\ {\mathrm{\η}}_{l+1}(2k+1)\end{array}\right]$

(0≤k≤N/2-1)

记信道响应矩阵为H(k)，即

$H\left(k\right)=\left[\begin{array}{cccc}{H}_1\left(k\right)& H_2\left(k\right)& H_3\left(k\right)& H_4\left(k\right)\\ {{-H}_2}^{*}\left(k\right)& {H_1}^{*}\left(k\right)& {{-H}_4}^{*}\left(k\right)& {H_3}^{*}\left(k\right)\\ {{-H}_3}^{*}\left(k\right)& {{-H}_4}^{*}\left(k\right)& {H_1}^{*}\left(k\right)& {H_2}^{*}\left(k\right)\\ H_4\left(k\right)& {-H}_3\left(k\right)& {-H}_2\left(k\right)& H_1\left(k\right)\end{array}\right]$ (0≤k≤N/2-1)

将接收信号表达式的两端均左乘H(k)的Hermit变换矩阵H^H(k)，可得最终的信号估计值X′(2k，l)、X′(2k+1，l)、X′(2k，l+1)和X′(2k+1，l+1)为：

$\left[\begin{array}{c}{X}^{\′}(2k,l)\\ X^{{\′}^{*}}(2k+1,l)\\ X^{{\′}^{*}}(2k,l+1)\\ X^{\′}(2k+1,l+1)\end{array}\right]\left[\begin{array}{cccc}a\left(k\right)& 0& 0& b\left(k\right)\\ 0& a\left(k\right)& -b\left(k\right)& 0\\ 0& -b\left(k\right)& a\left(k\right)& 0\\ b\left(k\right)& 0& 0& a\left(k\right)\end{array}\right]\left[\begin{array}{c}X(2k,l)\\ X^{*}(2k+1,l)\\ X^{*}(2k,l+1)\\ X(2k+1,l+1)\end{array}\right]+\left[\begin{array}{c}{\mathrm{\η}}_l^{\′}\left(2k\right)\\ {\mathrm{\η}}_l^{\′*}(2k+1)\\ {{\mathrm{\η}}_{l+1}^{\′}}^{*}\left(2k\right)\\ {\mathrm{\η}}_{l+1}^{\′}(2k+1)\end{array}\right]---(*)$

(0≤k≤N/2-1)

其中， $a\left(k\right)=\underset{i=1}{\overset{4}{\mathrm{\Σ}}}{\left|H_i\left(k\right)\right|}^2,$ b(k)＝2Re(H₁(k)H₄ ^*(k)-H₂(k)H₃ ^*(k))，η′_l仍然是复值加性白高斯噪声(AWGN)向量。

可见，在接收端进行最大似然(ML)译码时，可以将判决式分为两部分：(X(2k，l)，X(2k+1，l+1))和(X(2k+1，l)，X(2k，l+1))，这样可以使得译码的运算量大大减小。从上面的分析中还可以看到，采用本发明所提出的发射分集结构虽然分集增益会有所减小，但保持编码效率为1。

根据本发明提出的时频域联合发射分集方法结构(图4)，利用DFT变换的运算性质，在上面的分析中可以看到，四个发射链路的时域信号x_Tx1(n，l)、x_Tx2(n，l)、x_Tx3(n，l)、x_Tx4(n，l)都可以由缓存中的数据x₁ ^to(n，l)、x₁ ^te(n，l)、x₁ ^to(n，l+1)和x₁ ^te(n，l+1)通过简单的复数乘法和加法计算出来，也即在一个OFDM帧的时间段内平均做一次N点DFT运算以及3N/2次复数乘法和3N次复数加法。而在传统的STC-OFDM(“Lee K and Williams D.Aspace-time coded transmitterdiversity technique for frequency selective fading channels.in Proc.IEEE Sensor Array andMultichannel Signal Processing Workshop.Cambridge，MA：IEEE，2000.149-152”)和SFC-OFDM(“Lee K and Williams D.A space-frequency transmitter diversity technique for OFDM systems.inProc.IEEE GLOBECOM’00.San Francisco，CA：IEEE，2000，1473-1477”)系统中，在一个OFDM帧的时间段内平均做四次N点DFT运算。因此，当OFDM系统的子载波数N足够大时，本发明提出的方法的运算量约为原有方法的1/4。特别地，以TDS-OFDM系统为例，N＝3780，应用此方法运算量可以减少约70％。

一般地，在传统的STC-OFDM系统中，即将连续的k个OFDM符号按照一定的编码格式G所示进行编码。此时，为了保持传输矩阵G正交(或准正交)的性质，必须假设信道是准静态的，即在相邻的连续k个OFDM符号时间内信道保持不变。这个假设在快衰落信道下会产生很大误差，因而像在DTTB的移动接收等情况下并不适用。而在传统的SFC-OFDM系统中，将相邻的k个子载波上的数据编码。这样，为了保持传输矩阵G正交(或准正交)的性质，需要假设相邻的k个子载波上的频率响应相同。同样由于信道假设带来的误差，在频率选择性衰落信道中并不适用。而实际存在的信道环境大多是双选择性(时间和频率)的，因此STC-OFDM和SFC-OFDM系统在实际应用中都会带来较大误差。本发明提出的基于STBC的时频域联合发射分集方法只须假设相邻的两个时间段以及相邻的两个子载波之间的信道响应近似相同，因此更加适用于实际的时间、频率双选择性信道。

基于上述分析，对本发明所提出的TDS-OFDM系统中基于空时分组编码(STBC)的一种时频域联合发射分集方法进行了计算机仿真，仿真的系统结构与图5相同。在仿真中我们采用表1和表2所示的两种信道模型A和B。其中模型A的多径延时较短，而且回波强度较弱；模型B中引入了具有长延时的强多径(多径6)。

表1分集传输信道仿真模型A

表2分集传输信道仿真模型B

仿真中采用QPSK星座图和保护间隔为数据长度1/9的3780点TDS-OFDM系统，信道编码采用2/3码率的卷积码。并且假设接收机可以得到每条发射链路的准确信道信息，各个发射天线的信号经过的信道不相关。图6～9给出了不分集、传统的STC-OFDM、SFC-OFDM和本发明所提出的分集方法在不同信道情况下的误比特率(BER)仿真结果。为了使仿真结果具有可比性，图中的信噪比(Signal to Noise Ratio，SNR)是以接收天线为准的，也就是说，发射分集方法中每个发射天线的功率仅为单发射天线方法(不分集)中发射功率的1/4。

图6为短时延慢衰落信道，采用信道模型A，不加多普勒频移效应。三种发射分集系统相比不分集系统都有很明显的增益，系统性能近似相同。

图7为短时延快衰落信道下，仍采用信道模型A，最大多普勒频移f_d＝50Hz。由于信道假设误差的影响，STC-OFDM系统已不能提供分集增益，它的性能甚至比不分级系统还要差。而SFC-OFDM系统和本发明所提的分集方案表现出对多普勒频移的不敏感，仍然保持了较高增益。

图8为长时延慢衰落信道，采用信道模型B，不加多普勒频移效应。可以看出，由于信道假设误差的影响，SFC-OFDM系统已不能提供分集增益。而SFC-OFDM系统和本发明所提的分集方案仍可以提供较高增益。

图9所示信道既有长延时和强多径，又有快衰落，采用信道模型B，最大多普勒频移f_d＝50Hz。同样由于信道假设误差的影响，STC-OFDM和SFC-OFDM系统都已无法提供分集增益，而本发明所提的分集方案仍可以获得较高增益。

本发明所提出的TDS-OFDM系统中基于空时分组编码(STBC)的时频域联合发射分集方法分集系统工作的鲁棒性强，支持“软失败(Soft Failure)”，即如果部分传输链路工作不正常，不用对原有系统进行任何修改，接收端仍然可以正常译码，并且误码性能至少不低于单发射机系统下的情况。这是由于，如果某个发射链路无法正常工作，那么在最后的判决公式(*)中，只须将相应的信道响应矩阵中的元素置零即可，仍可以正常进行ML译码。在本发明中，如果任一个发射链路出现故障，那么其余三个发射链路仍可以组成准正交STBC结构，从而获得分集增益。当只有两个发射链路工作正常时，可以分为三种情况：1)Tx1和Tx3(或Tx2和Tx4)，可以组成一个2天线的正交STC-OFDM结构；2)Tx1和Tx2(或Tx3和Tx4)，可以组成一个2天线的正交SFC-OFDM结构；3)Tx1和Tx4(或Tx2和Tx3)。在1)和2)情况下，系统仍能获得分集增益。图10给出了系统中存在不同数目的正常工作的发射机的情况下的计算机仿真结果，仿真时采用信道模型A，不加多普勒频移效应。仿真结果很好的验证了上述分析结果。

附图说明

图1接收分集方案框图

图2为TDS-OFDM系统分级帧结构。

图3为TDS-OFDM系统信号帧结构图。

图4为本发明提出的时频域联合发射分集方法框图。

图5为本发明提出的发射分集方法的应用系统结构框图。

图6为本发明对信道模型A的对比仿真结果(f_d＝0Hz)。

图7为本发明对信道模型A的对比仿真结果(f_d＝50Hz)。

图8为本发明对信道模型B的对比仿真结果(f_d＝0Hz)。

图9为本发明对信道模型B的对比仿真结果(f_d＝50Hz)。

图10为本发明提出的分集系统中存在不同数目的正常工作的发射机的情况下的对比仿真结果。

具体实施方式

见图4。本系统中的频域输入序列是经过调制映射(星座图映射)后的复数信号，其可以是非编码信号，也可以是经过编码后的信号。频域输入序列首先按照子载波序号分为奇数子序列和偶数子序列，它们的长度均为N/2。然后将连续两帧输入数据的共四个N/2长的子序列分别作N/2点反离散Fourier变换，其结果x₁ ^to(n，l)、x₁ ^te(n，l)、x₁ ^to(n，l+1)和x₁ ^te(n，l+1)存入到缓存中存入缓存中。利用反离散Fourier变换的运算性质，分下列四种情况对缓存中的数据进行简单的复数乘法和加法计算，分别得到四个发射链路的时域信号：

1)对于天线Tx1：

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}1}(n,l)=[x_1^{\mathrm{te}}(n,l)+x_1^{\mathrm{to}}(n,l)W_N^{-n}]/2\\ x_{\mathrm{Tx}1}(n+N/2,l)=[x_1^{\mathrm{te}}(n,l)-x_1^{\mathrm{to}}(n,l)W_N^{-n}]/2\\ x_{\mathrm{Tx}1}(n,l+1)=[x_1^{\mathrm{te}}(n,l+1)+x_1^{\mathrm{to}}(n,l+1)W_N^{-n}]/2\\ x_{\mathrm{Tx}1}(n+N/2,l+1)=[x_1^{\mathrm{te}}(n,l+1)-x_1^{\mathrm{to}}(n,l+1)W_N^{-n}]/2\end{array}\right.$ (0≤n≤N/2-1)

2)对于天线Tx2，先经过空频编码：

$\left\{\begin{array}{c}{x}_2^{\mathrm{te}}(n,l)={x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l)\\ x_2^{\mathrm{to}}(n,l)={-x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l)\\ x_2^{\mathrm{te}}(n,l+1)={x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l+1)\\ x_2^{\mathrm{to}}(n,l+1)={-x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l+1)\end{array}\right.$ (0≤n≤N/2-1)

然后计算得到时域信号：

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}2}(n,l)=[x_2^{\mathrm{te}}(n,l)+x_2^{\mathrm{to}}(n,l)W_N^{-n}]/2\\ x_{\mathrm{Tx}2}(n+N/2,l)=[x_2^{\mathrm{te}}(n,l)-x_2^{\mathrm{to}}(n,l)W_N^{-n}]/2\\ x_{\mathrm{Tx}2}(n,l+1)=[x_2^{\mathrm{te}}(n,l+1)+x_2^{\mathrm{to}}(n,l+1)W_N^{-n}]/2\\ x_{\mathrm{Tx}2}(n+N/2,l+1)=[x_2^{\mathrm{te}}(n,l+1)-x_2^{\mathrm{to}}(n,l+1)W_N^{-n}]/2\end{array}\right.$ (0≤n≤N/2-1)

3)对于天线Tx3，先经过空时编码：

$\left\{\begin{array}{c}{x}_3^{\mathrm{te}}(n,l)={x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l+1)\\ x_3^{\mathrm{to}}(n,l)={x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l+1)\\ x_3^{\mathrm{te}}(n,l+1)=-{x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l)\\ x_3^{\mathrm{to}}(n,l+1)={-x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l)\end{array}\right.$ (0≤n≤N/2-1)

然后计算得到时域信号：

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}3}(n,l)=[x_3^{\mathrm{te}}(n,l)+x_3^{\mathrm{to}}(n,l)W_N^{-n}]/2\\ x_{\mathrm{Tx}3}(n+N/2,l)=[x_3^{\mathrm{te}}(n,l)-x_3^{\mathrm{to}}(n,l)W_N^{-n}]/2\\ x_{\mathrm{Tx}3}(n,l+1)=[x_3^{\mathrm{te}}(n,l+1)+x_3^{\mathrm{to}}(n,l+1)W_N^{-n}]/2\\ x_{\mathrm{Tx}3}(n+N/2,l+1)=[x_3^{\mathrm{te}}(n,l+1)-x_3^{\mathrm{to}}(n,l+1)W_N^{-n}]/2\end{array}\right.$ (0≤n≤N/2-1)

4)对于天线Tx4，将经过空时编码后的数据再经过空频编码：

$\left\{\begin{array}{c}{x}_4^{\mathrm{te}}(n,l)={x_3^{\mathrm{to}}}^{*}({(-n)}_{N/2},l)\\ x_4^{\mathrm{to}}(n,l)={-x_3^{\mathrm{te}}}^{*}({(-n)}_{N/2},l)\\ x_4^{\mathrm{te}}(n,l+1)={x_3^{\mathrm{to}}}^{*}({(-n)}_{N/2},l+1)\\ x_4^{\mathrm{to}}(n,l+1)={-x_3^{\mathrm{te}}}^{*}({(-n)}_{N/2},l+1)\end{array}\right.$ (0≤n≤N/2-1)

然后计算得到时域信号：

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}4}(n,l)=[x_4^{\mathrm{te}}(n,l)+x_4^{\mathrm{to}}(n,l)W_N^{-n}]/2\\ x_{\mathrm{Tx}4}(n+N/2,l)=[x_4^{\mathrm{te}}(n,l)-x_4^{\mathrm{to}}(n,l)W_N^{-n}]/2\\ x_{\mathrm{Tx}4}(n,l+1)=[x_4^{\mathrm{te}}(n,l+1)+x_4^{\mathrm{to}}(n,l+1)W_N^{-n}]/2\\ x_{\mathrm{Tx}4}(n+N/2,l+1)=[x_4^{\mathrm{te}}(n,l+1)-x_4^{\mathrm{to}}(n,l+1)W_N^{-n}]/2\end{array}\right.$ (0≤n≤N/2-1)

按照TDS-OFDM系统的信道帧结构，在四个发射链路的TDS-OFDM保护间隔内分别插入不同的PN序列作为帧头，将帧头PN序列和上述步骤中得到的帧体分别组成四个发射链路各自完整的信号帧。然后再将完整的TDS-OFDM信号进行成形滤波和数模变换处理，然后经过频率上变换和功放等前端处理，最后分别通过四个天线在预定的频道带宽中发射出去，完成发射天线分集。

在接收端，假设在相邻的两个时间段以及相邻的两个子载波之间的信道响应近似相同，经过OFDM解调后，接收信号在第2k、2k+1个子载波上的样值可以表示为：

$\left[\begin{array}{c}{R}_l\left(2k\right)\\ {R_l}^{*}(2k+1)\\ {R_{l+1}}^{*}\left(2k\right)\\ R_{l+1}(2k+1)\end{array}\right]=\left[\begin{array}{cccc}{H}_1\left(k\right)& H_2\left(k\right)& H_3\left(k\right)& H_4\left(k\right)\\ {{-H}_2}^{*}\left(k\right)& {H_1}^{*}\left(k\right)& {{-H}_4}^{*}\left(k\right)& {H_3}^{*}\left(k\right)\\ {{-H}_3}^{*}\left(k\right)& {{-H}_4}^{*}\left(k\right)& {H_1}^{*}\left(k\right)& {H_2}^{*}\left(k\right)\\ H_4\left(k\right)& {-H}_3\left(k\right)& {-H}_2\left(k\right)& H_1\left(k\right)\end{array}\right]\left[\begin{array}{c}X(2k,l)\\ X^{*}(2k+1,l)\\ X^{*}(2k,l+1)\\ X(2k+1,l+1)\end{array}\right]+\left[\begin{array}{c}{\mathrm{\η}}_l\left(2k\right)\\ {{\mathrm{\η}}_l}^{*}(2k+1)\\ {{\mathrm{\η}}_{l+1}}^{*}\left(2k\right)\\ {\mathrm{\η}}_{l+1}(2k+1)\end{array}\right]$

(0≤k≤N/2-1)

记信道响应矩阵为H(k)，即

$H\left(k\right)=\left[\begin{array}{cccc}{H}_1\left(k\right)& H_2\left(k\right)& H_3\left(k\right)& H_4\left(k\right)\\ {{-H}_2}^{*}\left(k\right)& {H_1}^{*}\left(k\right)& {{-H}_4}^{*}\left(k\right)& {H_3}^{*}\left(k\right)\\ {{-H}_3}^{*}\left(k\right)& {{-H}_4}^{*}\left(k\right)& {H_1}^{*}\left(k\right)& {H_2}^{*}\left(k\right)\\ H_4\left(k\right)& {-H}_3\left(k\right)& {-H}_2\left(k\right)& H_1\left(k\right)\end{array}\right]$ (0≤k≤N/2-1)

将接收信号表达式的两端均左乘H(k)的Hermit变换矩阵H^H(k)，可得最终的信号估计值X′(2k，l)、X′(2k+1，l)、X′(2k，l+1)和X′(2k+1，l+1)为：

$\left[\begin{array}{c}{X}^{\′}(2k,l)\\ X^{{\′}^{*}}(2k+1,l)\\ X^{{\′}^{*}}(2k,l+1)\\ X^{\′}(2k+1,l+1)\end{array}\right]\left[\begin{array}{cccc}a\left(k\right)& 0& 0& b\left(k\right)\\ 0& a\left(k\right)& -b\left(k\right)& 0\\ 0& -b\left(k\right)& a\left(k\right)& 0\\ b\left(k\right)& 0& 0& a\left(k\right)\end{array}\right]\left[\begin{array}{c}X(2k,l)\\ X^{*}(2k+1,l)\\ X^{*}(2k,l+1)\\ X(2k+1,l+1)\end{array}\right]+\left[\begin{array}{c}{\mathrm{\η}}_l^{\′}\left(2k\right)\\ {\mathrm{\η}}_l^{\′*}(2k+1)\\ {{\mathrm{\η}}_{l+1}^{\′}}^{*}\left(2k\right)\\ {\mathrm{\η}}_{l+1}^{\′}(2k+1)\end{array}\right]$

(0≤k≤N/2-1)

其中， $a\left(k\right)=\underset{i=1}{\overset{4}{\mathrm{\Σ}}}{\left|H_i\left(k\right)\right|}^2,$ b(k)＝2Re(H₁(k)H₄ ^*(k)-H₂(k)H₃ ^*(k))，η′_l仍然是复值加性白高斯噪声(AWGN)向量。在进行最大似然(ML)译码时，将判决式分为两部分：(X(2k，l)，X(2k+1，l+1))和(X(2k+1，l)，X(2k，l+1))，这样可以使得译码的运算量大大减小。

上面本发明的一般性具体实施例进行了说明，但本发明并不限制于上述实施例，在不脱离本申请的权利要求的精神和范围情况下，本领域的技术人员可作出各种修改或改型。

Claims (3)

1、时域同步正交频分复用系统的发射分集方法，其特征在于，所述方法是一种基于空时分组编码的时频域联合发射分集方法，它在专用数字集成电路中是按照以下步骤依次实现的：

步骤1.记频域输入序列为X(k，l)，其中k表示子载波序号，0≤k≤N-1，N为OFDM系统中的子载波数，l表示信号帧序号，将X(k，l)按照子载波序号分为奇数子序列X_o(k，l)和偶数子序列X_e(k，l)，它们的长度均为N/2；

步骤2.将X_o(k，l)和X_e(k，l)分别作N/2点反离散傅里叶变换，得到的时域序列为x₁ ^to(n，l)和x₁ ^te(n，l)；

步骤3.将连续两帧输入数据进行反离散傅里叶变换后的结果x₁ ^to(n，l)、x₁ ^te(n，l)、x₁ ^to(n，l+1)和x₁ ^te(n，l+1)存入到缓存中；

步骤4.然后将缓存中的数据按照下面所述四种情况进行不同运算，分别得到用于四个天线发射所需的时域信号：

(a)对于第一个发射天线Tx1，时域信号x_Tx1(n，l)、x_Tx1(n，l+1)为

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}1}(n,l)=[x_1^{\mathrm{te}}(n,l)+x_1^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}1}(n+N/2,l)=[x_1^{\mathrm{te}}(n,l)-x_1^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}1}(n,l+1)=[x_1^{\mathrm{te}}(n,l+1)+x_1^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\\ x_{\mathrm{Tx}1}(n+N/2,l+1)=[x_1^{\mathrm{te}}(n,l+1)-x_1^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\end{array}\right.$ 0≤n≤N/2-1；

其中， $W_N^{-n}=e^{j\frac{2\mathrm{\π}}{N}n},$ N为OFDM系统中的子载波数；

(b)对于第二个发射天线Tx2，先将缓存中的数据经过空频编码后得到

$\left\{\begin{array}{c}{x}_2^{\mathrm{te}}(n,l)={x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l),\\ x_2^{\mathrm{to}}(n,l)=-{x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l),\\ x_2^{\mathrm{te}}(n,l+1)={x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l+1),\\ x_2^{\mathrm{to}}(n,l+1)=-{x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l+1),\end{array}\right.$ 0≤n≤N/2-1；

其中，*表示复数共轭运算，(n)_N/2表示对n取模N/2运算，然后得到时域信号x_Tx2(n，l)、x_Tx2(n，l+1)为

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}2}(n,l)=[x_2^{\mathrm{te}}(n,l)+x_2^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}2}(n+N/2,l)=[x_2^{\mathrm{te}}(n,l)-x_2^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}2}(n,l+1)=[x_2^{\mathrm{te}}(n,l+1)+x_2^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\\ x_{\mathrm{Tx}2}(n+N/2,l+1)=[x_2^{\mathrm{te}}(n,l+1)-x_2^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\end{array}\right.$ 0≤n≤N/2-1；

(c)对于第三个发射天线Tx3，先将缓存中的数据经过空时编码后得到

$\left\{\begin{array}{c}{x}_3^{\mathrm{te}}(n,l)={x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l+1),\\ x_3^{\mathrm{to}}(n,l)={x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l+1),\\ x_3^{\mathrm{te}}(n,l+1)={-x_1^{\mathrm{te}}}^{*}({(-n)}_{N/2},l),\\ x_3^{\mathrm{to}}(n,l+1)=-{x_1^{\mathrm{to}}}^{*}({(-n)}_{N/2},l),\end{array}\right.$ 0≤n≤N/2-1；

然后得到时域信号x_Tx3(n，l)、x_Tx3(n，l+1)为

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}3}(n,l)=[x_3^{\mathrm{te}}(n,l)+x_3^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}3}(n+N/2,l)=[x_3^{\mathrm{te}}(n,l)-x_3^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}3}(n,l+1)=[x_3^{\mathrm{te}}(n,l+1)+x_3^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\\ x_{\mathrm{Tx}3}(n+N/2,l+1)=[x_3^{\mathrm{te}}(n,l+1)-x_3^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\end{array}\right.$ 0≤n≤N/2-1；

(d)对于第四个发射天线Tx4，缓存中的数据先经过空时编码，得到如(c)中所示的结果x₃ ^to(n，l)、x₃ ^te(n，l)、x₃ ^to(n，l+1)和x₃ ^te(n，l+1)，然后再经过空频编码得到

$\left\{\begin{array}{c}{x}_4^{\mathrm{te}}(n,l)={x_3^{\mathrm{to}}}^{*}({(-n)}_{N/2},l),\\ x_4^{\mathrm{to}}(n,l)=-{x_3^{\mathrm{te}}}^{*}({(-n)}_{N/2},l),\\ x_4^{\mathrm{te}}(n,l+1)={x_3^{\mathrm{to}}}^{*}({(-n)}_{N/2},l+1),\\ x_4^{\mathrm{to}}(n,l+1)=-{x_3^{\mathrm{te}}}^{*}({(-n)}_{N/2},l+1),\end{array}\right.$ 0≤n≤N/2-1；

最后得到时域信号x_Tx4(n，l)、x_Tx4(n，l+1)为

$\left\{\begin{array}{c}{x}_{\mathrm{Tx}4}(n,l)=[x_4^{\mathrm{te}}(n,l)+x_4^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}4}(n+N/2,l)=[x_4^{\mathrm{te}}(n,l)-x_4^{\mathrm{to}}(n,l)W_N^{-n}]/2,\\ x_{\mathrm{Tx}4}(n,l+1)=[x_4^{\mathrm{te}}(n,l+1)+x_4^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\\ x_{\mathrm{Tx}4}(n+N/2,l+1)=[x_4^{\mathrm{te}}(n,l+1)-x_4^{\mathrm{to}}(n,l+1)W_N^{-n}]/2,\end{array}\right.$ 0≤n≤N/2-1；

这四个发射天线发出的信号组成准正交设计的空时分组编码结构；

步骤5.根据时域同步正交频分复用系统信号帧帧头的长度420或945，生成四个不同的相应长度的PN序列；

步骤6.按照时域同步正交频分复用系统的信道帧结构，在发射天线Tx1、Tx2、Tx3和Tx4四个链路的TDS-OFDM保护间隔内分别插入上述不同的PN序列作为帧头，将帧头PN序列和步骤(4)得到的帧体x_Tx1(n，l)、x_Tx2(n，l)、x_Tx3(n，l)、x_Tx4(n，l)分别组成四个发射链路各自完整的信号帧；

步骤7.将上述完整的TDS-OFDM信号进行成形滤波和数模变换处理，然后经过包含频率上变换和功放在内的前端处理，最后分别通过天线Tx1、Tx2、Tx3和Tx4在预定的频道带宽中发射出去，完成发射天线分集。

2、根据权利要求1所述的时域同步正交频分复用系统的发射分集方法，其特征在于，当任意一个发射天线的链路出现故障时，其余三个发射天线的链路仍组成准正交设计的空时分组编码结构。

3、根据权利要求1所述的时域同步正交频分复用系统的发射分集方法，其特征在于，当只有两个发射天线的链路工作正常时，则有：

(1)Tx1和Tx3，或者Tx2和Tx4，分别组成一个2天线的正交空时编码的OFDM结构；

(2)Tx1和Tx2，或者Tx3和Tx4，分别组成一个2天线的正交空频编码的OFDM结构。

Images