在numpy或matlab中从满秩的非方阵中获得可逆方阵

import numpy as np

# 创建一个满秩的非方阵
A = np.array([[1, 2, 3],
              [4, 5, 6]])

# 计算矩阵的伪逆
A_pinv = np.linalg.pinv(A)

# 检查是否为可逆方阵
is_invertible = np.allclose(A, np.dot(A, np.dot(A_pinv, A)))

print("可逆方阵：")
print(A_pinv)
print("是否为可逆方阵：", is_invertible)

% 创建一个满秩的非方阵
A = [1, 2, 3; 4, 5, 6];

% 计算矩阵的伪逆
A_pinv = pinv(A);

% 检查是否为可逆方阵
is_invertible = isequal(A, A * A_pinv * A);

disp("可逆方阵：");
disp(A_pinv);
disp("是否为可逆方阵：" + is_invertible);