从“被砍掉的频谱“到无码间串扰：升余弦滚降滤波器的完全解读

从"被砍掉的频谱"到无码间串扰：升余弦滚降滤波器的完全解读

前言

在数字通信中，发送滤波器的频谱在奈奎斯特频率 fNf_NfN 之前就开始衰减，这并非"砍掉"有用频率，而是为了实现无码间串扰（ISI）的必然设计。本文将结合时域与频域条件，严格推导升余弦滚降滤波器的数学表达式，并验证其满足无ISI准则，同时区分升余弦与根升余弦的概念，解释滚降消除ISI的物理本质。

---

1. 码间串扰（ISI）的定义

在基带传输中，每个符号由一个脉冲表示。接收端在每个符号的中心时刻采样，理想情况下采样值仅与当前符号有关。然而，实际信道或滤波器是带限的，每个脉冲的响应会延伸至相邻符号的采样时刻。若该"拖尾"在相邻符号的采样点处不为零，则叠加后会导致采样值偏离真实值，此即码间串扰。

---

2. 无ISI的时域条件

若想消除码间串扰，即相邻码元的前一个码元的波形到达后一个码元抽样判决的时刻已经衰减到0，如下图示的波形就行满足要求。

其对应的系统冲激响应 h(t)h(t)h(t) 应满足：

