简介面向通信专业本科、硕士生及教研人员的Matlab仿真资源包围绕多径衰落信道模型及多普勒频移效应展开能够帮助读者掌握无线信道仿真基础流程。压缩包共16个文件包含7个可运行的MATLAB脚本、8张运行结果JPG图像和1份Word说明文档整体大小仅382KB脚本覆盖参数配置、多径叠加、信道响应计算到结果绘制等环节便于对照图片检查每一步输出。文档对多径传输机理和仿真设计做了梳理适合零基础或初学通信仿真的读者按步骤复现。目前已有204人学习浏览使用Matlab 2019a即可直接运行也可在课程设计、毕业设计或课堂教学中作为基础参考。读者通过修改时延、增益、多普勒频移等参数能进一步观察不同信道条件下的衰落变化比如幅度波动与频率选择性衰落为后续通信系统级仿真打下基础。1. 多径衰落信道matlab仿真从多径传输直通多普勒频移的脚本化拆解通信多径衰落信道matlab仿真这套资源主题落在多径传输与多普勒频移两条主线上。从u1.m到u7.m七个脚本把信道脉冲响应、瑞利包络、多普勒功率谱、误码率曲线依次跑通并附六张运行结果图供对照。适合本科、硕士做通信原理课程设计把教材里的相干带宽、相干时间变成可调参数的仿真图也适合工作两三年的算法工程师在做系统级链路仿真前确认自己对抽头延迟线和最大多普勒频移的理解没有偏差。重点不是代码本身多复杂而是能从这些图中读出延迟扩展怎么毁掉频响、多普勒频移怎么把星座点推成圆环这才是后续改参数与扩展仿真的起点。2. 多普勒频移与延迟扩展多径衰落信道仿真的两类基础参数2.1 先判定频繁选择性还是时间选择性再决定改哪个脚本变量多径衰落仿真的第一步不是写代码而是确认信道落在哪个区间。延迟扩展决定频率选择性多普勒频移决定时间选择性两者来自同一份信道脉冲响应但在脚本里由完全独立的变量控制。拆这类教学资源时我会先把参数和脚本的对应关系列出来改任何系数前先对着表查一次能省掉大量无效调试时间。参数与信道建模的关系对应脚本位置相干带宽 Bc ≈ 1/最大延迟扩展Bc大于信号带宽是平坦衰落否则是频率选择性衰落u1.m、u2.m 的路径延迟差最大多普勒频移 fd v·fc/cfd决定相干时间TcTc小于符号周期时信道在符号内快速变化u3.m、u5.m 的 fd 值可分辨路径数 L路径越密包络越逼近理论瑞利分布u1.m 抽头数移动速度 v先换成 fd再做后续所有仿真u5.m 的输入参数作业里最常见的坑是设个顺眼的 fd 就开始跑完全不管信号带宽和采样率是否匹配。比如载频 2GHz、移动速度 30m/s 时 fd 约 200Hz若符号率 1ksps一个符号持续 1ms而相干时间约 2.1ms此时按慢衰落近似是合理的。速度提到 120m/s 后 fd 变成 800HzTc 降到 0.53ms低于符号周期信道在一个符号内就明显变化此时 u6.m 的误码率曲线会出现明显下压不去的错误底限要么改仿真符号率要么引入分集或均衡否则曲线没有任何参考价值。注意判断块衰落是否成立的快捷指标是 fd 乘符号周期 Ts超过 0.01 就应该在符号内部对信道系数做插值更新这套资源里的 Jakes 叠加实现天然支持这种更新只需要把时间轴步长调细。2.2 抽头延迟线的离散化与采样延迟换算射频上的多径现象下变频到复基带后表现为一个多抽头时变滤波器每个抽头对应一条可分辨延迟路径抽头系数是复高斯过程。教学级仿真不必上太多数学常见做法是先把抽头延迟换算成整数倍采样周期再对基带信号做延迟加权叠加。在老版脚本里这一步用的是 rayleighchan 函数matlab2019a 之后更推荐 comm.RayleighChannel两者参数含义接近但 PathDelays 的单位和 AveragePathGains 的功率定义并不完全一致从旧脚本迁移过来最容易在这里翻车。% 抽头延迟线初始化把时间延迟换算为采样点延迟 fs 1e6; % 系统采样率 1MHz L 3; % 抽头数量 3 tau [0 2e-6 5e-6]; % 路径时延单位秒 avg_power [0 -5 -10]; % 平均路径功率单位dB amp 10.^(avg_power / 20); % dB转线性幅度 tap_delay round(tau * fs); % 时延转整数采样点 for k 1:L fprintf(Tap%d: delay%d samples, amplitude%.4f\n, ... k, tap_delay(k), amp(k)); end这段代码的逻辑很直接先定义采样率再用 round 把小数时延转成整数采样偏移最后输出每个抽头的实际延迟和幅度。参数上最容易出错的是 dB 与线性幅度的换算路径功率 -5dB 对应的线性值是 10^(-5/20)≈0.56不是 0.315第二是当 tau 不能被 1/fs 整除时直接 round 会损失时延精度工程做法是进一步提高仿真采样率或者使用多相滤波器做分数延迟后者在资源里没有体现但这类课程设计后续要接 OFDM 链路时是绕不开的。3. u1.m~u7.m脚本拆解多径脉冲响应、多普勒功率谱与BPSK误码率仿真3.1 七个脚本的分工映射这套资源的七个脚本是按信道建模的认知顺序拆开的每个文件独立可跑输出对应一张或多张运行结果图。下面是我按文件名和结果图整理的对照关系跑的时候建议按编号顺序执行方便验证中间结果脚本仿真内容对应结果图核心变量u1.m多径信道冲激响应路径延迟与增益运行结果.jpgtau、路径功率u2.m不同延迟扩展下的信道频率响应运行结果2.jpg路径时延差、fsu3.m瑞利衰落包络与理论分布对比运行结果3.jpgfd、采样时间u4.m瑞利与莱斯衰落包络对比运行结果4.jpgK因子u5.m多普勒功率谱的U形曲线运行结果5.jpg、运行结果5·.jpgfd、fsu6.mBPSK在平坦衰落与频率选择性衰落下的BER运行结果6.jpgEb/N0、蒙特卡洛符号数u7.m补充对比或辅助计算脚本—视前序脚本而定要说明的是u7.m 没有独立结果图多数情况下它承担的是辅助函数或批量参数扫描功能也可能是绘制多径叠加效果时的公共代码段。matlab2019a 下直接按文件名顺序执行即可如果某个脚本提示找不到函数大概率是当前路径没切到解压目录或者脚本间存在隐式依赖检查 cd 和 addpath 即可。3.2 用Jakes正弦叠加法生成瑞利衰落包络u3.m 最常见的实现方式是 Jakes 正弦叠加用有限个不同频率的复正弦叠加逼近瑞利过程。这个方法比直接调用通信工具箱函数更直观能看到每条路径的多普勒频移是怎么贡献到最终包络上的。% u3.m 关键片段Jakes模型生成瑞利衰落包络并做统计验证 N_path 16; % 振荡器数量8~16即可逼近瑞利分布 fd 20; % 最大多普勒频移单位Hz fs 1000; % 采样频率单位Hz T 2; % 仿真时长单位s t (0:1/fs:T-1/fs).; alpha (2*pi*(0:N_path-1))/(N_path1); % 均匀分布的到达角 h zeros(size(t)); for k 1:N_path phi 2*pi*rand; % 随机初始相位避免同相叠加出尖峰 h h exp(1j*(2*pi*fd*cos(alpha(k))*t phi)); end h h / sqrt(N_path 1); % 归一化总平均功率为1 env abs(h); figure; histogram(env, 60, Normalization, pdf); hold on; x 0:0.02:max(env); plot(x, x .* exp(-x.^2/2), r-, LineWidth, 1.5); xlabel(包络幅度); ylabel(概率密度); legend(仿真直方图, 理论瑞利分布);逻辑说明for 循环里每条路径的瞬时频率是 fd·cos(alpha)alpha 在 0 到 2π 间均匀分布cos 映射后一部分路径的多普勒频率接近 0一部分接近 ±fd合成功率谱就呈现经典的 Jakes U 形。求和后的实虚部是两个独立近似高斯过程模值近似瑞利分布直方图与理论曲线的重合度受 N_path 和仿真时长的共同影响。参数上 N_path 低于 8 时直方图会明显锯齿化fd 提高到 100Hz 后如果 fs 仍只有 1000Hz每个正弦周期内只有 10 个采样点包络峰值会出现欠采样纹路这就是很多同学发现直方图尾巴翘起来的原因。3.3 u2.m与u6.m从频选信道响应到BPSK误码率验证u2.m 对应的是频率选择性衰落验证核心是画出两径信道的频响凹陷。延迟差为 5us 的两径信道频响上相邻凹陷间隔约为 1/5us200kHz这个周期可以通过改变 tau(2) 的值直观看到。% u2.m 关键片段两径信道的幅度响应 fs 1e6; f linspace(-fs/2, fs/2, 4096); tau [0 5e-6]; % 路径时延差5us gain_db [0 -3]; % 第二径衰减3dB amp 10.^(gain_db/20); H amp(1)*exp(-1j*2*pi*f*tau(1)) amp(2)*exp(-1j*2*pi*f*tau(2)); plot(f/1e3, 20*log10(abs(H))); xlabel(频率/kHz); ylabel(幅度响应/dB); title(两径信道频响); grid on;这段代码用复数指数叠加计算信道频域响应两个指数项分别代表两条路径的相位旋转。第二径衰减设为 -3dB 意味着该径功率是主径的一半频响零点不会完全陷到负无穷若两个径增益相同零点就会深到 -60dB 以下仿真图中会看到明显的周期性深谷。u6.m 的背景是把 BPSK 信号送进平坦瑞利信道统计误码率% u6.m 关键片段平坦瑞利衰落下的BPSK误码率 M 1e5; % 蒙特卡洛符号数 EbN0dB 0:2:20; EbN0lin 10.^(EbN0dB/10); ber_sim zeros(size(EbN0dB)); for idx 1:length(EbN0dB) data randi([0 1], M, 1); tx 2*data - 1; % BPSK: 0映射-1, 1映射1 h (randn(M,1) 1j*randn(M,1))/sqrt(2); % 复高斯信道系数 n (randn(M,1) 1j*randn(M,1))/sqrt(2); % 复噪声基底 noise n / sqrt(2*EbN0lin(idx)); rx tx .* h noise; rx_eq rx ./ h; % 理想信道均衡需要完整CSI bits_hat real(rx_eq) 0; ber_sim(idx) mean(bits_hat ~ data); end ber_theory 0.5*(1 - sqrt(EbN0lin ./ (1 EbN0lin))); semilogy(EbN0dB, ber_sim, o); hold on; semilogy(EbN0dB, ber_theory, -); grid on; legend(仿真, 理论); xlabel(Eb/N0 (dB)); ylabel(BER);参数说明M 是蒙特卡洛符号数低于 1e4 时高信噪比区域曲线会剧烈抖动h 的实虚部各用独立高斯变量并用 sqrt(2) 归一化保证 E[|h|^2]1噪声功率按 Eb/N0 换算复噪声基底再除 sqrt(2EbN0lin) 是因为 BPSK 每个符号携带 1bit。理论曲线用的是平坦瑞利信道闭式解 0.5(1-sqrt(gamma/(1gamma)))资源里的运行结果6.jpg 应该能看到仿真散点与理论曲线基本贴合如果高信噪比处明显偏离多半是噪声功率算错或信道功率没有归一化。4. 运行结果图判读从脉冲响应与多普勒谱反推信道参数4.1 运行结果.jpg从冲激响应直读延迟扩展与相干带宽运行结果.jpg 画的是多径信道冲激响应横轴通常是对时延或采样点数纵轴是各路径增益。图上能看到 3 到 4 根离散竖线每根代表一条可分辨路径。读图第一步是确认最大延迟差 Δτ然后用 1/Δτ 估算相干带宽。假设 Δτ2us则 Bc 约 500kHz如果发送信号带宽是 2MHz就能判断该信道是频率选择性衰落需要均衡器或 OFDM 子载波划分如果信号带宽只有 100kHz则按平坦衰落仿真即可。不少情况下 u1.m 输出的路径增益是负数 dB对应线性幅度小于 1这时要留意总信道功率是否归一化。若各径线性功率之和明显偏离 1后续 u6.m 的 BER 曲线会整体平移几个 dB看起来像性能变好或变差但实际上是能量标定错误。因此每次改路径参数后先在命令行里执行 sum(10.^(gain_db/10))确认总功率接近 1 再继续。4.2 运行结果5.jpgU形多普勒谱与相干时间的可视化验证运行结果5.jpg 和 5·.jpg 对应多普勒功率谱理论形态是中心频率处最高、两侧在 ±fd 位置急剧抬升的双尖峰 U 形曲线。实际仿真中因为采样点数和窗函数影响谱线不可能无限尖但边缘位置仍然清晰可辨。读图时直接量谱的频带宽度半宽处对应的频率数值就是最大多普勒频移 fd。如果在这个图上看到频谱边缘出现明显的翻折回卷说明仿真采样率 fs 小于 2fd不满足奈奎斯特条件必须把 fs 加倍重新运行。这个图也是验证相干时间最方便的手段。fd20Hz 时相干时间约 21msfd200Hz 时降到约 2.1ms频带随之展宽。课程答辩时如果能指着这张图说出“多普勒谱宽度与相干时间成反比测到半功率带宽就能回推移动速度”通常比单纯贴代码效果好得多。4.3 运行结果6.jpg与异常曲线排查表运行结果6.jpg 在资源中显示的是 BPSK 在两类信道下的误码率对比。正常情况应看到三条曲线AWGN 参考曲线下降最快平坦瑞利曲线下降缓慢频率选择性衰落曲线介于两者之间或误码率更高具体取决于均衡策略。若频选信道未加均衡它的误码率曲线会明显高于平坦瑞利且高信噪比时出现“地板”抬升。异常现象可能原因排查位置多普勒谱边缘翻折fs 小于 2fdu5.m 中提高 fs包络直方图与理论瑞利分布偏差大路径数太少或初始相位固定N_path 调大到16以上检查 rand 是否启动频响凹陷间隔与理论不符时延差 tau 与 fs 换算错误u2.m 中确认 round 后的采样点偏移BER 高信噪比区不掉信道功率未归一化或噪声功率重复叠加h 除以 sqrt(2) 的归一化处BER 曲线比理论低均衡器用了未来信息或误用信道信息u6.m 中检查 rx_eq 是否引入了非因果量这张表的最后一行值得展开仿真里做理想均衡时如果直接把 h 的当前符号值用于当前符号解调是最保守的慢衰落假设计算合理。但如果用整个序列的统计量或前后符号均值去替代当前系数就引入非因果信息误码率会人为降低答辩时很容易被问住。按资源原本的脚本结构u6.m 采用逐符号独立信道系数不存在这个问题但当你把它改造成块衰落仿真时必须注意信道更新间隔。5. 把多径衰落仿真改成工程级代码的三个收尾技巧5.1 用随机种子与功率归一化锁定可复现结果课程设计报告里最大的问题是图和代码对不上第二次运行结果飘忽。原因是 rand 和 randi 没有固定种子。在 u1.m 开头增加一行 rng(2019)所有脚本后续随机序列即可复现matlab2019a 中建议把 rng 放在每个脚本顶部而不是只在主脚本设置因为这个资源的七个脚本相互独立直接运行任一脚本都可能作为入口。5.2 把速度、载频、多普勒频移封装成参数结构体直接在脚本里改 fd 很容易漏掉其他依赖项我习惯在最前面加一个参数结构体把物理量与仿真量一次性换算% 参数统一入口修改这里后面所有脚本跟着变 params.fc 2e9; % 载频 2GHz params.v 30; % 移动速度 30m/s params.fd params.v * params.fc / 3e8; % 最大多普勒频移 params.fs 10*params.fd; % 采样率至少10倍于多普勒这样从 u3.m 到 u5.m 都引用 params 结构体里的 fd既不必逐个脚本改还能避免速度与载频不匹配造成的物理失真。后续若要把仿真下沉到 Simulink 或导出数据给深度学习模型做信道预测这个结构体也可作为接口直接透传。5.3 做 fd 批量扫描输出相干时间与BER的关系表最后一个技巧是把单点仿真改成参数扫描一次性观察多普勒对系统的影响。简单做法是循环 fd 取值在每个 fd 下调用 u3.m 的核心函数计算包络自相关再记录 0.5 相关点的宽度作为相干时间估计值% check_coherence.m 辅助脚本示例 fd_list [20 50 100 200 400]; for i 1:length(fd_list) fd fd_list(i); % 调用Jakes模型生成信道包络代码与u3.m一致 % 这里省略具体实现直接输出Tc估计 tc_est 0.423 / fd; fprintf(fd%4dHz理论Tc%.3fms\n, fd, tc_est*1000); end输出结果与理论公式对照能快速确认仿真时间长度是否足够覆盖信道变化周期如果仿真时长比 Tc 短很多信道统计特性就无法体现这和“仿真发散”是另一类常见问题。做这种批量扫描时注意单次蒙特卡洛符号数不要缩减否则曲线抖动会被误判为信道特性差异。把这段辅助逻辑写在 u7.m 里正好让整套资源从教学演示升级为可做参数研究的工具链。本文还有配套的精品资源点击获取 
长河君 · 2026/9/11 18:27:43 · 阅读 1.2k 多径衰落信道Matlab仿真：从多普勒频移到瑞利衰落的脚本化实现
简介面向通信专业本科、硕士生及教研人员的Matlab仿真资源包围绕多径衰落信道模型及多普勒频移效应展开能够帮助读者掌握无线信道仿真基础流程。压缩包共16个文件包含7个可运行的MATLAB脚本、8张运行结果JPG图像和1份Word说明文档整体大小仅382KB脚本覆盖参数配置、多径叠加、信道响应计算到结果绘制等环节便于对照图片检查每一步输出。文档对多径传输机理和仿真设计做了梳理适合零基础或初学通信仿真的读者按步骤复现。目前已有204人学习浏览使用Matlab 2019a即可直接运行也可在课程设计、毕业设计或课堂教学中作为基础参考。读者通过修改时延、增益、多普勒频移等参数能进一步观察不同信道条件下的衰落变化比如幅度波动与频率选择性衰落为后续通信系统级仿真打下基础。1. 多径衰落信道matlab仿真从多径传输直通多普勒频移的脚本化拆解通信多径衰落信道matlab仿真这套资源主题落在多径传输与多普勒频移两条主线上。从u1.m到u7.m七个脚本把信道脉冲响应、瑞利包络、多普勒功率谱、误码率曲线依次跑通并附六张运行结果图供对照。适合本科、硕士做通信原理课程设计把教材里的相干带宽、相干时间变成可调参数的仿真图也适合工作两三年的算法工程师在做系统级链路仿真前确认自己对抽头延迟线和最大多普勒频移的理解没有偏差。重点不是代码本身多复杂而是能从这些图中读出延迟扩展怎么毁掉频响、多普勒频移怎么把星座点推成圆环这才是后续改参数与扩展仿真的起点。2. 多普勒频移与延迟扩展多径衰落信道仿真的两类基础参数2.1 先判定频繁选择性还是时间选择性再决定改哪个脚本变量多径衰落仿真的第一步不是写代码而是确认信道落在哪个区间。延迟扩展决定频率选择性多普勒频移决定时间选择性两者来自同一份信道脉冲响应但在脚本里由完全独立的变量控制。拆这类教学资源时我会先把参数和脚本的对应关系列出来改任何系数前先对着表查一次能省掉大量无效调试时间。参数与信道建模的关系对应脚本位置相干带宽 Bc ≈ 1/最大延迟扩展Bc大于信号带宽是平坦衰落否则是频率选择性衰落u1.m、u2.m 的路径延迟差最大多普勒频移 fd v·fc/cfd决定相干时间TcTc小于符号周期时信道在符号内快速变化u3.m、u5.m 的 fd 值可分辨路径数 L路径越密包络越逼近理论瑞利分布u1.m 抽头数移动速度 v先换成 fd再做后续所有仿真u5.m 的输入参数作业里最常见的坑是设个顺眼的 fd 就开始跑完全不管信号带宽和采样率是否匹配。比如载频 2GHz、移动速度 30m/s 时 fd 约 200Hz若符号率 1ksps一个符号持续 1ms而相干时间约 2.1ms此时按慢衰落近似是合理的。速度提到 120m/s 后 fd 变成 800HzTc 降到 0.53ms低于符号周期信道在一个符号内就明显变化此时 u6.m 的误码率曲线会出现明显下压不去的错误底限要么改仿真符号率要么引入分集或均衡否则曲线没有任何参考价值。注意判断块衰落是否成立的快捷指标是 fd 乘符号周期 Ts超过 0.01 就应该在符号内部对信道系数做插值更新这套资源里的 Jakes 叠加实现天然支持这种更新只需要把时间轴步长调细。2.2 抽头延迟线的离散化与采样延迟换算射频上的多径现象下变频到复基带后表现为一个多抽头时变滤波器每个抽头对应一条可分辨延迟路径抽头系数是复高斯过程。教学级仿真不必上太多数学常见做法是先把抽头延迟换算成整数倍采样周期再对基带信号做延迟加权叠加。在老版脚本里这一步用的是 rayleighchan 函数matlab2019a 之后更推荐 comm.RayleighChannel两者参数含义接近但 PathDelays 的单位和 AveragePathGains 的功率定义并不完全一致从旧脚本迁移过来最容易在这里翻车。% 抽头延迟线初始化把时间延迟换算为采样点延迟 fs 1e6; % 系统采样率 1MHz L 3; % 抽头数量 3 tau [0 2e-6 5e-6]; % 路径时延单位秒 avg_power [0 -5 -10]; % 平均路径功率单位dB amp 10.^(avg_power / 20); % dB转线性幅度 tap_delay round(tau * fs); % 时延转整数采样点 for k 1:L fprintf(Tap%d: delay%d samples, amplitude%.4f\n, ... k, tap_delay(k), amp(k)); end这段代码的逻辑很直接先定义采样率再用 round 把小数时延转成整数采样偏移最后输出每个抽头的实际延迟和幅度。参数上最容易出错的是 dB 与线性幅度的换算路径功率 -5dB 对应的线性值是 10^(-5/20)≈0.56不是 0.315第二是当 tau 不能被 1/fs 整除时直接 round 会损失时延精度工程做法是进一步提高仿真采样率或者使用多相滤波器做分数延迟后者在资源里没有体现但这类课程设计后续要接 OFDM 链路时是绕不开的。3. u1.m~u7.m脚本拆解多径脉冲响应、多普勒功率谱与BPSK误码率仿真3.1 七个脚本的分工映射这套资源的七个脚本是按信道建模的认知顺序拆开的每个文件独立可跑输出对应一张或多张运行结果图。下面是我按文件名和结果图整理的对照关系跑的时候建议按编号顺序执行方便验证中间结果脚本仿真内容对应结果图核心变量u1.m多径信道冲激响应路径延迟与增益运行结果.jpgtau、路径功率u2.m不同延迟扩展下的信道频率响应运行结果2.jpg路径时延差、fsu3.m瑞利衰落包络与理论分布对比运行结果3.jpgfd、采样时间u4.m瑞利与莱斯衰落包络对比运行结果4.jpgK因子u5.m多普勒功率谱的U形曲线运行结果5.jpg、运行结果5·.jpgfd、fsu6.mBPSK在平坦衰落与频率选择性衰落下的BER运行结果6.jpgEb/N0、蒙特卡洛符号数u7.m补充对比或辅助计算脚本—视前序脚本而定要说明的是u7.m 没有独立结果图多数情况下它承担的是辅助函数或批量参数扫描功能也可能是绘制多径叠加效果时的公共代码段。matlab2019a 下直接按文件名顺序执行即可如果某个脚本提示找不到函数大概率是当前路径没切到解压目录或者脚本间存在隐式依赖检查 cd 和 addpath 即可。3.2 用Jakes正弦叠加法生成瑞利衰落包络u3.m 最常见的实现方式是 Jakes 正弦叠加用有限个不同频率的复正弦叠加逼近瑞利过程。这个方法比直接调用通信工具箱函数更直观能看到每条路径的多普勒频移是怎么贡献到最终包络上的。% u3.m 关键片段Jakes模型生成瑞利衰落包络并做统计验证 N_path 16; % 振荡器数量8~16即可逼近瑞利分布 fd 20; % 最大多普勒频移单位Hz fs 1000; % 采样频率单位Hz T 2; % 仿真时长单位s t (0:1/fs:T-1/fs).; alpha (2*pi*(0:N_path-1))/(N_path1); % 均匀分布的到达角 h zeros(size(t)); for k 1:N_path phi 2*pi*rand; % 随机初始相位避免同相叠加出尖峰 h h exp(1j*(2*pi*fd*cos(alpha(k))*t phi)); end h h / sqrt(N_path 1); % 归一化总平均功率为1 env abs(h); figure; histogram(env, 60, Normalization, pdf); hold on; x 0:0.02:max(env); plot(x, x .* exp(-x.^2/2), r-, LineWidth, 1.5); xlabel(包络幅度); ylabel(概率密度); legend(仿真直方图, 理论瑞利分布);逻辑说明for 循环里每条路径的瞬时频率是 fd·cos(alpha)alpha 在 0 到 2π 间均匀分布cos 映射后一部分路径的多普勒频率接近 0一部分接近 ±fd合成功率谱就呈现经典的 Jakes U 形。求和后的实虚部是两个独立近似高斯过程模值近似瑞利分布直方图与理论曲线的重合度受 N_path 和仿真时长的共同影响。参数上 N_path 低于 8 时直方图会明显锯齿化fd 提高到 100Hz 后如果 fs 仍只有 1000Hz每个正弦周期内只有 10 个采样点包络峰值会出现欠采样纹路这就是很多同学发现直方图尾巴翘起来的原因。3.3 u2.m与u6.m从频选信道响应到BPSK误码率验证u2.m 对应的是频率选择性衰落验证核心是画出两径信道的频响凹陷。延迟差为 5us 的两径信道频响上相邻凹陷间隔约为 1/5us200kHz这个周期可以通过改变 tau(2) 的值直观看到。% u2.m 关键片段两径信道的幅度响应 fs 1e6; f linspace(-fs/2, fs/2, 4096); tau [0 5e-6]; % 路径时延差5us gain_db [0 -3]; % 第二径衰减3dB amp 10.^(gain_db/20); H amp(1)*exp(-1j*2*pi*f*tau(1)) amp(2)*exp(-1j*2*pi*f*tau(2)); plot(f/1e3, 20*log10(abs(H))); xlabel(频率/kHz); ylabel(幅度响应/dB); title(两径信道频响); grid on;这段代码用复数指数叠加计算信道频域响应两个指数项分别代表两条路径的相位旋转。第二径衰减设为 -3dB 意味着该径功率是主径的一半频响零点不会完全陷到负无穷若两个径增益相同零点就会深到 -60dB 以下仿真图中会看到明显的周期性深谷。u6.m 的背景是把 BPSK 信号送进平坦瑞利信道统计误码率% u6.m 关键片段平坦瑞利衰落下的BPSK误码率 M 1e5; % 蒙特卡洛符号数 EbN0dB 0:2:20; EbN0lin 10.^(EbN0dB/10); ber_sim zeros(size(EbN0dB)); for idx 1:length(EbN0dB) data randi([0 1], M, 1); tx 2*data - 1; % BPSK: 0映射-1, 1映射1 h (randn(M,1) 1j*randn(M,1))/sqrt(2); % 复高斯信道系数 n (randn(M,1) 1j*randn(M,1))/sqrt(2); % 复噪声基底 noise n / sqrt(2*EbN0lin(idx)); rx tx .* h noise; rx_eq rx ./ h; % 理想信道均衡需要完整CSI bits_hat real(rx_eq) 0; ber_sim(idx) mean(bits_hat ~ data); end ber_theory 0.5*(1 - sqrt(EbN0lin ./ (1 EbN0lin))); semilogy(EbN0dB, ber_sim, o); hold on; semilogy(EbN0dB, ber_theory, -); grid on; legend(仿真, 理论); xlabel(Eb/N0 (dB)); ylabel(BER);参数说明M 是蒙特卡洛符号数低于 1e4 时高信噪比区域曲线会剧烈抖动h 的实虚部各用独立高斯变量并用 sqrt(2) 归一化保证 E[|h|^2]1噪声功率按 Eb/N0 换算复噪声基底再除 sqrt(2EbN0lin) 是因为 BPSK 每个符号携带 1bit。理论曲线用的是平坦瑞利信道闭式解 0.5(1-sqrt(gamma/(1gamma)))资源里的运行结果6.jpg 应该能看到仿真散点与理论曲线基本贴合如果高信噪比处明显偏离多半是噪声功率算错或信道功率没有归一化。4. 运行结果图判读从脉冲响应与多普勒谱反推信道参数4.1 运行结果.jpg从冲激响应直读延迟扩展与相干带宽运行结果.jpg 画的是多径信道冲激响应横轴通常是对时延或采样点数纵轴是各路径增益。图上能看到 3 到 4 根离散竖线每根代表一条可分辨路径。读图第一步是确认最大延迟差 Δτ然后用 1/Δτ 估算相干带宽。假设 Δτ2us则 Bc 约 500kHz如果发送信号带宽是 2MHz就能判断该信道是频率选择性衰落需要均衡器或 OFDM 子载波划分如果信号带宽只有 100kHz则按平坦衰落仿真即可。不少情况下 u1.m 输出的路径增益是负数 dB对应线性幅度小于 1这时要留意总信道功率是否归一化。若各径线性功率之和明显偏离 1后续 u6.m 的 BER 曲线会整体平移几个 dB看起来像性能变好或变差但实际上是能量标定错误。因此每次改路径参数后先在命令行里执行 sum(10.^(gain_db/10))确认总功率接近 1 再继续。4.2 运行结果5.jpgU形多普勒谱与相干时间的可视化验证运行结果5.jpg 和 5·.jpg 对应多普勒功率谱理论形态是中心频率处最高、两侧在 ±fd 位置急剧抬升的双尖峰 U 形曲线。实际仿真中因为采样点数和窗函数影响谱线不可能无限尖但边缘位置仍然清晰可辨。读图时直接量谱的频带宽度半宽处对应的频率数值就是最大多普勒频移 fd。如果在这个图上看到频谱边缘出现明显的翻折回卷说明仿真采样率 fs 小于 2fd不满足奈奎斯特条件必须把 fs 加倍重新运行。这个图也是验证相干时间最方便的手段。fd20Hz 时相干时间约 21msfd200Hz 时降到约 2.1ms频带随之展宽。课程答辩时如果能指着这张图说出“多普勒谱宽度与相干时间成反比测到半功率带宽就能回推移动速度”通常比单纯贴代码效果好得多。4.3 运行结果6.jpg与异常曲线排查表运行结果6.jpg 在资源中显示的是 BPSK 在两类信道下的误码率对比。正常情况应看到三条曲线AWGN 参考曲线下降最快平坦瑞利曲线下降缓慢频率选择性衰落曲线介于两者之间或误码率更高具体取决于均衡策略。若频选信道未加均衡它的误码率曲线会明显高于平坦瑞利且高信噪比时出现“地板”抬升。异常现象可能原因排查位置多普勒谱边缘翻折fs 小于 2fdu5.m 中提高 fs包络直方图与理论瑞利分布偏差大路径数太少或初始相位固定N_path 调大到16以上检查 rand 是否启动频响凹陷间隔与理论不符时延差 tau 与 fs 换算错误u2.m 中确认 round 后的采样点偏移BER 高信噪比区不掉信道功率未归一化或噪声功率重复叠加h 除以 sqrt(2) 的归一化处BER 曲线比理论低均衡器用了未来信息或误用信道信息u6.m 中检查 rx_eq 是否引入了非因果量这张表的最后一行值得展开仿真里做理想均衡时如果直接把 h 的当前符号值用于当前符号解调是最保守的慢衰落假设计算合理。但如果用整个序列的统计量或前后符号均值去替代当前系数就引入非因果信息误码率会人为降低答辩时很容易被问住。按资源原本的脚本结构u6.m 采用逐符号独立信道系数不存在这个问题但当你把它改造成块衰落仿真时必须注意信道更新间隔。5. 把多径衰落仿真改成工程级代码的三个收尾技巧5.1 用随机种子与功率归一化锁定可复现结果课程设计报告里最大的问题是图和代码对不上第二次运行结果飘忽。原因是 rand 和 randi 没有固定种子。在 u1.m 开头增加一行 rng(2019)所有脚本后续随机序列即可复现matlab2019a 中建议把 rng 放在每个脚本顶部而不是只在主脚本设置因为这个资源的七个脚本相互独立直接运行任一脚本都可能作为入口。5.2 把速度、载频、多普勒频移封装成参数结构体直接在脚本里改 fd 很容易漏掉其他依赖项我习惯在最前面加一个参数结构体把物理量与仿真量一次性换算% 参数统一入口修改这里后面所有脚本跟着变 params.fc 2e9; % 载频 2GHz params.v 30; % 移动速度 30m/s params.fd params.v * params.fc / 3e8; % 最大多普勒频移 params.fs 10*params.fd; % 采样率至少10倍于多普勒这样从 u3.m 到 u5.m 都引用 params 结构体里的 fd既不必逐个脚本改还能避免速度与载频不匹配造成的物理失真。后续若要把仿真下沉到 Simulink 或导出数据给深度学习模型做信道预测这个结构体也可作为接口直接透传。5.3 做 fd 批量扫描输出相干时间与BER的关系表最后一个技巧是把单点仿真改成参数扫描一次性观察多普勒对系统的影响。简单做法是循环 fd 取值在每个 fd 下调用 u3.m 的核心函数计算包络自相关再记录 0.5 相关点的宽度作为相干时间估计值% check_coherence.m 辅助脚本示例 fd_list [20 50 100 200 400]; for i 1:length(fd_list) fd fd_list(i); % 调用Jakes模型生成信道包络代码与u3.m一致 % 这里省略具体实现直接输出Tc估计 tc_est 0.423 / fd; fprintf(fd%4dHz理论Tc%.3fms\n, fd, tc_est*1000); end输出结果与理论公式对照能快速确认仿真时间长度是否足够覆盖信道变化周期如果仿真时长比 Tc 短很多信道统计特性就无法体现这和“仿真发散”是另一类常见问题。做这种批量扫描时注意单次蒙特卡洛符号数不要缩减否则曲线抖动会被误判为信道特性差异。把这段辅助逻辑写在 u7.m 里正好让整套资源从教学演示升级为可做参数研究的工具链。本文还有配套的精品资源点击获取