解线性方程组的直接方法实验报告.doc

解线性方程组的直接方法实验报告
解线性方程组的直接方法实验报告
实验目的：
1、 通过该课题的实验，体会模块化结构程序设计方法 的优点；
2、 运用所学的计算方法，解决各类线性方程组的直接 算法；
3、 提高分析和解决问题的能力，解线性方程组的直接方法实验报告
解线性方程组的直接方法实验报告
实验目的：
1、 通过该课题的实验，体会模块化结构程序设计方法 的优点；
2、 运用所学的计算方法，解决各类线性方程组的直接 算法；
3、 提高分析和解决问题的能力，做到学以致用；
4、 通过三对角形线性方程组的解法，体会稀疏线性方 程组解法的特点。
：
实验代码：
ttinclude stdio. h
ttinclude math, h
ttinclude
using namespace std;
//Gauss 法
void Izy(double **a, double *b, int n)
{
int i, j, k;
double 1, x[10], temp;
for (k=0;kn-1;k++)
if (j==k)
temp=fabs(a[j][k]);
else if(tempfabs(a[j][k]))
temp=fabs(a[j][k]);
i=j；
}
}
if(temp==O)
{
cout无解 ；
return;
}
else
{
for(j=k;jn;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+l;in;i++)
{
l=a[i][k]/a[k][k];
for (j=k;jn;j++)
a[i] [j]=a[i] [j]-l*a[k] [j];
b[i]=b[i]-l*b[k];
}
}
if (a[n-l][n-l]==0)
(
cout无解 ；
return;
}
x[n-l]=b[n-l]/a[n-l][n-1];
for (i=n-2;i=0;i一)
{
temp=0;
for(j=i+l;jn;j++)
temp=temp+a[i][j]*x[j];
x[i] = (b[i]-temp)/a[i] [i];
printf(x%d=%lf ,i+l,x[i]);
printf();
}
}
〃平方根法
void pfg(double **a, double *b, int n)
{
int i, k, m;
double x[8], y[8], temp;
for (k=0;kn;k++)
{
temp=0;
for (m=0;mk;m++)
temp=temp+pow(a[k][m], 2);
if (a[k][k]temp)
return;
a[k] [k]=pow( (a[k] [k]-temp), 1. 0/2. 0);
for (i=k+l;in;i++)
{
temp=0;
for (m=0;mk;m++)
temp=temp+a[i][m]*a[k][m];
a[i] [k]