编程实现杜利特尔分解法.doc

《编程实现杜利特尔分解法.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、北京信息科技大学实验报告实验项目线性方程组求解班级通信1003班姓名/学号李叶林2010010422实验日期2011-12-29成绩一、实验内容利用杜利特尔分解法求解线性方程组二、实验原理一、程序#include#include#defineN3usingnamespacestd;intmain(){doubleA[N+1][N+1]={{0,0,0,0},{0,2,1,1},{0,1,3,2},{0,1,2,2}};doubleL[N+1][N+1]={0};doubleU[N+1][N+1]={0};doubleb[N

2、+1]={0,4,6,5};doubley[N+1];doublex[N+1];inti,j,k,p;for(j=1;j<=N;j++)U[1][j]=A[1][j];L[1][1]=1;for(i=2;i<=N;i++){L[i][1]=A[i][1]/U[1][1];L[i][i]=1;}for(k=2;k<=N;k++){for(j=k;j<=N;j++){doublel=0;for(p=1;p<=k-1;p++)l+=L[k][p]*U[p][j];U[k][j]=A[k][j]-l;}for(i=k+1;i<=N;i++){doublel=0;for(

3、p=1;p<=k-1;p++)l+=L[i][p]*U[p][k];L[i][k]=(A[i][k]-l)/U[k][k];}}y[1]=b[1];for(k=2;k<=N;k++){doublel=0;for(j=1;j<=k-1;j++)l+=L[k][j]*y[j];y[k]=b[k]-l;}x[N]=y[N]/U[N][N];for(k=N-1;k>=1;k--){doublel=0;for(j=k+1;j<=N;j++)l+=U[k][j]*x[j];x[k]=(y[k]-l)/U[k][k];}printf("向量y为:");for(i=1;i<=N

4、;i++)printf("%.1lft",y[i]);printf("");cout<<"L："<

5、="<