24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2437  |  回复: 1

dangyuluo

新虫 (初入文坛)

[求助] 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(1cesize-1))];
s = (-1)^(max(size(A))+1);
B = s*polyvalm(p,A);
回复此楼

» 猜你喜欢

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

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

libralibra

至尊木虫 (著名写手)

骠骑将军

【答案】应助回帖

感谢参与,应助指数 +1
在a_inv == a_adj / det(a)之前,加入代码,显示要比较的2个矩阵到command window看看想等不相等,语句末尾不要加分号
a_inv
a_adj/det(a)
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
2楼2014-04-08 17:54:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 dangyuluo 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿2110,化学学硕310分,本科重点双非求调剂 +13 努力奋斗112 2026-04-08 13/650 2026-04-09 00:27 by Evan_Liu
[考研] 求调剂,262机械专硕 +5 嗯yyl 2026-04-08 5/250 2026-04-08 23:27 by 6781022
[考研] 生物学328分求调剂 +9 闪电kkl 2026-04-08 10/500 2026-04-08 21:42 by liuhuiying09
[考研] 一志愿211,化学310分,本科重点双非,求调剂 +13 努力奋斗112 2026-04-08 13/650 2026-04-08 21:17 by 学员tURuqU
[考研] 368化学求调剂 +12 wwwwabcde 2026-04-07 13/650 2026-04-08 20:48 by LiuXingYu星星
[考研] 11408 325分 +3 jgtxuxgkx 2026-04-07 3/150 2026-04-07 23:10 by lbsjt
[考研] 生物工程求调剂 +13 喜欢还是不甘心 2026-04-05 13/650 2026-04-07 16:55 by Ecowxq666!
[考研] 285求调剂 +15 哦呦呼o 2026-04-04 17/850 2026-04-06 23:02 by chenzhimin
[考研] 327考研调剂推荐 +6 呜呜呜呜呢 2026-04-06 6/300 2026-04-06 21:39 by 啵啵啵0119
[考研] 307求调剂 +3 所念及所望 2026-04-06 3/150 2026-04-06 17:30 by 土木硕士招生
[考研] 第一志愿东南大学物理313,有科研竞赛获奖经历,希望物理复试调剂 +3 马内橙 2026-04-05 3/150 2026-04-06 10:32 by 蓝云思雨
[考研] 复试调剂 +5 asdasdassda 2026-04-05 5/250 2026-04-06 09:32 by dongzh2009
[考研] 求调剂,一志愿郑州大学材料与化工专硕,英二数二342分,求老师收留 +19 v12abo 2026-04-02 21/1050 2026-04-06 09:29 by 蓝云思雨
[考研] 085600调剂 +9 东照照照 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
[考研] 341求调剂 +3 学无止境,冲 2026-04-05 3/150 2026-04-05 09:40 by lbsjt
[考研] 可跨专业调剂 +3 周的得地 2026-04-04 6/300 2026-04-04 22:21 by barlinike
[考研] 环境285分,过六级,求调剂 +10 xhr12 2026-04-02 10/500 2026-04-04 21:53 by bn53987
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 303求调剂 +3 一色清羽 2026-04-02 4/200 2026-04-03 10:22 by 蓝云思雨
[考研] 298求B区调剂 +4 zzz,,r 2026-04-02 5/250 2026-04-02 12:17 by 土木硕士招生
信息提示
请填处理意见