最好使用Sage或其他合适的工具。在

以下只是简单的非专家尝试,但是旋转高斯消去法应该给出可逆性的确切结果:import random

from scipy.linalg import toeplitz

import numpy as np

def is_invertible_F2(a):

"""

Determine invertibility by Gaussian elimination

"""

a = np.array(a, dtype=np.bool_)

n = a.shape[0]

for i in range(n):

pivots = np.where(a[i:,i])[0]

if len(pivots) == 0:

return False

# swap pivot

piv = i + pivots[0]

row = a[piv,i:].copy()

a[piv,i:] = a[i,i:]

a[i,i:] = row

# eliminate

a[i+1:,i:] -= a[i+1:,i,None]*row[None,:]

return True

n = 10

column = [random.choice([0,1]) for x in xrange(n)]

row = [column[0]]+[random.choice([0,1]) for x in xrange(n-1)]

matrix = toeplitz(column, row)

print(is_invertible_F2(matrix))

print(int(np.round(np.linalg.det(matrix))) % 2)

请注意,np.bool_仅在有限的意义上与F_2相似——Fୱ2中的二进制运算+是bool的-,一元运算-是{}。不过,乘法是一样的。在

^{pr2}$

上面的高斯消去法只使用这些运算,所以它是有效的。在

matlab检验矩阵是否可逆,有限域上矩阵是否可逆的检验相关推荐

  1. matlab 矩阵列乘系数,matlab 给某一列乘上一个系数

    矩阵M是一个 mxn 的矩阵,现在要给M矩阵的第一列都要乘上10,使其第一列扩大10倍,那肿么做呢? 我第一时间用的是: M(:,1) = M(:,1)*10; //错误的 但是这个错了,结果是不对的 ...

  2. matlab判断矩阵不可约,用Matlab计算二元域GF(2)上的不可约多项式

    1 二元域 GF(2) 上的不可约多项式 二元域 GF(2)={0,1} 上的运算规则如下: 加法:+ 0 1 0 0 1 1 1 0 乘法:⋅ 0 1 0 0 0 1 0 1 二元域 GF(2) 上 ...

  3. 特征值与特征向量、特征方程、特征多项式、矩阵相似、相似变换、矩阵对角化、奇异值分解(Singular Value Decomposition)手算加MATLAB

    特征值(Eigenvalue)与特征向量(Eigenvector)   定义:若AAA为n×nn\times nn×n的矩阵,xxx为非零向量,若存在数λ\lambdaλ使得Ax=λxAx=\lamb ...

  4. 【机器学习的数学基础】(六)矩阵分解(Matrix Decomposition)(上)

    文章目录 4 矩阵分解(Matrix Decompositions)(上) 4.1 行列式与迹 4.2 特征值和特征向量 4 矩阵分解(Matrix Decompositions)(上) 在第2章和第 ...

  5. 可逆矩阵的秩等于矩阵的阶数_矩阵论一些总结点

    两周前考完了矩阵论,本来想把整本书的要点都扔上来,感觉没这个必要,就放一些总结性的东西吧!(排版比较简陋,大家凑活看吧 对角化的充要条件 矩阵有 n 个特征向量 矩阵的 s 重特征值对应有 s 个特征 ...

  6. Matlab 基础知识——矩阵操作及运算(矩阵、数组区别)

    看论文时,经常看到矩阵,但在记忆里又看到数组.那么问题来了,矩阵和数组分别是什么?二者有什么区别?看论文时,经常看到矩阵,但在记忆里又看到数组.那么问题来了,矩阵和数组分别是什么?二者有什么区别? 在 ...

  7. matlab矩阵检索、嵌套,矩阵操作笔记

    matlab矩阵检索.嵌套,矩阵操作笔记 学习是一个输入和输出的过程,遂老衲决定多做笔记,同时多举例子,因为举例子是从抽象到具象的过程 矩阵点乘就是逐元素操作,矩阵直接乘起来就是按照矩阵的乘法,左行右 ...

  8. matlab因子载荷矩阵正交旋转,因素分析中的矩阵旋转

    因素分析中的矩阵旋转 因素分析法 因素分析是一种统计技术,目的是从众多的可观测的"变量"中,概括和推论少数"因素".用最少数的"因素"来概括 ...

  9. Matlab中创建、串联和扩展矩阵

    目录 构建数据矩阵 专用矩阵函数 串联矩阵 生成数值序列 扩展矩阵 空数组 最基本的 MATLAB® 数据结构体是矩阵.矩阵是按行和列排列的数据元素的二维矩形数组.元素可以是数字.逻辑值(true 或 ...

  10. matlab矩阵指定行最大值,求Matlab程序:在2行矩阵中,如何求第1行最大值和第2行相应的最大,请问,如何在MATLAB中编写程序实现求两矩阵A*B,A.*...

    导航:网站首页 > 求Matlab程序:在2行矩阵中,如何求第1行最大值和第2行相应的最大,请问,如何在MATLAB中编写程序实现求两矩阵A*B,A.* 求Matlab程序:在2行矩阵中,如何求 ...

最新文章

  1. 2018-4-15摘录笔记,《网络表征学习前沿与实践》 崔鹏以及《网络表征学习中的基本问题初探》 王啸 崔鹏 朱文武
  2. 给.net程序打内存补丁-转
  3. python小程序-python学习—几个简单小程序
  4. Fedora 16 硬盘安装体会
  5. iOS后台运行的设置方法
  6. 最长非下降子序列(O(nlogn))(offer收割)
  7. 批量裁剪GIS数据(包含GDB,MDB,Shp)
  8. python导出excel文件数字签名_Python使用RSA+MD5实现数字签名
  9. 枚举函数enumerate
  10. vue elementui表格数据
  11. 中学办公室计算机管理制度,中学多媒体室使用管理制度条例
  12. Springboot物业管理系统
  13. 操作系统android9.0,三星公布了升级Android9.0操作系统的时间表
  14. 数字信号处理——DFT
  15. 【Sass/SCSS】预加载器中的“轩辕剑”
  16. List<? extends T>和List<? super T>之间的区别
  17. P5960 【模板】差分约束算法
  18. 机器学习--KNN算法应用,iris鸢尾花数据集的分类
  19. 如何使用git 生成密钥?
  20. 计算机英语教学设计反思,小学五年级英语教学设计与反思

热门文章

  1. 如何在小方框上打对号 小方框内打对勾 word 方框打对勾
  2. unreal-教程-第十八章:UI跟随的基本原理(c++)
  3. 计算机数据库基础知识填空题,数据库练习题(基础)
  4. LaTex 之 各类括号总结
  5. 复权、前复权和后复权
  6. maven使用modules、parent标签时遇到的问题
  7. adb 连接某个wifi_adb wifi 测试(无需root)
  8. html5_滑条等其他标签
  9. hc05刷hid固件_Arduino使用HC05蓝牙模块与手机连接
  10. CAD梦想画图中的“打断命令”