计算方法 上机操作(C语言).docx

#include stdlib.h#include math.h#include stdio.h#define N 4int agaus(double aNN,double bN,int n);/高斯全主元法求deltaX，具体见程序中注释void GetFun(double xN,double bN);/F(X),结果放在b里void GetJacobi(double xN,double JNN);/计算F(X)，结果放在J里double GetDeltaX(double JNN,double bN,double deltaXN);/求deltaX主程序，调用agausvoid UpdateX(double xN,double deltaXN);/xk+1=xk+deltaXvoid SolveNonlinear(double xN,double eps,int k);/非线性方程组主程序,迭代终止条件是最大迭代次数或者绝对值范数（街区距离）main() int i,k;double eps;double x4=1,1.0,1.0,1.0;double bN;eps=0.0000001; k=100;SolveNonlinear(x,eps,k);GetFun(x,b);printf(n);printf(n);for (i=0; ik)|(NormDeltaXeps)break;void GetFun(x,y)/F(X),结果放在b里double x,y;y0=x0*x0+x1*x1+x2*x2+x3*x3-1.0;y1=2.0*x0*x0+x1*x1-4.0*x2+2*x3*x3;y2=3.0*x0*x0-4.0*x1+x2*x2+3*x3*x3;y3=4.0*x0+x1*x1+2*x2*x2+5*x3*x3-3.0;return;void GetJacobi(x,J)/计算F(X)，结果放在J里double xN,JNN;J00=2*x0;J01=2*x1;J02=2*x2;J03=2*x3;J10=4*x0;J11=2*x1;J12=-4;J13=4*x3;J20=6*x0;J21=-4;J22=2*x2;J23=6*x3; J30=4;J31=2*x1;J32=4*x2;J33=10*x3;return;int agaus(a,b,n)/高斯全主元法解方程，结果放在b里int n;double aNN,bN; int *js,l,k,j,i,is;double d,t; js=malloc(n*sizeof(int);/记载列交换标号l=1;for (k=0;k=n-2;k+) d=0.0;for (i=k;i=n-1;i+)/求全主元for (j=k;jd)d=t; jsk=j;/记住列is=i;/记住行if (d+1.0=1.0)l=0;else /*if (jsk!=k)/列交换 for (i=0;i=n-1;i+)t=aik; aik=aijsk; aijsk=t;*/if (is!=k) for (j=k;j=n-1;j+)/行交换 t=akj; akj=aisj; aisj=t;t=bk; bk=bis; bis=t;/b交换if (l=0)/最大主元为零不能计算出正确结果free(js); printf(failn);return(0);/以下为消元过程d=akk;akk=1;/这一步是否要无关紧要，只要在回代时不除以该项即可。/当前行除以akkfor (j=k+1;j=n-1;j+)akj=akj/d;bk=bk/d;/for (i=k+1;i=n-1;i+)/消元for (j=k+1;j=0;i-) t=0.0;for (j=i+1;j=0;k-)if (jsk!=k) t=bk; bk=bjsk; bjsk=t;*/free(js);return(1);double GetDeltaX(double JNN,double bN,double deltaXN)/求deltaX主程序，调用agausint nResultFlag,i;double NormDeltaX;nResultFlag=agaus(J,b,N);if(nResultFlag)NormDeltaX=0;for (i=0;iN;i+)NormDeltaX=NormDeltaX+fabs(bi);deltaXi=bi;return NormD