用Python动态演示CDMA正交码片原理：从数学抽象到可视化理解

每次翻开通信原理教材看到CDMA那堆公式就头疼？明明知道正交码片序列很重要，却始终搞不懂为什么内积运算能区分不同信号？今天咱们换个学习方式——用Python代码把抽象概念变成可运行的动态演示。不需要死记硬背，跟着我边写代码边理解，你会发现正交性原来如此直观！

1. 为什么码片序列用-1和+1表示

通信教科书里总把0写成-1，这可不是数学家们的恶作剧。让我们用NumPy做个简单实验：

python复制import numpy as np

# 传统二进制表示
station_a_binary = [0, 0, 0, 1, 1, 0, 1, 0]  
# 转换为±1表示
station_a_polar = [-1 if bit == 0 else 1 for bit in station_a_binary]

print(f"二进制表示: {station_a_binary}")
print(f"极性表示: {station_a_polar}")

运行后会看到：

code复制二进制表示: [0, 0, 0, 1, 1, 0, 1, 0]
极性表示: [-1, -1, -1, 1, 1, -1, 1, -1]

关键优势：

• 乘法运算等价于XOR逻辑：1 * -1 = -1对应1 XOR 0 = 1
• 规格化内积计算更直观：直接相加后求平均
• 零均值特性：方便信号处理中的直流分量消除

提示：在无线通信中，-1和+1可以对应载波相位的180°反转（BPSK调制），这是物理层实现的数学基础

2. 正交码片的Python化验证

正交性是CDMA的核心魔法，让我们用代码拆解这个黑盒子：

python复制def normalized_inner_product(vec1, vec2):
    return np.dot(vec1, vec2) / len(vec1)

# 定义四个正交码片序列
code_A = np.array([-1, -1, -1, 1, 1, -1, 1, 1])
code_B = np.array([-1, -1, 1, -1, 1, 1, 1, -1])
code_C = np.array([-1, 1,