我在 Octave 中的 MatLab 有問題我需要制作一個程式，該程式使用反向替換方法（ https://i.stack.imgur.com/ETPxn.png ）解決三角線性方程組我有一個問題我的代碼，八度顯示此錯誤：
error: =: nonconformant arguments (op1 is 1x1, op2 is 2x2)
error: called from
backsub at line 6 column 10
prog31 at line 3 column 1
run at line 93 column 5
這是我的代碼：
Backsub 函式 (backsub.m)
function X = backsub(A,B)
n = length(B);
X = zeros(n, 1);
X(n) = B(n)./A(n,n);
for k = n-1:-1:1
X(k) = (B(k)-A(k,k 1:n) X(k 1:n))./A(k,k);
end
主要功能（prog31.m）
A = [3 -2 1 -1; 0 4 -1 2; 0 0 2 3; 0 0 0 5];
B = [8; -3; 11; 15];
X = backsub(A, B);
disp(X);
uj5u.com熱心網友回復：
您在-loop 中輸入錯誤，因此出現錯誤*。這是修復： for
function X = backsub(A,B)
n = length(B);
X = zeros(n, 1);
X(n) = B(n)/A(n,n); % remove . because A(n,n) is scalar
for k = flip(1:n-1)
X(k) = (B(k)-(A(k,k 1:n)*X(k 1:n)))/A(k,k); % remove . and replace with *
end % end the for loop
end
附帶說明：for-loop 中的乘積是點積，因為A(k,k 1:n)是行向量，X(k 1:n)而是適當大小的列向量。如果您想逐字使用提供的公式（使用 sigma 和），它將是
X(k) = (B(k)-sum(A(k,k 1:n)'.*X(k 1:n)))/A(k,k);
轉載請註明出處，本文鏈接：https://www.uj5u.com/ruanti/456636.html
上一篇：matlab中向量的連續元素