NTSC色度副载波相位抖动与VHS磁化衰减的数学建模

问题背景：模拟信号瑕疵的物理根源

NTSC 制式在 1953 年引入彩色电视时，为了兼容黑白接收机，采用了频分复用策略：亮度信号 (Y) 占据低频段，色度信号 (I、Q) 调制在约 3.58MHz 的副载波上。这种设计带来了固有的技术妥协 —— 色度与亮度通道的时序对齐误差、副载波相位抖动、以及磁性存储介质的磁化衰减，共同构成了 VHS 录像带独特的视觉特征。

与简单的 LUT 滤镜不同，ntsc-rs 这类工具的核心价值在于从物理层面还原这些瑕疵的生成机制。本文聚焦于两个关键技术点：色度副载波相位抖动的数学描述，以及磁带磁化衰减的时域模型。

色度副载波相位抖动的数学建模

副载波调制基础

NTSC 色度信号通过正交幅度调制 (QAM) 承载于副载波：

C(t) = I(t)·cos(ωt) + Q(t)·sin(ωt)

其中 ω = 2π・3.579545MHz 为副载波角频率。理想情况下，接收端通过同步解调分离 I、Q 分量。然而实际系统中存在多种相位误差源。

相位抖动的来源分类

1. 时基误差 (Time Base Error)

VHS 磁头与磁带之间的机械速度波动导致水平同步脉冲位置抖动。这种时基抖动可建模为：

φ_jitter(t) = φ_0 + φ_1·sin(2πf_drum·t) + n(t)

• φ_0：静态相位偏移（磁头方位角误差）
• φ_1：磁鼓旋转引起的周期性分量（VHS 磁鼓转速约 30Hz）
• n (t)：高频随机噪声（磁带表面粗糙度、电机转速波动）

2. 色度 - 亮度延迟误差 (Y/C Delay)

VHS 信号处理链中，色度通道经过降频转换（626.9kHz）后再升频，与亮度通道的群延迟差异导致：

Δτ_typical = 100-300ns（设备相关）

在像素层面，这表现为颜色边缘相对于亮度边缘的水平偏移，形成 "color bleeding" 现象。

3. 相位噪声的频谱特性

实际测量表明，VHS 播放机的相位噪声呈现 1/f 特性：

S_φ(f) = S_0·(1 + f_c/f)  （f < f_c时以1/f主导）

其中 f_c 约为数百 Hz，对应机械系统的惯性截止频率。

相位误差对解调的影响

当存在相位误差 φ_e 时，解调后的 I、Q 分量发生串扰：

I' = I·cos(φ_e) - Q·sin(φ_e)
Q' = I·sin(φ_e) + Q·cos(φ_e)

对于小角度近似 (φ_e << 1)：

I' ≈ I - Q·φ_e
Q' ≈ Q + I·φ_e

这意味着 I、Q 通道间的串扰与相位误差成正比。当 φ_e 随时间变化时，产生色度噪声和色调漂移。

VHS 磁带磁化衰减模型

磁性记录的物理基础

VHS 使用氧化铁或二氧化铬磁粉涂覆的聚酯带基。记录时，视频磁头产生交变磁场，使磁粉沿磁道方向磁化。理想磁化强度 M 与记录电流的关系遵循磁滞回线。

磁化衰减机制

磁带老化导致的信号劣化主要源于三种机制：

1. 热衰减（磁各向异性能弛豫）

单畴颗粒的磁化方向随时间热波动发生反转，遵循 Arrhenius 型弛豫：

M(t)/M_0 = exp(-t/τ)
τ = τ_0·exp(KV/kT)

其中 K 为磁各向异性常数，V 为颗粒体积，τ_0 约为 10^-9 秒。对于 VHS 磁带，KV/kT 通常在 25-40 范围，对应 τ 从数年到数十年。

2. 退磁场引起的自退磁

相邻磁化区域的退磁场导致短波长信号衰减更严重。对于波长 λ 的信号，衰减因子：

A(λ) = exp(-2πd/λ)

d 为磁性层厚度（VHS 约 4-5μm）。这意味着高频分量（对应图像细节）比低频衰减更快。

3. 磁头 - 磁带接触劣化

反复播放导致磁粉脱落、带基变形，引入附加噪声。这种非线性损伤可用随机脉冲噪声建模：

n_dropout(x,y,t) = Σ A_i·δ(x-x_i, y-y_i, t-t_i)

综合衰减模型

将上述机制整合，磁带传输函数可表示为：

H_tape(f,t) = H_0(f)·exp(-t/τ_eq(f))·(1 + n_dropout)

其中 τ_eq (f) 为频率相关的等效弛豫时间，高频分量 τ_eq 更小。

工程实现参数清单

基于上述模型，在 ntsc-rs 或类似系统中实现时，建议采用以下参数化方案：

相位抖动参数

磁化衰减参数

关键实现要点

1. YIQ 域处理

所有色度相关操作应在 YIQ 空间进行：

// RGB → YIQ转换矩阵
let yiq = rgb_to_yiq(rgb);
// 应用相位抖动到I、Q通道
let (i, q) = apply_phase_jitter(yiq.i, yiq.q, phase_error);
// 应用磁化衰减（主要影响高频）
let y_filtered = apply_hf_rolloff(yiq.y, hf_loss);
// YIQ → RGB转换
let rgb_out = yiq_to_rgb(YIQ { y: y_filtered, i, q });

2. 相位抖动的帧间连续性

为避免闪烁，相位误差应在帧间保持连续性：

// 使用一阶AR模型保持时序相关性
phase_error[t] = α·phase_error[t-1] + (1-α)·noise[t]

建议 α 取值 0.8-0.95，对应约 5-20 帧的时间常数。

3. 磁化衰减的空间分布

真实磁带的老化不均匀，可引入空间变化：

// 基于噪声纹理的局部老化系数
let local_aging = base_aging + noise_texture(x, y) · aging_variance;

验证与调优策略

主观质量评估：使用 SMPTE 彩条测试图，观察以下特征：

• 白黄青绿品红蓝黑条的边缘是否出现彩虹效应（验证相位误差）
• 高饱和度区域的色度噪声水平
• 黑色区域的亮度噪声颗粒度

客观指标：

• 色度 - 亮度互相关峰值位置（验证延迟误差）
• 高频响应衰减曲线斜率

参数联动关系：

• 增加 phase_drift_amp 会放大 chroma_delay 的视觉影响
• tape_age_years 应与 hf_loss_db 联动设置，典型比例约 0.3-0.5 dB / 年

总结

NTSC 色度副载波相位抖动与 VHS 磁化衰减的精确建模，需要将模拟信号的物理机制转化为可计算的数学描述。核心在于理解：相位误差导致的 I/Q 串扰遵循三角函数关系，而磁化衰减呈现指数型频率相关特性。通过 YIQ 域的分通道处理、帧间连续的相位噪声生成、以及频率相关的增益衰减，可以在现代数字系统中高保真地重现模拟时代的视觉特征。

参考来源

• ntsc-rs - 开源 VHS 效果模拟器
• zhuker/ntsc - NTSC 视频模拟器 Python 实现
• Marginal Futility - NTSC Part 1 - NTSC 信号处理技术分析

systems