CSI-RS 信道估计仿真与评估平台

面向 NR CSI-RS 信道估计的工程化对标实现

构建可复现、可对照、可报告的端到端估计链路

【CSI-RS 对标实现】【NR 信道估计】【全链路仿真】【自动化报告】

📌 为什么选择

该工程聚焦 CSI-RS 信道估计的完整闭环实现。它把配置、建模、估计与可视化统一到单一工程中。这样可以在同一输入条件下稳定复现实验结论。

🎯 核心价值

🔬 学术研究价值

工程实现与标准条款和公式推导保持同向映射。

• 导频构造机理验证
• CDM 去扩行为分析
• 插值去噪链路评估
• 完美参考对照实验

💼 工程应用价值

工程结构强调可运行、可维护和可回归。

• 统一参数配置入口
• 模块职责清晰划分
• 脚本级一键运行
• 报告级结果沉淀

⚡ 技术亮点

🌊 CSI-RS 估计链路：传统分散脚本 vs 本工程闭环实现

| 特性      | 传统方案       | 本方案                                           |

| ----------- | ---------------- | -------------------------------------------------- |

| 参数管理  | 多脚本分散配置 | nr_csirs_example_config 集中配置            |

| 导频处理  | 局部临时实现   | nr_csi_rs 与 nr_csi_rs_indices 独立实现 |

| OFDM 链路 | 仅验证局部步骤 | nr_ofdm_info/modulate/demodulate 闭环       |

| 信道估计  | 单一 LS 视角   | LS + CDM + 插值 + CIR 去噪                       |

| 结果沉淀  | 图与数据分离   | 统一结果结构体 + 自动报告输出                    |

📊 性能指标（实测数据）

数据来源： reports/nr_csirs_advanced_demo_2026-03-23_103625/高级演示图报告.md

| 场景               | 基线               | 本方案              | 结论             |

| -------------------- | -------------------- | --------------------- | ------------------ |

| 全局 NMSE          | 工程目标 ≤ -30 dB | -35.09 dB           | 达到并优于目标   |

| 误差 P90           | 工程目标 ≤ 0.02   | 0.01081             | 误差尾部可控     |

| 有限样本率         | 工程目标 100%      | 100.00%             | 估计样本全部有效 |

| 逐符号 NMSE 稳定性 | 目标波动小于 1 dB  | [-35.13, -35.02] dB | 跨符号稳定       |

🎯 同步与估计一体化能力

流程在同步、估计与可视化之间保持单次贯通。输出字段可直接用于误差诊断和对照分析。

| 参数     | 配置                             | 性能                       |

| ---------- | ---------------------------------- | ---------------------------- |

| 定时估计 | 互相关峰值法 + 参考 CSI-RS       | 稳定输出 TimingOffset |

| 信道估计 | CDMLengths=[2,1]+ CIR 去噪 | 全局 NMSE -35.09 dB        |

| 报告生成 | 6 图自动导出 + 指标汇总          | 一次运行形成完整报告       |

🖥️ 运行环境

工程默认面向 MATLAB 仿真工作流。运行依赖保持最小化，并与脚本入口保持一致。

• 语言：MATLAB（建议 R2025b）
• 依赖：Communications Toolbox、Signal Processing Toolbox
• 硬件：Windows x64 （建议内存 16 GB 及以上）

📁 项目结构

csirs_ce/
 ├── src/                          # 核心算法源码
 │   ├── core/                     # 配置、栅格与通用工具
 │   ├── csirs/                    # CSI-RS 符号与索引
 │   └── channel/                  # 信道、同步与估计
 ├── examples/                     # 示例与报告脚本
 │   ├── run_nr_channel_estimation_using_csirs.m         # 端到端示例入口
 │   └── run_nr_channel_estimation_advanced_report.m     # 高级图与报告入口

📄 文档体系

文档体系由算法与代码两条主线构成。两类文档共同支撑复现与二次分析。

📘 算法文档

docs/算法文档.md：覆盖 CSI-RS 标准语义、估计原理与公式推导。

📒 代码文档

docs/代码文档.md：覆盖模块职责、接口定义、主流程与数据结构。

