443解三對角方程組的matlab函數(shù)文件(編輯修改稿)
2024-11-17 16:34 本頁面
【文章內(nèi)容簡介】 for i=j+1:n L(i,j)=U(i,j)/ U(j,j)。 第四章方程組的直接解法 U(i,i)=U(I,:)L(I,j)*u(j,:)。 end end function[L,U]=doolittle(A) [n,m]=size(A)。 U=zeros(n,n),L=eye(n)。 for k=1:n U(1,k)=A(1,k)。L(k,1)=A(k,1)/U(1,1)。 第四章方程組的直接解法 end for k=2:n U(k,k)=A(k,k)L(k,1:k1)*U(1:k1,k)。 for j=k+1:n U(k,j)=A(k,j)L(k,1:k1)*U(1:k1,j)。 L(j,k)=(A(j,k)L(j,1:k1)*U(1:k1,k))/U(k,k)。 end
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號-1