数值计算基础实验报告模板

《数值计算基础实验报告模板》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、WuhanUniversityofScience&Technology〈〈数值计算基础》实验报告专业：软件工程班级：11级02班学号：_201113138064&名王伟实验名称解线性方程组的直接法实验日期4.12实验目的(10分)实验内容一(35分)实验内容二(35分)实验小结(20分)总分【实验目的】1.理解直接法求解线性方程组的两种方法高斯消元法与三角分解法；2.掌握高斯消元元的思想及方法；3.掌握矩阵的LU分解的两种方法Doolittle和Cmunt分解及使用该方法解线性方程组;4.掌握对正定矩阵Cholesky分解及使用该方法解线性方程组；5.掌握解三对角方

2、程组的追赶法。【实验内容一】一、问题用Doolittle法解方程组二、算法描述(1)

3、ujdoolite函数传入矩阵。(2)计ffU的第一•行，再计算L的第一列，(/(U)=6Z(iy),j=l，2,3...n;/(zl)=6z(zl)/f/(ll),i:2,3...n。(3)计算U的第r行n♦j)=》(咖⑹=IAr印⑹+UM，jkr人=1*=1U、rj、=乞l、rk^J、kj、，j=r，••”tr，r=2,…，n.A=l⑷计算L的第r列6Z(/r)=[/(//:)(7(々r)=^l{ik)U(kr)^-l(ir)U(rri>rk=k=々=1U(rr)(5)输

4、出单位上三角矩阵和单位K三角矩阵。doublesuml，sum2;double*L，*M;L=creat(n);M二creat(n);三、关键源代码voiddoolite(double*p，intn,doubleb[J){inti，j，k;doublex[N],y[N];for(i=0;i

5、)sum1+=*(L+i*n+k)*(*(M+k*n+j))；)*(M+i*n+j)=*(p+i*n+j)-suml;}for(j=i+l;j

6、"单位上三角矩阵为Ann);四、运行结果for(j=0;j

7、j];y[n=bfi]-suml;}printf("Y的值为for(i=

8、0;i=0;i—){sum2=0;for(j=i+l;j

9、13孩矩阵为，2.0010.000.00-3.00-3.00-4.00-12.0013.001.002.003-00-4.004.0014.009.00-13.00单位上三角矩阵为：2.000010.00000.0000一3■00000.0醐11.0000—12.00000.0醐0.0000-0.2727一0.1818e.@0000.00000.0000—4-0000单位下二角矩阵为■■1.00000.00000.00000.0000-1.50001.00000.00000.00000.5000-0.27271.00000.00002.0000-0.5455-9