Can you please right the matlab code for this question? Finish the function. Out
ID: 3403938 • Letter: C
Question
Can you please right the matlab code for this question?
Finish the function. Output should follow these rules (note it uses functions from above):
if diag_dom(A) = 1, use Jacobi iteration to solve the system (Ax = b) and output x (out = x)
else, if the iteration matrix, M, for the Jacobi method is such that will_converge(M)= 1, solve using Jacobi Method, output x
else, return an empty matrix (i. e. []). function out = iterative_solve(A,b,initial_x,tol)
Code:
%A is a square matrix, b a column vector with equal rows
%If possible solves Ax= b using Jacobi Method starting with
%initial_x, ending when norm(A*x_i -b,inf)
Explanation / Answer
clear all;
if diag_dom(A) = 1
out = iterative_solve(A,b,intial_x,tol);
end
elseif will_converge(M)=1
out = iterative_solve (A,b,intial_x,tol);
end
else
out = 0;
end
%function for Jacobian Method
function out = iterative_solve(A,b,initial_x,tol)
%Matrix A is split into the three matrices L, U and D such that A=D-L-U
D = diag(diag(A));
L = tril(-A,-1);
U = triu(-A,1);
%transition matrix and constant vector used for iterations
% As Ax = b, or (D-L-U)x = b, or Dx = (L+U)x + b, or x = Inv(D)*(L+U)x + Inv(D)*b
% transition matrix, T = Inv (D)*(L+U) & constant matrix, c = Inv(D)*b
T = inv(D)*(L+U);
c = inv(D)*b;
x = intial_x;
while 1
x(:,k+1) = T*x(:,k) + c;
if norm(A*x(:,k+1)-b,inf)
out = x(:,k+1);
break
end
end