解线性方程组的直接方法实验报告（实验报告）.docx

解线性方程组的直接方法实验报告
　　解线性方程组的直接方法实验报告
　　：
　　1、通过该课题的实验，体会模块化结构程序设计方法的优点；
　　2、运用所学的计算方法，解决各类线性方程组的直接
解线性方程组的直接方法实验报告
　　解线性方程组的直接方法实验报告
　　：
　　1、通过该课题的实验，体会模块化结构程序设计方法的优点；
　　2、运用所学的计算方法，解决各类线性方程组的直接算法；
　　3、提高分析和解决问题的能力，做到学以致用；
　　4、通过三对角形线性方程组的解法，体会稀疏线性方程组解法的特点。
　　：
　　实验代码：
　　#include ""
　　#include ""
　　#include
　　using namespace std;
　　//Gauss法
　　void lzy(double a,double *b,int n)
　　{
　　int i,j,k;
　　double l,x[10],temp;
　　for(k=0;k　　{
　　for(j=k,i=k;j　　{

　　if(j==k)
　　temp=fabs(a[j][k]);
　　else if(temp　　{
　　temp=fabs(a[j][k]);
　　i=j;
　　}
　　}
　　if(temp==0)
　　{
　　cout　　return;
　　}
　　else
　　{
　　for(j=k;j　　{
　　temp=a[k][j];
　　a[k][j]=a[i][j];
　　a[i][j]=temp;
　　}
　　temp=b[k];
　　b[k]=b[i];
　　b[i]=temp;
　　}

　　for(i=k+1;i　　{
　　l=a[i][k]/a[k][k];
　　for(j=k;j　　a[i][j]=a[i][j]-l*a[k][j];
　　b[i]=b[i]-l*b[k];
　　}
　　}
　　if(a[n-1][n-1]==0)
　　{
　　cout　　return;
　　}
　　x[n-1]=b[n-1]/a[n-1][n-1];
　　for(i=n-2;i>=0;i--)
　　{
　　temp=0;
　　for(j=i+1;j　　temp=temp+a[i][j][j];
　　x[i]=(b[i]-temp)/a[i][i];
　　}
　　for(i=0;i　　{
　　printf("x%d=%lf ",i+1,x[i]);
　　printf(" ");
　　}
　　}

　　//平方根法
　　void pfg(double a,double *b,int n)
　　{
　　int i,k,m;
　　double x[8],y[8],temp;
　　for(k=0;k　　{
　　temp=0;
　　for(m=0;m　　temp=temp+pow(a[k][m],2);
　　if(a