简介：面向通信原理学习与 Simulink 仿真实训场景，这一压缩包提供了一套完整的 2ASK（二进制幅度键控）调制解调仿真方案。RAR 包内共有 2 个文件，包含 1 个可运行 M 脚本和 1 个 Simulink 模型，整体体积仅 22KB，轻量易用，能够直接打开并运行查看仿真结果。模型以模块化方式串联信号源、2ASK 调制器、加性高斯白噪声信道、解调器与误码率统计环节，可直观观察原始二进制序列、已调波形、噪声叠加后波形以及最终解调恢复的数据。目前已有 1624 人学习使用，适合通信工程、电子信息类学生进行课程设计、毕业设计或课外实验。使用者可自行调整载波频率、信噪比等参数，对比不同条件下的误码率变化，从而深入理解 2ASK 对噪声的敏感程度及幅度键控的基本原理；同时，这套小体积模型也可作为继续学习 PSK、QAM 等调制技术的可复用仿真基础。
1. 2ASK 仿真在 Simulink 里卡住的从来不是调制本身
2ASK 的调制原理一句话就能说完：二进制 1 时输出载波，二进制 0 时输出 0，也就是“通断键控”。但把这句话搬进 Simulink，很多人第一个晚上是这样度过的：Scope 里看不到正弦包络，报错说采样时间冲突，BER 不是恒为 0.5 就是低到不真实。问题几乎都出在模型层面，而不是通信原理层面：采样率设多少、载波频率填 rad/s 还是 Hz、低通滤波器群延迟怎么补偿、AWGN 模块里的 SNR 和 EbN0 差多少，这四个点任何一个没对齐，整条链路的误码率都是假的。
这篇文章就是围绕这四个点展开的。先讲清楚 2ASK 在频域和波形上的两个基本事实，再给出一条通带模型下可复现的调制解调链路，然后用误码率曲线和排查表把参数怎么调说透。内容适合三类人：通信方向做课设或毕设的学生、要把 2ASK 当作跳板去理解 QPSK 和 OFDM 的从业者，以及想在 Simulink 里做算法链路验证的工程师。
2. 2ASK 调制原理与 Simulink 建模的两条路线
2.1 从数学表达式看 2ASK：为什么必须用单极性 NRZ
2ASK 的时域表达式只有一行：
s(t) = A · m(t) · cos(2π·fc·t)其中m(t)是单极性不归零码，取值 0 或 1。这里最容易犯的理解错误是：把m(t)看成 ±1 的双极性信号。如果 m(t) 取 ±1，乘上载波之后相位会在 0 和 π 之间跳变，那已经不是 ASK，而是 BPSK。单极性 NRZ 的 0 电平把载波“关断”，这正是 OOK 名字的由来，也决定了它和 PSK 在误码性能上的本质差异。
从频域看，单极性 NRZ 的频谱包含一个直流量和以 Rb 为第一个零点的连续谱，调制后相当于把这段谱搬到 ±fc 两侧，双边主瓣总宽约2·Rb。所以 2ASK 的射频带宽通常取两倍比特率，频带利用率不高，换来的是接收端可以用包络检波这种极其简单的结构。频谱宽度直接决定了后面通带仿真里的滤波器阶数和采样率，这是后续所有参数设置的源头。
2.2 通带模型与基带等效：波形可看性和误码率可算性
在 Simulink 里做 2ASK 仿真有两条路线：通带模型和基带等效模型。通带模型保留真实的余弦载波，Scope 里能直接看到正弦包络，信息直观，适合做波形演示、报告截图和硬件接口验证。基带等效模型把载波频率归零，信号变成复包络，仿真步长只受基带带宽约束，计算量小一个数量级以上，适合批量跑误码率曲线。
| 对比维度 | 通带模型 | 基带等效模型 |
|---|---|---|
| 信号表达 | 实信号，载波可见 | 复包络，载波被搬移到 0 频 |
| 采样率要求 | Fs 必须覆盖 fc 和 2·Rb | Fs 只需覆盖基带带宽，通常 4 到 8 倍 Rb |
| 常见用途 | 波形验证、课设演示、射频链路联调 | EbN0 扫描、BER 统计、算法对比 |
| 典型模块 | Sine Wave + Product | Comm Toolbox 的 ASK 基带调制模块 |
两条路线并不冲突。我一般建议先用通带模型把波形和链路走通，确认滤波器、判决和对齐没问题，再决定是否切换到基带等效模型去做大批量仿真。本文第 3 章搭的是通带模型，第 5 章的批量扫描也直接在同一个通带模型上跑，因为通带模型更容易定位“BER 为什么不对”这类问题。
2.3 基础模块与 Comm Toolbox 模块的选型建议
Simulink 里实现 2ASK 有两种模块选择。一种是纯基础模块：Bernoulli Binary Generator、Sine Wave、Product、Discrete FIR Filter，所有参数都摊在明面上；另一种是 Communications Toolbox 提供的 ASK 基带调制解调模块，一条线拉完，内部把采样率和星座映射都封装好了。
封装模块省事，但会把“每个比特重复 N 个样本”这类关键操作藏起来，一旦 BER 曲线对不上理论值，很难判断是噪声设置问题还是信号帧格式问题。基础模块虽然连线多，但每一步都可以用 Scope 和 To Workspace 单独检查。对于想真正搞懂 2ASK 仿真过程的人，先搭一遍基础模块链路是值得的。下面这段代码对应基础模块里“单极性 NRZ”信号的生成逻辑，理解它之后，Simulink 里 Repeat 模块的作用也就清楚了：
Rb = 1000; % 比特率 1 kbps N = 100; % 每比特采样点数 Fs = Rb * N; % 模型采样率 data = randi([0 1], 30, 1); % 30 个随机比特 x = repelem(data, N); % 每个比特重复 N 点，形成单极性 NRZ t = (0:length(x)-1) / Fs; plot(t * 1000, x); % 横轴单位 ms，查看前几个比特repelem(data, N)把每个比特持续 N 个采样点，形成矩形脉冲序列，这和 Simulink 里 Bernoulli Binary Generator 后接 Repeat 模块做的事情完全一致。理解了这段代码，后面调试采样率冲突和波形失真时就有对照基准。
3. 在 Simulink 里搭出一条 2ASK 通带收发链路
3.1 发射端搭建：数据源、NRZ 成形与载波相乘
发射端参数先定下来：比特率Rb = 1000，每比特采样数N = 100，模型采样率Fs = Rb * N = 100 kHz，载波频率fc = 10 kHz。这样每个载波周期有 10 个采样点，Scope 里能看出正弦形状，又不至于让仿真时间过长。
| 模块 | 所在库 | 关键参数设置 |
|---|---|---|
| Bernoulli Binary Generator | Communications Toolbox / Comm Sources | Probability of a zero = 0.5；Sample time = 1/Rb；Output data type = double |
| Repeat | DSP System Toolbox / Signal Operations | Repetition factor = N |
| Sine Wave | Simulink / Sources | Sine type = Time based；Frequency = 2pifc；Sample time = 1/Fs |
| Product | Simulink / Math Operations | 两个输入，默认乘法 |
| Scope | Simulink / Sinks | Number of input ports = 2，同时观察 NRZ 和调制输出 |
连线顺序是 Bernoulli Binary Generator 输出每个比特一个样本，Repeat 把它扩展成每比特 N 个样本，再和离散正弦波相乘。这里有一个高频错误：Sine Wave 模块的频率参数单位是 rad/s，必须填2*pi*fc，直接填 fc 会让载波周期缩短 2π 倍，Scope 里看起来就像一团乱码。
如果机器上没有 Communications Toolbox，Bernoulli 模块不可用，可以用 MATLAB 脚本生成时间序列后从 Workspace 导入：
rng(7); numBits = 5000; data = randi([0 1], numBits, 1); x = repelem(data, N); x_ts = timeseries(x, (0:length(x)-1)' / Fs);在 Simulink 里放一个 From Workspace 模块，参数填x_ts，输出就是单极性 NRZ 波形。这条替代路径用到的都是 MATLAB 基础功能，和后续接收端的连接方式完全一致。
3.2 接收端搭建：相干解调中的滤波器、阈值与延迟补偿
接收端采用相干解调，整个信号链是：AWGN Channel 加噪 → 乘以本地载波 → 低通滤波 → 阈值判决 → 与参考序列对齐统计误码。本地载波在仿真中直接复用发射端的 Sine Wave，相位设为 0，这相当于假设载波同步已经完成。
相乘之后，信号变成两项：m(t)/2和中心频率在2*fc的高频分量。低通滤波器要把第二项滤掉，保留第一项。滤波器用 FIR 而不是 IIR，原因是 FIR 线性相位，群延迟恒定为(Ntap-1)/2个样本，方便用 Integer Delay 精确补偿；IIR 的群延迟在通带内有波动，判决点不好找。
Ntap = 64; % 滤波器阶数 63 h = fir1(Ntap - 1, 1/N); % 归一化截止 0.5*Rb，约 1/N参数说明：1/N是归一化截止频率。真实截止频率为0.5 * Rb，除以半采样率Fs/2后得到(0.5*Rb) / (Rb*N/2) = 1/N。取 0.5 倍 Rb 是因为单极性 NRZ 的频谱主瓣从 0 延伸到 Rb，切掉后半段能滤掉大量噪声，且不会对判决造成明显码间干扰。Ntap 取 64，群延迟为(64-1)/2 = 31.5个样本，Integer Delay 补偿 31 或 32 都可以。
滤波器系数h在 MATLAB 工作区生成后，写到 Discrete FIR Filter 模块的 Numerator 参数里。低通输出大约是0.5 * m(t)，判决阈值取 0.25，用 Relational Operator 模块做大于等于比较即可，输出为 0/1 的判决序列。这里不要用 Relay 模块，Relay 默认带滞回特性，参数设置不当会在阈值附近产生毛刺。
3.3 误码率统计：Error Rate Calculation 的接法与对齐
误码率统计是整个仿真里最容易被“蒙混”过去的部分。把判决输出和发射端 NRZ 信号分别拉进 Error Rate Calculation 模块，如果不做延迟补偿，滤波器造成的 31 个样本延迟会让比较的两个序列错位。最直接的后果是 BER 恒定为 0.5，因为两个随机序列在时间错位时统计出来就是一半一致一半不一致。
对齐方法是：在发射端 Repeat 模块的输出后面接一个 Integer Delay 模块，延迟值设为 31，用它作为参考序列进入 Error Rate Calculation。接收端的判决序列本身已经包含了滤波器延迟，两端同时处于同一时间基准上。这样不用做符号定时恢复，因为每个比特被重复了 N 个样本，只要延迟误差远小于 N，逐样本比较的结果和逐比特比较几乎一致。
Error Rate Calculation 模块参数里有一项 Output data，选 Port 后模块输出一个向量：[BER, 错误比特数, 比较比特数]。把这个输出接一个 Display 模块，模型就能实时显示当前信噪比下的误码率。这个模块允许前几个数据被丢弃，但在这里没有意义，因为我们的序列已经在模块外手动对齐了。
4. 仿真参数怎么调：误码率曲线、信噪比换算与三个必查项
4.1 相干与非相干 2ASK 的理论误码率与仿真曲线怎么对齐
相干 2ASK 的理论误码率公式和 BPSK 形式上一致：
Pe = Q(sqrt(Eb/N0))非相干包络检波的近似公式为：
Pe ≈ 0.5 * exp(-Eb / (2*N0))注意这里用的是平均比特能量。OOK 的两个符号不是等能量信号，发送 1 时能量为A²·Tb/2，发送 0 时能量为 0，平均 Eb 是A²·Tb/4。代入判决距离后，相干检测恰好得到Q(sqrt(Eb/N0))而不是Q(sqrt(2*Eb/N0))。很多人的仿真曲线比理论曲线好 3 dB，就是因为在公式里多算了一个 2。
在 MATLAB 里把理论曲线画出来，用于和仿真结果叠图：
EbN0_dB = 0:2:14; EbN0_lin = 10.^(EbN0_dB/10); Pe_coherent = qfunc(sqrt(EbN0_lin)); % 相干 OOK Pe_nocoherent = 0.5 * exp(-EbN0_lin/2); % 非相干近似 semilogy(EbN0_dB, Pe_coherent, '-', EbN0_dB, Pe_nocoherent, '--');观察这两条曲线的差距可以发现，非相干检测在高信噪比下大约比相干检测差 1 dB 左右。如果你的 Simulink 模型里包络检波路线跑出来的 BER 反而比相干曲线还好，那一定是噪声设置或者判决阈值出了问题，而不是包络检波更优。
4.2 AWGN 模块里的 SNR 与 EbN0 之间差了一个每比特采样数
这是整个 2ASK 仿真里最隐蔽的一个参数陷阱。AWGN Channel 模块通常以 SNR 为参数，而理论误码率公式用的是 EbN0。两者在离散系统中的关系为：
SNR_dB = EbN0_dB - 10 * log10(N)其中 N 是每比特采样数，也就是模型采样率 Fs 除以比特率 Rb。推导过程很直接：OOK 平均信号功率约为A²/4，噪声功率为N0 * Fs，所以 SNR 等于 EbN0 除以每比特采样数 N。在仿真脚本里设置 AWGN 模块时，正确的做法是：
SNR_dB = EbN0_dB - 10 * log10(N); set_param([model '/AWGN Channel'], 'SNR', num2str(SNR_dB));如果直接拿 EbN0 填进 SNR 参数，仿真曲线会比理论曲线右移10*log10(N)dB。N 取 100 时这个偏差是 20 dB，整条 BER 曲线形状对、位置完全不对。反过来，如果看到曲线左移，就要检查是不是在计算 EbN0 时把平均功率当成了峰值功率。
4.3 三个最常见的“BER 假死”原因和排查顺序
| 现象 | 可能原因 | 排查与处理 |
|---|---|---|
| BER 恒定 0.5 | 参考序列与判决序列未对齐，或判决阈值被噪声淹没 | 用 Integer Delay 补偿滤波器群延迟；把 AWGN 暂时设为 0 确认链路本身无误 |
| 波形没有正弦包络 | 载波频率填错单位，或每比特采样数太少 | 检查 Sine Wave 的 Frequency 是否为2*pi*fc；把 N 提高到 100 以上 |
| 仿真曲线整体右移 10~20 dB | SNR 与 EbN0 换算错误 | 核对 AWGN 模块里填的是不是EbN0_dB - 10*log10(N) |
排查顺序有讲究。先把 AWGN Channel 的 SNR 设成一个很大的负数或直接断开启用，看 BER 是否为 0。这一步不过关，后面所有噪声相关分析都没有意义。然后固定一个 SNR，把发射端 NRZ、接收端判决前信号、参考序列三个信号同时画到一个 Scope 里，肉眼看延迟和阈值位置。最后才去核对误码率公式。
判断阈值也可以微调。理论最优阈值在A/4附近，但滤波器通带增益不完全是 1 时，实际信号幅度会偏离 0.5。一个更稳的做法是先用 Scope 量出判决前信号在发送 1 时的稳态电平，取该电平的一半作为阈值。
5. 用 MATLAB 脚本批量扫描 EbN0，省掉手动点模块
5.1 批量仿真脚本：set_param 加 sim 循环
手动在 AWGN 模块里改 SNR、点一次运行、抄一次 BER，两三个点还能忍，跑到 8 个 EbN0 点就变成了纯体力活。批量扫描的标准做法是用set_param在每次仿真前改写模块参数，再用sim驱动模型运行。
EbN0_dB = 0:2:14; N = 100; % 与模型里的 Repeat 因子一致 model = 'ask2_sim'; load_system(model); ber_sim = zeros(1, numel(EbN0_dB)); for k = 1:numel(EbN0_dB) SNR_dB = EbN0_dB(k) - 10 * log10(N); set_param([model '/AWGN Channel'], 'SNR', num2str(SNR_dB)); simOut = sim(model, 'StopTime', '5'); % 5 秒仿真，对应 5000 比特 ber = simOut.get('ber'); % 从 To Workspace 读取 ber_sim(k) = ber(1); % 取 BER，即输出向量第一个元素 end theory = qfunc(sqrt(10.^(EbN0_dB/10))); semilogy(EbN0_dB, max(ber_sim, 1e-6), 'o-', EbN0_dB, theory, '-'); legend('Simulink 仿真', '相干 OOK 理论');脚本逻辑说明：每次循环先由 EbN0 换算出 SNR，写入 AWGN Channel 模块，再调用sim运行整个模型。StopTime设 5 秒对应 5000 个比特，这个数据量在 BER 低于1e-3时勉强够看，曲线尾部会有些抖动，属于正常现象。max(ber_sim, 1e-6)是为了半对数坐标上不出现 0 值，避免画图时丢点。
值得注意的一个细节是：sim每次调用都会对模型做完整的初始化，因此 Error Rate Calculation 模块的内部计数器每次都会从零开始，不需要额外接 Reset 信号。要让每次仿真使用相同的数据序列，可以在脚本里固定随机种子，或者用第 3.1 节的timeseries方式把预生成数据从 Workspace 输入，这样不同 SNR 点之间只差噪声，信号序列完全一致，曲线抖动更小。
如果安装了 Parallel Computing Toolbox，可以把for换成parsim，配合Simulink.SimulationInput数组并行跑不同 SNR 点。但并行模式下每次仿真的模型状态是隔离的，需要在循环外先构造好全部 SimulationInput 对象再统一启动。先把串行脚本跑通，再改并行，定位问题会省力得多。
本文还有配套的精品资源，点击获取
