版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

dangyuluo

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
帖子: 4
在線: 4.9小時(shí)
蟲號(hào): 3098025
注冊(cè): 2014-03-28
專業(yè): 工程熱力學(xué)

[求助] Matlab對(duì)兩個(gè)矩陣做==相等比較運(yùn)算結(jié)果不正確已有1人參與

是這樣，這個(gè)問題我遇到過好幾次，一直沒有解決。下面是問題的樣例

我先創(chuàng)建一個(gè)rand矩陣
a = rand(5);

求其逆矩陣a_inv
a_inv = inv(a);
再求其伴隨矩陣a_adj //adj函數(shù)在文末附上
a_adj = adj(a);

按照線性代數(shù)的結(jié)論，a_inv = a_adj / |a|

但是在程序里面我輸入：
a_inv == a_adj / det(a)

返回的是一個(gè)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

我覺得是精度的問題，哪位大神能幫忙解決這個(gè)基礎(chǔ)的問題？

附 adj函數(shù):

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);

回復(fù)此樓

» 猜你喜歡

290求調(diào)劑已經(jīng)有3人回復(fù)
一志愿西南交通專碩材料355 本科雙非求調(diào)劑已經(jīng)有4人回復(fù)
295復(fù)試調(diào)劑已經(jīng)有6人回復(fù)
279分求調(diào)劑一志愿211 已經(jīng)有8人回復(fù)
工科材料085601 279求調(diào)劑已經(jīng)有8人回復(fù)
317求調(diào)劑已經(jīng)有8人回復(fù)
一志愿南昌大學(xué)，327分，材料與化工085600 已經(jīng)有5人回復(fù)
274求調(diào)劑已經(jīng)有7人回復(fù)
317求調(diào)劑已經(jīng)有9人回復(fù)
招收調(diào)劑碩士已經(jīng)有11人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

matlab如何將得到的結(jié)果帶到下一個(gè)方程中已經(jīng)有12人回復(fù)
matlab中fmincon的fun是矩陣運(yùn)算怎么辦？已經(jīng)有5人回復(fù)
用matlab組集剛度矩陣時(shí)，總剛度矩陣對(duì)應(yīng)的行列式不為零已經(jīng)有5人回復(fù)
matlab簡(jiǎn)單矩陣運(yùn)算問題請(qǐng)教已經(jīng)有7人回復(fù)
matlab中std2函數(shù)計(jì)算矩陣的標(biāo)準(zhǔn)差時(shí)，分母是除以了n-1還是n？已經(jīng)有3人回復(fù)
Fortran 求解復(fù)數(shù)矩陣SVD時(shí)，用到cgesvd，結(jié)果與matlab不同，為什么？已經(jīng)有15人回復(fù)
Matlab矩陣計(jì)算已經(jīng)有3人回復(fù)
matlab求逆矩陣是不是有問題？已經(jīng)有13人回復(fù)
matlab 超出矩陣維度問題已經(jīng)有12人回復(fù)
請(qǐng)大俠幫忙看一下我這個(gè)matlab小程序結(jié)果為什么不返回矩陣已經(jīng)有7人回復(fù)
matlab循環(huán)計(jì)算轉(zhuǎn)為矩陣 30金幣求助已經(jīng)有7人回復(fù)
matlab中如何隨機(jī)不重復(fù)選取矩陣中的元素已經(jīng)有4人回復(fù)
MATLAB FEM 剛度矩陣已經(jīng)有12人回復(fù)
如何用Matlab畫矩陣的散點(diǎn)圖？已經(jīng)有3人回復(fù)
matlab 有限元計(jì)算擴(kuò)散問題，建立整體矩陣好慢。大家?guī)兔纯创a 已經(jīng)有8人回復(fù)
matlab中系數(shù)是矩陣的積分已經(jīng)有5人回復(fù)
matlab線性矩陣不等式求解做優(yōu)化問題已經(jīng)有10人回復(fù)
matlab 不規(guī)則矩陣求平均值已經(jīng)有6人回復(fù)

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

libralibra

至尊木蟲 (著名寫手)

驃騎將軍

程序強(qiáng)帖: 40
應(yīng)助: 817 (博后)
金幣: 12914.1
紅花: 64
帖子: 2238
在線: 287.3小時(shí)
蟲號(hào): 696514
注冊(cè): 2009-02-05
專業(yè): 計(jì)算機(jī)軟件

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

在a_inv == a_adj / det(a)之前,加入代碼,顯示要比較的2個(gè)矩陣到command window看看想等不相等,語句末尾不要加分號(hào)
a_inv
a_adj/det(a)

贊一下

回復(fù)此樓

matlab/VB/python/c++/Java寫程序請(qǐng)發(fā)QQ郵件:790404545@qq.com

2樓2014-04-08 17:54:07

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 dangyuluo 的主題更新