简介这套资料是一份关于语音信号频带传输通信系统仿真的课程设计文档基于PCM编码与PSK调制技术使用MATLAB/Simulink平台完成。内容覆盖PCM采样、量化与编码PSK调制解调高斯白噪声信道建模以及完整收发端系统流程适合通信工程、电子信息类专业学生用于课程设计、仿真实践或毕业设计参考。文档共1个docx文件压缩包大小约3.29MB包含26页设计报告从引言、设计原理到模块参数设置与仿真结果分析均有展开结构完整。资源已有359人学习下载说明其对同类课题具参考价值。读者可直接借助其中的系统框图、模块配置说明和结果对比思路快速搭建自己的语音频带传输仿真模型减少调试摸索时间。1. 从PCM编码到PSK调制这套频带传输仿真到底在做什么你在设计一个语音对讲模块时话筒送出的模拟信号并不能直接发射必须先经过PCM编码转成比特流再通过PSK调制把比特流搬到合适的载波频段上这就是“语音信号频带传输”的基本路径。仿真这件事不是为了把流程图跑一遍而是要把每一级信号的波形、频谱、误码率都量化出来回答几个实际问题采样率取多少量化位数够不够载波频率会不会和符号率打架信道噪声大到什么程度还能听懂语音。本文按“先数字化、再调制、最后跑通全链路”的顺序用MATLAB把这条链路从零搭起来并给出参数调整和验证方法。适合正在做通信课程设计、以及想从“看公式”转向“看波形”的从业者。本文不依赖任何硬件平台所有结果都能在普通电脑上复现。2. PCM编码仿真把模拟语音变成可传输的比特流2.1 语音数字化中的采样、量化、编码三步拆解语音信号是带限模拟信号典型电话语音频带是300Hz到3.4kHz。根据奈奎斯特定理采样率至少是最高频率的两倍实际通常取8kHz。采样只完成了时间离散化幅度仍然是连续的需要量化到有限个电平。PCM的核心就是“采样-量化-编码”三步采样确定时间点量化把连续幅度映射到2^n个离散电平编码把每个电平写成n位二进制码字。这里n通常取8再加上8kHz采样率得到64kbps的速率也就是标准电话PCM的速率。量化误差是PCM最直接的失真源。均匀量化在信号幅度分布不均匀时浪费比特位因此实际系统常采用A律或μ律压扩对小信号用更短的量化间隔。在仿真中我们既可以先忽略压扩用均匀量化验证链路也可以用MATLAB的compand函数模拟压扩效果。量化位数n直接决定信噪比每增加1位量化信噪比约提升6dB——这是选型时判断“8位还是12位”的基本依据。真正做语音仿真时还要注意语音信号的峰值因子较大归一化方式会影响量化精度一般用max(abs(x))做归一化而不是用均方根值。2.2 用MATLAB实现PCM编码从采样到输出比特流的最小代码下面这段代码生成一段测试音完成采样、量化和二进制编码。它不读取外部语音文件便于你首先确认算法本身是否正确。Fs 8000; % 采样率 8kHz t 0:1/Fs:0.08; % 80ms 时长 x sin(2*pi*800*t) 0.5*sin(2*pi*2000*t); % 双音信号 nBits 8; % 量化位数 L 2^nBits; % 量化电平数 xMax max(abs(x)); % 归一化 xNorm x / xMax; % 归一化到 [-1, 1] % 均匀量化把 [-1,1] 映射到 [0, L-1] xq round((xNorm 1) * (L-1) / 2); xq max(min(xq, L-1), 0); % 防止舍入越界 % 十进制量化电平转二进制比特流 bits de2bi(xq, nBits, left-msb); % 每行一个采样点左对齐 bitStream reshape(bits., 1, []); % 按时间顺序串行输出逻辑说明先归一化信号幅度保证量化器满量程利用xq是量化索引范围0到L-1de2bi函数把每个索引转换为nBits位二进制向量left-msb表示最高有效位在前最后reshape成一行连续的比特流。整个过程保持了采样点顺序后续PSK调制直接把bitStream作为输入。参数说明Fs8000对应窄带语音下限若信号包含更高频率必须按实际最高频率提高采样率nBits8是常规语音质量追求更高保真可换成16但输出比特率翻倍这里使用均匀量化若改为μ律压扩只需在量化前对xNorm用compand处理。2.3 PCM仿真落地的3个必调参数第一是采样率。仿真中看到频谱混叠首先检查采样率是否大于最高频率的2倍若信号是通用的4kHz音频8kHz刚好够但留一点余量更稳比如取16kHz。第二是量化位数。位数的效果要结合回放或SNR来看建议分别跑8、12、16位对比误码率或量化信噪比你会看到每多1位约6dB的改善。第三是压扩选择。A律和μ律的差异主要在低幅度区间的信噪比如果测试信号中包含极小声段使用μ律/律能让重建波形在示波器上更贴近原信号。下表给出了PCM仿真中需要权衡的核心参数作为初始配置参考参数典型取值影响对象采样率 Fs8kHz / 16kHz最高可恢复频率决定频谱占用量化位数 nBits8 / 12 / 16输出比特率量化信噪比约6n dB压扩律无、A律、μ律小信号量化信噪比语音主观质量这里要注意PCM仿真的“坑”常出在量化边界当信号幅度刚好等于最大值时round可能产生L越界后必须clip。上段的max(min())就是干这个的。另一个坑是de2bi的列顺序默认最低有效位在前不加left-msb会导致比特顺序错乱解调后语音完全不可懂。3. PSK调制仿真把比特流搬上载波实现频带传输3.1 为什么频带传输要选PSK而不是ASK或FSKPCM输出的比特流是基带信号频率集中在低频段无法直接通过带通型信道比如无线频道、电缆频分复用传输。频带传输的核心是选择一个载波频率把比特信息调制到载波的某个参数上。ASK调制幅度FSK调制频率PSK调制相位。PSK在相同信噪比下误码率最低功率效率高但相位模糊是它的固有缺陷。在仿真中我们常先做BPSK每符号1比特验证链路后再扩展QPSK每符号2比特以提升频谱效率。选择PSK也是因为它很适合理想信道仿真接收端用相干解调即可达到理论误码率不需要复杂的信道估计就能对照公式验证结果。调制过程的等效含义是把一串比特映射成一串复平面上带幅相的点再通过脉冲成型滤波器限制带宽最后乘以载波完成频谱搬移。仿真中“搬移”的体现是发射信号频谱在载波频率两侧对称出现直接用periodogram看即可。3.2 BPSK调制仿真从比特流到已调信号的最小实现首先把前面PCM产生的bitStream当输入。BPSK将比特0/1映射为±1每符号只携带1比特。为了模拟真实系统需要过采样和脉冲成型再与载波相乘。% 系统参数 sps 16; % 每个符号的采样点数 Rb 64000; % 比特率 64kbps来自 8kHz * 8bit Rsym Rb; % BPSK符号率 比特率 Fs Rsym * sps; % 仿真采样率 1.024MHz % 比特流 - 符号 symbols 2 * bitStream - 1; % 0 - -1, 1 - 1 % 上采样 根升余弦脉冲成型 symbolsUp upsample(symbols, sps); beta 0.35; rrc rcosdesign(beta, 6, sps, sqrt); txBaseband filter(rrc, 1, symbolsUp); % 载波调制 fc 200e3; % 载波频率 200kHz t (0:length(txBaseband)-1) / Fs; carrier cos(2*pi*fc*t); txPassband txBaseband .* carrier;逻辑说明upsample在每两个符号间插入sps-1个零相当于升采样rcosdesign生成根升余弦滤波器这里滚降系数0.35滤波器长度为6个符号周期filter完成成型滤波得到连续基带波形。最后乘以余弦载波频谱被搬到fc两侧。这里的仿真采样率Fs必须大于2*(fc信号带宽)否则会产生混叠上例中fc200kHz信号带宽约(1beta)*Rsym/2≈43kHzFs1.024MHz满足要求。参数说明sps影响时域波形的平滑度和仿真负载一般取8~32过小则波形不平滑过大则运行慢。beta决定频带宽度取值0.15~0.5越小越省带宽但脉冲的时域衰减越慢对定时误差越敏感。fc的选择要兼顾仿真采样率太大会导致Fs成倍增大太小则基带和频带频谱容易重叠。输出txPassband就是实际的频带传输信号。下表汇总了PSK调制仿真中的关键参数及其选择逻辑方便你在批量实验时快速定位参数常见范围对系统的影响sps8 ~ 32精度与仿真速度平衡过小区间干扰大beta0.15 ~ 0.5带宽占用与定时鲁棒性fc与Fs满足Fs2(fcB)频谱位置过高增加仿真成本3.3 接收端相干解调本地载波、低通滤波与匹配滤波相干解调要求接收端产生与发射端同频同相的本地载波。仿真中最简单的做法是直接发射载波的副本等于假设理想载波同步。用本地载波乘以接收信号产生基带和2倍载频分量低通滤除高频分量再通过匹配滤波器与发射端相同的根升余弦滤波器并下采样判决。% 假设 rxPassband 是接收信号发射端参数复用 rxBaseband rxPassband .* cos(2*pi*fc*t); % 混频到基带 % 低通滤波用发射端成型滤波器匹配滤波这里只做一次 rxFiltered filter(rrc, 1, rxBaseband); % 下采样在每个符号中心采样 samplingPoints 1 sps/2 : sps : length(rxFiltered); rxSymbols rxFiltered(samplingPoints); % 判决BPSK 0 判为1否则0 rxBits rxSymbols 0;这段代码里混频后没有单独低通滤波器而是直接用了匹配滤波器。由于根升余弦滤波器的频响能压制混频产生的2fc分量通常在仿真中省去一级额外低通前提是采样率足够高、fc远离基带。samplingPoints从sps/2开始取是为了对准脉冲成型的峰值点。如果定时偏移眼图会闭合。需要提醒的是真正系统存在相位偏差、频偏和定时误差。仿真时可以人为给载波加一个相移观察解调星座图的旋转这是后面升级到非理想同步的切入点。4. 全链路仿真加入信道噪声测误码率并调整参数4.1 给频带信号加高斯白噪声用SNR反推噪声功率信道仿真的基本模型是加性高斯白噪声AWGN。给频带信号加噪需要先计算信号平均功率再根据目标信噪比折算噪声功率。这里常用两种方式直接以采样点为单位计算snr或用通信系统常用的Eb/N0。误码率曲线习惯横轴用Eb/N0每比特能量与噪声功率谱密度之比BPSK的理论误码率为Q(sqrt(2*Eb/N0))。代码中可以用awgn函数一步到位EbN0_dB 6; % 目标 Eb/N0单位dB snr_dB EbN0_dB 10*log10(2); % BPSK 时 Es2Eb符号SNR % 注意awgn 默认 snr 按信号功率与噪声功率比 rxPassband awgn(txPassband, snr_dB, measured);awgn的measured表示它以输入信号的实际功率为基准计算噪声。snr_dB换算关系对BPSK每个符号有两个比特能量因此 Es/N0 Eb/N0 10log10(2) dB。若做QPSK每个符号2比特Es/N0 Eb/N0 10log10(2) 同样成立因为每符号能量是2Eb注意区分。如果不想用awgn也可以手动生成噪声noise sqrt(noisePower)*randn(size(txPassband))其中noisePower signalPower / 10^(snr_dB/10)。手动方式便于控制噪声种子和功率适合批量仿真。4.2 搭建完整收发链路从原始语音到重建语音的闭环现在把第2章PCM、第3章调制解调和第4章加噪串成一个脚本。流程是读取/生成语音 → PCM编码 → BPSK调制 → 过信道 → 相干解调 → PCM解码 → 播放/对比。下面给出主流程的关键片段% ---- 发射端 ---- % 假设 x 是归一化语音信号长度 N [x, Fs_audio] audioread(test.wav); % 也可用自己合成的信号 x x(:,1)./max(abs(x)); % 单声道归一化 % PCM编码得到 bitStream调用前面函数 % 注意 bitStream 长度为 8*length(x)因为 nBits8 % ---- 调制 ---- [txPassband, params] psk_modulate(bitStream, sps, fc, Fs); % ---- 信道 ---- rxPassband awgn(txPassband, snr_dB, measured); % ---- 解调 ---- rxBits psk_demodulate(rxPassband, params); % ---- PCM解码 ---- % 比特流转回量化索引 rxBitsMat reshape(rxBits, nBits, []).; xq_recv bi2de(rxBitsMat, left-msb); xr xq_recv * (2/(L-1)) - 1; % 还原到 [-1,1] xr xr * xMax; % 反归一化 % 评估失真 SNR_rec snr(x(1:length(xr)), xr - x(1:length(xr))); sound(xr, Fs_audio); % 听一下这里psk_modulate和psk_demodulate是自定义函数封装第3节代码。逻辑说明reshape把串行比特流还原成每行nBitsbi2de是de2bi的反函数注意发射端如果做过压扩这里也要做逆压扩。snr函数计算重建信号相对噪声的清晰度。需要特别强调的是比特流长度必须能被nBits整除而且信道噪声引起的误码会使重建样本出现跳变用耳朵听会有“咔哒”声这就是误码在语音表现上的典型效果。4.3 误码率曲线与参数调整的4个方向先跑一组不同Eb/N0下的BER绘图与理论曲线对比EbN0_dB_list 0:2:10; berSim zeros(size(EbN0_dB_list)); for i 1:length(EbN0_dB_list) % 重复上面链路统计误码率 errs sum(bitStream ~ rxBits); berSim(i) errs / length(bitStream); end theoBER qfunc(sqrt(2*10.^(EbN0_dB_list/10))); semilogy(EbN0_dB_list, berSim, o, EbN0_dB_list, theoBER, -);运行后如果仿真BER偏离理论值优先检查以下四个参数。第一sps过小导致码间串扰特征是在高信噪比时BER曲线出现地板效应试调大到16或32。第二fc与采样率比例失调若Fs/fc小于4混频后2fc分量可能混进基带误差在低信噪比时会被淹没但高信噪比时曲线不收敛。第三beta过小比如0.1但定时采样点偏离峰值会造成ISI检查眼图是否张开。第四PCM解码时比特顺序错位即发射端和接收端reshape的行列顺序不一致此时BER约为0.5语音变成类似白噪声。用下面表格能快速定位问题原因现象可能参数问题调整方向BER高信噪比时地板sps过小或beta过小增大sps到16以上beta调到0.35左右频谱形状异常fc与Fs不匹配提高Fs使Fs≥4(fcB)解调符号旋转载波未同步检查接收端载波频率/相位语音完全不可懂比特顺序或reshape错误检查de2bi/bi2de的列顺序5. 验证仿真可信度的3个技巧星座图、眼图和频谱搬移检查5.1 用星座图和眼图快速看调制质量在解调器匹配滤波后、判决前把rxSymbols画成散点图就能得到星座图。BPSK理想状态是两个点QPSK是四个。当噪声增大时散点围绕理想点扩散如果出现整体旋转则说明接收载波不同步。眼图则是把连续基带波形按符号周期分段叠加绘制观察眼睛张开程度。在信号加噪前后分别运行下面代码eyediagram(rxBaseband(1:length(rxBaseband)), sps); constellation rxSymbols(:); scatterplot(constellation);眼睛张开越小说明ISI和噪声越严重。若眼图中心模糊优先调整定时采样点位置。5.2 用频谱图验证“频带传输”的真实搬移频带传输是否成功的直接证据是发射信号频谱以载波频率为中心。用periodogram绘制功率谱密度[psd, f] periodogram(txPassband, hann(length(txPassband)), [], Fs); plot(f/1e3, 10*log10(psd)); xlabel(频率 (kHz)); ylabel(PSD (dB));在f200kHz处应有明显峰值两侧展宽约43kHz。如果峰值出现在基带处说明txBaseband直接输出了而不是经过载波调制如果频谱左右不对称多半是滤波器系数不对称或调制过程维度不一致。5.3 把单个仿真封装成函数批量扫参避免“仿真发散”很多人在调整参数时反复改全局变量导致结果不稳定。建议把从PCM编码到解调的整条链路封装成一个函数ber run_psk_pcm_sim(EbN0_dB, sps, fc, beta, nBits)。这样批量扫参时只需循环调用。如果仿真数值发散先检查滤波器是否因阶数过长产生NaN局部采样点是否过少以及接收端载波频率和发射端是否严格一致。仿真发散不是随机事件通常是采样率低于2倍载波频率或滤波器不稳定导致的。这些验证手段加在一起就能确认你搭的不是“看起来像通信”的黑盒而是一套可解释、可复现的语音信号频带传输仿真链路。本文还有配套的精品资源点击获取 MATLAB语音信号PCM编码与PSK调制频带传输仿真指南
简介这套资料是一份关于语音信号频带传输通信系统仿真的课程设计文档基于PCM编码与PSK调制技术使用MATLAB/Simulink平台完成。内容覆盖PCM采样、量化与编码PSK调制解调高斯白噪声信道建模以及完整收发端系统流程适合通信工程、电子信息类专业学生用于课程设计、仿真实践或毕业设计参考。文档共1个docx文件压缩包大小约3.29MB包含26页设计报告从引言、设计原理到模块参数设置与仿真结果分析均有展开结构完整。资源已有359人学习下载说明其对同类课题具参考价值。读者可直接借助其中的系统框图、模块配置说明和结果对比思路快速搭建自己的语音频带传输仿真模型减少调试摸索时间。1. 从PCM编码到PSK调制这套频带传输仿真到底在做什么你在设计一个语音对讲模块时话筒送出的模拟信号并不能直接发射必须先经过PCM编码转成比特流再通过PSK调制把比特流搬到合适的载波频段上这就是“语音信号频带传输”的基本路径。仿真这件事不是为了把流程图跑一遍而是要把每一级信号的波形、频谱、误码率都量化出来回答几个实际问题采样率取多少量化位数够不够载波频率会不会和符号率打架信道噪声大到什么程度还能听懂语音。本文按“先数字化、再调制、最后跑通全链路”的顺序用MATLAB把这条链路从零搭起来并给出参数调整和验证方法。适合正在做通信课程设计、以及想从“看公式”转向“看波形”的从业者。本文不依赖任何硬件平台所有结果都能在普通电脑上复现。2. PCM编码仿真把模拟语音变成可传输的比特流2.1 语音数字化中的采样、量化、编码三步拆解语音信号是带限模拟信号典型电话语音频带是300Hz到3.4kHz。根据奈奎斯特定理采样率至少是最高频率的两倍实际通常取8kHz。采样只完成了时间离散化幅度仍然是连续的需要量化到有限个电平。PCM的核心就是“采样-量化-编码”三步采样确定时间点量化把连续幅度映射到2^n个离散电平编码把每个电平写成n位二进制码字。这里n通常取8再加上8kHz采样率得到64kbps的速率也就是标准电话PCM的速率。量化误差是PCM最直接的失真源。均匀量化在信号幅度分布不均匀时浪费比特位因此实际系统常采用A律或μ律压扩对小信号用更短的量化间隔。在仿真中我们既可以先忽略压扩用均匀量化验证链路也可以用MATLAB的compand函数模拟压扩效果。量化位数n直接决定信噪比每增加1位量化信噪比约提升6dB——这是选型时判断“8位还是12位”的基本依据。真正做语音仿真时还要注意语音信号的峰值因子较大归一化方式会影响量化精度一般用max(abs(x))做归一化而不是用均方根值。2.2 用MATLAB实现PCM编码从采样到输出比特流的最小代码下面这段代码生成一段测试音完成采样、量化和二进制编码。它不读取外部语音文件便于你首先确认算法本身是否正确。Fs 8000; % 采样率 8kHz t 0:1/Fs:0.08; % 80ms 时长 x sin(2*pi*800*t) 0.5*sin(2*pi*2000*t); % 双音信号 nBits 8; % 量化位数 L 2^nBits; % 量化电平数 xMax max(abs(x)); % 归一化 xNorm x / xMax; % 归一化到 [-1, 1] % 均匀量化把 [-1,1] 映射到 [0, L-1] xq round((xNorm 1) * (L-1) / 2); xq max(min(xq, L-1), 0); % 防止舍入越界 % 十进制量化电平转二进制比特流 bits de2bi(xq, nBits, left-msb); % 每行一个采样点左对齐 bitStream reshape(bits., 1, []); % 按时间顺序串行输出逻辑说明先归一化信号幅度保证量化器满量程利用xq是量化索引范围0到L-1de2bi函数把每个索引转换为nBits位二进制向量left-msb表示最高有效位在前最后reshape成一行连续的比特流。整个过程保持了采样点顺序后续PSK调制直接把bitStream作为输入。参数说明Fs8000对应窄带语音下限若信号包含更高频率必须按实际最高频率提高采样率nBits8是常规语音质量追求更高保真可换成16但输出比特率翻倍这里使用均匀量化若改为μ律压扩只需在量化前对xNorm用compand处理。2.3 PCM仿真落地的3个必调参数第一是采样率。仿真中看到频谱混叠首先检查采样率是否大于最高频率的2倍若信号是通用的4kHz音频8kHz刚好够但留一点余量更稳比如取16kHz。第二是量化位数。位数的效果要结合回放或SNR来看建议分别跑8、12、16位对比误码率或量化信噪比你会看到每多1位约6dB的改善。第三是压扩选择。A律和μ律的差异主要在低幅度区间的信噪比如果测试信号中包含极小声段使用μ律/律能让重建波形在示波器上更贴近原信号。下表给出了PCM仿真中需要权衡的核心参数作为初始配置参考参数典型取值影响对象采样率 Fs8kHz / 16kHz最高可恢复频率决定频谱占用量化位数 nBits8 / 12 / 16输出比特率量化信噪比约6n dB压扩律无、A律、μ律小信号量化信噪比语音主观质量这里要注意PCM仿真的“坑”常出在量化边界当信号幅度刚好等于最大值时round可能产生L越界后必须clip。上段的max(min())就是干这个的。另一个坑是de2bi的列顺序默认最低有效位在前不加left-msb会导致比特顺序错乱解调后语音完全不可懂。3. PSK调制仿真把比特流搬上载波实现频带传输3.1 为什么频带传输要选PSK而不是ASK或FSKPCM输出的比特流是基带信号频率集中在低频段无法直接通过带通型信道比如无线频道、电缆频分复用传输。频带传输的核心是选择一个载波频率把比特信息调制到载波的某个参数上。ASK调制幅度FSK调制频率PSK调制相位。PSK在相同信噪比下误码率最低功率效率高但相位模糊是它的固有缺陷。在仿真中我们常先做BPSK每符号1比特验证链路后再扩展QPSK每符号2比特以提升频谱效率。选择PSK也是因为它很适合理想信道仿真接收端用相干解调即可达到理论误码率不需要复杂的信道估计就能对照公式验证结果。调制过程的等效含义是把一串比特映射成一串复平面上带幅相的点再通过脉冲成型滤波器限制带宽最后乘以载波完成频谱搬移。仿真中“搬移”的体现是发射信号频谱在载波频率两侧对称出现直接用periodogram看即可。3.2 BPSK调制仿真从比特流到已调信号的最小实现首先把前面PCM产生的bitStream当输入。BPSK将比特0/1映射为±1每符号只携带1比特。为了模拟真实系统需要过采样和脉冲成型再与载波相乘。% 系统参数 sps 16; % 每个符号的采样点数 Rb 64000; % 比特率 64kbps来自 8kHz * 8bit Rsym Rb; % BPSK符号率 比特率 Fs Rsym * sps; % 仿真采样率 1.024MHz % 比特流 - 符号 symbols 2 * bitStream - 1; % 0 - -1, 1 - 1 % 上采样 根升余弦脉冲成型 symbolsUp upsample(symbols, sps); beta 0.35; rrc rcosdesign(beta, 6, sps, sqrt); txBaseband filter(rrc, 1, symbolsUp); % 载波调制 fc 200e3; % 载波频率 200kHz t (0:length(txBaseband)-1) / Fs; carrier cos(2*pi*fc*t); txPassband txBaseband .* carrier;逻辑说明upsample在每两个符号间插入sps-1个零相当于升采样rcosdesign生成根升余弦滤波器这里滚降系数0.35滤波器长度为6个符号周期filter完成成型滤波得到连续基带波形。最后乘以余弦载波频谱被搬到fc两侧。这里的仿真采样率Fs必须大于2*(fc信号带宽)否则会产生混叠上例中fc200kHz信号带宽约(1beta)*Rsym/2≈43kHzFs1.024MHz满足要求。参数说明sps影响时域波形的平滑度和仿真负载一般取8~32过小则波形不平滑过大则运行慢。beta决定频带宽度取值0.15~0.5越小越省带宽但脉冲的时域衰减越慢对定时误差越敏感。fc的选择要兼顾仿真采样率太大会导致Fs成倍增大太小则基带和频带频谱容易重叠。输出txPassband就是实际的频带传输信号。下表汇总了PSK调制仿真中的关键参数及其选择逻辑方便你在批量实验时快速定位参数常见范围对系统的影响sps8 ~ 32精度与仿真速度平衡过小区间干扰大beta0.15 ~ 0.5带宽占用与定时鲁棒性fc与Fs满足Fs2(fcB)频谱位置过高增加仿真成本3.3 接收端相干解调本地载波、低通滤波与匹配滤波相干解调要求接收端产生与发射端同频同相的本地载波。仿真中最简单的做法是直接发射载波的副本等于假设理想载波同步。用本地载波乘以接收信号产生基带和2倍载频分量低通滤除高频分量再通过匹配滤波器与发射端相同的根升余弦滤波器并下采样判决。% 假设 rxPassband 是接收信号发射端参数复用 rxBaseband rxPassband .* cos(2*pi*fc*t); % 混频到基带 % 低通滤波用发射端成型滤波器匹配滤波这里只做一次 rxFiltered filter(rrc, 1, rxBaseband); % 下采样在每个符号中心采样 samplingPoints 1 sps/2 : sps : length(rxFiltered); rxSymbols rxFiltered(samplingPoints); % 判决BPSK 0 判为1否则0 rxBits rxSymbols 0;这段代码里混频后没有单独低通滤波器而是直接用了匹配滤波器。由于根升余弦滤波器的频响能压制混频产生的2fc分量通常在仿真中省去一级额外低通前提是采样率足够高、fc远离基带。samplingPoints从sps/2开始取是为了对准脉冲成型的峰值点。如果定时偏移眼图会闭合。需要提醒的是真正系统存在相位偏差、频偏和定时误差。仿真时可以人为给载波加一个相移观察解调星座图的旋转这是后面升级到非理想同步的切入点。4. 全链路仿真加入信道噪声测误码率并调整参数4.1 给频带信号加高斯白噪声用SNR反推噪声功率信道仿真的基本模型是加性高斯白噪声AWGN。给频带信号加噪需要先计算信号平均功率再根据目标信噪比折算噪声功率。这里常用两种方式直接以采样点为单位计算snr或用通信系统常用的Eb/N0。误码率曲线习惯横轴用Eb/N0每比特能量与噪声功率谱密度之比BPSK的理论误码率为Q(sqrt(2*Eb/N0))。代码中可以用awgn函数一步到位EbN0_dB 6; % 目标 Eb/N0单位dB snr_dB EbN0_dB 10*log10(2); % BPSK 时 Es2Eb符号SNR % 注意awgn 默认 snr 按信号功率与噪声功率比 rxPassband awgn(txPassband, snr_dB, measured);awgn的measured表示它以输入信号的实际功率为基准计算噪声。snr_dB换算关系对BPSK每个符号有两个比特能量因此 Es/N0 Eb/N0 10log10(2) dB。若做QPSK每个符号2比特Es/N0 Eb/N0 10log10(2) 同样成立因为每符号能量是2Eb注意区分。如果不想用awgn也可以手动生成噪声noise sqrt(noisePower)*randn(size(txPassband))其中noisePower signalPower / 10^(snr_dB/10)。手动方式便于控制噪声种子和功率适合批量仿真。4.2 搭建完整收发链路从原始语音到重建语音的闭环现在把第2章PCM、第3章调制解调和第4章加噪串成一个脚本。流程是读取/生成语音 → PCM编码 → BPSK调制 → 过信道 → 相干解调 → PCM解码 → 播放/对比。下面给出主流程的关键片段% ---- 发射端 ---- % 假设 x 是归一化语音信号长度 N [x, Fs_audio] audioread(test.wav); % 也可用自己合成的信号 x x(:,1)./max(abs(x)); % 单声道归一化 % PCM编码得到 bitStream调用前面函数 % 注意 bitStream 长度为 8*length(x)因为 nBits8 % ---- 调制 ---- [txPassband, params] psk_modulate(bitStream, sps, fc, Fs); % ---- 信道 ---- rxPassband awgn(txPassband, snr_dB, measured); % ---- 解调 ---- rxBits psk_demodulate(rxPassband, params); % ---- PCM解码 ---- % 比特流转回量化索引 rxBitsMat reshape(rxBits, nBits, []).; xq_recv bi2de(rxBitsMat, left-msb); xr xq_recv * (2/(L-1)) - 1; % 还原到 [-1,1] xr xr * xMax; % 反归一化 % 评估失真 SNR_rec snr(x(1:length(xr)), xr - x(1:length(xr))); sound(xr, Fs_audio); % 听一下这里psk_modulate和psk_demodulate是自定义函数封装第3节代码。逻辑说明reshape把串行比特流还原成每行nBitsbi2de是de2bi的反函数注意发射端如果做过压扩这里也要做逆压扩。snr函数计算重建信号相对噪声的清晰度。需要特别强调的是比特流长度必须能被nBits整除而且信道噪声引起的误码会使重建样本出现跳变用耳朵听会有“咔哒”声这就是误码在语音表现上的典型效果。4.3 误码率曲线与参数调整的4个方向先跑一组不同Eb/N0下的BER绘图与理论曲线对比EbN0_dB_list 0:2:10; berSim zeros(size(EbN0_dB_list)); for i 1:length(EbN0_dB_list) % 重复上面链路统计误码率 errs sum(bitStream ~ rxBits); berSim(i) errs / length(bitStream); end theoBER qfunc(sqrt(2*10.^(EbN0_dB_list/10))); semilogy(EbN0_dB_list, berSim, o, EbN0_dB_list, theoBER, -);运行后如果仿真BER偏离理论值优先检查以下四个参数。第一sps过小导致码间串扰特征是在高信噪比时BER曲线出现地板效应试调大到16或32。第二fc与采样率比例失调若Fs/fc小于4混频后2fc分量可能混进基带误差在低信噪比时会被淹没但高信噪比时曲线不收敛。第三beta过小比如0.1但定时采样点偏离峰值会造成ISI检查眼图是否张开。第四PCM解码时比特顺序错位即发射端和接收端reshape的行列顺序不一致此时BER约为0.5语音变成类似白噪声。用下面表格能快速定位问题原因现象可能参数问题调整方向BER高信噪比时地板sps过小或beta过小增大sps到16以上beta调到0.35左右频谱形状异常fc与Fs不匹配提高Fs使Fs≥4(fcB)解调符号旋转载波未同步检查接收端载波频率/相位语音完全不可懂比特顺序或reshape错误检查de2bi/bi2de的列顺序5. 验证仿真可信度的3个技巧星座图、眼图和频谱搬移检查5.1 用星座图和眼图快速看调制质量在解调器匹配滤波后、判决前把rxSymbols画成散点图就能得到星座图。BPSK理想状态是两个点QPSK是四个。当噪声增大时散点围绕理想点扩散如果出现整体旋转则说明接收载波不同步。眼图则是把连续基带波形按符号周期分段叠加绘制观察眼睛张开程度。在信号加噪前后分别运行下面代码eyediagram(rxBaseband(1:length(rxBaseband)), sps); constellation rxSymbols(:); scatterplot(constellation);眼睛张开越小说明ISI和噪声越严重。若眼图中心模糊优先调整定时采样点位置。5.2 用频谱图验证“频带传输”的真实搬移频带传输是否成功的直接证据是发射信号频谱以载波频率为中心。用periodogram绘制功率谱密度[psd, f] periodogram(txPassband, hann(length(txPassband)), [], Fs); plot(f/1e3, 10*log10(psd)); xlabel(频率 (kHz)); ylabel(PSD (dB));在f200kHz处应有明显峰值两侧展宽约43kHz。如果峰值出现在基带处说明txBaseband直接输出了而不是经过载波调制如果频谱左右不对称多半是滤波器系数不对称或调制过程维度不一致。5.3 把单个仿真封装成函数批量扫参避免“仿真发散”很多人在调整参数时反复改全局变量导致结果不稳定。建议把从PCM编码到解调的整条链路封装成一个函数ber run_psk_pcm_sim(EbN0_dB, sps, fc, beta, nBits)。这样批量扫参时只需循环调用。如果仿真数值发散先检查滤波器是否因阶数过长产生NaN局部采样点是否过少以及接收端载波频率和发射端是否严格一致。仿真发散不是随机事件通常是采样率低于2倍载波频率或滤波器不稳定导致的。这些验证手段加在一起就能确认你搭的不是“看起来像通信”的黑盒而是一套可解释、可复现的语音信号频带传输仿真链路。本文还有配套的精品资源点击获取