雅可比迭代c语言实验报告,Jacobi迭代和Gauss迭代 c语言实现

/*

*TestMain.cpp

*功能：利用Jacobi迭代和Gauss迭代求解方程组

*时间：2013.4.1

*/

#include "Jacobi.h"

#include "Gauss.h"

int main ()

{

printf("Jacobi迭代的结果：\n");

Jacobi(3 , 0.0001);

printf("*****************************************************\n");

printf("Gauss迭代的结果：\n");

Gauss(3,0.0001);

return 0;

}

/*

*Gauss.cpp

*功能：高斯迭代

*时间：2013.4.1

*/

#include "Gauss.h"

void Gauss(int n , double E)

{

int i; // 用来循环

int k=0; //记录迭代次数

double maxE;

//设置初始值

double *a = (double *)malloc(sizeof(double)*(n+1));

for(i=1;i<=n;i++)

{

a[i]=0;

printf("%d :X0=%f ,",k,a[i]);

}

printf("\n");

//存放新的解

double *temp = (double *)malloc(sizeof(double)*(n+1));

do

{

CalculateG(a,temp);

k++;

printf("%d :",k);

for(i=1;i<=n;i++)

{

printf("X%d=%f ,",i,temp[i]);

}

printf("\n");

swapX(temp,a,n);

maxE = fabs(temp[1]-a[1]);

for(i=2;i<=n;i++)

{

if(maxE

maxE = fabs(temp[i]-a[i]);

}

}while(maxE>=E);

return ;

}

//计算x

void CalculateG(double *a,double *temp)

{

temp[1] = -2.4-0.4*a[2]-0.2*a[3];

temp[2] = 5+0.25*temp[1]+0.5*a[3];

temp[3] = 0.3-0.2*temp[1]+0.3*temp[2];

return ;

}

/*

*Jacobi.cpp

*功能：Jacobi迭代

*时间：2013.4.1

*/

#include "Jacobi.h"

void Jacobi(int n , double E)

{

int i; // 用来循环

int k=0; //记录迭代次数

double maxE;

//设置初始值

double *a = (double *)malloc(sizeof(double)*(n+1));

for(i=1;i<=n;i++)

{

a[i]=0;

printf("%d :X0=%f ,",k,a[i]);

}

printf("\n");

//存放新的解

double *temp = (double *)malloc(sizeof(double)*(n+1));

do

{

CalculateJ(a,temp);

k++;

printf("%d :",k);

for(i=1;i<=n;i++)

{

printf("X%d=%f ,",i,temp[i]);

}

printf("\n");

swapX(temp,a,n);

maxE = fabs(temp[1]-a[1]);

for(i=2;i<=n;i++)

{

if(maxE

maxE = fabs(temp[i]-a[i]);

}

}while(maxE>=E);

return ;

}

//计算x一课本3.1为例

void CalculateJ(double *a,double *temp)

{

temp[1] = -2.4-0.4*a[2]-0.2*a[3];

temp[2] = 5+0.25*a[1]+0.5*a[3];

temp[3] = 0.3-0.2*a[1]+0.3*a[2];

return ;

}

/*

*common.cpp

*功能:Jacobi和Gauss迭代公用函数

*/

#include "common.h"

//交换解

void swapX(double *temp,double *a,int n)

{

int i;

double tmp;

for(i=1;i<=n;i++)

{

tmp = temp[i];

temp[i] = a[i];

a[i] = tmp;

}

return ;

}

相关头文件:

/*

*Jacobi.h

*/

#include "common.h"

void Jacobi(int n , double E);

void CalculateJ(double *a,double *temp);

/*

*Gauss.h

*/

#include "common.h"

void Gauss(int n,double E);

void CalculateG(double *a,double *temp);

/*

*common.h

*/

#include

#include

#include

void swapX(double *temp,double *a,int n);