尧图网站建设 
简介本资源是一份面向通信工程专业本科生及无线通信入门学习者的Alamouti空时编码仿真实践材料聚焦MIMO系统中经典的2发1收、2发2收等典型场景帮助读者理解分集增益原理与编码矩阵设计逻辑。压缩包共8个文件含5个MATLAB源码文件.m与3张关键仿真结果图.jpg其中Alamoudi_RUN.m为主运行脚本rey.m实现瑞利信道建模Alamoudi_2T1R.m与Alamoudi_2T2R.m分别对应不同天线配置的编解码流程配套图片直观展示误码率曲线与信号星座图便于验证理论性能。资源体积仅69KB轻量易部署代码结构清晰、注释完整涵盖信道建模、符号映射、Alamouti编码矩阵构造、最大比合并MRC解码等核心环节可直接运行复现经典文献结论。目前已有180人学习下载适合课堂实验、课程设计或自学巩固空时编码基础概念。1. Alamouti 编码不是“加个天线就抗衰落”——它用确定性正交结构在2×1 MIMO中实现满分分集增益很多人第一次看到 Alamouti 编码以为只是“两根发射天线发不同信号”甚至误以为它靠功率叠加提升信噪比。实际上Alamouti 的核心价值在于在不增加接收天线、不牺牲带宽、不依赖信道估计的前提下仅靠发射端的符号配对与时间-空间正交映射就让单接收天线系统获得满分2阶分集增益。这意味着在瑞利衰落信道下误码率随 SNR 的下降速度比 SISO 快一倍斜率翻倍且解码复杂度仅为线性——不需要矩阵求逆连 MMSE 都不用。本资源包正是用 MATLAB 实现了这一经典方案的全链路验证从 QPSK 符号生成、Alamouti 矩阵构造、AWGN 信道建模到接收端的共轭转置合并与硬判决全部可运行、可调试、可修改调制阶数或天线配置。适合通信原理课设、MIMO 入门实验、或作为 OFDM 系统中空时编码模块的 baseline 参考。2. 为什么 Alamouti 能用 2T1R 实现满分分集从正交设计到接收合并的数学闭环2.1 Alamouti 矩阵的构造逻辑时间-空间正交性是分集增益的唯一来源Alamouti 编码的本质是将两个复数符号 $s_1$ 和 $s_2$ 映射为一个 $2 \times 2$ 的发射矩阵$$ \mathbf{X} \begin{bmatrix} s_1 s_2 \ -s_2^* s_1^* \end{bmatrix} $$该矩阵满足关键性质$\mathbf{X}^H \mathbf{X} (|s_1|^2 |s_2|^2) \mathbf{I}_2$即列向量严格正交且能量守恒。这意味着第1时刻天线1发 $s_1$天线2发 $s_2$第2时刻天线1发 $-s_2^$天线2发 $s_1^$。这种设计使两个符号在接收端被“天然分离”——即使信道系数 $h_1, h_2$ 完全未知且随时间变化只要在两个时刻内保持不变块衰落假设就能通过共轭转置操作消除符号间干扰。提示Alamoudi_2T1R.m中第 32–35 行正是该矩阵的 MATLAB 实现注意conj()的使用位置和负号方向错一位就会破坏正交性导致分集增益归零。2.2 接收端合并无需信道估计的线性处理直接还原符号能量设接收信号为 $\mathbf{y} [y_1, y_2]^T$其中$y_1 h_1 s_1 h_2 s_2 n_1$ 时刻1$y_2 -h_1 s_2^* h_2 s_1^* n_2$ 时刻2接收端构造合并向量 $$ \tilde{\mathbf{y}} \begin{bmatrix} h_1^* h_2 \ -h_2^* h_1 \end{bmatrix} \begin{bmatrix} y_1 \ y_2 \end{bmatrix}(|h_1|^2 |h_2|^2) \begin{bmatrix} s_1 \ s_2 \end{bmatrix} \begin{bmatrix} h_1^* n_1 h_2 n_2 \ -h_2^* n_1 h_1 n_2 \end{bmatrix} $$可见$s_1$ 和 $s_2$ 被独立放大 $(|h_1|^2 |h_2|^2)$ 倍而噪声项仍为高斯白噪声等效 SNR 提升为原信道增益之和。这正是满分分集增益的数学体现错误概率 $P_e \propto \left( \frac{1}{\text{SNR}} \right)^2$。2.2.1 MATLAB 中的合并实现rey.m的关键三步打开rey.m定位到第 48–52 行其核心逻辑如下% 假设 h1, h2 为已知信道仿真中可设为随机复数 y1 h1 * s1 h2 * s2 noise1; y2 -conj(h2) * s1 conj(h1) * s2 noise2; % 注意此处是标准接收模型与发射矩阵共轭转置匹配 % 合并构造合并矩阵并应用 y_tilde1 conj(h1) * y1 h2 * y2; y_tilde2 -conj(h2) * y1 h1 * y2; % 硬判决QPSK s1_hat sign(real(y_tilde1)) 1i * sign(imag(y_tilde1)); s2_hat sign(real(y_tilde2)) 1i * sign(imag(y_tilde2));这段代码的关键在于y2的构造必须与 Alamouti 发射矩阵严格对应-conj(h2)和conj(h1)否则合并后无法消除干扰项合并系数不是任意选的而是发射矩阵 $\mathbf{X}$ 的共轭转置 $\mathbf{X}^H$sign()函数仅适用于 BPSK/QPSK 等星座点位于单位圆上的调制若改用 16-QAM此处需替换为最近邻判决min(abs(y_tilde - constellation), [], 2)。2.3 仿真参数表控制变量法验证分集阶数为确认是否真达到 2 阶分集需固定其他参数仅改变 SNR 并统计误码率。下表为Alamoudi_RUN.m中推荐的基准配置可直接复制修改参数名默认值说明修改建议M4QPSK 调制阶数sqrt(M)为每维电平数改为2测试 BPSK观察斜率是否一致N10000总发送符号数≥5000 才能稳定统计 $10^{-3}$ 量级误码率SNR_dB0:2:20SNR 扫描范围dB必须覆盖 10–20 dB 区间以观察斜率拐点h1,h2randn1i*randn瑞利信道系数复高斯每次循环重新生成模拟快衰落noise_var10^(-SNR_dB/10)AWGN 方差由 SNR 自动计算勿手动赋值注意Alamouti_T2R1.jpg和Alamouti_T2R2.jpg是该配置下生成的 BER 曲线图横轴为 SNR纵轴为 log10(BER)。你会发现 2T1R 曲线斜率明显陡于 SISO对比图中虚线而 2T2R 斜率更陡——这正是分集阶数 $L N_t \times N_r$ 的直观体现。3. 四种典型场景的代码复现从单接收起步逐步扩展至双收双发3.1 2T1R 基础版Alamoudi_2T1R.m的逐行执行与断点调试这是理解 Alamouti 最小可行单元。运行前请确保工作路径包含该文件并在命令行输入 Alamoudi_2T1R程序将自动执行以下流程生成 $10^4$ 个 QPSK 符号s randsrc(1, N, [11i,1-1i,-11i,-1-1i])按 Alamouti 规则分组每2个符号构成1个 $2\times2$ 块加入瑞利信道与 AWGN接收端合并并判决绘制 BER-SNR 曲线保存为Alamouti_T2R1.jpg。3.1.1 关键断点设置建议提升调试效率在第 28 行X [s1, s2; -conj(s2), conj(s1)];处设断点检查X是否为 2×2 复矩阵且满足 $\det(\mathbf{X}) |s_1|^2 |s_2|^2$在第 45 行y_tilde1 conj(h1)*y1 h2*y2;处查看y_tilde1实部/虚部分布应近似为原始s1的缩放噪声若 BER 高于理论值如 10 dB 时 BER 1e-3立即检查noise_var是否按10^(-SNR/10)计算而非10^(SNR/10)常见笔误。3.2 2T2R 扩展版Alamoudi_2T2R.m中的信道矩阵与合并权重升级当接收端有2根天线时系统变为 2×2 MIMO信道变为 $2\times2$ 矩阵 $\mathbf{H} [h_{11}, h_{12}; h_{21}, h_{22}]$。此时 Alamouti 合并需对每个接收天线单独进行再能量叠加% 对接收天线1做合并 y1_tilde1 conj(h11)*y1_1 h12*y2_1; % y1_1, y2_1 为天线1在两时刻的接收 y1_tilde2 -conj(h12)*y1_1 h11*y2_1; % 对接收天线2做合并 y2_tilde1 conj(h21)*y1_2 h22*y2_2; y2_tilde2 -conj(h22)*y1_2 h21*y2_2; % 能量合并最大比合并 MRC s1_hat y1_tilde1 y2_tilde1; s2_hat y1_tilde2 y2_tilde2;Alamoudi_2T2R.m第 67–75 行即为此逻辑。注意y1_1表示天线1在时刻1的接收y1_2表示天线1在时刻2的接收命名易混淆需对照注释合并后不再直接sign()因能量叠加使幅度增大需先归一化s1_hat s1_hat / (abs(h11)^2 abs(h12)^2 abs(h21)^2 abs(h22)^2)该配置下理论分集阶数为 $2 \times 2 4$BER 曲线斜率应为 2T1R 的两倍见Alamouti_T2R1_T2R2.jpg中实线 vs 虚线。3.3 多天线配置统一框架Refresh.m的模块化重构思路Refresh.m并非独立仿真脚本而是提供了一套可配置的 Alamouti 封装函数。其核心价值在于将发射天线数Nt、接收天线数Nr、调制方式mod_type抽象为输入参数避免重复写矩阵构造逻辑。例如% 构造 Alamouti 码字支持 Nt2 固定 function X alamouti_encode(s1, s2, Nt) if Nt 2 X [s1, s2; -conj(s2), conj(s1)]; else error(Alamouti only supports 2 transmit antennas); end end % 通用接收合并支持任意 Nr function [s1_hat, s2_hat] alamouti_decode(y, H, Nr) s1_hat zeros(1, size(y,2)); s2_hat zeros(1, size(y,2)); for i 1:Nr h1 H(i,1); h2 H(i,2); y1 y(i,1); y2 y(i,2); s1_hat(i) conj(h1)*y1 h2*y2; s2_hat(i) -conj(h2)*y1 h1*y2; end s1_hat sum(s1_hat); s2_hat sum(s2_hat); % MRC 合并 end提示若需扩展至 4T1R不能直接套用 Alamouti它仅定义于 2T应切换至空时格码STTC或空频编码SFCRefresh.m中未实现但其函数接口设计已预留扩展槽位如encode_method参数。3.4 信道建模进阶从理想瑞利到实际相关信道的替换路径当前所有.m文件默认使用独立同分布瑞利信道h (randn 1i*randn)/sqrt(2)。但在实际部署中天线间距过小会导致信道相关。要模拟相关信道只需替换h1,h2的生成方式% 引入相关系数 rho 0.7强相关 rho 0.7; h_corr sqrt(rho) * h1 sqrt(1-rho) * (randn 1i*randn)/sqrt(2); % 此时 h1 与 h_corr 相关系数 ≈ rho在Alamoudi_RUN.m中将第 38 行h1 ...和第 39 行h2 ...替换为上述逻辑重新运行即可观察相关性对分集增益的削弱效应——BER 曲线斜率将趋近于 1 阶类似 SISO。4. 误码率曲线验证与参数敏感性分析识别真实分集阶数的三步法4.1 对数坐标下的斜率提取用最小二乘拟合确认分集阶数理论分集阶数 $L$ 体现在 BER 曲线的渐近斜率上$\log_{10} \text{BER} \approx -L \cdot \log_{10} \text{SNR} C$。要实证是否达到满分 $L2$需在高 SNR 区间≥15 dB提取斜率% 假设 snr_db 和 ber 为已知向量来自 Alamoudi_RUN.m 输出 idx snr_db 15; % 取高 SNR 区间 X [ones(sum(idx),1), log10(10.^(snr_db(idx)/10))]; % 设计矩阵 [1, log10(SNR)] coeff X \ log10(ber(idx)); % 最小二乘拟合 L_est -coeff(2); % 斜率即为估计分集阶数 fprintf(Estimated diversity order: %.2f\n, L_est);运行此段代码于Alamouti_T2R1.jpg对应数据输出应接近2.00若为1.2则说明存在实现错误如噪声方差计算错误或合并系数符号颠倒。4.2 关键参数敏感性表格哪些改动会直接杀死分集增益参数错误设置示例对 BER 的影响诊断方法发射矩阵负号X [s1, s2; conj(s2), conj(s1)]漏负号分集增益消失BER 与 SISO 重合检查det(X)是否恒等于abs(s1)^2 abs(s2)^2接收合并共轭位置y_tilde1 h1*y1 h2*y2未共轭 h1合并后残留干扰项BER 下降缓慢观察y_tilde1的实部分布是否围绕real(s1)集中噪声方差noise_var 10^(SNR_dB/10)指数符号反所有 SNR 点 BER 偏高 10–100 倍对比var(noise)与noise_var数值是否相等符号分组长度s1 s(1:2:end); s2 s(2:2:end)未按块分组符号配对错位完全失效在X构造前disp([s1(1:3); s2(1:3)])查看是否成对4.3 一个实用技巧用Alamouti_T2R1_T2R2.jpg快速定位性能瓶颈该图同时绘制了 2T1R蓝线与 2T2R红线的 BER 曲线。若你修改代码后发现蓝线变平缓但红线仍陡峭→ 问题出在 2T1R 的接收合并逻辑如rey.m中h1,h2赋值错误两条线同步变平缓→ 问题出在公共模块如符号生成、噪声添加、判决方式蓝线正常红线异常→ 问题出在 2T2R 的信道矩阵应用或 MRC 合并如Alamoudi_2T2R.m第 72 行未对y2_tilde1归一化。提示直接双击Alamouti_T2R1_T2R2.jpg查看图像MATLAB 会自动打开 figure 窗口。右键点击曲线 → “Data Cursor”可精确读取任意 SNR 点的 BER 值与理论值 $0.5 \cdot \text{erfc}(\sqrt{\text{SNR}})$2T1R或 $0.5 \cdot \text{erfc}(\sqrt{2 \cdot \text{SNR}})$2T2R比对误差 0.5 dB 即需排查。本文还有配套的精品资源点击获取 Alamouti编码原理与MATLAB实现：2T1R满分分集增益详解
简介本资源是一份面向通信工程专业本科生及无线通信入门学习者的Alamouti空时编码仿真实践材料聚焦MIMO系统中经典的2发1收、2发2收等典型场景帮助读者理解分集增益原理与编码矩阵设计逻辑。压缩包共8个文件含5个MATLAB源码文件.m与3张关键仿真结果图.jpg其中Alamoudi_RUN.m为主运行脚本rey.m实现瑞利信道建模Alamoudi_2T1R.m与Alamoudi_2T2R.m分别对应不同天线配置的编解码流程配套图片直观展示误码率曲线与信号星座图便于验证理论性能。资源体积仅69KB轻量易部署代码结构清晰、注释完整涵盖信道建模、符号映射、Alamouti编码矩阵构造、最大比合并MRC解码等核心环节可直接运行复现经典文献结论。目前已有180人学习下载适合课堂实验、课程设计或自学巩固空时编码基础概念。1. Alamouti 编码不是“加个天线就抗衰落”——它用确定性正交结构在2×1 MIMO中实现满分分集增益很多人第一次看到 Alamouti 编码以为只是“两根发射天线发不同信号”甚至误以为它靠功率叠加提升信噪比。实际上Alamouti 的核心价值在于在不增加接收天线、不牺牲带宽、不依赖信道估计的前提下仅靠发射端的符号配对与时间-空间正交映射就让单接收天线系统获得满分2阶分集增益。这意味着在瑞利衰落信道下误码率随 SNR 的下降速度比 SISO 快一倍斜率翻倍且解码复杂度仅为线性——不需要矩阵求逆连 MMSE 都不用。本资源包正是用 MATLAB 实现了这一经典方案的全链路验证从 QPSK 符号生成、Alamouti 矩阵构造、AWGN 信道建模到接收端的共轭转置合并与硬判决全部可运行、可调试、可修改调制阶数或天线配置。适合通信原理课设、MIMO 入门实验、或作为 OFDM 系统中空时编码模块的 baseline 参考。2. 为什么 Alamouti 能用 2T1R 实现满分分集从正交设计到接收合并的数学闭环2.1 Alamouti 矩阵的构造逻辑时间-空间正交性是分集增益的唯一来源Alamouti 编码的本质是将两个复数符号 $s_1$ 和 $s_2$ 映射为一个 $2 \times 2$ 的发射矩阵$$ \mathbf{X} \begin{bmatrix} s_1 s_2 \ -s_2^* s_1^* \end{bmatrix} $$该矩阵满足关键性质$\mathbf{X}^H \mathbf{X} (|s_1|^2 |s_2|^2) \mathbf{I}_2$即列向量严格正交且能量守恒。这意味着第1时刻天线1发 $s_1$天线2发 $s_2$第2时刻天线1发 $-s_2^$天线2发 $s_1^$。这种设计使两个符号在接收端被“天然分离”——即使信道系数 $h_1, h_2$ 完全未知且随时间变化只要在两个时刻内保持不变块衰落假设就能通过共轭转置操作消除符号间干扰。提示Alamoudi_2T1R.m中第 32–35 行正是该矩阵的 MATLAB 实现注意conj()的使用位置和负号方向错一位就会破坏正交性导致分集增益归零。2.2 接收端合并无需信道估计的线性处理直接还原符号能量设接收信号为 $\mathbf{y} [y_1, y_2]^T$其中$y_1 h_1 s_1 h_2 s_2 n_1$ 时刻1$y_2 -h_1 s_2^* h_2 s_1^* n_2$ 时刻2接收端构造合并向量 $$ \tilde{\mathbf{y}} \begin{bmatrix} h_1^* h_2 \ -h_2^* h_1 \end{bmatrix} \begin{bmatrix} y_1 \ y_2 \end{bmatrix}(|h_1|^2 |h_2|^2) \begin{bmatrix} s_1 \ s_2 \end{bmatrix} \begin{bmatrix} h_1^* n_1 h_2 n_2 \ -h_2^* n_1 h_1 n_2 \end{bmatrix} $$可见$s_1$ 和 $s_2$ 被独立放大 $(|h_1|^2 |h_2|^2)$ 倍而噪声项仍为高斯白噪声等效 SNR 提升为原信道增益之和。这正是满分分集增益的数学体现错误概率 $P_e \propto \left( \frac{1}{\text{SNR}} \right)^2$。2.2.1 MATLAB 中的合并实现rey.m的关键三步打开rey.m定位到第 48–52 行其核心逻辑如下% 假设 h1, h2 为已知信道仿真中可设为随机复数 y1 h1 * s1 h2 * s2 noise1; y2 -conj(h2) * s1 conj(h1) * s2 noise2; % 注意此处是标准接收模型与发射矩阵共轭转置匹配 % 合并构造合并矩阵并应用 y_tilde1 conj(h1) * y1 h2 * y2; y_tilde2 -conj(h2) * y1 h1 * y2; % 硬判决QPSK s1_hat sign(real(y_tilde1)) 1i * sign(imag(y_tilde1)); s2_hat sign(real(y_tilde2)) 1i * sign(imag(y_tilde2));这段代码的关键在于y2的构造必须与 Alamouti 发射矩阵严格对应-conj(h2)和conj(h1)否则合并后无法消除干扰项合并系数不是任意选的而是发射矩阵 $\mathbf{X}$ 的共轭转置 $\mathbf{X}^H$sign()函数仅适用于 BPSK/QPSK 等星座点位于单位圆上的调制若改用 16-QAM此处需替换为最近邻判决min(abs(y_tilde - constellation), [], 2)。2.3 仿真参数表控制变量法验证分集阶数为确认是否真达到 2 阶分集需固定其他参数仅改变 SNR 并统计误码率。下表为Alamoudi_RUN.m中推荐的基准配置可直接复制修改参数名默认值说明修改建议M4QPSK 调制阶数sqrt(M)为每维电平数改为2测试 BPSK观察斜率是否一致N10000总发送符号数≥5000 才能稳定统计 $10^{-3}$ 量级误码率SNR_dB0:2:20SNR 扫描范围dB必须覆盖 10–20 dB 区间以观察斜率拐点h1,h2randn1i*randn瑞利信道系数复高斯每次循环重新生成模拟快衰落noise_var10^(-SNR_dB/10)AWGN 方差由 SNR 自动计算勿手动赋值注意Alamouti_T2R1.jpg和Alamouti_T2R2.jpg是该配置下生成的 BER 曲线图横轴为 SNR纵轴为 log10(BER)。你会发现 2T1R 曲线斜率明显陡于 SISO对比图中虚线而 2T2R 斜率更陡——这正是分集阶数 $L N_t \times N_r$ 的直观体现。3. 四种典型场景的代码复现从单接收起步逐步扩展至双收双发3.1 2T1R 基础版Alamoudi_2T1R.m的逐行执行与断点调试这是理解 Alamouti 最小可行单元。运行前请确保工作路径包含该文件并在命令行输入 Alamoudi_2T1R程序将自动执行以下流程生成 $10^4$ 个 QPSK 符号s randsrc(1, N, [11i,1-1i,-11i,-1-1i])按 Alamouti 规则分组每2个符号构成1个 $2\times2$ 块加入瑞利信道与 AWGN接收端合并并判决绘制 BER-SNR 曲线保存为Alamouti_T2R1.jpg。3.1.1 关键断点设置建议提升调试效率在第 28 行X [s1, s2; -conj(s2), conj(s1)];处设断点检查X是否为 2×2 复矩阵且满足 $\det(\mathbf{X}) |s_1|^2 |s_2|^2$在第 45 行y_tilde1 conj(h1)*y1 h2*y2;处查看y_tilde1实部/虚部分布应近似为原始s1的缩放噪声若 BER 高于理论值如 10 dB 时 BER 1e-3立即检查noise_var是否按10^(-SNR/10)计算而非10^(SNR/10)常见笔误。3.2 2T2R 扩展版Alamoudi_2T2R.m中的信道矩阵与合并权重升级当接收端有2根天线时系统变为 2×2 MIMO信道变为 $2\times2$ 矩阵 $\mathbf{H} [h_{11}, h_{12}; h_{21}, h_{22}]$。此时 Alamouti 合并需对每个接收天线单独进行再能量叠加% 对接收天线1做合并 y1_tilde1 conj(h11)*y1_1 h12*y2_1; % y1_1, y2_1 为天线1在两时刻的接收 y1_tilde2 -conj(h12)*y1_1 h11*y2_1; % 对接收天线2做合并 y2_tilde1 conj(h21)*y1_2 h22*y2_2; y2_tilde2 -conj(h22)*y1_2 h21*y2_2; % 能量合并最大比合并 MRC s1_hat y1_tilde1 y2_tilde1; s2_hat y1_tilde2 y2_tilde2;Alamoudi_2T2R.m第 67–75 行即为此逻辑。注意y1_1表示天线1在时刻1的接收y1_2表示天线1在时刻2的接收命名易混淆需对照注释合并后不再直接sign()因能量叠加使幅度增大需先归一化s1_hat s1_hat / (abs(h11)^2 abs(h12)^2 abs(h21)^2 abs(h22)^2)该配置下理论分集阶数为 $2 \times 2 4$BER 曲线斜率应为 2T1R 的两倍见Alamouti_T2R1_T2R2.jpg中实线 vs 虚线。3.3 多天线配置统一框架Refresh.m的模块化重构思路Refresh.m并非独立仿真脚本而是提供了一套可配置的 Alamouti 封装函数。其核心价值在于将发射天线数Nt、接收天线数Nr、调制方式mod_type抽象为输入参数避免重复写矩阵构造逻辑。例如% 构造 Alamouti 码字支持 Nt2 固定 function X alamouti_encode(s1, s2, Nt) if Nt 2 X [s1, s2; -conj(s2), conj(s1)]; else error(Alamouti only supports 2 transmit antennas); end end % 通用接收合并支持任意 Nr function [s1_hat, s2_hat] alamouti_decode(y, H, Nr) s1_hat zeros(1, size(y,2)); s2_hat zeros(1, size(y,2)); for i 1:Nr h1 H(i,1); h2 H(i,2); y1 y(i,1); y2 y(i,2); s1_hat(i) conj(h1)*y1 h2*y2; s2_hat(i) -conj(h2)*y1 h1*y2; end s1_hat sum(s1_hat); s2_hat sum(s2_hat); % MRC 合并 end提示若需扩展至 4T1R不能直接套用 Alamouti它仅定义于 2T应切换至空时格码STTC或空频编码SFCRefresh.m中未实现但其函数接口设计已预留扩展槽位如encode_method参数。3.4 信道建模进阶从理想瑞利到实际相关信道的替换路径当前所有.m文件默认使用独立同分布瑞利信道h (randn 1i*randn)/sqrt(2)。但在实际部署中天线间距过小会导致信道相关。要模拟相关信道只需替换h1,h2的生成方式% 引入相关系数 rho 0.7强相关 rho 0.7; h_corr sqrt(rho) * h1 sqrt(1-rho) * (randn 1i*randn)/sqrt(2); % 此时 h1 与 h_corr 相关系数 ≈ rho在Alamoudi_RUN.m中将第 38 行h1 ...和第 39 行h2 ...替换为上述逻辑重新运行即可观察相关性对分集增益的削弱效应——BER 曲线斜率将趋近于 1 阶类似 SISO。4. 误码率曲线验证与参数敏感性分析识别真实分集阶数的三步法4.1 对数坐标下的斜率提取用最小二乘拟合确认分集阶数理论分集阶数 $L$ 体现在 BER 曲线的渐近斜率上$\log_{10} \text{BER} \approx -L \cdot \log_{10} \text{SNR} C$。要实证是否达到满分 $L2$需在高 SNR 区间≥15 dB提取斜率% 假设 snr_db 和 ber 为已知向量来自 Alamoudi_RUN.m 输出 idx snr_db 15; % 取高 SNR 区间 X [ones(sum(idx),1), log10(10.^(snr_db(idx)/10))]; % 设计矩阵 [1, log10(SNR)] coeff X \ log10(ber(idx)); % 最小二乘拟合 L_est -coeff(2); % 斜率即为估计分集阶数 fprintf(Estimated diversity order: %.2f\n, L_est);运行此段代码于Alamouti_T2R1.jpg对应数据输出应接近2.00若为1.2则说明存在实现错误如噪声方差计算错误或合并系数符号颠倒。4.2 关键参数敏感性表格哪些改动会直接杀死分集增益参数错误设置示例对 BER 的影响诊断方法发射矩阵负号X [s1, s2; conj(s2), conj(s1)]漏负号分集增益消失BER 与 SISO 重合检查det(X)是否恒等于abs(s1)^2 abs(s2)^2接收合并共轭位置y_tilde1 h1*y1 h2*y2未共轭 h1合并后残留干扰项BER 下降缓慢观察y_tilde1的实部分布是否围绕real(s1)集中噪声方差noise_var 10^(SNR_dB/10)指数符号反所有 SNR 点 BER 偏高 10–100 倍对比var(noise)与noise_var数值是否相等符号分组长度s1 s(1:2:end); s2 s(2:2:end)未按块分组符号配对错位完全失效在X构造前disp([s1(1:3); s2(1:3)])查看是否成对4.3 一个实用技巧用Alamouti_T2R1_T2R2.jpg快速定位性能瓶颈该图同时绘制了 2T1R蓝线与 2T2R红线的 BER 曲线。若你修改代码后发现蓝线变平缓但红线仍陡峭→ 问题出在 2T1R 的接收合并逻辑如rey.m中h1,h2赋值错误两条线同步变平缓→ 问题出在公共模块如符号生成、噪声添加、判决方式蓝线正常红线异常→ 问题出在 2T2R 的信道矩阵应用或 MRC 合并如Alamoudi_2T2R.m第 72 行未对y2_tilde1归一化。提示直接双击Alamouti_T2R1_T2R2.jpg查看图像MATLAB 会自动打开 figure 窗口。右键点击曲线 → “Data Cursor”可精确读取任意 SNR 点的 BER 值与理论值 $0.5 \cdot \text{erfc}(\sqrt{\text{SNR}})$2T1R或 $0.5 \cdot \text{erfc}(\sqrt{2 \cdot \text{SNR}})$2T2R比对误差 0.5 dB 即需排查。本文还有配套的精品资源点击获取