互补判断矩阵p构造一个G python代码

以下是构造互补判断矩阵 p 并生成 G 矩阵的 Python 代码：

import numpy as np

n = 4  # 矩阵维度
p = np.zeros((n, n))  # 初始化互补判断矩阵 p

# 构造互补判断矩阵 p
for i in range(n):
    for j in range(n):
        if i == j:
            p[i][j] = 0
        elif i < j:
            p[i][j] = 1
            p[j][i] = -1

# 生成 G 矩阵
G = np.kron(np.eye(n), p) + np.kron(p.T, np.eye(n))

print("互补判断矩阵 p：\n", p)
print("生成的 G 矩阵：\n", G)

代码中，我们首先定义了矩阵的维度 n，并初始化了一个 n*n 的零矩阵 p。然后，按照互补判断矩阵的定义，对 p 进行构造。

接着，我们利用 `np.kron` 函数生成 Kronecker 积矩阵，将 p 矩阵与单位矩阵相乘，生成左半部分的 G 矩阵；将 p 的转置矩阵与单位矩阵相乘，生成右半部分的 G 矩阵，再将两个 G 矩阵相加即可得到最终的 G 矩阵。

最后，我们通过 `print` 函数输出生成的互补判断矩阵 p 和 G 矩阵。