我目前用的办法是用
magic(n)生成随机的n阶矩阵
然后用det()检查是否非奇异矩阵
但这种办法效率太低了。
大家有好点的办法吗?
对了,v2ex能加一个Matlab板块吗?
magic(n)生成随机的n阶矩阵
然后用det()检查是否非奇异矩阵
但这种办法效率太低了。
大家有好点的办法吗?
对了,v2ex能加一个Matlab板块吗?
1
yangzh Mar 14, 2015
|
2
meteor2013 OP |
3
h4x3rotab Mar 14, 2015
非奇异矩阵 <=> 行列式不为0 <=> 可以高斯-约当消元得到单位矩阵
反过来,拿一个单位矩阵,随机做线性变换,是不是就能得到一个非奇异矩阵? 脑洞出来的,不知道对不对 |
4
ruoyu0088 Mar 14, 2015 你确定magic函数是生成随机矩阵的吗?
|
5
jokester Mar 14, 2015
re 4f
然后用det判定singuarity有点那个 rank应该会快些 |
6
facat Mar 14, 2015
我以前试过用rand,生成的矩阵几乎不奇异。
或者楼主先找出一个非奇异矩阵,然后在此基础上做初等变换,依据线性代数的知识,得到的矩阵仍然是非奇异的。 |