h(kTB)={1,k=00,k为非零整数 h(kT_B) = \begin{cases} 1, & k=0 \\ 0, & k \text{为非零整数} \end{cases} h(kTB)={1,0,k=0k为非零整数

即在本符号采样时刻取最大值，在其他符号采样时刻为零（过零点条件），如下图所示。

---

3. 无ISI的频域条件：奈奎斯特第一准则

时域过零点条件等价于频域条件：

∑i=−∞∞H(ω+2πiTB)=TB,∣ω∣≤πTB \sum_{i=-\infty}^{\infty} H\left(\omega + \frac{2\pi i}{T_B}\right) = T_B, \quad |\omega| \leq \frac{\pi}{T_B} i=−∞∑∞H(ω+TB2πi)=TB,∣ω∣≤TBπ

该式表明：将系统频率响应 H(ω)H(\omega)H(ω) 以 2π/TB2\pi/T_B2π/TB 为周期重叠相加，在基带 ∣ω∣≤π/TB|\omega| \leq \pi/T_B∣ω∣≤π/TB 内应为常数。这意味着 H(ω)H(\omega)H(ω) 必须关于 fN=1/(2TB)f_N = 1/(2T_B)fN=1/(2TB) 奇对称滚降。具体图示如下：

---

4. 理想低通传输特性

理想低通滤波器满足奈奎斯特准则，其系统函数如图2左图所示：

H(ω)={TB,∣ω∣≤πTB0,∣ω∣>πTB H(\omega) = \begin{cases} T_B, & |\omega| \leq \frac{\pi}{T_B} \\ 0, & |\omega| > \frac{\pi}{T_B} \end{cases} H(ω)={TB,0,∣ω∣≤TBπ∣ω∣>TBπ

其冲激响应为 sinc(πt/TB)\text{sinc}(\pi t/T_B)sinc(πt/TB)，满足 h(kTB)=0 (k≠0)h(kT_B)=0\ (k\neq0)h(kTB)=0 (k=0)。

• 优点 ：频带利用率最高 η=2\eta = 2η=2 Baud/Hz。
• 缺点 ：
  1. 物理不可实现（陡峭截止）。
  2. 时域拖尾按 1/t1/t1/t 衰减，缓慢且对定时误差敏感。

---

5. 升余弦滚降滤波器

为克服理想低通的缺陷，引入平滑滚降的频谱，即升余弦滚降滤波器。

5.1 频域表达式

设符号周期为 TTT，滚降系数 α∈[0,1]\alpha \in [0,1]α∈[0,1]，奈奎斯特频率 fN=1/(2T)f_N = 1/(2T)fN=1/(2T)。升余弦频谱为：

H(f)={T,0≤∣f∣≤1−α2TT2[1+cos⁡(πTα(∣f∣−1−α2T))],1−α2T<∣f∣≤1+α2T0,∣f∣>1+α2T H(f) = \begin{cases} T, & 0 \le |f| \le \dfrac{1-\alpha}{2T} \\[6pt] \dfrac{T}{2}\left[1 + \cos\left(\dfrac{\pi T}{\alpha}\left(|f| - \dfrac{1-\alpha}{2T}\right)\right)\right], & \dfrac{1-\alpha}{2T} < |f| \le \dfrac{1+\alpha}{2T} \\[6pt] 0, & |f| > \dfrac{1+\alpha}{2T} \end{cases} H(f)=⎩ ⎨ ⎧T,2T[1+cos(απT(∣f∣−2T1−α))],0,0≤∣f∣≤2T1−α2T1−α<∣f∣≤2T1+α∣f∣>2T1+α

在滚降段，表达式为 T2[1+cos⁡(⋅)]\frac{T}{2}[1+\cos(\cdot)]2T[1+cos(⋅)]，即抬升的余弦函数，故得名"升余弦"。

5.2 时域表达式

对上述频谱作傅里叶逆变换，得到时域冲激响应：

h(t)=sinc⁡(πtT)⋅cos⁡(παt/T)1−(2αt/T)2 h(t) = \operatorname{sinc}\left(\frac{\pi t}{T}\right) \cdot \frac{\cos(\pi \alpha t / T)}{1 - (2\alpha t / T)^2} h(t)=sinc(Tπt)⋅1−(2αt/T)2cos(παt/T)

其中 sinc⁡(x)=sin⁡(x)/x\operatorname{sinc}(x) = \sin(x)/xsinc(x)=sin(x)/x。

---

6. 验证升余弦满足无ISI条件

6.1 时域验证

• 当 t=kTt = kTt=kT（k≠0k \neq 0k=0 整数），sinc⁡(πk)=0\operatorname{sinc}(\pi k)=0sinc(πk)=0，分子 cos⁡(παk)\cos(\pi \alpha k)cos(παk) 有限，分母 1−(2αk)21-(2\alpha k)^21−(2αk)2 非零（奇异点处极限仍为0），故 h(kT)=0h(kT)=0h(kT)=0。
• 当 k=0k=0k=0，lim⁡t→0h(t)=1\lim_{t\to0}h(t)=1limt→0h(t)=1（归一化后）。
  因此时域过零点条件成立。

6.2 频域验证

在 fN=1/(2T)f_N = 1/(2T)fN=1/(2T) 处，升余弦谱满足奇对称：H(fN+Δf)+H(fN−Δf)=TH(f_N + \Delta f) + H(f_N - \Delta f) = TH(fN+Δf)+H(fN−Δf)=T（常数）。

• 平坦段 0≤∣f∣≤(1−α)/(2T)0 \le |f| \le (1-\alpha)/(2T)0≤∣f∣≤(1−α)/(2T)，幅值为 TTT。
• 滚降段 (1−α)/(2T)<∣f∣≤(1+α)/(2T)(1-\alpha)/(2T) < |f| \le (1+\alpha)/(2T)(1−α)/(2T)<∣f∣≤(1+α)/(2T)，由余弦抬升的互补性保证和值为常数。
• 超出部分为零。
  故而满足奈奎斯特第一准则，实现无ISI传输。

---

7. 升余弦与根升余弦的区别

特性	升余弦 (RC)	根升余弦 (RRC)
定义	总系统响应	HRC(f)\sqrt{H_{\text{RC}}(f)}HRC(f)
时域过零点	在 kTkTkT 处为零	不为零（单独使用有ISI）
发送端单独使用	可，接收端只采样	不可，需接收端匹配RRC
总响应	HRC(f)H_{\text{RC}}(f)HRC(f)	HRRC(f)⋅HRRC(f)=HRC(f)H_{\text{RRC}}(f) \cdot H_{\text{RRC}}(f) = H_{\text{RC}}(f)HRRC(f)⋅HRRC(f)=HRC(f)
抗噪声性能	非匹配，较差	匹配滤波，最优
实际系统	少见	主流（QAM、OFDM等）

系统冲激响应 h(t)h(t)h(t) 总是指整个传输链路的总响应，在实际采用根升余弦成形的系统中，总响应即为升余弦。

---

8. 代码演示

通过MATLAB代码可以直观对比：矩形脉冲经理想低通滤波后，采样点受相邻符号干扰而偏离（有ISI）；升余弦总脉冲的采样点则精确等于原始符号（无ISI）。

clear; clc; close all;

%% 参数
T = 1;                    % 符号周期
sps = 100;                % 采样率
t = -3:1/sps:6;           % 时间轴（覆盖三个符号及前后拖尾）
fc = 0.5/T;               % 理想低通截止频率 = 符号率/2

%% 发送三个符号：+1, -1, +1
a = [1, -1, 1];
t0 = 0; t1 = T; t2 = 2*T;   % 符号发送时刻

%% 1) 矩形脉冲 → 通过理想低通（产生ISI）
% 矩形脉冲（持续时间为T）
rect_pulse = @(tau) double(abs(tau) < T/2);
% 构造矩形脉冲串
s_rect = a(1)*rect_pulse(t - t0) + a(2)*rect_pulse(t - t1) + a(3)*rect_pulse(t - t2);
% 理想低通滤波器的冲激响应 (sinc)
h_lpf = 2*fc * sinc(2*fc * t);   % 带宽 fc，采样率足够高
% 卷积（滤波）
s_lpf = conv(s_rect, h_lpf, 'same') / sps;   % 归一化
% 归一化：使得单个矩形脉冲通过后的峰值约为1
s_lpf = s_lpf / max(abs(s_lpf));   % 简单归一化，方便观察

% 采样时刻（符号周期中心）
sample_times = [t0, t1, t2];
samples_lpf = interp1(t, s_lpf, sample_times, 'linear');

%% 2) 升余弦滚降（无ISI）
alpha = 0.5;                  % 滚降系数
span = 6;                     % 脉冲长度（符号数）
h_rc = rcosdesign(alpha, span, sps, 'normal');
% 构造三个移位脉冲之和
t_h = (-span/2 : 1/sps : span/2);
pulse0 = interp1(t_h, h_rc, t - t0, 'linear', 0);
pulse1 = interp1(t_h, h_rc, t - t1, 'linear', 0);
pulse2 = interp1(t_h, h_rc, t - t2, 'linear', 0);
s_rc = a(1)*pulse0 + a(2)*pulse1 + a(3)*pulse2;
% 采样点
samples_rc = interp1(t, s_rc, sample_times, 'linear');

%% 绘图对比
figure('Position', [100 100 1200 500]);

% 左图：矩形+理想低通（有ISI）
subplot(1,2,1);
plot(t, s_lpf, 'b', 'LineWidth', 1.5); hold on;
stem(sample_times, samples_lpf, 'ro', 'LineWidth', 2);
stem(sample_times, a, 'k^', 'LineWidth', 1.5, 'MarkerSize', 8);
grid on; title('矩形脉冲 → 理想低通 (严重 ISI)');
xlabel('时间 t / T'); ylabel('幅度');
legend('滤波后波形', '接收采样值', '原始符号值', 'Location', 'best');
text(sample_times(2)+0.1, samples_lpf(2), sprintf('采样 = %.2f', samples_lpf(2)), 'FontSize', 11);
ylim([-1.5 1.8]);

% 右图：升余弦滚降（无ISI）
subplot(1,2,2);
plot(t, s_rc, 'b', 'LineWidth', 1.5); hold on;
stem(sample_times, samples_rc, 'ro', 'LineWidth', 2);
stem(sample_times, a, 'k^', 'LineWidth', 1.5, 'MarkerSize', 8);
grid on; title(sprintf('升余弦滚降 (α = %.1f) → 无 ISI', alpha));
xlabel('时间 t / T'); ylabel('幅度');
legend('发送波形', '接收采样值', '原始符号值', 'Location', 'best');
text(sample_times(2)+0.1, samples_rc(2), sprintf('采样 = %.2f', samples_rc(2)), 'FontSize', 11);
ylim([-1.5 1.8]);

% 显示数值差异
fprintf('理想低通（矩形脉冲）第二个符号采样值 = %.4f  (期望 -1)\n', samples_lpf(2));
fprintf('升余弦滚降第二个符号采样值     = %.4f  (期望 -1)\n', samples_rc(2));

---

9. 常见误区澄清

• 误区1 ：滚降"砍掉"了有用频谱。
  事实 ：滚降将能量从 0∼fN0 \sim f_N0∼fN 重新分配到 0∼(1+α)fN0 \sim (1+\alpha)f_N0∼(1+α)fN，实现频域奇对称，时域拖尾快速归零。

• 误区2 ：升余弦优于根升余弦。
  事实：两者总响应等价，但根升余弦配合匹配滤波可获得最佳信噪比，实用中更常见。

• 误区3 ：平滑波形仅为美观。
  事实 ：平滑带来实际益处：对定时误差不敏感、物理可实现、拖尾衰减快（α>0\alpha>0α>0 时按 1/t31/t^31/t3 衰减）。

---

10. 总结与补充

• 升余弦的时域与频域表达式严格满足无ISI的时域过零点条件和频域奈奎斯特第一准则。

• 系统总响应（升余弦）是设计目标，根升余弦是收发分配合法实现方式。

• 滚降消除ISI的本质：频域奇对称 ⇔ 时域拖尾在整数倍符号周期处为零。

• 可以通过调整滚降系数 α\alphaα 可观察时域波形从 sinc 状剧烈振荡到平滑过零点的变化，加深理解。

📢 如果本文对你有帮助，欢迎关注「手搓物理层 」

后台回复关键字 "ISI" 或者获取本文完整PDF + 配套MATLAB代码包。

文中图例均引用于樊昌信《通信原理》