简介面向无线通信与信号处理初学者的空时分组码STBC编码译码实现示例围绕STBC编译码流程与误码率性能评估展开适合用于课程实验、仿真入门或算法验证前的快速理解。压缩包共3个文件均为.m脚本整体约1KB体量轻便便于在MATLAB环境中直接查看和运行其中包含编译码实现与误码率曲线绘制相关脚本可帮助读者把抽象的分组码原理对应到可执行代码与可视化结果。目前已有2322人学习下载说明其入门参考价值较受关注。通过这份材料读者可以了解STBC从编码、译码到BER曲线对比的基本实现思路并据此复现实验、观察不同条件下的性能变化为后续扩展仿真与排错提供起点。1. 从一条对不上的 BER 曲线说起STBC 到底在发射端做了什么把空时分组码STBC的 BER 曲线跑出来形状对了、位置却整体偏 3 dB是这套仿真里最常见的翻车点。STBC 解决的是发射端拿不到信道信息时怎么用多根天线换来分集增益——两个时隙发两个符号接收端一次共轭合并原本要做二维联合搜索的检测就退化成两个标量判决。它适合基站不动、终端乱跑的控制信道和上行链路也适合当 MIMO 入门的第一份能跑的代码码矩阵只有两行两列每一步都能手算验证。后面把编码矩阵、合并公式、Eb/N0 归一化和五类排查记录拆开配一份能直接跑的 Python 链路。2. Alamouti 编码矩阵两个时隙、两根天线怎么排符号Alamouti 码之所以被写进教科书不是因为它的增益最高而是因为它在“发射端完全不知道信道”这个约束下用最少的维度换到了满分集。理解它只需要盯住一个 2×2 的复数矩阵但真正写代码时第一个会绊住人的问题是索引怎么排行是天线的还是时间是天线的负号和共轭加在哪一侧。2.1 时隙与天线两个维度先把索引约定固定下来Alamouti 码块占两个连续的符号周期跨越两根发射天线。约定矩阵的行是时间、列是天线一个码块写成天线 1天线 2时隙 t0s1s2时隙 t1-s2*s1*两个时隙送出两个独立符号所以码率为 1。这一点常被忽略很多入门资料先讲重复编码两个时隙重复同一个符号码率 1/2分集阶数是 2 但频谱效率砍半Alamouti 用同样的两个时隙装了两个符号分集和速率同时保住这才是它当年被广泛采用的原因。功率分配是第二个必须提前定死的事。总发射功率受限的系统中两根天线各自只承担一半能量所以发送符号要乘 1/√Nt。这一步在仿真里极容易被漏掉漏掉的直接后果是曲线看着更“漂亮”但和理论值对不上而且换一根天线的配置时结果不可比。发射端完全不需要信道信息——这正是 STBC 相对波束赋形、预编码这类方案的定位差异。波束赋形要靠下行信道估计反馈STBC 靠的是码字本身的正交结构接收端用信道做合并就行。选型时可以这么判断如果发射端能稳定拿到信道状态信息空分复用加预编码的吞吐更高如果拿不到或者只想在恶劣信道下保一个稳健的低速通道STBC 是更省心的选择。2.2 正交性从哪来Gram 矩阵必须是标量乘单位阵把上面的矩阵记作 S两行分别是 $\mathbf{r}_0(s_1, s_2)$ 和 $\mathbf{r}_1(-s_2^, s_1^)$。计算行内积$\langle \mathbf{r}_0, \mathbf{r}_1 \rangle s_1 \cdot \overline{(-s_2^)} s_2 \cdot \overline{(s_1^)} -s_1 s_2 s_2 s_1 0$两行正交且各自范数都是 $|s_1|^2|s_2|^2$。于是 $S^H S (|s_1|^2|s_2|^2) I_2$。这个等式是整个方案的支点正因为 Gram 矩阵是标量乘单位阵接收端的两路信号在被合并之后符号 $s_1$ 和 $s_2$ 之间不会互相串扰二维联合最大似然搜索才被拆成两个独立的一维搜索。负号和共轭不是随便加的它必须与“第二个时隙把两个符号的动位置换一下”配合。写成 $(s_1, s_2)$ 和 $(s_2^, -s_1^)$ 也成立只是天线列顺序反了本质等价。但写成 $(s_1, s_2)$ 和 $(s_2, s_1)$ 就完蛋——两行内积变成 $2s_1 s_2$正交性没了译码端会看到一对互相干扰的量BER 曲线在高 SNR 处直接出现平台。正交结构的存在性有硬上限这点决定了 STBC 的适用边界。Tarokh 等人证明过复星座下码率为 1 的正交空时分组码只对两根发射天线存在三根和四根天线最高只能做到 3/4 码率需要往码矩阵里填零五个时隙发三个或四个符号。天线数再往上正交性只能靠牺牲码率换所以大规模天线系统里基本见不到 OSTBC走的是空分复用那条路。做链路选型时记住一句STBC 的分集增益随天线数涨但码率和检测复杂度也一起涨2×2 是甜点。2.3 写一个可复用的编码函数并做形状自检下面这段把码块结构固化成函数返回三维数组(时隙, 块索引, 天线索引)。这种排法在后面加接收天线维度时最省事广播不会出错。import numpy as np def alamouti_encode(sym, nt2, normalizeTrue): sym: (N,) 复数符号序列N 必须为偶数 返回 tx: (2, N//2, nt) —— [时隙, 块索引, 天线索引] sym np.asarray(sym, dtypenp.complex128) if sym.size % 2 ! 0: raise ValueError(符号数必须为偶数否则凑不满一个 Alamouti 块) if nt ! 2: raise ValueError(Alamouti 只对 2 根发射天线成立) s1 sym[0::2] # 每个块的第一个符号 s2 sym[1::2] # 每个块的第二个符号 tx np.empty((2, s1.size, 2), dtypenp.complex128) tx[0, :, 0] s1 # t0, 天线1 tx[0, :, 1] s2 # t0, 天线2 tx[1, :, 0] -np.conj(s2) # t1, 天线1 tx[1, :, 1] np.conj(s1) # t1, 天线2 return tx / np.sqrt(2) if normalize else tx逻辑上有三处值得单独说。sym[0::2]和sym[1::2]把序列拆成奇偶两路这对应“一个码块吃两个符号”的节奏后面译码端按同样的规则拼回去就能对齐tx的第一维是时隙而不是天线因为接收端求和是在时隙上做接收采样、在天线上做信道相乘这个顺序和物理过程一致最后除以√2是总功率归一把两根天线的辐射总能量拉回 1如果整条链路统一用“每根天线满功率”的口径把normalizeFalse打开即可但全链路必须一致。写完先别急着跑 BER做一个形状和正交性自检几行代码就能挡住一大半低级错误def check_orthogonality(sym, tol1e-10): tx alamouti_encode(sym, normalizeFalse) # (2, nblk, 2) S tx[:, 0, :] # 取第一个块的 2x2 gram S.conj().T S target (np.abs(sym[0])**2 np.abs(sym[1])**2) * np.eye(2) assert np.allclose(gram, target, atoltol), f正交性破坏:\n{gram} check_orthogonality(np.array([12j, -0.50.3j]))S.conj().T S就是 Gram 矩阵如果它不是对角阵说明共轭或负号放错了位置如果它是对角阵但两个对角元素不相等说明两个时隙的符号位置没对称。这一步跑通发射端的活儿就基本干净了。3. 接收端合并与 ML 译码把两路信号压成一个标量发射端把两个符号摊在了两根天线、两个时隙上接收端拿到的是一堆混在一起的观测。合并做的事就是用信道系数当权重把属于 $s_1$ 的那部分能量从两个时隙里全部抽出来同时把 $s_2$ 的干扰压掉。这一章的每一步都能手推推完再写代码就不会对着公式猜。3.1 平坦瑞利信道与准静态假设先把信道模型钉死。用平坦衰落、准静态的复数高斯信道$h_{n} \sim \mathcal{CN}(0,1)$且在一个码块的两个时隙内保持不变。准静态这个假设不是形式上的偷懒它直接决定了合并公式能不能成立——推导里默认 $h_1$、$h_2$ 在两个时隙里是同一组值如果它们每个时隙都重新生成合并出来的量就会带上一个无法消除的交叉项。工程上这个假设合不合理看相干时间。两个符号周期的跨度通常远小于信道相干时间所以准静态是站得住的但在高速移动场景下多普勒扩展会让相干时间掉到几个符号量级这时候 STBC 会出现误码平台靠加天线也压不下去。这类平台不是实现 bug是模型失效排查时要能区分开。另外这里的信道是“发射天线 n 到接收天线 r”的复增益多个接收天线之间假设独立同分布。真做硬件验证时天线间相关性强分集增益会打折仿真里可以加一个相关系数矩阵乘进去看退化幅度。3.2 共轭合并公式的推导与代码落地记接收天线上两个时隙的观测为 $r_0, r_1$噪声为 $n_0, n_1$$r_0 h_1 s_1 h_2 s_2 n_0$ $r_1 -h_1 s_2^* h_2 s_1^* n_1$要凑出只含 $s_1$ 的量思路是把 $r_0$ 乘 $\overline{h_1}$、把 $r_1$ 取共轭后乘 $h_2$$\overline{h_1} r_0 h_2 \overline{r_1} (|h_1|^2|h_2|^2) s_1 \overline{h_1} n_0 h_2 \overline{n_1}$交叉项正好抵消。同理要凑 $s_2$ 就用 $\overline{h_2} r_0 - h_1 \overline{r_1}$得到 $(|h_1|^2|h_2|^2) s_2$ 加噪声。合并增益记为 $g |h_1|^2|h_2|^2$它是一个实数正标量这一条极其关键——正因为它没有相位判决时不需要做相位补偿直接比距离就行。def alamouti_combine(rx, h): rx: (2, nblk, nr) 两个时隙的接收样本 h : (nblk, 2, nr) 信道nr 为接收天线数 返回 y1, y2: (nblk,) 合并后的两个软输出 g : (nblk,) 合并增益实数 r0, r1 rx[0], rx[1] # (nblk, nr) h0, h1 h[:, 0, :], h[:, 1, :] # (nblk, nr) y1 np.sum(np.conj(h0) * r0 h1 * np.conj(r1), axis1) y2 np.sum(np.conj(h1) * r0 - h0 * np.conj(r1), axis1) g np.sum(np.abs(h0) ** 2 np.abs(h1) ** p a hrefhttps://download.csdn.net/download/yuanclkk/10411105 stylecolor:#ec7500;font-size:14px; 本文还有配套的精品资源点击获取 /a img altmenu-r.4af5f7ec.gif srchttps://csdnimg.cn/release/wenkucmsfe/public/img/menu-r.4af5f7ec.gif stylewidth:16px;margin-left:4px;vertical-align:text-bottom;cursor:text; /p
简介面向无线通信与信号处理初学者的空时分组码STBC编码译码实现示例围绕STBC编译码流程与误码率性能评估展开适合用于课程实验、仿真入门或算法验证前的快速理解。压缩包共3个文件均为.m脚本整体约1KB体量轻便便于在MATLAB环境中直接查看和运行其中包含编译码实现与误码率曲线绘制相关脚本可帮助读者把抽象的分组码原理对应到可执行代码与可视化结果。目前已有2322人学习下载说明其入门参考价值较受关注。通过这份材料读者可以了解STBC从编码、译码到BER曲线对比的基本实现思路并据此复现实验、观察不同条件下的性能变化为后续扩展仿真与排错提供起点。1. 从一条对不上的 BER 曲线说起STBC 到底在发射端做了什么把空时分组码STBC的 BER 曲线跑出来形状对了、位置却整体偏 3 dB是这套仿真里最常见的翻车点。STBC 解决的是发射端拿不到信道信息时怎么用多根天线换来分集增益——两个时隙发两个符号接收端一次共轭合并原本要做二维联合搜索的检测就退化成两个标量判决。它适合基站不动、终端乱跑的控制信道和上行链路也适合当 MIMO 入门的第一份能跑的代码码矩阵只有两行两列每一步都能手算验证。后面把编码矩阵、合并公式、Eb/N0 归一化和五类排查记录拆开配一份能直接跑的 Python 链路。2. Alamouti 编码矩阵两个时隙、两根天线怎么排符号Alamouti 码之所以被写进教科书不是因为它的增益最高而是因为它在“发射端完全不知道信道”这个约束下用最少的维度换到了满分集。理解它只需要盯住一个 2×2 的复数矩阵但真正写代码时第一个会绊住人的问题是索引怎么排行是天线的还是时间是天线的负号和共轭加在哪一侧。2.1 时隙与天线两个维度先把索引约定固定下来Alamouti 码块占两个连续的符号周期跨越两根发射天线。约定矩阵的行是时间、列是天线一个码块写成天线 1天线 2时隙 t0s1s2时隙 t1-s2*s1*两个时隙送出两个独立符号所以码率为 1。这一点常被忽略很多入门资料先讲重复编码两个时隙重复同一个符号码率 1/2分集阶数是 2 但频谱效率砍半Alamouti 用同样的两个时隙装了两个符号分集和速率同时保住这才是它当年被广泛采用的原因。功率分配是第二个必须提前定死的事。总发射功率受限的系统中两根天线各自只承担一半能量所以发送符号要乘 1/√Nt。这一步在仿真里极容易被漏掉漏掉的直接后果是曲线看着更“漂亮”但和理论值对不上而且换一根天线的配置时结果不可比。发射端完全不需要信道信息——这正是 STBC 相对波束赋形、预编码这类方案的定位差异。波束赋形要靠下行信道估计反馈STBC 靠的是码字本身的正交结构接收端用信道做合并就行。选型时可以这么判断如果发射端能稳定拿到信道状态信息空分复用加预编码的吞吐更高如果拿不到或者只想在恶劣信道下保一个稳健的低速通道STBC 是更省心的选择。2.2 正交性从哪来Gram 矩阵必须是标量乘单位阵把上面的矩阵记作 S两行分别是 $\mathbf{r}_0(s_1, s_2)$ 和 $\mathbf{r}_1(-s_2^, s_1^)$。计算行内积$\langle \mathbf{r}_0, \mathbf{r}_1 \rangle s_1 \cdot \overline{(-s_2^)} s_2 \cdot \overline{(s_1^)} -s_1 s_2 s_2 s_1 0$两行正交且各自范数都是 $|s_1|^2|s_2|^2$。于是 $S^H S (|s_1|^2|s_2|^2) I_2$。这个等式是整个方案的支点正因为 Gram 矩阵是标量乘单位阵接收端的两路信号在被合并之后符号 $s_1$ 和 $s_2$ 之间不会互相串扰二维联合最大似然搜索才被拆成两个独立的一维搜索。负号和共轭不是随便加的它必须与“第二个时隙把两个符号的动位置换一下”配合。写成 $(s_1, s_2)$ 和 $(s_2^, -s_1^)$ 也成立只是天线列顺序反了本质等价。但写成 $(s_1, s_2)$ 和 $(s_2, s_1)$ 就完蛋——两行内积变成 $2s_1 s_2$正交性没了译码端会看到一对互相干扰的量BER 曲线在高 SNR 处直接出现平台。正交结构的存在性有硬上限这点决定了 STBC 的适用边界。Tarokh 等人证明过复星座下码率为 1 的正交空时分组码只对两根发射天线存在三根和四根天线最高只能做到 3/4 码率需要往码矩阵里填零五个时隙发三个或四个符号。天线数再往上正交性只能靠牺牲码率换所以大规模天线系统里基本见不到 OSTBC走的是空分复用那条路。做链路选型时记住一句STBC 的分集增益随天线数涨但码率和检测复杂度也一起涨2×2 是甜点。2.3 写一个可复用的编码函数并做形状自检下面这段把码块结构固化成函数返回三维数组(时隙, 块索引, 天线索引)。这种排法在后面加接收天线维度时最省事广播不会出错。import numpy as np def alamouti_encode(sym, nt2, normalizeTrue): sym: (N,) 复数符号序列N 必须为偶数 返回 tx: (2, N//2, nt) —— [时隙, 块索引, 天线索引] sym np.asarray(sym, dtypenp.complex128) if sym.size % 2 ! 0: raise ValueError(符号数必须为偶数否则凑不满一个 Alamouti 块) if nt ! 2: raise ValueError(Alamouti 只对 2 根发射天线成立) s1 sym[0::2] # 每个块的第一个符号 s2 sym[1::2] # 每个块的第二个符号 tx np.empty((2, s1.size, 2), dtypenp.complex128) tx[0, :, 0] s1 # t0, 天线1 tx[0, :, 1] s2 # t0, 天线2 tx[1, :, 0] -np.conj(s2) # t1, 天线1 tx[1, :, 1] np.conj(s1) # t1, 天线2 return tx / np.sqrt(2) if normalize else tx逻辑上有三处值得单独说。sym[0::2]和sym[1::2]把序列拆成奇偶两路这对应“一个码块吃两个符号”的节奏后面译码端按同样的规则拼回去就能对齐tx的第一维是时隙而不是天线因为接收端求和是在时隙上做接收采样、在天线上做信道相乘这个顺序和物理过程一致最后除以√2是总功率归一把两根天线的辐射总能量拉回 1如果整条链路统一用“每根天线满功率”的口径把normalizeFalse打开即可但全链路必须一致。写完先别急着跑 BER做一个形状和正交性自检几行代码就能挡住一大半低级错误def check_orthogonality(sym, tol1e-10): tx alamouti_encode(sym, normalizeFalse) # (2, nblk, 2) S tx[:, 0, :] # 取第一个块的 2x2 gram S.conj().T S target (np.abs(sym[0])**2 np.abs(sym[1])**2) * np.eye(2) assert np.allclose(gram, target, atoltol), f正交性破坏:\n{gram} check_orthogonality(np.array([12j, -0.50.3j]))S.conj().T S就是 Gram 矩阵如果它不是对角阵说明共轭或负号放错了位置如果它是对角阵但两个对角元素不相等说明两个时隙的符号位置没对称。这一步跑通发射端的活儿就基本干净了。3. 接收端合并与 ML 译码把两路信号压成一个标量发射端把两个符号摊在了两根天线、两个时隙上接收端拿到的是一堆混在一起的观测。合并做的事就是用信道系数当权重把属于 $s_1$ 的那部分能量从两个时隙里全部抽出来同时把 $s_2$ 的干扰压掉。这一章的每一步都能手推推完再写代码就不会对着公式猜。3.1 平坦瑞利信道与准静态假设先把信道模型钉死。用平坦衰落、准静态的复数高斯信道$h_{n} \sim \mathcal{CN}(0,1)$且在一个码块的两个时隙内保持不变。准静态这个假设不是形式上的偷懒它直接决定了合并公式能不能成立——推导里默认 $h_1$、$h_2$ 在两个时隙里是同一组值如果它们每个时隙都重新生成合并出来的量就会带上一个无法消除的交叉项。工程上这个假设合不合理看相干时间。两个符号周期的跨度通常远小于信道相干时间所以准静态是站得住的但在高速移动场景下多普勒扩展会让相干时间掉到几个符号量级这时候 STBC 会出现误码平台靠加天线也压不下去。这类平台不是实现 bug是模型失效排查时要能区分开。另外这里的信道是“发射天线 n 到接收天线 r”的复增益多个接收天线之间假设独立同分布。真做硬件验证时天线间相关性强分集增益会打折仿真里可以加一个相关系数矩阵乘进去看退化幅度。3.2 共轭合并公式的推导与代码落地记接收天线上两个时隙的观测为 $r_0, r_1$噪声为 $n_0, n_1$$r_0 h_1 s_1 h_2 s_2 n_0$ $r_1 -h_1 s_2^* h_2 s_1^* n_1$要凑出只含 $s_1$ 的量思路是把 $r_0$ 乘 $\overline{h_1}$、把 $r_1$ 取共轭后乘 $h_2$$\overline{h_1} r_0 h_2 \overline{r_1} (|h_1|^2|h_2|^2) s_1 \overline{h_1} n_0 h_2 \overline{n_1}$交叉项正好抵消。同理要凑 $s_2$ 就用 $\overline{h_2} r_0 - h_1 \overline{r_1}$得到 $(|h_1|^2|h_2|^2) s_2$ 加噪声。合并增益记为 $g |h_1|^2|h_2|^2$它是一个实数正标量这一条极其关键——正因为它没有相位判决时不需要做相位补偿直接比距离就行。def alamouti_combine(rx, h): rx: (2, nblk, nr) 两个时隙的接收样本 h : (nblk, 2, nr) 信道nr 为接收天线数 返回 y1, y2: (nblk,) 合并后的两个软输出 g : (nblk,) 合并增益实数 r0, r1 rx[0], rx[1] # (nblk, nr) h0, h1 h[:, 0, :], h[:, 1, :] # (nblk, nr) y1 np.sum(np.conj(h0) * r0 h1 * np.conj(r1), axis1) y2 np.sum(np.conj(h1) * r0 - h0 * np.conj(r1), axis1) g np.sum(np.abs(h0) ** 2 np.abs(h1) ** p a hrefhttps://download.csdn.net/download/yuanclkk/10411105 stylecolor:#ec7500;font-size:14px; 本文还有配套的精品资源点击获取 /a img altmenu-r.4af5f7ec.gif srchttps://csdnimg.cn/release/wenkucmsfe/public/img/menu-r.4af5f7ec.gif stylewidth:16px;margin-left:4px;vertical-align:text-bottom;cursor:text; /p
