手把手实战 MATLAB 通信系统仿真 06

第 6 章 抽样、量化与 PCM本章目标通过 MATLAB 直观观察采样混叠现象理解采样定理的工程约束手工实现均匀量化器中间型量化掌握量化限幅、分区、索引映射全过程完成量化索引与 PCM 二进制比特流的编码/解码闭环建立模拟信号数字化ADC的完整认知。6.1采样定理不止是「Fs2fmax」很多教材简写采样定理为采样率大于信号最高频率的2倍。但在仿真与工程中这句话并不完整。完整表述带限信号的采样率必须严格大于2倍最高频率才能保证频谱不重叠、无混叠可无失真恢复原始信号。工程补充关键细节理想无陡峭滤波器场景下不能将采样率刚好卡在 $$F_s2f_\text{max}$$会造成频谱混叠实际硬件抗混叠滤波器存在过渡带必须预留保护带宽因此工程采样率通常远高于2倍最高频率混叠是不可逆失真一旦发生采样混叠高频信号会折叠到低频无法通过后续滤波恢复。MATLAB 仿真的核心价值肉眼观察混叠现象——同一个高频连续信号在欠采样条件下会完全等效为一个低频离散信号。代码观察混叠f0 9000; % 原始信号 9kHz Fs 12000; % 采样率 12kHz不满足奈奎斯特 N 120; n (0:N-1).; % 对9kHz高频信号欠采样 x cos(2*pi*f0*n/Fs); figure; stem(n,x,filled); grid on; xlabel(n); ylabel(x[n]); title(9 kHz 信号以 12 kHz 采样发生混叠); % 混叠折叠频率Fs-f0 3kHz xAlias cos(2*pi*3000*n/Fs); fprintf(max difference %.3e\n,max(abs(x-xAlias)));实验现象与原理输出差值趋近于 0说明欠采样后的 9kHz 序列与 3kHz 低频序列完全一致。离散时间频率折叠公式。这证明计算机看到的离散波形无法区分原始高频和混叠后的低频混叠失真彻底不可逆。6.2 手工实现中间型均匀量化化是模拟信号数字化的核心步骤将连续幅度的信号映射为有限个离散电平。本章实现中间型均匀量化每段量化区间取中点作为重构值是最标准、最常用的量化模型。均匀量化核心参数B量化比特数量化总层数量化步长量化完整流程限幅 → 区间划分 → 计算索引 → 中点重构。代码B 位均匀量化器B 4; L 2^B; xmin -1; xmax 1; Delta (xmax-xmin)/L; % 量化步长 % 生成超出量化范围的输入观察饱和失真 x linspace(-1.2,1.2,1000).; % 1. 限幅超出量程部分强制截断 xClip min(max(x,xmin),xmax-eps); % 2. 计算量化索引 0 ~ L-1 idx floor((xClip-xmin)/Delta); % 3. 中点重构中间型均匀量化 xq xmin (idx0.5)*Delta; figure; plot(x,x,--); hold on; plot(x,xq,LineWidth,1.2); grid on; xlabel(Input Amplitude); ylabel(Quantized Output); title(4bit 均匀量化特性曲线);实验观察输入信号在 [-1,1] 范围内输出呈现标准阶梯量化波形输入超出 ±1 量程后输出不再变化出现量化饱和失真过载失真采用区间中点重构量化误差范围严格控制在。6.3 从量化索引到 PCM 比特流PCM脉冲编码调制的本质将模拟信号量化后的十进制索引转为固定长度二进制比特流传输接收端再反向恢复索引与量化电平。这是数字语音、数字通信最基础的编码框架核心是固定比特长度、无损索引映射。代码PCM 编码与解码【基础 MATLAB】B 4; % 一组量化索引示例0~15 合法范围 idx [0; 1; 5; 10; 15]; % 编码十进制索引 → B位二进制矩阵 → 一维比特流 bitMat dec2bin(idx,B) - 0; % 转为数值0/1矩阵 bitStream reshape(bitMat.,[],1);% 串行比特流 % 解码一维比特流 → 重构比特矩阵 → 恢复索引 bitMatRx reshape(bitStream,B,[]).; idxHat bin2dec(char(bitMatRx 0)); % 校验编解码无损 assert(isequal(idx,idxHat)); disp(PCM编解码闭环验证成功无失真);关键核心细节dec2bin自动补齐高位保证每个索引严格占用B 个比特reshape 转置必须成对反向否则比特顺序错乱、解码错误assert 无报错证明 PCM 编解码是无损索引传输。6.4 量化信噪比的直观理解在满幅、均匀分布、无过载、均匀量化的理想条件下量化比特数每增加 1bit量化信噪比大约提升 6dB。⚠️ 重要误区纠正6dB/bit 是理想理论值并非通用公式实际场景中语音信号非均匀分布、信号过载失真、非均匀量化、压扩处理都会打破 6dB/bit 的规律比特位数提升可减小量化步长、降低量化噪声是提升量化精度的唯一核心手段。本章小结采样不足会产生不可逆频率混叠工程采样率必须预留滤波过渡带不能仅满足 2 倍频率下限。均匀量化四步流程限幅 → 区间分割 → 索引计算 → 中点重构。信号超量程会出现过载饱和失真是量化失真的重要来源。PCM 本质是量化索引的固定长度二进制序列化传输可实现无损编解码闭环。6dB/bit 量化信噪比规律仅适用于理想均匀量化场景不可通用。动手练习1. 把 B 从 4 改到 8比较量化误差均方值。2. 设计一个输入幅度经常超过 ±1 的信号观察过载失真。3. 把 PCM 比特流随机翻转 1 bit观察恢复的量化索引可能发生多大变化。