| 查看: 2404 | 回复: 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); |
» 猜你喜欢
有没有人能给点建议
已经有5人回复
假如你的研究生提出不合理要求
已经有12人回复
实验室接单子
已经有7人回复
全日制(定向)博士
已经有5人回复
萌生出自己或许不适合搞科研的想法,现在跑or等等看?
已经有4人回复
Materials Today Chemistry审稿周期
已经有4人回复
参与限项
已经有3人回复
对氯苯硼酸纯化
已经有3人回复
所感
已经有4人回复
要不要辞职读博?
已经有7人回复
» 本主题相关价值贴推荐,对您同样有帮助:
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












cesize-1))];
回复此楼