💻 核心代码展示

🔥 CSI-RS 参考信号生成模块

该模块负责导频符号生成与资源配置解析。实现覆盖多行配置、密度模式与周期调度。

function [sym, info] = nr_csi_rs(carrier, csirs_cfg)
 % NR_CSI_RS 5G NR 信道状态信息参考信号生成
 carrier = nr_utils_to_carrier_struct(carrier);
 
 nslot = carrier.NSlot;
 nframe = nr_utils_get_field(carrier, 'NFrame', 0);
 scs = nr_utils_get_field(carrier, 'SubcarrierSpacing', 15);
 
 if strcmpi(carrier.CyclicPrefix, 'extended')
     symPerSlot = 12;
 else
     symPerSlot = 14;
 end
 
 rowNumber = nr_utils_get_field(csirs_cfg, 'RowNumber', 3);
 if numel(rowNumber) > 1
     [sym, info] = generate_multi_resource(carrier, csirs_cfg, symPerSlot, 10 * scs / 15);
     return;
 end

🌟 OFDM 调制主链模块

该模块负责栅格维度校验、相位预补偿与调制主过程。其输出直接进入信道建模链路。

%  获取 OFDM 参数
 info = nr_ofdm_info(carrier, varargin{:});
 
 %  验证栅格维度
 K = carrier.NSizeGrid * 12;
 assert(size(grid, 1) == K, ...
     'nr_ofdm_modulate:invalidGrid', ...
     '栅格行数 (%d) 必须等于 NSizeGrid*12 (%d)', size(grid, 1), K);
 
 N = size(grid, 2);
 P = size(grid, 3);
 nfft = info.Nfft;
 
 %  获取当前时隙的 CP 长度和符号相位
 nSlot = carrier.NSlot;
 [cp_lengths, sym_phases] = get_slot_params(info, nSlot, N);
 
 %  相位预补偿
 if any(sym_phases ~= 0)
     phase_matrix = reshape(exp(1j * sym_phases(1:N)), 1, N);
     grid = grid .* phase_matrix;
 end

🚀 实用信道估计核心模块

该模块执行 LS 估计、CDM 去扩和扩展网格写回。后续再衔接插值、去噪与噪声方差估计。

%  LS 估计与 CDM 解扩（扩展子载波网格）
 eRB = 4;
 eK = eRB * 12;
 Kext = K + eK;
 
 Hext = complex(zeros(Kext, N, R, P, 'like', rxGrid));
 Hn = complex(nan(Kext, N, R, P, 'like', rxGrid));
 
 for p = 1:P
     [ksubs, nsubs, thisPort] = get_port_subscripts(K, N, P, refInd, p);
     refSymThisPort = refSym(thisPort);
     un = unique(nsubs).';
 
     for uni = 1:numel(un)
         n = un(uni);
         mask = (nsubs == n);
         k = ksubs(mask);
         if isempty(k)
             continue;
         end
 
         refS = refSymThisPort(mask);
         for r = 1:R
             hLS = rxGrid(k, n, r) ./ refS;
             if fdCDM > 1
                 hLS = despread_fd_cdm(hLS, fdCDM);
             end
             Hext(k + eK/2, n, r, p) = hLS;
         end
     end
 end

🎬 一键运行

run('startup_csirs_ce.m');
 run('examples/run_nr_channel_estimation_using_csirs.m');
 run('examples/run_nr_channel_estimation_advanced_report.m');

结果预览

在 SNR=50 dB、 Seed=2026、 2Tx/2Rx 配置下，最新报告给出全局 NMSE 为 -35.09 dB。  图示说明：单次运行会输出资源映射、信道图谱、误差统计、时频趋势与 MIMO 指标图。

📸 演示图片预览

🛒 获取方式

本文代码仅为核心片段，完整版工程已整理好。 关注公众号 【 3GPP 仿真实验室】进行获取。

📚 参考文献

1. 3GPP TS 38.211, NR; Physical channels and modulation.
2. 3GPP TR 38.901, Study on channel model for frequencies from 0.5 to 100 GHz.
3. 3GPP TS 38.214, NR; Physical layer procedures for data.