伊莉討論區

標題: 請問要如何用cofactor matrix求解矩陣A的inverse (限使用 MATLAB撰寫程式) [打印本頁]

---

作者: ~夜神月~    時間: 2016-10-12 10:55 AM     標題: 請問要如何用cofactor matrix求解矩陣A的inverse (限使用 MATLAB撰寫程式)

A =[5 4 2 1;7 1 6 2;4 2 3 1;5 5 6 8];
不可以用inv和det及反除 我目前卡在餘因子的地方?


有位大大寫這樣
function C = cofactor(A, i, j)

% cofactor Matrix of cofactors.
%
% C = cofactor(A) returns the matrix of cofactors of A.
% If A is invertible, then inv(A) = C' / det(A).
%
% C = cofactor(A, i, j) returns the cofactor of
% row i, column j of A.

if nargin == 3
% Remove row i and column j to produce the minor.
M = A;
M(i, = [];
M(:,j) = [];
C = (-1)^(i+j) * det(M);
else
[n,n] = size(A);
for i = 1:n
for j = 1:n
C(i,j) = cofactor(A, i, j);
end
end
end
但nargin == 3 以及[n,n] = size(A); 這邊 我看不太懂 還望各位大大解答

另外這是我自己寫的
A=[5 4 2 1;7 1 6 2;4 2 3 1;5 5 6 8];
i=1:4;
j=1:4;
b=A(i,j);
M=zeros(4,4);
ans=0;
while i<4;
M=A;
M(i, = [];
M(:, j) = [];;
ans=ans+b*-1^(i+j)*det(M)
end

不知為何ans還是0
請各位幫忙解答一下 拜託 拜託

那如果有哪位大大的寫法可以PO上來給在下參考一下嗎?

---

歡迎光臨 伊莉討論區 (http://www21.eyny.com/)

Powered by Discuz!