简介：面向通信工程、电子信息类专业的学生与课程设计使用者，这套MATLAB源码程序基于matlab2021a，用于瑞利信道仿真与多径衰落建模，可帮助理解信道参数设置及仿真流程。包内共3个文件，包括2个m源码文件和1个avi操作录像，压缩包整体约553KB。源码在main.m中预设了信号长度、最大多普勒频移、载波频率、多径时延与功率等关键参数，并通过正弦叠加信号作为输入来观察信道响应对传输信号的影响；Rayl.m则完成瑞利信道核心建模。操作录像使用Windows Media Player播放，详细演示了程序运行、结果输出以及MATLAB左侧当前文件夹路径的配置方法，可有效避免因路径设置错误导致的运行失败。目前已有647人学习下载，适合需要快速上手信道仿真实验、参考源码参数配置或核对仿真结果的初学者。
1. 瑞利信道仿真MATLAB源码程序的适用场景与版本选型
瑞利信道仿真MATLAB源码程序是无线通信链路设计里最常被拿出来反复改的一套代码，它用来模拟基站与移动终端之间没有视距路径时的多径衰落。通信专业研究生做物理层算法，射频测试工程师搭基带验证平台，或者做车联网、WiFi、LoRa 这类移动场景的工程师，都会用到这套仿真。使用版本 matlab2021a 是个稳妥选择：系统对象接口稳定，Communications Toolbox 和 Signal Processing Toolbox 功能完整，内存占用比 R2023b 以后版本更保守，而且大部分高校实验课还停留在这个环境上。标题里的“操作录像”是个很实用的附加资源，能让新接触仿真的人按真实操作路径跑通，而不是对着静态 README 猜。
2. 瑞利信道仿真的核心模型与MATLAB 2021a实现
2.1 瑞利衰落模型的统计特性与仿真原理
在移动通信里，信号从发射机到接收机往往经过建筑物、地面、车辆等反射体，形成多条路径。当这些反射分量在接收端叠加时，如果每条路径的相位随机并且反射体数量足够多，根据中心极限定理，接收信号的同相分量和正交分量都会接近高斯分布。用复数表示信道增益 h = x + i*y，那么 x 和 y 都是零均值、方差相等的高斯随机变量，包络 r = |h| 就服从瑞利分布，相位服从均匀分布。这是瑞利信道仿真最基本的统计依据。
很多刚拿到源码程序的人会把瑞利信道误当成一个固定衰减系数，这是错的。真正的瑞利信道应该是一段随时间变化的复数序列，包络在深衰落点可能快速跌到接近 0，而相位也会随机跳变。这种时间相关性主要由多普勒扩展决定，所以仿真必须引入多普勒频谱形状。MATLAB 2021a 的 comm.RayleighChannel 系统对象内部使用滤波器法实现多普勒谱，默认是经典 Jakes 谱的离散近似。相比老式 rayleighchan 函数，这个系统对象更适合按帧处理数据，也更容易和 OFDM 或 MIMO 的系统级仿真集成。
2.2 用MATLAB 2021a生成瑞利信道的最小可运行代码
这里给出一段可以直接复制运行的最小实现，目的是先建立从信道对象到输出波形的最小闭环：
% raleigh_demo.m % 适用MATLAB 2021a，生成单径瑞利信道并绘制包络与功率谱 fs = 1000; % 采样率，单位Hz fd = 10; % 最大多普勒频移，单位Hz Ts = 1/fs; % 采样周期 t = (0:fs*2-1)*Ts; % 2秒时间向量 rayChan = comm.RayleighChannel(... 'SampleRate', fs, ... 'MaximumDopplerShift', fd, ... 'RandomStream', 'mt19937ar with seed', ... 'Seed', 23); % 固定随机种子，保证结果可复现 [rxSig, pathGains] = rayChan(ones(length(t),1)); figure; subplot(2,1,1); plot(t, abs(rxSig)); xlabel('时间 (s)'); ylabel('包络幅度'); title('瑞利信道包络'); subplot(2,1,2); [psd, f] = pwelch(rxSig, [], [], [], fs); plot(f, 10*log10(psd)); xlabel('频率 (Hz)'); ylabel('功率谱密度 (dB/Hz)'); title('多普勒谱');代码先把采样率和多普勒参数提出来，再创建瑞利信道对象。输入一个全 1 列向量，是为了让信道本身不被调制信息干扰，这样输出 rxSig 能直接反映信道包络。pathGains 是真实的信道增益复数序列，它才是后续误码率计算要用的量。pwelch 计算功率谱密度，可以直观看到多普勒谱的形状是否符合 Jakes 谱特征。
参数说明：
- SampleRate 设置信道工作的采样率，单位必须是 Hz，不能小于最大多普勒频移的两倍，否则多普勒谱会发生混叠；
- MaximumDopplerShift 由移动速度和工作频段决定，典型值 5 到 100 Hz；
- RandomStream 和 Seed 配合使用能复现同一组随机衰落，做蒙特卡洛仿真时建议每个信噪比点都换不同 Seed，避免相关性导致误码率曲线抖动；
- pathGains 在单径信道里是一个复数列向量，在多径信道里则是矩阵，列数代表路径数。
2.3 关键参数：多普勒频移、采样率、路径时延
| 参数 | 含义 | 典型值 | 影响 |
|---|---|---|---|
| MaximumDopplerShift | 最大多普勒频移 | 1-100 Hz | 决定包络的起伏速度，影响深度衰落持续时间 |
| SampleRate | 信道采样率 | 1000 Hz 到 20 MHz | 限制多普勒谱范围与时延分辨率 |
| PathDelays | 各路径相对时延 | [0 1e-6 2e-6] 秒 | 决定频率选择性衰落特性 |
| AveragePathGains | 各路径平均功率 | [0 -3 -6] dB | 决定功率延迟分布与总输出功率 |
这组参数是调整仿真行为最常用的入口。很多源码程序默认只设单径，也就是把 PathDelays 设为 0，AveragePathGains 设为 0，输出的是平坦瑞利衰落。但实际 OFDM 系统更关心多径，因为多径时延会导致不同子载波上增益不同。需要特别留意，SampleRate 是后续离散基带信号的基础，不是符号率。如果 OFDM 系统采样率是 20 MHz，而源码程序里仍用 1000 Hz，那么时延 0.1us 会在仿真中被错误量化，甚至直接被忽略。
3. 源码程序的文件结构、操作录像与运行方法
3.1 源码文件组织与模块划分
网上和课题组里流传的瑞利信道仿真MATLAB源码程序，文件结构大多围绕主脚本、参数配置、信道生成、误码率计算和画图来组织。我比较推荐下面这种布局：
rayleigh_sim/ ├── main.m % 主入口，依次调用配置、仿真、绘图 ├── config.m % 所有可调参数集中在这里 ├── create_rayleigh_channel.m % 封装信道对象创建 ├── ber_sim.m % 误码率蒙特卡洛仿真 ├── plot_results.m % 画图与保存结果 ├── README.md % 版本要求和运行说明 └── video/ └── demo.mp4 % 操作录像main.m 是总入口，运行时先执行 config.m 把参数写进工作区，然后调用 create_rayleigh_channel.m 生成信道对象，接着用 ber_sim.m 做循环仿真，最后调用 plot_results.m 画图。把参数单独放一个文件的好处是，做参数扫描时不用反复修改主脚本。操作录像演示的就是这一套流程，你会在录像里看到作者如何修改 config.m 里的 fd，然后重新运行 main.m 观察曲线变化。
如果你的源码程序里没有 config.m，而是把所有参数堆在主脚本开头，那说明这个程序偏实验性质，结构不够工程化。建议你在读完源码后自己拆一次，把参数提取出去，后续跑多组对比会方便很多。
3.2 在MATLAB 2021a中运行仿真并查看操作录像
拿到源码程序后，先在 MATLAB 2021a 命令行里执行：
% 切换到源码程序根目录，路径按实际位置修改 cd('D:\rayleigh_sim'); % 清理工作区和图形窗口，避免上次运行残留数据干扰 clear variables; close all; clc; % 加载配置参数到工作区 run('config.m'); % 运行主仿真脚本 run('main.m');这段命令展示了最安全的启动顺序。clear variables 和 close all 可以清掉旧变量和旧图窗，尤其当你在同一个工作区里反复尝试不同 fd 时，残留的 pathGains 会让某个中间变量没被覆盖，最后画出的曲线张冠李戴。config.m 使用 run 而不是 function 调用，是因为配置脚本通常只是给变量赋值，不需要返回值。
操作录像建议在运行 main.m 之前先完整看一遍。录像里通常包含三个关键信息：第一，作者在哪个位置修改 fd；第二，运行后命令行窗口出现什么输出提示；第三，绘制的图形窗口里哪些曲线是真实的，哪些是辅助参考。比如有的源码程序会在控制台打印“EbNo = 10 dB, BER = 3.2e-3”，录像里能看到每个循环的进度，这对你判断程序有没有卡死有很大帮助。
3.3 常见运行错误与版本兼容性问题
| 错误信息 | 原因 | 解决办法 |
|---|---|---|
| Unrecognized function or variable 'comm.RayleighChannel' | 缺少 Communications Toolbox | 执行 ver 检查工具箱，安装对应工具包 |
| Invalid maximum Doppler shift | fd 超过 SampleRate/2 | 提高 SampleRate 或降低 fd |
| Error using step; input signal length must be... | 系统对象在有状态情况下被重复调用 | 使用 release(rayChan) 重置状态 |
| RayChan output size mismatch | 多径输入输出尺寸理解错误 | 检查输入信号列数是否等于 PathDelays 长度 |
从旧版本迁移到 MATLAB 2021a 时，最典型的问题是老代码里用了 rayleighchan 函数，而新版本推荐使用 comm.RayleighChannel。前者是旧通信工具箱接口，后者是系统对象，两者参数结构完全不同。如果源码程序里出现 “rayleighchan” 关键词，说明核心代码还没更新，而操作录像却用 2021a 运行，这时你需要把信道创建部分全部替换成系统对象写法。
还有一种隐蔽错误是随机种子固定后，多个信噪比点复用同一段衰落序列，这会让 BER 曲线出现错误的相关性。排查的方法是修改 Seed 或直接删除 Seed 参数，看曲线是否发生变化。如果曲线完全不变，说明程序内部可能用了全局变量缓存信道。这种问题在源码程序里很难一眼看出来，通常要配合断点调试才能定位。
4. 仿真参数对瑞利衰落影响的分析与调优
4.1 多普勒频移对衰落速率的影响
最大多普勒频移 fd 的计算公式是 fd = v / lambda，其中 v 是移动速度，lambda 是载波波长。2.4 GHz 下载波波长约为 0.125 米，所以步行速度 1.4 m/s 对应的 fd 大约 11 Hz，而车速 30 m/s 对应 fd 约为 240 Hz。在 MATLAB 2021a 中设置不同 fd，包络曲线会呈现明显的时间尺度差异。
% doppler_test.m fdList = [5 20 50]; t = (0:1999)/1000; for k = 1:length(fdList) rayChan = comm.RayleighChannel(... 'SampleRate', 1000, ... 'MaximumDopplerShift', fdList(k), ... 'RandomStream', 'mt19937ar with seed', ... 'Seed', 1); [~, pathGains] = rayChan(ones(2000,1)); subplot(3,1,k); plot(t, abs(pathGains)); title(['fd = ', num2str(fdList(k)), ' Hz']); end运行这段代码后，你会看到 fd=50 Hz 的包络在 2 秒内经历了十几次深衰落，而 fd=5 Hz 时包络变化很慢。这直接影响系统的交织器设计，如果你的符号周期是 1 ms，那么 fd=10 Hz 时信道相干时间大约 50 ms，符号间衰落相关性很高。做性能仿真时不能只选一个固定 fd，而应把 fd 作为扫描参数，才能看出系统在不同移动速度下的鲁棒性。
多普勒频谱的形状也是验证仿真器是否正确的重要手段。经典 Jakes 谱在 [−fd, fd] 范围内两端高、中间平坦，接近 U 形。用 pwelch 观察时，如果频谱能量在超出 fd 的位置仍然很明显，说明 SampleRate 不够，或者多普勒整形滤波器阶数太低。MATLAB 默认实现已经能保证谱形，但当你自定义采样率和 fd 组合时，还是建议把频谱画出来检查一次。
4.2 路径时延与频率选择性衰落的调整
多径时延扩展是判断信道是平坦衰落还是频率选择性衰落的核心指标。在 OFDM 仿真中，如果循环前缀长度大于最大时延扩展，可以近似认为每个子载波经历独立平坦衰落。常用方法是给 comm.RayleighChannel 设置 PathDelays 和 AveragePathGains：
% multipath_config.m fs = 20e6; % OFDM基带采样率 tdelays = [0 0.1e-6 0.3e-6]; % 路径时延，单位秒 tgains = [0 -3 -6]; % 平均路径增益，单位dB rayChan = comm.RayleighChannel(... 'SampleRate', fs, ... 'PathDelays', tdelays, ... 'AveragePathGains', tgains, ... 'MaximumDopplerShift', 100);这段配置描述了一条主径和两条延迟分别为 0.1us 和 0.3us 的反射径。时延差 0.3us 对应 6 个采样周期，在频率轴上会产生间隔约 3.3 MHz 的周期性凹陷，所以不同 OFDM 子载波信噪比差异会很大。要验证频率选择性，可以输入一个线性调频信号，然后观察频域增益曲线。
需要注意 PathDelays 的解析度受采样率限制。如果两个时延之间的间隔小于 1/fs，MATLAB 会默认把它们量化到最近的采样点，结果就等效成两条路径合并。因此，当你需要仿真 0.05us 时延差时，fs 至少要到 20 MHz，而不是只提高 OFDM 的采样率。很多源码程序在带宽参数修改后，完全忘掉同步修改 PathDelays，导致时延差和实际物理场景对不上。
4.3 平均路径增益与功率归一化
在多径信道仿真里，每一条路径都有平均功率，AveragePathGains 的单位是 dB。如果不做归一化，信道输出的总功率会随路径数目增加而增加，导致后续噪声叠加时 Eb/N0 计算偏差。常见做法是在创建信道对象前先归一化：
% normalize_gain.m xdB = [0 -3 -6 -9]; xLin = 10.^(xdB/10); xLinNorm = xLin / sum(xLin); newdB = 10*log10(xLinNorm); rayChan = comm.RayleighChannel(... 'SampleRate', 20e6, ... 'PathDelays', [0 0.1e-6 0.2e-6 0.3e-6], ... 'AveragePathGains', newdB);这段代码把各径平均功率之和缩放为 1。归一化后，信道在所有随机实现上的平均输出功率约等于输入功率，但单次时隙里的瞬时功率仍会波动，因为瑞利衰落的包络本身就服从瑞利分布，方差不为零。很多新手看到某一次仿真里输出功率是 1.5 就认为是 bug，实际上这属于正常波动，统计意义下才是 1。
功率归一化直接关系到误码率曲线的信噪比标定。如果你不归一化就按 N0 加噪声，等效于改变了每个比特的能量，会让仿真结果与理论公式产生系统性偏移。源码程序里如果包含归一化代码，通常会在注释里写明“Normalize total path power to unity”，没有的话建议自己补上。
5. 用误码率曲线验证瑞利信道仿真源码程序
5.1 单径瑞利信道下的理论误码率
判断一个瑞利信道仿真MATLAB源码程序是否正确，最有力的证据是让误码率曲线和理论公式重合。平坦瑞利衰落信道下，BPSK 的平均误码率公式为：
Pe = 0.5 * (1 - sqrt(EbN0_bar / (1 + EbN0_bar)))
其中 EbN0_bar 是平均比特信噪比。这里要注意，这个公式假设接收端已知信道增益并能做理想均衡，或者采用差分检测。如果不做均衡，误码率会明显变差，因为衰落造成深度低功率点让符号严重损坏。
和 AWGN 信道不同，瑞利信道误码率曲线在信噪比升高时会逐渐变缓，不会像 AWGN 那样陡峭下降。这是因为包络出现深衰落的概率不会随噪声降低而减少，形成误码率的“地板效应”。在你的源码程序里看到这种形状，反而是瑞利信道工作正常的信号。
5.2 用MATLAB 2021a实现BPSK误码率仿真
下面这段代码是完整的 BPSK 在单径瑞利信道下的蒙特卡洛仿真，可以直接替换源码程序里的误码率部分：
% ber_rayleigh_bpsk.m EbNo = 0:2:20; % 平均比特信噪比，dB numBits = 1e6; % 每个信噪比点的比特数 ber = zeros(size(EbNo)); data = randi([0 1], numBits, 1); % 固定数据序列 mod = 2*data - 1; % BPSK映射: 0->-1, 1->1 for idx = 1:length(EbNo) % 每次创建独立信道对象，重置随机种子 rayChan = comm.RayleighChannel(... 'SampleRate', 1000, ... 'MaximumDopplerShift', 20, ... 'RandomStream', 'mt19937ar with seed', ... 'Seed', 100+idx); [~, pathGains] = rayChan(mod); % 信道增益与发送信号同长度 % 计算噪声功率 Es = mean(abs(mod).^2); EbN0Lin = 10^(EbNo(idx)/10); N0 = Es / EbN0Lin; noise = sqrt(N0/2) * (randn(numBits,1) + 1j*randn(numBits,1)); % 过信道并加噪 rx = mod .* pathGains + noise; % 理想信道均衡 est = rx ./ pathGains; bits = real(est) > 0; ber(idx) = mean(bits ~= data); end figure; semilogy(EbNo, ber, 'o-'); hold on; % 理论曲线 EbN0Lin = 10.^(EbNo/10); theory = 0.5 * (1 - sqrt(EbN0Lin ./ (1 + EbN0Lin))); semilogy(EbNo, theory, 'r-'); legend('仿真', '理论'); xlabel('Eb/N0 (dB)'); ylabel('BER'); grid on;这段代码的关键点有三个。第一，每个信噪比点使用不同的 Seed，避免所有 Eb/N0 点上的衰落序列完全一样；第二，噪声功率用 Es/N0 计算，其中 N0 = Es / EbN0Lin，因为 BPSK 一个符号携带一个比特；第三，均衡采用复数除法，而不是 abs(pathGains)，否则相位信息会丢失，误码率在下行阶段无法收敛到理论值。
运行后如果仿真曲线在 14 dB 以上仍然高于理论值，那么请检查 pathGains 是否在加噪声之前被归一化到单位平均功率。如果 pathGains 的平均功率是 0.7，那么实际 SNR 比预期低 1.5 dB，曲线会右移。这个问题在源码程序里是非常常见的。
5.3 用操作录像配合快速定位曲线不匹配
操作录像的价值在误码率验证阶段体现得最充分。视频里通常会有暂停点，显示星座图或误码率中间结果。你也可以自己加入下面这段快速观察代码：
% scatter_view.m if idx == 5 % 在 EbNo=8 dB 处查看星座图 scatterplot(est(1:1000)); title('Eb/N0=8 dB 均衡后星座'); end星座图能直观看出衰落和噪声叠加后的云团形状。理想情况下，两个星座点云团集中在 ±1 附近，扩散半径正比于噪声标准差。如果云团出现旋转或尺度发散，说明均衡时除的是不正确的复数，建议打印出 pathGains(1:5) 的实部和虚部做检查。
如果你把数据分析法和录像结合，就能快速判断误码率偏差是来自信道生成、噪声功率还是均衡逻辑。这样也不用把整个源码程序重新读一遍，通常十几分钟的定位就能找到问题。我通常建议在脚本里记录每次运行的信道平均功率和噪声方差，作为最基础的调试日志。
6. 将瑞利信道仿真源码程序沉淀为可复用模块
6.1 封装成带参数校验的函数
把源码程序里的信道创建逻辑封装成函数，是提升复用率最直接的技巧。函数签名可以设计成：
function [rxSig, pathGains] = rayleigh_channel_sim(fs, fd, pathDelays, pathGainsDb, inputSignal)函数内部使用 inputParser 做参数检查，比如 fs 必须为正数，fd 小于 fs/2，pathDelays 和 pathGainsDb 长度一致。这样后续在 BER 仿真、误码率扫描或者多天线仿真里都可以直接调用，而不是复制粘贴。封装后的模块也能放进自己的 matlab 工具箱文件夹，方便其他项目引用。
6.2 用脚本批量跑参数扫描并自动保存图片
当你需要输出多个 fd 条件下的误码率曲线时，脚本化扫描是最高效的方案。
fdList = [5 10 20 50]; for i = 1:length(fdList) [rxSig, pathGains] = rayleigh_channel_sim(1000, fdList(i), 0, 0, ones(2000,1)); figure; plot(abs(pathGains)); title(sprintf('fd=%d Hz', fdList(i))); saveas(gcf, sprintf('fd_%dHz.png', fdList(i))); close(gcf); end批量保存图片时一定要在循环里显式 figure 和 close，否则同一窗口的句柄会被反复覆盖，脚本到最后只保留最后一张图。这也是操作录像里容易被忽略的细节，如果你在录制的代码里看到 close(gcf) 被注释掉，会导致生成多张图时比较麻烦。
6.3 用 README 记录版本和操作录像的元数据
使用版本 matlab2021a 不仅要在 README 里写清楚，还应该把操作录像的录制环境一并记录，例如 MATLAB 版本、操作系统、工具箱版本号。这样可以避免你换到 R2022b 时，因为工具箱接口差异而对源码程序产生不信任。操作录像文件中最好包含一个简单的时间戳和版本号，放在文件名里是推荐做法，例如 demo_R2021a_20240115.mp4。
当你需要维护或更新源码程序时，第一步不是读代码，而是打开操作录像，对照 README 里的版本说明确认环境。第二步才是运行 main.m 验证基本路径，最后再由浅入深修改参数。这套流程能让你把别人写的瑞利信道仿真MATLAB源码程序，快速变成自己手里的标准模块，并且在后续仿真任务里直接调用。
本文还有配套的精品资源，点击获取
