数值计算实验2.docVIP

本科实验报告 课程名称： 计算机数值方法B 实验项目： 线性方程组的直接解法 实验地点： 逸夫楼302 专业班级： 学号： 学生姓名： 指导教师： 2013年 4月 13日 一、实验目的和要求 （1）了解线性方程组常见的直接解法，如Guass消元法、LU分解法、追赶法。 （2）加深对线性方程组求解方法的认识，掌握算法。 （3）会进行误差分析，并能对不同方法进行比较。 二、实验内容和原理 合理利用Gauss消元法、LU分解法或追赶法求解下列方程组： ① ② ③ ④（n=5,10,100,…） 1、高斯分解法： ⑴将原方程组化为三角形方阵的方程组： lik=aik/akk aij= aij- lik* akj k=1,2,…,n-1 i=k+1,k+2, …,n j=k+1,k+2, …,n+1 ⑵由回代过程求得原方程组的解： xn= ann+1/ ann xk=( akn+1-∑akj xj)/ akk (k=n-1,n-2, …,2,1) 2、LU分解法： 将系数矩阵A转化为A=L*U, L为单位下三角矩阵，U为普通上三角矩阵，然后通过解方程组l*y=b,u*x=y,来求解x。 3、追赶法: 用来求对角方程组；将系数矩阵A转化为A=L*U, L为普通下n-1对角矩阵，U为单位上n-1对角矩阵，然后通过解方程组l*y=b,u*x=y,来求解x。 三、主要仪器设备 惠普ProBook6460b VC++运行环境 四、操作方法与实验步骤 （1）Guass消元法 完全主元素消元法 #include stdafx.h #includemath.h #includeiostream using namespace std; float a[100][101]; float x[10]; int N; void output() { for(int i=1;i=N;i++) { for(int j=1;j=N+1;j++) { couta[i][j] ; } coutendl; } } void input() { cout请输入矩阵阶数:endl; cinN; cout请输入矩阵各项:endl; for(int i=1;i=N;i++) for(int j=1;j=N+1;j++) { cina[i][j]; } coutendl; } void main() { int z[10]; int maxi,maxj,i; input(); for(int i=1;i=N;i++) z[i]=i; for(int k=1;kN;k++) { maxi=k;maxj=k;float maxv=abs(a[k][k]); for(i=k;i=N;i++) for(int j=k;j=N;j++) if(abs(a[i][j])maxv) { maxv=abs(a[i][j]);maxi=i;maxj=j; } if(maxi!=k) { for(int j=1;j=N+1;j++) { float t=a[k][j];a[k][j]=a[maxi][j];a[maxi][j]=t; } } if(maxj!=k) { for(i=1;i=N;i++) { float t=a[i][k];a[i][k]=a[i][maxj];a[i][maxj]=t; } int t=z[k];z[k]=z[maxj];z[maxj]=t; } for(int i=k+1;i=N;i++) { float l=a[i][k]/a[k][k]; for(int j=k;j=N+1;j++) { a[i][j]+=-l*a[k][j]; } } } for(i=N;i0;i--) { float s=0; for(int j=i+1;j=N;j++) { s+=a[i][j]*x[z[j]]; } x[z[i]]=(a[i][N+1]-