从生成器矩阵生成Paritycheck矩阵

Question

（numpy）中是否有函数或经过良好测试的函数来计算奇偶校验矩阵（https://en.wikipedia.org/wiki/Parity-check_matrix）

从发电机矩阵

附言

我在这个网站上没有找到解决方案

Answer 1

若我的理解是正确的，那个么奇偶校验矩阵就是生成器的零空间 模2中的矩阵。在scipy中有解决方案，但是这个函数 给非整数零空间。您可以使用sympy，但对于大型应用程序，它可能会很慢 矩阵

"""
>>> np.set_string_function(str)
>>> h
[[0 1 1 1 1 0 0]
 [1 0 1 1 0 1 0]
 [1 1 0 1 0 0 1]]
>>> (g @ h.T) % 2
[[0 0 0]
 [0 0 0]
 [0 0 0]
 [0 0 0]]
"""
import sympy
import numpy as np

g = np.array([[1, 1, 1, 0, 0, 0, 0],
       [1, 0, 0, 1, 1, 0, 0],
       [0, 1, 0, 1, 0, 1, 0],
       [1, 1, 0, 1, 0, 0, 1]])
h = np.array(sympy.Matrix(g).nullspace()) % 2

其中h是奇偶校验矩阵