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

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

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

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

📌 为什么选择

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

🎯 核心价值

🔬 学术研究价值

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

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

💼 工程应用价值

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

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

⚡ 技术亮点

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

📊 性能指标（实测数据）

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

🎯 同步与估计一体化能力

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

🖥️ 运行环境

工程默认面向 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 指标图。

📸 演示图片预览

🛒 获取方式

📚 参考文献

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.