
利用MATLAB解题,希望带解释。我理解不了。希望有人可以帮助我理解,附带MATLAB的运行图和过程。
举个栗子 高斯消元法解方程
列主元消去法

解方程组

②LU分解法解方程组
function agui_lu(a,b)
%求可逆矩阵a的Doolittle分解,1返回单位下三角矩阵,u返回上三角矩阵n=length(a);
n=length(a);
u=zeros(n,n);
l=eye(n,n);
u(1,:)=a(1,:);
l(2:n,1)=a(2:n,1)/u(1,1);
for k=2:n
u(k,k:n)=a(k,k:n)-l(k,1:k-1)*u(1:k-1,k:n);
l(k+1:n,k)=(a(k+1:n,k)-l(k+1:n,1:k-1)*u(1:k-1,k))/u(k,k);
end
l
u
%解Ly=b
y=zeros(n,1);
y(1)=b(1);
for k=2:n
y(k)=b(k)-l(k,1:k-1)*y(1:k-1);
end
%解Ux=y
x=zeros(n,1);
x(n)=y(n)/u(n,n);
for k=(n-1):-1:1
x(k)=(y(k) -u(k,(k+1):n)*x((k+1):n))/u(k,k);
end
x
运行截图

方程你自己列嘛