简介这份资源面向通信原理学习者与C语言开发者聚焦BPSK调制解调技术的完整实现帮助读者在AWGN信道下理解数字通信链路并评估误码率性能。压缩包共8个文件以6个cpp源文件为主另有1个c文件与1个bgi图形驱动文件整体约8KB涵盖调制解调核心算法、噪声生成、误差函数计算及数据签名等辅助模块结构紧凑便于逐文件研读。已有726人学习下载说明其在同类基础实验中具有一定参考价值。读者可从中获得BPSK相位映射、高斯白噪声模拟、误码率统计等关键环节的C语言实现思路并借助误差函数与图形驱动模块完成性能曲线绘制适合作为课程设计、通信仿真入门或算法移植的参考素材。1. 从零手搓 BPSK为什么我劝你先用 C 语言把调制器写一遍如果你正在做数字通信的课设、准备电赛的数据链模块或者单纯想搞明白“一个比特到底是怎么变成一段正弦波”的那这份基于 C 语言的 BPSK 实现资源大概率能帮你省下不少翻教材的时间。BPSK二进制相移键控是所有数字调制里最朴素的一种一个比特对应一个相位0 和 1 分别映射到 0 度和 180 度。听起来简单到不像话但真让你用 C 语言从比特流一路写到采样点输出中间关于相位连续性、符号速率、载波频率、采样率这几件事怎么摆平才是真正卡人的地方。这份资源给的不是一个调库调参的 Demo而是一套能编译、能跑、能拿示波器或者音频接口看波形的完整链路适合刚学完《通信原理》但没动过手的学生也适合想快速验证解调算法、需要一个干净基线的工程师。我拿到之后第一件事就是把它拆开看它到底怎么处理相位跳变和定时关系——这两处要是含糊后面加噪声、加频偏全是玄学。2. 拆开 BPSK 调制器比特映射、载波生成与采样率怎么定2.1 比特到相位的映射逻辑与差分编码的取舍BPSK 最核心的一步就是把 0/1 映射成相位。常见做法有两种一种直接映射比特 0 对应相位 0比特 1 对应相位 π另一种先做差分编码用“相位是否变化”来携带信息。这份资源用的是直接映射代码里一个if就写完了但这里有个容易被忽略的坑直接映射在接收端需要知道绝对相位参考一旦信道引入固定相偏解调就会整体翻转。差分编码能绕开这个问题代价是误码会成对出现。我一般建议课设阶段先用直接映射把链路跑通等要加信道仿真了再换成差分改动量其实很小就在映射表那里动几行。/* 直接映射bit 为 0 输出 1为 1 输出 -1 */ float bpsk_map(int bit) { return (bit 0) ? 1.0f : -1.0f; } /* 差分编码版本用当前比特与前一比特的异或决定相位变化 */ int diff_encode(int bit, int *prev) { int out bit ^ (*prev); *prev bit; return out; }上面两个函数是整条链路的起点。bpsk_map返回的是归一化幅度后面乘上载波就是实际波形。diff_encode里*prev需要调用方维护初始值一般设 0。注意差分编码后仍然要过bpsk_map别直接把异或结果当相位用。2.2 载波频率、符号速率与采样率的三角关系这三个参数定不好后面全是血泪。符号速率 Rs 决定一个比特占多长时间载波频率 fc 决定正弦波多密采样率 fs 决定你一个符号能采多少个点。工程上一般要求 fs 至少是 fc 的 4 倍以上否则波形看起来像锯齿解调时积分判决会偏。这份资源里默认给了一组参数Rs 1000 bpsfc 10 kHzfs 100 kHz也就是一个符号 100 个采样点一个载波周期 10 个采样点。这个比例比较舒服既不会让数组太大也能看清相位跳变。#define SYMBOL_RATE 1000.0f /* 符号速率单位 bps */ #define CARRIER_FREQ 10000.0f /* 载波频率单位 Hz */ #define SAMPLE_RATE 100000.0f /* 采样率单位 Hz */ #define SAMPLES_PER_SYMBOL (int)(SAMPLE_RATE / SYMBOL_RATE) /* 生成一个符号内的 BPSK 波形phase 为当前符号对应的相位 */ void gen_symbol(float *buf, int symbol_val, float *phase) { float w 2.0f * M_PI * CARRIER_FREQ / SAMPLE_RATE; for (int i 0; i SAMPLES_PER_SYMBOL; i) { buf[i] symbol_val * cosf(*phase); *phase w; if (*phase 2.0f * M_PI) *phase - 2.0f * M_PI; } }gen_symbol里用*phase累加而不是每次重新算cos(2πfct)是为了保证符号之间相位连续。如果每个符号都从 0 开始算符号边界会出现人为的相位跳变频谱会鼓包。symbol_val就是bpsk_map的输出1 或 -1。SAMPLES_PER_SYMBOL必须是整数所以 fs 要选成 Rs 的整数倍否则会累积定时误差。2.3 用文件做比特源和波形输出的最小闭环资源里没有搞复杂的实时采集而是用文本文件读比特、用二进制文件写浮点采样点这样在 PC 上就能跑完整流程。输入文件一行一个 0 或 1输出文件是 float32 的 IQ 或者实信号。我一般会再加一个 Python 脚本把输出文件画出来确认眼图张开、相位跳变干净。/* 从文本文件读比特生成 BPSK 采样点写入二进制文件 */ int main(void) { FILE *fin fopen(bits.txt, r); FILE *fout fopen(bpsk.bin, wb); if (!fin || !fout) return -1; float phase 0.0f; float buf[SAMPLES_PER_SYMBOL]; int bit; while (fscanf(fin, %d, bit) 1) { int sym bpsk_map(bit); gen_symbol(buf, sym, phase); fwrite(buf, sizeof(float), SAMPLES_PER_SYMBOL, fout); } fclose(fin); fclose(fout); return 0; }这段main就是最小可跑闭环。bits.txt里随便写几十个 0/1编译后运行bpsk.bin就是基带实信号。用 Python 的numpy.fromfile读出来matplotlib一画符号边界和载波周期一目了然。注意fscanf的返回值判断别用while(!feof(fin))那个会多读一次。3. 解调端怎么接相干积分、判决门限与定时同步的简化处理3.1 相干解调的核心乘同频载波再积分调制是把比特搬上载波解调就是把它搬下来。BPSK 相干解调的标准动作是接收信号乘以同频同相的载波然后在一个符号周期内积分积分结果大于 0 判 0小于 0 判 1。这份资源里解调部分和调制共用同一组参数载波相位从 0 开始相当于理想同步。实际工程里载波恢复和定时恢复是两大块独立模块但课设阶段先把这个理想模型跑通再逐步加频偏和相偏思路会清晰很多。/* 对一个符号周期的采样点做相干积分并判决 */ int bpsk_demod(const float *buf, float *phase) { float w 2.0f * M_PI * CARRIER_FREQ / SAMPLE_RATE; float acc 0.0f; for (int i 0; i SAMPLES_PER_SYMBOL; i) { acc buf[i] * cosf(*phase); *phase w; if (*phase 2.0f * M_PI) *phase - 2.0f * M_PI; } return (acc 0.0f) ? 0 : 1; }acc就是积分结果这里用累加代替连续积分。phase的推进必须和调制端完全一致否则积分结果会随相偏余弦衰减。判决门限设 0 是因为 BPSK 两个星座点对称如果信道有直流偏置门限要相应调整。3.2 定时误差从哪来符号起始点对不齐会怎样上面解调代码默认每个符号的起始采样点就是数组第 0 个但实际接收时你根本不知道符号边界在哪。如果起始点偏了哪怕几个采样点积分区间就跨到了相邻符号判决结果会随机出错。这份资源里为了简化假设定时完美但我在复现时故意把接收数组偏移了 5 个点误码率直接从 0 跳到 0.3 以上。所以如果你要拿它做信道仿真定时同步这一环不能省常见做法是用早迟门或者平方环先估计符号边界再对齐积分区间。/* 简单演示接收数组偏移 offset 个采样点后再解调 */ int demod_with_offset(const float *rx, int offset, int nbits) { float phase 0.0f; int err 0; for (int k 0; k nbits; k) { int bit bpsk_demod(rx offset k * SAMPLES_PER_SYMBOL, phase); /* 这里假设已知发送比特实际用误码统计 */ if (bit ! tx_bits[k]) err; } return err; }offset就是人为引入的定时偏差。你可以写个循环从 0 扫到SAMPLES_PER_SYMBOL-1看误码数怎么随偏移变化一般会在符号边界附近出现一个平坦区那就是最佳采样时刻。这个扫描过程能帮你直观理解为什么定时同步模块要输出一个“误差信号”去调整采样相位。3.3 加高斯白噪声Eb/N0 怎么折算成采样点上的方差要验证解调器性能得加噪声。BPSK 在 AWGN 下的理论误码率是Q(sqrt(2Eb/N0))但代码里加的是采样点级别的高斯噪声方差要按采样率折算。这份资源里给了一个add_noise函数思路是先算符号能量再根据目标 Eb/N0 算噪声功率。我见过不少人直接把 Eb/N0 当 SNR 用结果曲线和理论对不上翻车就翻在这里。/* 向采样点添加高斯白噪声target_ebn0 单位 dB */ void add_noise(float *buf, int n, float target_ebn0) { float ebn0 powf(10.0f, target_ebn0 / 10.0f); float eb 1.0f; /* 每个符号能量幅度为 1 时 */ float n0 eb / ebn0; float sigma sqrtf(n0 * SAMPLE_RATE / (2.0f * SYMBOL_RATE)); for (int i 0; i n; i) { float u1 (float)rand() / RAND_MAX; float u2 (float)rand() / RAND_MAX; float g sqrtf(-2.0f * logf(u1)) * cosf(2.0f * M_PI * u2); buf[i] sigma * g; } }sigma的表达式是关键n0 * fs / (2 * Rs)再开根号。因为一个符号内有fs/Rs个采样点噪声功率要按带宽比例分配。u1不能为 0实际用的时候加个极小值保护。跑完加噪再解调统计误码和erfc算出来的理论值对比如果差得远先检查sigma对不对。4. 避坑与排查BPSK 从仿真到波形输出的五个常见翻车点4.1 波形看起来像调幅而不是调相现象用音频软件播放bpsk.bin听到的是忽大忽小的声音示波器上看包络起伏明显。原因符号映射后没有做相位连续累加每个符号都从cos(0)重新开始符号边界出现幅度跳变。解决像gen_symbol那样维护一个全局phase变量跨符号累加保证相位连续。如果用的是查表法表索引也要跨符号递增。4.2 解调误码率始终在 0.5 附近现象不管加不加噪声解出来的比特和发送比特一半对一半错。原因解调载波相位和调制载波相位差了 π导致积分结果整体反号。解决检查两边phase初始值是否一致或者把判决门限反过来。如果是差分编码检查prev初始值是否两端一致。4.3 采样率不是符号速率整数倍导致累积偏移现象跑长比特流时前面解调正确后面误码越来越多。原因SAMPLES_PER_SYMBOL用整数除法截断每个符号少采了零点几个点累积几千个符号后定时偏出一个符号。解决把SAMPLE_RATE改成SYMBOL_RATE的整数倍或者用浮点累加器记录每个符号的实际采样区间。4.4 噪声方差算错导致 Eb/N0 曲线整体平移现象误码率曲线形状对但和理论值差好几个 dB。原因把 Eb/N0 直接当 SNR 用忽略了采样率与符号速率的比值。解决按sigma sqrt(N0 * fs / (2 * Rs))折算其中N0 Eb / ebn0。验证方法不加噪声时误码为 0加噪声后曲线斜率应与理论一致。4.5 用feof控制读文件循环多读一个比特现象bits.txt里明明只有 100 个比特输出波形却有 101 个符号。原因while(!feof(fin))在最后一次读取后 feof 还没置位循环体多执行一次bit保持上一次的值。解决用while(fscanf(fin, %d, bit) 1)让读取函数的返回值控制循环。5. 进阶玩法把 BPSK 链路改成可配置的仿真框架跑通基础版之后我习惯把它改成一个带命令行参数的仿真框架这样换参数不用重新编译。用getopt解析-r符号速率、-f载波频率、-s采样率、-eEb/N0然后动态分配缓冲区。这样一套代码能同时跑课设演示和批量误码率扫描。/* 简易命令行参数解析编译后运行 ./bpsk -r 1000 -f 10000 -s 100000 -e 8 */ int main(int argc, char **argv) { float rs 1000.0f, fc 10000.0f, fs 100000.0f, ebn0 8.0f; int opt; while ((opt getopt(argc, argv, r:f:s:e:)) ! -1) { switch (opt) { case r: rs atof(optarg); break; case f: fc atof(optarg); break; case s: fs atof(optarg); break; case e: ebn0 atof(optarg); break; } } int sps (int)(fs / rs); /* 后续用 sps 替代宏动态分配 buf */ return 0; }改完之后sps不再是编译期常量gen_symbol和bpsk_demod都要改成接收sps参数。动态分配用malloc(sps * sizeof(float))记得在程序退出前free。这样你可以写个 shell 循环把 Eb/N0 从 0 扫到 12每个点跑一万个比特统计误码率最后用 Python 画出来和erfc理论曲线叠在一起。我每次改完参数都会先跑一个 Eb/N0 10 dB 的点看误码率是不是在 1e-3 量级如果差太多先查sigma和sps的整数性。从那以后我每次动采样率相关参数都强制走一遍这个单点验证省得后面批量跑完才发现曲线整体歪了。希望帮到你。本文还有配套的精品资源点击获取