作者：绳匠_ZZ0
从C语言数据到Matlab曲线，我终于看到了LDPC译码的"性能悬崖"！
一、前言：为什么要把数据可视化？
SNR(dB) SPA_BER MinSum_BER
0.0 1.050e-01 1.186e-01
1.0 5.200e-02 6.250e-02
2.0 1.750e-02 2.330e-02
3.0 3.200e-03 5.400e-03
4.0 4.500e-04 8.500e-04
5.0 5.000e-05 1.000e-04数字不够直观！SPA到底比Min-Sum强多少？关键增益区间在哪？用Matlab画出BER曲线，一切豁然开朗！
二、数据准备：C程序输出到文本
2.1 修改C程序输出
在仿真主函数中，将printf改为fprintf输出到文件：
c
FILE *fp = fopen("ber_results.txt", "w");
fprintf(fp, "SNR_dB\tSPA_BER\tMinSum_BER\n");
for (int snr_idx = 0; snr_idx < num_snr; snr_idx++) {
// 获取ber_spa, ber_min
fprintf(fp, "%.1f\t%e\t%e\n", snr_db, ber_spa, ber_min);
}
fclose(fp);生成文件格式：
SNR_dB SPA_BER MinSum_BER
0.0 1.050e-01 1.186e-01
1.0 5.200e-02 6.250e-02
...三、Matlab绘图实战
3.1 读取数据
matlab
data = readtable('ber_results.txt', 'FileType', 'text', 'Delimiter', '\t');
snr = data.SNR_dB;
ber_spa = data.SPA_BER;
ber_min = data.MinSum_BER;3.2 绘制半对数曲线
matlab
figure('Position', [100, 100, 800, 600]);
semilogy(snr, ber_spa, 'b-o', 'LineWidth', 2, 'MarkerSize', 8, 'MarkerFaceColor', 'b');
hold on;
semilogy(snr, ber_min, 'r-s', 'LineWidth', 2, 'MarkerSize', 8, 'MarkerFaceColor', 'r');
grid on;
xlabel('SNR (dB)');
ylabel('BER');
title('LDPC译码性能对比：SPA vs Min-Sum');
legend('SPA', 'Min-Sum', 'Location', 'southwest');
xlim([0, 8]); ylim([1e-6, 1e-0]);
text(3.5, 1e-3, 'SPA优于Min-Sum约0.5~1dB', 'FontSize', 10);
saveas(gcf, 'ber_comparison.png');效果示意图（纵轴对数刻度）：
BER
10^0 |●-------●------------------ Min-Sum
| ●
10^-1| ●
| ●
10^-2| ●
| ●
10^-3| ●
| ●
10^-4| ●----------- SPA
| ●
10^-5| ●
+----------------------------→ SNR(dB)
0 1 2 3 4 5 6 7四、关键发现
- 全区间优势：SPA曲线始终低于Min-Sum
- 最大增益区间：在BER=10⁻³时，SPA仅需3.2dB，Min-Sum需3.8dB（增益0.6dB）
- 低信噪比差距扩大：SNR=2dB时，SPA的BER低25%
- 瀑布效应：3-5dB区间出现典型LDPC"性能悬崖"，SPA的下降更陡峭
五、无Matlab解决方案
5.1 Python+matplotlib
python
import matplotlib.pyplot as plt
plt.semilogy(snr, ber_spa, 'b-o', label='SPA')
plt.semilogy(snr, ber_min, 'r-s', label='Min-Sum')
plt.grid(True); plt.xlabel('SNR(dB)'); plt.ylabel('BER')
plt.legend(); plt.savefig('ber_comparison.png')5.2 Excel方案
- 数据粘贴到Excel
- 插入→散点图→带平滑线的散点图
- 右键纵轴→设置坐标轴格式→对数刻度
进阶：Matlab直连C程序（MEX）
通过MEX接口将C译码器编译为Matlab可调用的二进制文件，实现：
- 避免中间文件读写
- 仿真速度提升10-100倍
- 实时动态绘制曲线
六、我的绘图心得
1. 半对数坐标是必须的 ：第一次我用plot(snr, ber)，结果所有点都挤在底部，根本看不清。后来才知道BER要用semilogy。
2. 数据点不够多时，曲线会不平滑：我只跑了5个SNR点，每个点5000帧。如果想得到光滑曲线，需要增加SNR采样点（比如步长0.5dB）和每点帧数（比如20000帧）。
3. 处理BER=0的情况 ：当BER=0时，semilogy会报错（log(0)无穷大）。所以我在数据中把0改为了一个很小的数，如1e-8，或者直接不画那个点（用NaN）。
matlab
ber_spa(ber_spa == 0) = NaN; % 使该点不显示4. 图例位置 ：我习惯把图例放在左下角（Location, southwest），因为曲线集中在左上区域。
七、结语
从C语言跑出数据，到Matlab画出曲线，我终于完整地验证了：SPA确实比Min-Sum更优，尤其在低信噪比区域。虽然计算量大了几倍，但在深空通信、卫星广播等对可靠性要求极高的场景中，这点代价是值得的。
曲线图不仅是给审稿人看的，更是给自己看的------它让你一眼就看到算法的优缺点，指引你下一步的优化方向。比如，看到SPA在3dB后性能陡降，你会思考：能不能在2~3dB区间进一步改进？归一化Min-Sum会不会是更好的折中？
如果你也跑出了自己的BER曲线，欢迎在评论区分享你的结果。完整代码（C + Matlab）已上传GitHub。
最后，安利一个在线工具：Desmos也可以画半对数图，不过需要手动输入数据。