第10课：Reed-Solomon码

Reed-Solomon码：字节级纠错

Reed-Solomon(RS)码是BCH码在GF(2^m)上的推广，以符号(symbol)为单位纠错，特别适合纠正突发错误。RS码广泛应用于CD/DVD、QR码、DVB、深空通信和5G NR。

RS码的数学结构

🔧 Verilog实现

// rs_encoder.v - RS(7,5)编码器
module rs_encoder_7_5 (
    input  wire       clk, rst_n,
    input  wire [2:0] symbol_in,
    input  wire       valid_in,
    output reg  [2:0] symbol_out,
    output reg        valid_out
);
    reg [2:0] reg1, reg2;
    wire [2:0] fb = symbol_in ^ reg2;
    wire [2:0] g1 = fb;  // alpha^3系数
    wire [2:0] g0 = fb ^ reg1 ^ reg2;  // alpha^0系数
    reg [2:0] cnt;
    always @(posedge clk or negedge rst_n) begin
        if (!rst_n) begin reg1<=0;reg2<=0;cnt<=0;valid_out<=0; end
        else if (valid_in) begin
            if (cnt &lt; 5) begin
                symbol_out &lt;= symbol_in;
                reg1 &lt;= g1; reg2 &lt;= reg1 ^ g0;
            end else begin
                symbol_out &lt;= reg2; reg2 &lt;= reg1; reg1 &lt;= 0;
            end
            cnt &lt;= cnt + 1;
            valid_out &lt;= 1;
        end else valid_out &lt;= 0;
    end
endmodule

🐍 Python仿真

#!/usr/bin/env python3
import numpy as np
import matplotlib.pyplot as plt
# RS code simulation using univariate polynomial over GF(2^3)
def rs_encode(msg, n=7, k=5):
    # Simplified RS(7,5) encoding
    g = [1, 3, 1]  # generator polynomial coefficients in GF(2^3)
    # Shift message and compute remainder
    codeword = list(msg) + [0]*(n-k)
    for i in range(k):
        coef = codeword[i]
        if coef != 0:
            for j in range(len(g)):
                codeword[i+j] ^= g[j] * coef  # simplified GF mult
    return codeword
print("RS(7,5) encoder ready")

RS码的编码增益与应用

RS码作为MDS(最大距离可分)码，具有最优的最小距离。RS(n,k)码的d_min=n-k+1达到了Singleton界的理论最大值。

RS码编码的硬件实现

系统级仿真与性能评估

完整的通信系统仿真需要考虑多个因素：信道模型、编码增益、同步误差、实现损耗等。以下Python代码提供了完整的系统级仿真框架。

#!/usr/bin/env python3
# 第10课系统级仿真
import numpy as np
import matplotlib.pyplot as plt
from scipy.special import erfc

def ber_theory_bpsk(snr_db):
    return 0.5 * erfc(np.sqrt(10**(snr_db/10)))

def simulate_system(mod_type='bpsk', coding_gain_db=0, num_bits=50000):
    np.random.seed(42)
    snr_range = np.arange(0, 20)
    ber_sim = []
    for snr_db in snr_range:
        effective_snr = snr_db + coding_gain_db
        snr_lin = 10**(effective_snr/10)
        bits = np.random.randint(0, 2, num_bits)
        symbols = 1 - 2*bits
        noise_std = 1.0 / np.sqrt(2*snr_lin)
        noise = noise_std * np.random.randn(len(symbols))
        rx = symbols + noise
        dec = (rx < 0).astype(int)
        ber = np.sum(bits != dec) / num_bits
        ber_sim.append(max(ber, 1e-7))
    return snr_range, ber_sim

snr, ber_u = simulate_system('bpsk', 0)
_, ber_coded = simulate_system('bpsk', 2)

fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(14, 6))
ax1.semilogy(snr, ber_u, 'c-o', markersize=3, label='未编码')
ax1.semilogy(snr, ber_coded, '#10b981-s', markersize=3, label='编码(+2dB)')
ax1.set_xlabel('Eb/N0 (dB)'); ax1.set_ylabel('BER')
ax1.set_title('第10课：BER仿真'); ax1.legend()
ax1.grid(True, alpha=0.3, which='both'); ax1.set_ylim(1e-7, 1)

snr_range2 = np.arange(0, 25)
throughput = [(1 - ber_theory_bpsk(s)) * 1e6 for s in snr_range2]
ax2.plot(snr_range2, np.array(throughput)/1e6, '#f59e0b', linewidth=2)
ax2.set_xlabel('SNR (dB)'); ax2.set_ylabel('吞吐率 (Mbps)')
ax2.set_title('吞吐率 vs SNR'); ax2.grid(True, alpha=0.3)
plt.tight_layout()
plt.savefig('/var/www/ttl/digital-comm/lesson10_sys.png', dpi=100,
            facecolor='#0f172a', edgecolor='none')
print("系统级仿真图已保存")

Verilog实现细节与优化