版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

dangyuluo

新虫 (初入文坛)

应助: 0 (幼儿园)
帖子: 4
在线: 4.9小时
虫号: 3098025
注册: 2014-03-28
专业: 工程热力学

[求助] 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);

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

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人回复

1楼 2014-04-08 16:50:05

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

程序强帖: 40
应助: 817 (博后)
金币: 12914.1
红花: 64
帖子: 2238
在线: 287.3小时
虫号: 696514
注册: 2009-02-05
专业: 计算机软件

【答案】应助回帖

感谢参与，应助指数 +1

在a_inv == a_adj / det(a)之前,加入代码,显示要比较的2个矩阵到command window看看想等不相等,语句末尾不要加分号
a_inv
a_adj/det(a)

赞一下