一、理论基础

1.码间干扰的产生机理

数字基带传输系统的总传输特性由发送滤波器\(G_T(\omega)\)、信道\(C(\omega)\)与接收滤波器\(G_R(\omega)\)级联构成，即 \(H(\omega) = G_T(\omega) \cdot C(\omega) \cdot G_R(\omega)\)。

根据时频对偶性，频带受限的系统，其时域冲激响应必然无限延展。单个码元经过限带系统后会产生拖尾，延伸至后续码元周期；多码元波形叠加时，相邻码元的拖尾分量会在当前码元的抽样判决时刻叠加，造成判决幅值偏差，这一现象即为码间干扰。

2.无码间干扰的判定条件

2.1 时域条件

假设抽样定时无偏差，无 ISI 的时域充要条件为：系统单位冲激响应\(h(t)\)仅在本码元抽样时刻有确定幅值，在其余所有码元抽样时刻取值恒为 0：

\(h(mT_B)= \begin{cases} 1, & m=0 \\ 0, & m为非零整数 \end{cases}\)

此时相邻码元的拖尾不会在抽样点产生串扰。

2.2 频域条件：奈奎斯特第一准则

时域条件经傅里叶变换映射到频域，即为奈奎斯特第一准则，表达式为：

\(\sum_{i=-\infty}^{\infty} H\left( \omega + \frac{2\pi i}{T_B} \right) = T_B \quad , \quad |\omega| \leq \frac{\pi}{T_B}\)

其物理意义可概括为分段、平移、叠加、常数：将传输特性\(H(\omega)\)以\(2\pi/T_B\)为间隔分段，全部平移至\([-\pi/T_B,\ \pi/T_B]\)主区间叠加，若结果为幅值恒定的理想低通特性，则码元速率\(R_B=1/T_B\)时系统无码间干扰。

3.典型无码间干扰基带传输特性

3.1 理想低通特性

理想低通是带宽最窄的无 ISI 特性，对应基带传输的性能理论极限：

• 传递函数： \(H(\omega)= \begin{cases} T_B, & |\omega| \leq \frac{\pi}{T_B} \\ 0, & |\omega| > \frac{\pi}{T_B} \end{cases}\)
• 冲激响应：\(h(t)=\text{Sa}\left( \frac{\pi t}{T_B} \right)\)
• 核心特性：频带利用率达 2 Baud/Hz 的理论上限；但时域拖尾按\(1/t\)衰减，对抽样定时误差敏感，且陡峭的过渡带物理上无法精确实现，仅作为理论基准。

3.2 升余弦滚降特性

工程主流方案，通过滚降因子\(\alpha\)（\(0\leq\alpha\leq1\)）平衡带宽与时域性能，\(\alpha=1\)为全滚晋升余弦：

• 传递函数： \(H(\omega)= \begin{cases} \frac{T_B}{2}\left( 1+\cos\frac{\omega T_B}{2} \right), & |\omega| \leq \frac{2\pi}{T_B} \\ 0, & |\omega| > \frac{2\pi}{T_B} \end{cases}\)
• 冲激响应：\(h(t) = \frac{\sin(\pi t/T_B)}{\pi t/T_B} \cdot \frac{\cos(\pi t/T_B)}{1-4t^2/T_B^2}\)
• 核心特性：总带宽\(B=f_N(1+\alpha)\)，\(\alpha\)越大带宽越高、频带利用率越低；时域拖尾按\(1/t^3\)快速衰减，抗定时误差能力强，过渡带平滑易实现，是实际系统的主流波形成型方案。

二、系统设计

1.系统功能

1.1 整体布局

系统主界面采用左右分栏结构：

• 左侧控制面板占窗口约 1/4 宽度，自上而下分为参数配置区、仿真结果文本区与 ISI 定量分析表格区，集成所有参数输入控件与功能按钮。
• 右侧波形展示区占窗口约 3/4 宽度，采用 “上排 3 幅、下排 2 幅” 的对称布局，共 5 个独立坐标轴，分别对应传输链路不同环节的波形可视化，所有子图均支持独立缩放与平移。

1.2 可配置参数

系统提供 4 项核心传输参数，可快速切换不同的基带传输场景：

