数值分析上机作业1.docxVIP

数值分析第一次上机作业一、三对角方程组的追赶法三对角方程组的一般形式为 （3-18）或表示为，。当系数矩阵A是严格对角占优时，即满足条件， ，A存在唯一的分解。三对角矩阵的分解有以下很规则的形式得到了的分解式以后，可将求解化为等价的求解以下两个更简单的方程组求解三对角方程组追赶法的计算步骤如下计算分解因子阵 ()（3-19）求解 , （3-20）求解, （3-21）求解三对角方程组追赶法的程序如下function x=zg(a,b,c,d)n=length(b);u(1)=b(1); %LU分解。for k=2:n if u(k-1)==0,D=0,return;endl(k)=a(k)/u(k-1); u(k)=b(k)-l(k)*c(k-1);endy(1)=d(1); % 追赶法求解之追过程，求解Ly=d。for k=2:n y(k)=d(k)-l(k)*y(k-1);endif u(n)==0,D=0,return;endx(n)=y(n)/u(n); % 追赶法求解之赶过程，求解Uz=yfor k=n-1:-1:1 if u(k)==0 ,D=0,break; end x(k)=(y(k)-c(k)*x(k+1))/u(k);end检查上面程序是否正确，并用此程序求解三对角方程组解：上述程序正确，在matlab中输入 b=[2 2 2 2]