1. 项目缘起为什么还要折腾部分响应系统在数字通信的仿真与教学领域MATLAB/Simulink 几乎是绕不开的工具。提到基带传输很多人会立刻想到经典的升余弦滚降滤波器它通过牺牲带宽来换取码间串扰的消除。但今天我想聊一个相对“古典”但思想极其精妙的技术第I类部分响应系统。你可能在教科书里见过它的原理框图知道它能达到 2Baud/Hz 的理论极限频带利用率并且能消除码间串扰。但当你真正打开 Simulink试图从零搭建一个能跑通、能直观看到其“可控码间串扰”特性的仿真模型时往往会遇到一堆教科书上没写的细节问题。比如预编码器到底该怎么实现为什么我的眼图看起来总是不对劲接收端的模2判决在 Simulink 里如何优雅地建模这些坑我当年都一一踩过。网络上关于 Simulink 仿真的资料很多但针对这种特定编码系统的完整、可复现的建模指南却很少。很多人止步于理论公式而无法通过仿真亲眼验证其“在码元速率等于奈奎斯特带宽时依然无码间串扰”这一反直觉的结论。因此我决定把这次完整的仿真建模过程、核心参数设置、关键模块的选型逻辑以及调试中遇到的典型问题和解法系统地梳理出来。无论你是正在完成课程设计的学生还是希望深入理解部分响应技术的工程师这篇基于 Simulink 的实战指南都能让你少走弯路真正把理论“跑”起来。2. 第I类部分响应系统核心原理再透视在动手搭建模型之前我们必须吃透它的工作原理。这不仅是仿真的前提更是后续调试时定位问题的“地图”。2.1 从“消除”到“控制”串扰的思想跃迁传统的奈奎斯特第一准则追求的是在每个采样时刻其他所有码元的响应叠加为零即彻底消除码间串扰。这需要系统的冲激响应在除本码元采样点外的所有其他码元采样点上过零。升余弦滚降特性就是这一思想的产物代价是带宽扩展。部分响应系统则采取了一种截然不同的策略它故意引入一种已知的、固定的码间串扰。更准确地说它让当前码元的响应与相邻的一个或多个码元的响应在采样时刻以确定的方式叠加。对于第I类部分响应其目标频率响应是理想的矩形带宽为W 1/(2T)T为码元周期而冲激响应是sin(x)/x的形式。直接实现这样的系统是非因果且不可实现的。但通过将两个时间上相差一个码元周期T的sin(x)/x响应相加我们得到了第I类部分响应的冲激响应h(t)h(t) \frac{\sin(\pi t / T)}{\pi t / T} \frac{\sin[\pi (t-T) / T]}{\pi (t-T) / T} \frac{\sin(\pi t / T)}{\pi t / T} - \frac{\sin(\pi t / T)}{\pi (t-T) / T}化简后h(t)在t0和tT的采样点上的值均为1在其他整数倍T的采样点上为0。这意味着当前码元的采样值会受到前一个码元采样值的“干扰”且干扰是“加性”的。如果我们发送的原始二进制序列是{a_k}(a_k ∈ {0, 1})那么经过这个系统后在采样时刻kT得到的输出c_k将是c_k a_k a_{k-1}这里c_k的取值可能是{0, 1, 2}。你看串扰被精确地控制为前一个码元的值。这就是“部分响应”中“部分”的含义——串扰没有完全消除但被限制和已知化了。2.2 预编码与模L判决破解串扰的钥匙直接传输c_k会遇到一个大问题在接收端当我们收到c_k后想恢复a_k需要做a_k c_k - a_{k-1}。这要求接收端必须已知a_{k-1}才能解出a_k。一旦某个码元判决错误错误会持续传播下去这叫差错传播。为了解决这个问题天才的“预编码”技术被引入。我们在发送端对原始序列{b_k}(b_k ∈ {0, 1}) 先进行一次变换得到发送序列{a_k}。对于第I类部分响应系统预编码规则为a_k b_k ⊕ a_{k-1} \quad (\text{模2加即异或})注意这里的a_k和上一节的a_k是同一个它是实际经过部分响应系统的输入序列。经过这个预编码后神奇的事情发生了部分响应系统输出的c_k a_k a_{k-1}。如果我们对c_k进行模2运算[c_k]_{mod2} [a_k a_{k-1}]_{mod2} a_k ⊕ a_{k-1}而根据预编码公式a_k b_k ⊕ a_{k-1}我们可以推导出a_k ⊕ a_{k-1} b_k。因此[c_k]_{mod2} b_k也就是说接收端只需要对采样值c_k直接进行模2判决即看c_k是偶数还是奇数就能无差错地恢复出原始信息b_k完全避免了差错传播c_k的取值是{0, 1, 2}对应关系为c_k为偶数(0或2)时判为0为奇数(1)时判为1。这个过程完全不需要知道前一个码元a_{k-1}是什么。注意这里有一个关键点容易混淆。部分响应系统h(t)作用的对象是预编码后的序列{a_k}。c_k是系统在采样点的输出。而我们要传输的原始信息是{b_k}。整个流程是b_k- (预编码) -a_k- (部分响应系统) -c_k- (信道) -\hat{c}_k- (模2判决) -\hat{b_k}。2.3 系统频带利用率分析第I类部分响应系统的频率响应H(f)是两个时移T的理想低通滤波器频率响应之和H(f) T \cdot [e^{-j\pi f T} e^{j\pi f T}] \cdot \text{rect}(f/(2W)) 2T \cos(\pi f T) \cdot \text{rect}(f/(2W))其中W 1/(2T)rect(f/(2W))在|f| W时为1否则为0。它的幅度谱是一个升余弦的形状在f W处为0。关键点在于它的带宽严格等于W 1/(2T) Hz。而码元速率是R_s 1/T Baud。因此其频带利用率为\eta R_s / W (1/T) / (1/(2T)) 2 \quad (\text{Baud/Hz})这达到了在无码间串扰这里是“无未知码间串扰”条件下二进制传输的极限频带利用率奈奎斯特极限。这是它最核心的优势。3. Simulink 仿真模型搭建全流程理论清晰后我们进入实战环节。我将按照信号流的方向逐一讲解每个模块的选型、参数配置和背后的考量。3.1 信源与预编码模块实现首先我们需要一个随机的二进制信源。在 Simulink 中有多个选择Bernoulli Binary Generator生成0/1序列可设置概率和初始种子。这是最直接的选择。Random Integer Generator设置M-ary number为2输出0和1。PN Sequence Generator生成伪随机序列相关性好常用于测试。这里我选择Bernoulli Binary Generator因为它最贴近“等概二进制信源”的概念。将其重命名为Binary Source采样时间Sample time设置为T比如1e-6(代表1us一个码元即1M Baud)。Probability of a zero设为0.5Initial seed可以任意设置比如12345。接下来是预编码器。根据公式a_k b_k ⊕ a_{k-1}这是一个带反馈的异或操作。在 Simulink 中实现反馈环路需要特别注意避免代数环。标准做法是引入单位延迟。搭建结构将信源输出接入一个XOR逻辑运算模块的一个输入端。XOR模块的输出即为a_k。将a_k连接到一个Unit Delay模块的输入。Unit Delay模块的输出即为a_{k-1}将其反馈到XOR模块的另一个输入端。Unit Delay模块的Initial condition需要设置。通常我们假设初始时刻a_{-1} 0。因此这里设为0。关键参数Unit Delay的Sample time必须设置为-1(继承) 或明确的T确保其延迟周期与信源码元周期一致。整个预编码子系统应封装成一个Subsystem便于管理。将其命名为Precoder。实操心得预编码模块的仿真步长必须与信源周期同步。如果使用变步长求解器Simulink 可能会在非码元周期时刻尝试计算导致反馈环路出现非整数延迟问题。一个稳妥的办法是在模型配置参数Solver中将Type设置为Fixed-stepFixed-step size设置为信源周期T或其整数分之一如T/10以提高波形分辨率。这能保证仿真在离散的时间点上精确推进。3.2 第I类部分响应滤波器设计与实现这是模型的核心。我们需要实现冲激响应为h(t)的系统。根据其定义它是两个理想低通滤波器冲激响应的叠加。但在实际仿真中我们通常直接在时域或频域逼近这个特性。方法一时域卷积法直接但低效理论上我们可以用Discrete Impulse Response滤波器直接输入h(t)的采样值作为系数。但h(t)是无限长的需要截断。计算系数也比较麻烦且这种方法对仿真步长非常敏感不推荐作为首选。方法二频域成型法推荐更贴近原理的方法是先产生冲激序列由a_k驱动然后通过一个具有第I类部分响应目标频率特性的滤波器。我们可以利用两个理想低通滤波器的并联来实现。生成发送波形预编码输出的a_k是离散的0/1序列。我们需要将其转换为连续时间波形。使用Zero-Order Hold模块采样时间设置为T。这样每个码元周期内电压值保持恒定。构建滤波器第I类部分响应H(f) 2T cos(πfT) * rect(f/(2W)),W1/(2T)。我们可以将其理解为一个截止频率为W的理想低通滤波器其输出与其延迟T后的输出相加。Simulink实现 a. 将 ZOH 输出的信号分为两路。 b. 一路直接进入一个Analog Filter Design或Discrete Filter模块配置为低通滤波器。关键点为了逼近“理想”特性我们通常使用阶数非常高的 FIR 滤波器或者使用Windowed FIR设计方法截止频率设为W。例如使用fir1(N, Wn)函数设计一个阶数N100甚至更高的 FIR 滤波器。在 Simulink 中可以使用Discrete FIR Filter模块并将系数Numerator设置为设计好的b。 c. 另一路经过一个Transport Delay模块延迟时间设置为T然后送入一个完全相同的低通滤波器。 d. 将两路滤波器的输出用Add模块相加。简化实现更常用实际上由于我们追求的是系统的净效果且仿真通常在离散时间进行我们可以直接设计一个具有目标频率响应H(f)的成型滤波器。在 MATLAB 中我们可以根据H(f)的公式直接计算其冲激响应的采样值作为 FIR 滤波器系数。% 示例设计一个用于仿真的第I类部分响应 FIR 滤波器 T 1e-6; % 码元周期 Rs 1/T; % 码元速率 Fs 10 * Rs; % 仿真采样率通常是码元速率的整数倍如10倍 W Rs / 2; % 奈奎斯特带宽 N_taps 101; % 滤波器阶数取奇数以便有中心抽头 t (-(N_taps-1)/2 : (N_taps-1)/2) / Fs; % 时间向量 h_ideal (sin(pi*t/T)./(pi*t/T)) (sin(pi*(t-T)/T)./(pi*(t-T)/T)); h_ideal(isnan(h_ideal)) 1 (t(isnan(h_ideal))0); % 处理 t0 和 tT 时的 NaN 值 % t0时第一项为1第二项为 -1/π * lim_{x-0} sin(x)/x ? 实际上需要计算极限。 % 更严谨的系数计算 n -(N_taps-1)/2 : (N_taps-1)/2; h zeros(size(n)); for i 1:length(n) if n(i) 0 h(i) 1; elseif n(i) 1 % 注意n 是索引偏移对应时间 n*Ts。需要映射到时间。 % 我们需要找到对应时间 T 的索引。如果 Fs 10*Rs, 则 T 10 个采样点。 idx_T round(T * Fs); % 延迟 T 对应的采样点数 if abs(n(i)) idx_T h(i) 1; end else h(i) (sin(pi*n(i)/idx_T)/(pi*n(i)/idx_T)) (sin(pi*(n(i)-idx_T)/idx_T)/(pi*(n(i)-idx_T)/idx_T)); end end h h / sum(abs(h)); % 可选归一化然后将计算好的h系数填入Discrete FIR Filter模块的Coefficients参数中。这种方法更直接性能也足够用于仿真。我将这个滤波器子系统封装起来命名为Partial Response Filter (Type I)。3.3 信道与噪声添加为了模拟真实环境我们添加一个简单的加性高斯白噪声信道。使用AWGN Channel模块。Mode设置为Signal to noise ratio (SNR)或Eb/No。对于基带仿真用SNR更直观。SNR (dB)这是一个重要的调试变量。可以从较高的值如30dB开始观察无噪声下的理想情况然后逐步降低如10dB5dB来测试系统抗噪性能。Input signal power (watts)需要根据发送信号的功率来设置。一个简单的方法是先不加噪声运行一次仿真用Display模块或工作空间变量查看发送滤波器输出信号的方差var(signal)将其值填入此处。或者设置为1然后通过调整发送信号幅度来归一化信号功率。3.4 接收端匹配滤波、采样与判决接收端的目标是从含噪波形r(t)中恢复出发送的原始信息b_k。匹配滤波器在最大信噪比准则下最佳接收滤波器是发送滤波器的匹配滤波器。由于我们的发送滤波器h(t)是实系数且对称第I类部分响应冲激响应是对称的所以匹配滤波器就是它本身。因此我们可以让接收信号通过一个与发送端完全相同的Partial Response Filter (Type I)。注意在实际系统中发送和接收滤波器通常各为sqrt(H(f))以实现匹配并满足总体响应但这里为简化我们使用相同的滤波器这会影响绝对性能但不影响原理验证。采样匹配滤波器输出后我们需要在最佳时刻进行采样。使用Zero-Order Hold模块的逆过程——采样。Simulink 中没有直接的“采样器”通常有两种方法使用Downsample模块将采样率降低到码元速率。例如如果仿真采样率Fs 10*Rs那么Downsample的Downsample factor设为10Sample offset需要调整以找到最佳采样相位通常是匹配滤波器输出眼图张开最大的时刻。使用Rate Transition模块设置不同的采样率。更精确的方法使用MATLAB Function模块或Interpreted MATLAB Fcn编写简单的采样脚本或者利用Pulse Generator触发一个Sample and Hold模块。对于教学仿真使用Downsample并仔细调整Sample offset是最简单的。模2判决采样后得到序列\hat{c}_k。根据原理我们需要进行模2判决判断\hat{c}_k最接近{0, 1, 2}中的哪个值然后输出对应的二进制位。实现使用一个MATLAB Function模块内部代码非常简单function y mod2_decision(u) % u: 输入采样值 % y: 判决输出的二进制比特 (0/1) levels [0, 1, 2]; % 可能的发送电平 [~, idx] min(abs(u - levels)); % 找到最接近的电平索引 y mod(levels(idx), 2); % 对该电平模2得到原始信息比特 end或者可以用Compare To Constant模块设置阈值。由于理想情况下c_k取值为0,1,2阈值可以设在0.5和1.5。小于0.5判为0大于1.5判为2两者之间判为1。然后再对结果0,1,2进行模2运算即判断奇偶性。这种方法更直观地体现了“门限判决”的思想。误码率计算将判决输出的\hat{b_k}与原始信源b_k进行比较。注意由于系统存在延迟滤波器延迟、匹配滤波器延迟等原始信源序列需要相应延迟后才能对齐比较。使用Error Rate Calculation模块设置好Receive delay通常为滤波器群延迟的2倍左右Computation delay设为0Computation mode设为Entire frame。将输出的误码率连接至Display模块观察。3.5 完整模型架构与信号流图将上述所有模块连接起来完整的 Simulink 模型顶层架构应如下所示[Bernoulli Binary] -- [Precoder] -- [ZOH] -- [Partial Response Filter] -- [AWGN Channel] | v [原始信源延迟对齐] -- [Delay] -- [Bernoulli Binary] [匹配滤波器] -- [AWGN Channel输出] | | | v v v [Error Rate Calc] -- (比较) -- [判决器] -- [采样器] -- [匹配滤波器输出]务必使用Scope模块在关键点观察波形预编码输出a_k、发送滤波器输出波形、匹配滤波器输出波形观察眼图、采样值\hat{c}_k、判决输出\hat{b_k}。4. 仿真参数配置与关键调试技巧模型搭好了但直接运行很可能得不到预期结果。以下是确保仿真成功的核心参数设置和调试步骤。4.1 求解器与步长设置这是 Simulink 动态系统仿真的基石设置不当会导致结果错误或仿真失败。进入Model Configuration Parameters(CtrlE)。Solver部分Type:强烈建议选择Fixed-step。对于这种数字通信系统仿真固定步长能保证采样时刻的精确性避免变步长带来的相位漂移和同步问题。Solver:discrete (no continuous states)。因为我们的模型本质是离散时间系统虽有连续波形但由离散信号驱动。Fixed-step size: 这是最重要的参数。它决定了仿真的时间分辨率。规则它必须是码元周期T的整数分之一。通常设置为T/10、T/20或T/50。更高的分辨率能画出更光滑的波形和更清晰的眼图但会降低仿真速度。例如T1e-6(1us)可以设Fixed-step size 1e-7(100ns)。确保所有离散模块如Unit Delay,Discrete FIR Filter的采样时间与此固定步长一致或为其整数倍。Data Import/Export部分可以取消勾选Save simulation output to workspace除非你需要详细数据做后期分析这能节省内存和仿真时间。4.2 各模块采样时间同步Simulink 中混合了连续和离散模块采样时间冲突是常见错误。Bernoulli Binary Generator、Unit Delay、Discrete FIR Filter这些是离散模块其采样时间应设置为-1(继承) 或明确的T。在固定步长仿真下继承的通常是固定步长值。但像Unit Delay在预编码中我们明确希望它延迟一个码元周期因此最好将其采样时间设置为T。如果固定步长是T/10那么Unit Delay的采样时间设为TSimulink 会自动处理多速率问题。Zero-Order Hold它的采样时间必须设置为码元周期T以确保每个码元保持T秒。AWGN Channel通常设置为-1(继承)它会根据输入信号的采样率来添加噪声。Downsample其Sample time设为-1Downsample factor设为Fs / Rs即仿真采样率与码元速率的比值。调试技巧在模型编译后按CtrlD可以点击菜单栏的Display - Sample Time - Colors。不同采样时间的模块和信号线会以不同颜色显示。检查是否有意外的红色固定步长不匹配或品红色连续与离散混合问题。确保信号流路径上的采样时间转换是合理且受控的。4.3 眼图观测与采样相位调整眼图是评估基带传输系统性能最直观的工具。在匹配滤波器的输出端连接一个Scope。设置Scope参数通常将Time span设置为2*T或3*T这样能显示2到3个码元的波形。在Scope的Configuration Properties中找到Display选项卡勾选Eye Diagram。将Eye diagram period设置为码元周期T。运行仿真。在无噪声或高 SNR 下你应该能看到一个清晰的眼图。第I类部分响应系统的眼图特征它应该有三个明显的电平开口对应c_k的0,1,2而不是传统二进制系统的两个。这是其核心特征。调整采样相位眼图的张开度最大处是最佳采样时刻。通过调整Downsample模块的Sample offset参数从0到Downsample factor-1之间尝试观察误码率的变化找到误码率最低的偏移量即为最佳采样相位。4.4 误码率性能测试逐步降低AWGN Channel模块中的SNR (dB)记录对应的误码率。可以将仿真结果与理论误码率曲线进行比较。对于第I类部分响应系统在采用匹配滤波和最佳采样后其误码性能比理想二进制系统如双极性不归零码约有2-3 dB的损失。这是因为它的三个电平判决距离更近抗噪声能力减弱。这是用性能换取带宽效率的典型体现。你可以通过绘制BER vs. SNR曲线来验证这一趋势。5. 常见问题排查与解决方案即使按照上述步骤你可能还是会遇到一些棘手的问题。以下是我在多次仿真中总结的“坑”和填坑方法。5.1 眼图模糊或不稳定现象眼图线条很粗重叠不整齐或者根本看不出三个电平。可能原因及解决仿真时间或瞬态效应滤波器尤其是高阶 FIR 滤波器在初始阶段有瞬态响应这不是稳态信号。解决在Error Rate Calculation中设置Computation delay忽略开头一段时间的输出。或者让仿真运行足够长的时间如10000个码元以上眼图会自动稳定。滤波器系数不正确这是最常见的原因。用于仿真的 FIR 滤波器系数未能准确逼近理想的第I类部分响应。解决将你设计的滤波器系数h绘制出来 (stem(h))看其波形是否近似于理论冲激响应h(t)在t0和tT处有峰值。检查滤波器的频率响应 (freqz(h,1)) 是否在f1/(2T)处为零。采样率不足固定步长Fixed-step size设置得太大导致波形分辨率太低无法刻画细节。解决减小固定步长例如从T/10提高到T/50。噪声过大SNR 设置太低。解决先将 SNR 设为很大如50dB观察无噪声下的眼图是否清晰。如果清晰再逐步降低 SNR。5.2 误码率始终为0.5或0.5左右现象无论 SNR 高低误码率都在0.5附近这意味着判决完全随机系统失效。可能原因及解决预编码与判决逻辑不匹配这是最可能的原因。检查预编码规则是否为a_k b_k ⊕ a_{k-1}。检查判决逻辑是否为对c_k进行模2判决奇偶判决。一个快速验证方法在无噪声下用Scope或To Workspace模块导出b_k、a_k、c_k、\hat{c}_k、\hat{b_k}几个序列的前20个值手动验证关系是否成立。采样时刻完全错误采样点正好落在眼图的闭合点即码间干扰最大的时刻。解决通过眼图观察并系统性地扫描Downsample的Sample offset寻找误码率最低的点。信号极性反转检查整个链路中是否有意外的Gain模块设置为-1或者加法器、减法器用错。确保发送和接收滤波器系数一致。5.3 仿真速度极慢现象模型很简单但仿真进度条蠕动缓慢。可能原因及解决使用了高阶连续滤波器如果采用Analog Filter Design等连续滤波器模块Simulink 需要使用小步长来求解微分方程极慢。解决全部改用离散滤波器 (Discrete FIR Filter)。仿真时间过长或步长过小码元数太多或固定步长太小。解决对于 BER 测试可以先跑较短序列如1e4个码元看趋势。步长在能分辨眼图的前提下不要设置得过小T/20通常足够。Scope 数据记录过多Scope默认记录所有仿真数据到工作区。解决在Scope配置中取消Logging选项卡下的Log data to workspace。或者使用Terminator模块终结不需要的信号线。5.4 代数环错误现象仿真无法启动报错提示代数环。可能原因及解决预编码反馈环路XOR和Unit Delay的环路是正常的不会产生代数环因为Unit Delay打破了代数环。如果错误检查Unit Delay的Sample time是否设置为-1或一个具体值而不是0。其他直接馈通环路检查模型中是否有其他没有延迟的反馈路径例如两个Gain模块直接首尾相连。确保反馈路径中至少有一个离散模块如Unit Delay,Memory,Discrete Filter。通过以上系统的建模、参数配置和调试你应该能够成功搭建并运行一个第I类部分响应系统的 Simulink 仿真模型并直观地验证其高频带利用率和独特的“三电平”眼图特性。这个过程中对细节的把握比如采样时间同步、滤波器系数设计、眼图观测和采样相位调整正是仿真工作从“知道”到“做到”的关键跨越。希望这份详细的指南能帮助你不仅完成一次仿真作业更能深刻理解部分响应系统的工作机制与仿真建模的精髓。