• 输入波形：包含理想冲激、50% 占空比方波、100% 占空比方波三类，分别对应理论分析模型与实际基带码型。
• 滤波器类型：提供理想低通、升余弦滚降两种经典无 ISI 传输特性，切换滤波器时自动联动滚降系数的输入状态。
• 滚降系数 α：取值范围为 (0,1]，控制升余弦滤波器过渡带的宽度，理想低通模式下自动置 0 并禁用编辑。
• 抽样间隔 TB：取值范围为 (0,2]，决定码元抽样时刻与奈奎斯特带宽，直接影响码间干扰的判定结果。

1.3 核心功能

1. 伪随机码元生成：内置 5 级线性反馈移位寄存器，生成周期为 31 的 PN 伪随机序列，并转换为双极性码元用于传输仿真，模拟真实的随机基带信号。
2. 全链路波形可视化：从时域到频域、从单脉冲到叠加信号，完整呈现基带传输的信号变换过程。
3. 码间干扰定性判定：通过抽样判决恢复输出码元序列，自动判定系统是否存在码间干扰。
4. ISI 定量分析：逐码元计算抽样点的理想值、实际值与干扰幅值，量化码间干扰的强度。

function updateAlphaState() handles = guidata(fig); filterType = get(handles.filterPopup, 'Value'); if filterType == 1 set(handles.alphaEdit, 'String', '0'); set(handles.alphaEdit, 'Enable', 'off'); set(handles.alphaEdit, 'BackgroundColor', disableColor); else set(handles.alphaEdit, 'String', '0.5'); set(handles.alphaEdit, 'Enable', 'on'); set(handles.alphaEdit, 'BackgroundColor', [0.18 0.18 0.18]); end end function alphaChanged(~, ~) handles = guidata(fig); filterType = get(handles.filterPopup, 'Value'); if filterType == 1 set(handles.alphaEdit, 'String', '0'); return; end alpha = str2double(get(handles.alphaEdit, 'String')); if isnan(alpha) || alpha <= 0 || alpha > 1 set(handles.alphaEdit, 'String', '0.5'); set(handles.resultBox, 'String', ... '参数错误：滚降系数α必须在(0,1]范围内，已重置为0.5'); end end function genPN31() handles = guidata(fig); reg = randi([0 1], 1, 5); while all(reg == 0) reg = randi([0 1], 1, 5); end pn = zeros(1, 31); for i = 1:31 pn(i) = reg(5); fb = mod(reg(5) + reg(2), 2); reg = [fb, reg(1:4)]; end handles.pn = pn; handles.symbols = 2*pn - 1; guidata(fig, handles); end

2.波形显示

波形显示是系统的核心可视化载体，5 个子图从不同维度呈现基带传输的物理过程，辅助理解码间干扰的产生与消除机理。

2.1 时域波形组

• 输入信号：根据选定的波形类型与码元序列，生成对应的发送端基带波形，图中标注码元分隔虚线与各码元的二进制数值，清晰对应每个码元的时间位置与幅值特性。
• 单个脉冲响应叠加：将每个码元单独经过系统后的冲激响应以不同颜色绘制，同时叠加得到总接收波形。通过该图可直观理解码间干扰的来源：单个码元的脉冲拖尾延伸至相邻码元周期，多码元叠加后形成抽样点上的干扰分量。
• 接收信号：展示所有码元冲激响应叠加后的最终接收波形，并以绿色标记标注每个抽样时刻的位置与幅值。抽样点的时间间隔由参数 TB 控制，通过观察抽样点幅值与理想码元幅值的偏差，可直观判断码间干扰的强弱。

2.2 频域与时域特性组

• 传递函数 H (ω)：频域展示系统的幅频特性曲线，同时叠加对应抽样速率的奈奎斯特理想低通门限。该图可直观验证奈奎斯特第一准则的几何意义：将传输特性以\(2\pi/T_B\)为间隔分段平移叠加后，若等效为理想低通，则系统无码间干扰。
• 单位冲激响应 h (t)：展示系统的时域单位冲激响应，并标注对应抽样间隔的抽样点。通过观察抽样点是否落在冲激响应的零点位置，可直接验证无 ISI 的时域条件，同时可直观对比不同滚降系数下冲激响应拖尾的衰减速度差异。

