实验2LU分解求线性方程组.docVIP

《数值分析》实验2 一．实验名称：LU分解法解线性方程组 二、实验目的： （1）学会对矩阵进行LU分解。解线性代数方程组。 在以上程序的基础上进一步编写出利用LU分解的方法求解线性方程组的程序。 五、实验结果 输出L、U及方程的解x。 #include #include int main int n 4,i,k,j,m n+1; float a[m][m],u[m][m],l[m][m],t,t1,x[m],y[m],b[m]; for i 1;i n;i++ for j 1;j n;j++ u[i][j] 0; l[i][j] 0; for i 1;i n;i++ l[i][i] 1; b[1] 2; b[2] 8; b[3] 4; b[4] 9; printf 请输入%d阶矩阵:\n,n ; for i 1;i n;i++ for j 1;j n;j++ scanf %f,a[i][j] ; for i 1;i n;i++ if i 1 l[1][1] 1; for k 1;k n;k++ u[i][k] a[i][k]; for k 2;k n;k++ l[k][i] a[k][i]/u[i][i]; if 2 ii n for k i;k n;k++ t 0; for j 1;j i-1;j++ t+ l[i][j]*u[j][k]; u[i][k] a[i][k]-t; for k i+1;k n;k++ t1 0; for j 1;j i-1;j++ t1+ l[k][j]*u[j][i]; l[k][i] a[k][i]-t1 /u[i][i]; y[1] b[1]; for i 2;i n;i++ t 0; for j 1;j i-1;j++ t+ l[i][j]*y[j]; y[i] b[i]-t; x[n] y[n]/u[n][n]; for i n-1;i 1;i-- t 0; for j n;j i+1;j-- t+ u[i][j]*x[j]; x[i] y[i]-t /u[i][i]; for i 1;i n;i++ for j 1;j n;j++ printf l[%d][%d] %0.4f,i,j,l[i][j] ; printf \n ; printf \n ; for i 1;i n;i++ for j 1;j n;j++ printf u[%d][%d] %0.4f,i,j,u[i][j] ; printf \n ; for i 1;i n;i++ printf x[%d] %0.4f,i,x[i] ;