简介面向通信原理课程实验的数字带通传输系统接收滤波器设计资料适合通信工程、电子信息类学生完成带通滤波器设计与图像传输验证。整套资源共15个文件包含8个Matlab脚本、4个txt数据文件、2个jpg和1个png图片素材压缩包仅543KB。其中Matlab脚本覆盖LP、BP滤波器设计、QPSK调制解调、信噪比分析及误码率统计等核心环节txt文件用于存储不同分辨率图像的传输数据图片素材则便于直观对比传输前后的视觉效果。另附实验代码说明文档可辅助理解各模块功能与参数设置。使用者能够对照完整链路学习接收滤波器特性对图像恢复质量的影响借助误码率分析量化系统性能还可以根据课程要求修改滤波器参数或调制方式重新仿真。已有368人学习下载适合作为通信原理实验参考或课程设计基础。1. 数字带通传输系统接收滤波器为什么硬件没动星座图还是歪的做通信原理课程设计时最容易掉进去的坑是以为接收滤波器只是“在解调前面随意加一个带通”。我拆过好几版这类设计包真正决定系统能不能在带外干扰和噪声下撑住的往往是这颗滤波器——中心频率、带宽、阶数、群延迟四个参数任何一个没对齐星座图和恢复图像都会翻车。这套数字带通传输系统接收滤波器设计资源里包含完整的滤波器系数设计脚本、图像传输链路仿真和误码率测试流程适合正在做通信课程设计、或者想把带通接收端理论落到仿真代码上的从业者拿来改参数就能跑出自己的结果。2. 滤波器选型与参数边界fir1 和 firpm 在接收端的真实差距2.1 设计指标先定这四个参数再谈算法接收端带通滤波器的作用不是把“高频”和“低频”分开而是在载波频点附近留出一道窗口让信号主瓣过去、把带外噪声和邻道干扰压掉。很多人上来就翻工具箱里的fir1我认为这是顺序反了应该先把四个参数定下来再决定用哪个函数。第一个是中心频率工程上要跟发送端载波频率严格对齐设计资源里通常把中心频率设为载波频率本身留出频偏余量第二个是通带带宽下限是符号速率决定的信号主瓣宽度上限受邻道约束第三个是阶数决定过渡带能做多陡代价是仿真变慢第四个是通带纹波和阻带衰减直接决定带内平坦度和对邻道的抑制深度。这四个参数之间的牵制关系可以用一张表概括。设计参数取小了取大了主要影响中心频率信号落在通带边缘带外干扰混入星座图旋转、误码率抬高通带带宽码间串扰、眼图闭合噪声功率上升BER 曲线整体劣化滤波器阶数过渡带平缓、选择性差仿真大幅变慢、边界效应明显实时性与精度权衡阻带衰减邻道信号压不掉阶数膨胀、群延迟变大带外选择性与资源开销我一般会在参数区里把这些写成可调变量而不是把数字直接写死在设计语句里否则后面换载波频率时容易漏改。资源里的脚本文件开头就是参数区通带边界是按载波频率的百分比算出来的这种做法在换频点测试时非常省事。2.2 窗函数法还是等波纹法选型逻辑与两段可直接套用的代码定完参数后设计函数的选择本质上是“过渡带、纹波、阶数”三者之间的交易。窗函数法实现简单、系数天然线性相位但过渡带宽度不可独立控制等波纹法也叫最优一致逼近法能在相同阶数下把阻带波动压得更均匀适合接收端这种对带外抑制要求高的场景。fs 200e3; % 采样率 200 kHz fc 40e3; % 载波频率 40 kHz Rs 10e3; % 符号速率 10 kHz % 边界频率通带 fx(fc-2k, fc2k)阻带压在 fc±5k 外 f_edge [0, 35e3, 38e3, 42e3, 45e3, fs/2] / (fs/2); a_mag [0, 0, 1, 1, 0, 0]; dev [0.001, 0.001, 0.01, 0.01, 0.001, 0.001]; % 等波纹设计阶数由 remezord 同族的 firpmord 估算 [n, f0, ao, w] firpmord(f_edge, [0 1 0], dev, fs); h_bpf firpm(n, f0, ao, w); freqz(h_bpf, 1, 1024, fs);这段代码里的f_edge是归一化边界频率数组注意firpmord要求边界频率严格递增而且通带两个边频要对称包住载波。dev数组对应每个频段的纹波值这里通带纹波设成 0.01阻带衰减对应 0.001换算成 dB 大约是 60dB 的带外抑制。设计完直接freqz看一眼幅频响应如果发现过渡带太宽就把通带边界往阻带边界方向收紧阶数会相应涨上去。对比窗函数法% 凯塞窗设计自动估算阶数 fc_n fc / (fs/2); bw_n 4e3 / (fs/2); % 通带宽度 ±2kHz win_bw 6e3 / (fs/2); % 过渡带宽度按 6kHz 估算 [n_ka, Wn, beta, ftype] kaiserord([fc_n-bw_n, fc_nbw_n, ... fc_n-win_bw, fc_nwin_bw], [0 1 0], [0.001 0.01 0.001], fs); h_win fir1(n_ka, Wn, ftype, kaiser(n_ka1, beta), noscale);两种方法在这个参数下跑出来的阶数有明显差异等波纹法通常给出三百多阶凯塞窗会略高一些。接收端选型我倾向于直接上等波纹带外抑制更均匀且通带边缘没有窗函数那种明显的拖尾。实际换到 60dB 阻带衰减时等波纹法的优势会进一步放大窗函数法的阶数已经涨到很难看的程度。2.3 群延迟和相位响应接收端最容易被忽略的指标滤波器设计包里很少有人主动看相位但这恰恰是接收端翻车的重灾区。带通滤波器会让不同频率分量产生不同的时延如果时延差太大符号波形会被拉畸变眼图开口变小抽样点偏移。FIR 滤波器只要系数对称群延迟就是常数(N-1)/2个采样点不会造成波形畸变IIR 滤波器虽然幅频响应做得很漂亮但群延迟在通带边缘会明显鼓包QPSK 信号解调后星座点发散得厉害。% 检查群延迟是否平坦 [gd, w] grpdelay(h_bpf, 1, 1024, fs); plot(w, gd); xlabel(Frequency (Hz)); ylabel(Group Delay (samples)); title(Filter Group Delay);这段代码用来打印群延迟曲线遇到设计结果时先看一眼通带范围内是否是一条水平直线。水平直线说明线性相位可以放心用如果有弯曲说明该滤波器不是线性相位解调端的符号采样点需要额外补偿。资源里所有接收滤波器都要求线性相位这也是为什么设计脚本里看不到任何butter或cheby1的原因。3. 图像传输链路搭建从调制到接收滤波的完整仿真脚本3.1 参数初始化采样率、载波频率和符号速率的匹配关系滤波器设计好后要把这颗滤波器真正塞进一条完整的数字带通传输链路里。这里最常见的问题是采样率乱选采样率过低滤波器过渡带归一化后占掉半个频带设计出来阶数极高采样率过高仿真数据量成倍增长跑一次 BER 曲线要等半天。% 链路基础参数 fs 200e3; % 采样率至少大于 2*(fcB/2) fc 40e3; % 载波频率 Rs 10e3; % 符号速率 sps fs / Rs; % 每符号采样点数这里等于 20 % 图像参数64x64 灰度图每个像素 8bit img_h 64; img_w 64; total_bits img_h * img_w * 8; % 仿真符号数按整帧图像折算 n_symbols ceil(total_bits / 1); % BPSK1bit 对应 1 个符号三个频率参数要满足两个约束载波频率不能超过采样率的一半符号速率对应的信号带宽不能和载波频率叠在一起。这里的组合是 40kHz 载波、10kHz 符号速率信号双边带占用 30~50kHz采样率 200kHz 留了四倍余量滤波器设计时归一化频率不会被挤到边缘。sps20意味着每个符号采样 20 个点这个值不是随便拍的。太小时滤波器对符号边缘的“刻画”精度不够太大时上方kron产生的时间和内存开销偏大。常见的做法是把 sps 控制在 10 到 40 之间资源里默认 20属于仿真速度和精度的折中值。3.2 发送端与信道图像怎么变成带通信号发送端要完成的任务是灰度图像 → 比特流 → 符号序列 → 连续时间波形。整个过程分三步先做比特映射再做符号映射最后上变频成带通信号。% 读图并转成比特流 img imread(gray_image.png); % 64x64 灰度图 img_vec double(img(:)); bits de2bi(img_vec, 8, left-msb); % 每个像素 8bit bits bits(:); % BPSK 调制 symbols pskmod(bits, 2); % 零阶保持把符号序列扩展成连续波形 tx_wave kron(symbols, ones(sps, 1)); % 上变频复数载波再取实部避免双边带重叠 t (0:length(tx_wave)-1). / fs; carrier exp(1j * 2 * pi * fc * t); tx real(tx_wave .* carrier); % 加噪声信道 snr_db 12; rx awgn(tx, snr_db, measured);de2bi的第三个参数left-msb保证每个字节的最高位在前否则接收端恢复出来的图像会有像素错位。pskmod(bits, 2)做 BPSK 调制映射结果是 1 和 -1 的实数序列。上变频时用复数载波的实部这样自然产生上下两个边带如果直接乘cos波形效果一样但复数写法在接收端做正交解调时更容易理解。信道这里先用awgn加白噪声measured表示按照发送信号实际功率计算信噪比避免因信号功率归一化问题导致信噪比虚高。这个参数在比较不同滤波器性能时很重要固定同样的snr_db才能对比出滤波器的真实增益差异。3.3 接收端滤波与正交解调滤波器放在哪一步很关键接收端链路是整套资源的核心滤波器和解调的衔接顺序决定能否正确恢复出图像。很多初学者把带通滤波器放在最前面这没问题但滤波完之后直接抽样判决那就是灾难——混频后的高频分量还在星座点全部带毛刺。% 接收端先过带通滤波器滤除带外噪声 filtered filter(h_bpf, 1, rx); % 下变频到基带得到 I/Q 两路 baseband filtered .* exp(-1j * 2 * pi * fc * t); % 基带低通去掉混频产生的 2*fc 高频分量 lp_coeff ones(sps, 1) / sps; % 移动平均滤波器窗口符号宽度 I_base real(baseband); Q_base imag(baseband); I_lp conv(I_base, lp_coeff, same); Q_lp conv(Q_base, lp_coeff, same); % 抽样判决每 sps 个点取一个符号 idx (n_filter_taps/2 1):sps:length(I_lp); I_sym I_lp(idx); Q_sym Q_lp(idx); rx_symbols sign(I_sym); % BPSK 只取 I 路这里有两处容易写错第一filter(h_bpf,1,rx)是零初始状态的卷积滤波输出序列开头有(N-1)/2个点的瞬态抽样起点要跳过这段第二下变频后的baseband里除了基带信号还有位于2*fc附近的高频镜像必须用低通滤波器去掉。资源里用长度为sps的移动平均作为简易低通原理上等效于对每个符号宽度内的能量求平均。第二路 Q 路信号在 BPSK 里不参与判决但仍然要完成同样的低通和抽样因为后续如果改成 QPSK 或 QAM两路都要用。这套代码的结构在资源里是独立函数换调制方式时只需要改pskmod的参数和最后的判决规则。3.4 结果评估误码率统计与恢复图像保存链路搭完后不能只看恢复图像像不像还要量化误码率。通信原理课程设计的评分也往往落在误码率对比曲线上所以资源里专门提供了一个统计函数。% 对齐发送符号和接收符号的长度 min_len min(length(symbols), length(rx_symbols)); tx_ref symbols(1:min_len); rx_ref rx_symbols(1:min_len); % 统计误码率BPSK 星座点为 1/-1 err sum(tx_ref ~ rx_ref); ber err / min_len; % 恢复图像 rx_bits double(rx_ref 0); rx_bits reshape(rx_bits, [], 8); img_rec bi2de(rx_bits, left-msb); img_rec reshape(img_rec, img_h, img_w); imwrite(uint8(img_rec), recovered_image.png);误码率统计前必须做长度对齐滤波器群延迟会让接收符号流比发送流短若干点错误统计不对齐的话第一帧的错位会产生连锁误码。资源里统计时强制min_len截断确保对比的是真正对齐的符号序列。恢复图像的reshape顺序要和发送端完全一致否则画面会出现横竖条纹。我在实际跑这套流程时发现误码率不是唯一指标恢复图像的视觉效果有时更直观。信噪比低于某个阈值后误码率会先涨图像上表现为一块一块的椒盐噪声而滤波器带宽过窄时图像不是噪点而是整体模糊两类现象要对着误码率一起看才能判断问题出在滤波器还是信道。4. 接收滤波器集成时的四个隐蔽坑现象、原因与解决4.1 参数错配类的坑中心频率偏移与带宽过窄数字带通滤波器集成到链路里后最容易出现的问题反而不是滤波器本身而是参数与链路不匹配。现象接收端频谱看起来正常星座图却一直在缓慢旋转误码率卡在某个值下不去怎么调信噪比都没用。 原因滤波器通带中心频率和实际载波之间存在固定频偏。设计中直接用了fc作为滤波器中心频率但发送端的上变频频率或接收端的下变频频率稍微有偏差通带边缘就把信号主瓣削掉了一部分。BPSK 信号被非对称滤波后星座点轨迹变成椭圆旋转是残余频偏叠加的结果。 解决先用fft估计接收信号的实际频谱峰值再把峰值频率替换到滤波器设计的fc位置或者在解调后加一个锁相环跟踪残余频偏。课程设计里最省事的办法是确认发送端乘的载波和接收端滤波器中心频率用的是同一个变量而不是各自写了近似数字。现象加了带通滤波器后频谱上白噪声确实压下去了但眼图反而比不过滤还差。 原因滤波器带宽收得太窄把信号主瓣的过渡能量也滤掉了。很多人设计滤波器时只想着“滤噪声”把通带边界压到符号速率附近实际通带至少要覆盖信号主瓣加过渡带通常取带宽大于等于符号速率的一到两倍。 解决回到第二章的参数表把带宽从Rs放宽到1.5*Rs甚至2*Rs重新设计滤波器。带宽和噪声的权衡不该靠压带宽解决而应该靠匹配滤波和均衡器解决。4.2 实现时序类的坑群延迟、边界瞬态与镜像边带现象接收端抽样点怎么移都有固定偏移误码率统计总是偏差。 原因线性相位 FIR 滤波器会给信号引入(N-1)/2个采样点的群延迟。这个延迟在高阶滤波器下可能达到一两百个采样点直接按1:sps:end抽样相当于整体提前判决。 解决抽样索引要从filter输出后的第(N-1)/2 1个点开始。资源里把这个偏移量做成了自动计算换滤波器阶数时不需要手动改。现象滤波器阶数很高时仿真一跑就是几分钟而且首尾几百个符号的错误率明显高于中间段。 原因高保真设计下滤波器阶数很容易冲到三百到五百阶filter卷积计算量随阶数线性增长同时滤波起始阶段的瞬态响应没有完全建立开头的若干个符号在错误环境下会被拉偏。 解决仿真数据长度增加 padding例如在符号序列前后各补一段循环前缀或在滤波后再截掉瞬态段。课程设计里没必要上太高阶数阻带衰减做到 50~60dB 就够用了再往上属于过度设计拿换阶数的时间调均衡器更划算。现象恢复图像出现重复的横条纹误码率却不高。 原因发送端把基带信号直接用实载波上变频时产生了以载波为对称轴的镜像边带接收端下变频后两条边带叠在一起符号间形成干扰。这是发送端结构造成的接收滤波器再努力也无法完全分开。 解决发送端使用复数载波exp(1j*2*pi*fc*t)后取实部这是希尔伯特变换意义上的单边带构造能保留所需边带结构。这类问题仿真代码里非常隐蔽因为频谱图上看起来只是正常的双边带恢复图像的条纹才是最明显的信号。这类坑在做图像传输仿真时几乎是必然遇到的我每次拆到类似资源都会先把这几条检查一遍再去动误码率曲线。5. 用 BER 曲线验证接收滤波器同一条链路跑完整个测试5.1 最小验收套件的组成整套资源验证到最后需要形成一套可判定的输出。我的习惯是固定五个检查项幅频响应曲线、群延迟曲线、无噪声星座图、BER 随 SNR 变化曲线、恢复图像。这五项全部通过滤波器和链路才算真正适配。幅频响应检查滤波器通带位置和阻带衰减群延迟曲线检查线性相位无噪声星座图最能暴露群延迟补偿是否到位——无噪声条件下星座点如果还是散的说明抽样点没对准。5.2 一套可复用的 BER 扫描脚本% 复用同一组滤波器系数跑 BER避免每次重新设计 snr_list 0:2:16; ber_list zeros(size(snr_list)); for k 1:length(snr_list) rng(42); % 固定种子确保对比公平 rx_cur awgn(tx, snr_list(k), measured); filtered filter(h_bpf, 1, rx_cur); baseband filtered .* exp(-1j * 2 * pi * fc * t); I_lp conv(real(baseband), lp_coeff, same); idx (n_filter_taps/2 1):sps:length(I_lp); rx_sym sign(I_lp(idx)); cur_err sum(tx_ref ~ rx_sym(1:length(tx_ref))); ber_list(k) cur_err / length(tx_ref); end semilogy(snr_list, ber_list, o-); hold on; semilogy(snr_list, qfunc(sqrt(2*10.^(snr_list/10))), --);跑 BER 扫描时有两个细节第一滤波系数必须在循环外设计好循环内只复用否则每个信噪比点的滤波器都不同曲线无法反映真实趋势第二rng固定随机种子确保每个信噪比点使用的是相同的噪声序列否则低误码率区域的抖动会被误认为规律性波动。扫描结果和理论曲线对比时BPSK 理论误码率是qfunc(sqrt(2*snr))snr 需要从 dB 转成线性值。如果实测曲线和理论曲线平行但整体偏移通常是滤波器引入的噪声带宽大于信号带宽如果曲线形状都不一样大概率是前面第四章里的某一个坑。从那以后我每次拿到这类接收滤波器设计资源都会强制自己先跑一遍幅频响应和群延迟再跑无噪声星座图最后才上 BER 扫描。这三步帮我避开了绝大多数翻车现场希望帮到你。本文还有配套的精品资源点击获取
