| 查看: 2493 | 回复: 1 | ||
[求助]
Matlab对两个矩阵做==相等比较运算结果不正确 已有1人参与
|
|
是这样，这个问题我遇到过好几次，一直没有解决。下面是问题的样例 我先创建一个rand矩阵 a = rand(5); 求其逆矩阵a_inv a_inv = inv(a); 再求其伴随矩阵a_adj //adj函数在文末附上 a_adj = adj(a); 按照线性代数的结论，a_inv = a_adj / |a| 但是在程序里面我输入： a_inv == a_adj / det(a) 返回的是一个0矩阵，按理来说是相等的啊 ans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 我觉得是精度的问题，哪位大神能帮忙解决这个基础的问题？ 附 adj函数: function B = adj(A) %求伴随矩阵 %ADJ Matrix adjoint. % ADJ(A) is the adjoint matrix of square matrix A. % It is computed using the Cayley-Hamilton Theorem. % The inverse of A is: INV(A) = ADJ(A)/det(A). % % Matrices that are not invertable still have an adjoint. %written by Paul Godfrey, April, 1998 %pjg@mlb.semi.harris.com ce = poly(eig(A)); cesize = max(size(ce)); p = [0 ce(1  cesize-1))];s = (-1)^(max(size(A))+1); B = s*polyvalm(p,A); |
回复此楼» 猜你喜欢
FileCode能看出啥？
已经有23人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有36人回复
-
关于Filecode分析方法
已经有9人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有10人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab如何将得到的结果带到下一个方程中
已经有12人回复
- matlab中fmincon的fun是矩阵运算怎么办？
已经有5人回复
- 用matlab组集刚度矩阵时，总刚度矩阵对应的行列式不为零
已经有5人回复
- matlab简单矩阵运算问题请教
已经有7人回复
- matlab中std2函数计算矩阵的标准差时，分母是除以了n-1还是n？
已经有3人回复
- Fortran 求解复数矩阵SVD时，用到cgesvd，结果与matlab不同，为什么？
已经有15人回复
- Matlab矩阵计算
已经有3人回复
- matlab求逆矩阵是不是有问题？
已经有13人回复
- matlab 超出矩阵维度问题
已经有12人回复
- 请大侠帮忙看一下我这个matlab小程序结果为什么不返回矩阵
已经有7人回复
- matlab循环计算转为矩阵 30金币求助
已经有7人回复
- matlab中如何随机不重复选取矩阵中的元素
已经有4人回复
- MATLAB FEM 刚度矩阵
已经有12人回复
- 如何用Matlab画矩阵的散点图？
已经有3人回复
- matlab 有限元计算扩散问题，建立整体矩阵好慢。大家帮忙看看代码
已经有8人回复
- matlab中系数是矩阵的积分
已经有5人回复
- matlab线性矩阵不等式求解做优化问题
已经有10人回复
- matlab 不规则矩阵 求平均值
已经有6人回复
libralibra
至尊木虫 (著名写手)
骠骑将军
- 程序强帖: 40
- 应助: 817 (博后)
- 金币: 12914.1
- 红花: 64
- 帖子: 2238
- 在线: 287.3小时
- 虫号: 696514
- 注册: 2009-02-05
- 专业: 计算机软件
2楼2014-04-08 17:54:07
