matlab实验课件6.pptVIP

Introduction to Scientific Computing-- A Matrix Vector Approach Using Matlab Written by Charles F.Van Loan 陈 文 斌 Wbchen@fudan.edu.cn 复旦大学 Chapter 6 Linear systems Triangular Problems Banded Problems Full Problems Analysis Computer L function [L,U] = GE(A) [n,n] = size(A); for k=1:n-1 A(k+1:n,k) = A(k+1:n,k)/A(k,k); A(k+1:n,k+1:n) = A(k+1:n,k+1:n) - A(k+1:n,k)*A(k,k+1:n); end L = eye(n,n) + tril(A,-1); U = triu(A); [L,U]= GE(A); y = LTriSol(L,b); x = UtriSol(U,y); ShowGe.m Stability ? Delta x(1) x(2) ---------------------------------------------------