function input_wave = genInputWave(t, symbols, N, Tb0, waveType) input_wave = zeros(size(t)); switch waveType case 1 for i = 1:N [~, idx] = min(abs(t - (i-0.5)*Tb0)); input_wave(idx) = symbols(i); end case 2 for i = 1:N t_start = (i-1)*Tb0; t_mid = (i-0.5)*Tb0; t_end = i*Tb0; if symbols(i) == 1 idx1 = (t >= t_start) & (t < t_mid); idx2 = (t >= t_mid) & (t < t_end); input_wave(idx1) = 1; input_wave(idx2) = -1; else idx = (t >= t_start) & (t < t_end); input_wave(idx) = -1; end end input_wave(end) = -1; case 3 for i = 1:N t_start = (i-1)*Tb0; t_end = i*Tb0; idx = (t >= t_start) & (t < t_end); input_wave(idx) = symbols(i); end input_wave(end) = symbols(end); end end function H = raisedCosineSpectrum(w, alpha, Tb0) wc = pi/Tb0; H = zeros(size(w)); for i = 1:length(w) wi = abs(w(i)); if wi <= wc*(1-alpha) H(i) = 1; elseif wi <= wc*(1+alpha) H(i) = 0.5 * (1 + cos(pi/(2*alpha*wc) * (wi - wc*(1-alpha)))); else H(i) = 0; end end end function [t_imp, h] = calcImpulseResponse(filterType, alpha, span, fs, Tb0) t_imp = linspace(-span/2*Tb0, span/2*Tb0, span*fs + 1); h = zeros(size(t_imp)); switch filterType case 1 h = sinc(t_imp / Tb0); case 2 for i = 1:length(t_imp) ti = t_imp(i); if ti == 0 h(i) = 1; elseif abs(ti) == Tb0/(2*alpha) && alpha > 0 h(i) = (pi/4) * sinc(1/(2*alpha)); else h(i) = sinc(ti/Tb0) * cos(pi*alpha*ti/Tb0) / (1 - (2*alpha*ti/Tb0)^2); end end end h = h / max(abs(h)); end

3.结果验证

结果验证模块从定性结论与定量数据两个层面，完成对无码间干扰传输条件的仿真验证。

3.1 定性结果输出

系统自动完成抽样判决流程，输出三项核心结论：

• 输入序列：展示当前仿真使用的前 10 位二进制伪随机码元。
• 输出序列：对接收波形按 TB 间隔抽样，以 0 为判决门限进行硬判决，恢复出二进制码元序列。
• ISI 判定：计算所有抽样点实际值与理想值的最大偏差，当偏差低于设定阈值时判定为无码间干扰，否则判定为存在码间干扰。

3.2 定量分析表格

系统以表格形式输出前 10 个码元的量化分析数据，包含四项指标：

• 码元序号：对应传输码元的位置编号。
• 理想值：发送码元对应的标准双极性幅值（±1）。
• 实际值：接收波形在该抽样时刻的实际幅值。
• ISI 值：实际值与理想值的绝对偏差，精准量化该码元受到的码间干扰强度。

function updateSimResult(received, t, symbols) handles = guidata(fig); pn = handles.pn; N = 10; sample_Tb = str2double(get(handles.tbEdit, 'String')); if isnan(sample_Tb) || sample_Tb <= 0 || sample_Tb > 2 sample_Tb = 1; end in_seq = sprintf('%d', pn(1)); for i = 2:10 in_seq = [in_seq sprintf(' %d', pn(i))]; end out_vals = zeros(1, N); for i = 1:N sample_time = (i-0.5)*sample_Tb; [~, idx] = min(abs(t - sample_time)); out_vals(i) = received(idx); end out_bits = out_vals > 0; out_seq = sprintf('%d', out_bits(1)); for i = 2:10 out_seq = [out_seq sprintf(' %d', out_bits(i))]; end isi_vals = abs(out_vals - symbols(1:N)); if max(isi_vals) < 1e-6 isi_status = '无'; else isi_status = '有'; end result_str = {['输入序列：' in_seq], ... ['输出序列：' out_seq], ... ['有无ISI：' isi_status]}; set(handles.resultBox, 'String', result_str); end function updateISITable(symbols, received, t, N, sample_Tb) handles = guidata(fig); tableData = cell(N, 4); for i = 1:N sample_time = (i-0.5)*sample_Tb; [~, idx] = min(abs(t - sample_time)); actual_val = received(idx); ideal_val = symbols(i); isi_val = abs(actual_val - ideal_val); tableData{i, 1} = sprintf('%d', i); tableData{i, 2} = sprintf('%.3f', ideal_val); tableData{i, 3} = sprintf('%.3f', actual_val); tableData{i, 4} = sprintf('%.3f', isi_val); end set(handles.isiTable, 'Data', tableData); end

三、效果展示