problem with back-substitution code

Question

Sienna Phillips 2020-8-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/582002-problem-with-back-substitution-code

评论： Shaima Al-Shalwi 2022-9-5

I want to make an algorithm for back substitution and am testing it on some U and some b. when I call the function though it says U is not defined. Can anyone see where I have made an error? I am very new to MATLAB!

%--INPUTS------------------------------------------------------------------

%

% U - upper triangular matrix

% b - column vector

%

%--OUTPUT------------------------------------------------------------------

%

% b - solution x of Ux=b, stored in b

%

function b=myBackwardSubstitution(U,b)

d=size(U,1);

U=[1 2 3; 4 5 6; 7 8 9];

b=[1;2;3];

for i=1:d

for j=i+1:d

b(i)=b(i)-L(i,j)*b(j);

end

b(i)=b(i)/L(i,i);

end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Alan Stevens 2020-8-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/582002-problem-with-back-substitution-code#answer_482576

在 MATLAB Online 中打开

Define U and b outside function to which you pass them.

Define L.

%--INPUTS------------------------------------------------------------------
%
% U          -     upper triangular matrix
% b          -     column vector
%
%--OUTPUT------------------------------------------------------------------
%
% b          -     solution x of Ux=b, stored in b
%
% Define U and b outside of function.
U=[1 2 3; 4 5 6; 7 8 9];
b=[1;2;3];
% Call function
b = myBackwardSubstitution(U,b);
function b=myBackwardSubstitution(U,b)
    d=size(U,1);
for i=1:d
    for j=i+1:d
        % You haven't defined L
        b(i)=b(i)-L(i,j)*b(j);
    end
    b(i)=b(i)/L(i,i);
end
end