計算方法試驗數(shù)據(jù)

1、計算方法實驗數(shù)據(jù)Exp 1#include #define M 4#define N 3#define V 3#define Y 2int main()(FILE *fp1;float AMN,BNV,CYM,DMN,EMN,FYM;int i,j,k;if(fp1=fopen(f1.txt,r)=NULL)(printf(errorn);return 0;)for(i=0; iM; i+)for(j=0; jV; j+)fscanf(fp1,%f,&Aij);FILE *fp2;if(fp2=fopen(f2.txt,r)=NULL)(printf(errorn);return 0;)for

2、(i=0; iN; i+)for(j=0; jV; j+)fscanf(fp2,%f,&Bij);FILE *r1;r1=fopen(result1.txt,w);for(i=0; iM; i+)(for(j=0; jN; j+)(Cij=0;for(k=0; kV; k+)Ci昨Cij+A皿k*Bkj;printf(%10.6f,Cij);fprintf(r1,%10.2f ,Cij);)printf(n);fprintf(r1,n);Exp 2#include #include float D(double x)return(sqrt(3+(1/x);/return(3*x+1)/(x*x

3、);return(x*x*x-1)/3);return(1/(x*x-3);)int main()double x1,x2,x3;float eps=0.000000005;scanf(%f,&x2);dox1=x2;x2=D(x1);printf(n);printf(x=%10.6fn”,x2);while(fabs(x1-x2)=eps);Exp 3.1#include #include int main()FILE *fp1;fp1=fopen(f1.txt,r);FILE *fp2;fp2=fopen(f2.txt,r);double A1515,M1515;double B15,s,

4、X15;int i,j,k=0,n=10;for(i=0; i=n-1; i+)for(j=0; j=n-1; j+)fscanf(fp1,%lf,&AijD;fscanf(fp2,%lf,&Bi);)fclose(fp1);fclose(fp2);for(k=0; k=n-2; k+)(for(i=k+1; i=n-1; i+)(Mik=Aik/Akk;for(j=k+1; j=0; i-)(s=0;for(j=i+1; j=n-1; j+)s=s+Aij*Xj;Xi=(Bi-s)/Aii;)for(i=0; i=n-1; i+)printf(%lfn,Xi);return 0;Exp 3.

5、2 #include #define M 8int main()FILE *f1;f1=fopen(1.dat,r);FILE *f2;f2=fopen(2.dat,r);double A1515,C15,D15,B1515,Y15,X15;double sum=0;int i,j,k;if (f1=NULL)printf(Open failedn);for(i=1; i=M; i+)for(j=1; j=i; j+)fscanf(f1,%lf,&Aij);Aji=Aij; fscanf(f2,%lf,&Ci);fclose(fl);fclose(f2);for(int i = 1; i =

6、M; i+)(for(int j = 1; j = M; j+)printf(%5.1lf, Aij);printf(n);)for(i=1;i=M;i+)printf(%5.1lf,Ci);printf(n);for(i=1; iM; i+)(for(j=1; j=i-1; j+)(sum=0;for(k=1; k=j-1; k+)sum=sum+(Dk*Bik*Bjk);Bij=(Aij-sum)/Dj;sum=0;for(k=1; k=i-1; k+)(sum=0;sum=sum+(Dk*Bik*Bik);)Di=Aii卜sum;)for(i=1; i=M; i+)(sum=0;for(

7、k=1; k=1; i-)(sum=0;for(k=i+1; k=M; k+)sum=sum+(Bki*Xk);Xi=Yi/Di-sum;for(i=1; i=M; i+)printf(X%d=%lfn,i,Xi);return 0;)Exp 4#include #include #define M 3#define MAX 100int main()FILE *f1;f1=fopen(data.txt,r);FILE *f2;f2=fopen(data1.txt,r);double AMM,BM,XM=0,0,0,丫M;double sum,eps=0.000000001,max;int i

8、,j,k;for(i=0; iM; i+)for(j=0; jM; j+)fscanf(f1,%lf,&Aij);)k=1;while(kMAX)(for(i=0; iM; i+)(sum=0;for(j=0; jM; j+)if(j-i)sum+=Aij*Xj;Yi=(Bi-sum)/Aii;)max=0;for(i=0; iM; i+)if(maxfabs(Xi-Yi)max=fabs(Xi卜Yi);if(maxeps)break;printf(nk=%d,k);for(i=0; iM; i+)printf(Y%d=%lf,i+1,Yi);k+;for(i=0; iM; i+)Xi=Yi；

9、)if(k=MAX)(printf(errorn);)printf(nk=%d,k);for(i=0; iM; i+)printf(Y%d=%lf,i,丫i);/*#include #define M 3#include math.hmain()(double aMM= 10,-1,-2,-1,10,-2,-1,-1,5;double bM = 7.2,8.3, 4.2 ;double s, max, eps = 0.000000001;int k, i, j, N = 100;k = 1;while( kN )for ( i=0; iM ; i+ )s = 0;for ( j=0; jM;

10、j+ )if( j-i )s += aij * xj;yi = ( bi-s) / aii;max = 0;for ( i=0; iM; i+ )if( max fabs( xi-yi)max = fabs( xi - yi);if ( max eps )break;printf( nk=%d, k );for( i=0; iM; i+ )printf( y%d=%lf,i, yi);k+;for( i=0; iM; i+ )xi = yi;)if(k = N)(printf( ERROR!n);)printf( nk=%d, k );for ( i=0; i M; i+ )printf( y

11、%d=%lf, i, y i );)*/Exp 5#include stdio.h#define M 5#define N 6float x= 0.4,0.55,0.65,0.80,0.95,1.05;float y= 0.41075,0.57815,0.69675,0.90,1.00,1.25382;float x2= 1,2,3,4,5,6,7;float y2= 0.368,0.135,0.050,0.018,0.007,0.002,0.001;float f1(float xx)int i,k;float t=1,j=0;float sum1=0,sum2=0;for(i=0; i=M

12、; i+)t=1;sum1=1;sum2=1;for(k=0; k=M; k+)if(i!=k)sum1*=xx-xk;sum2*=xi-xk;t=sum1/sum2;j+=yi*t;return j;float f2(float xx) (int i=0,k=0;float t=1,z=0;float sum1,sum2;for(i=0; i=N; i+)(t=1;sum1=1;sum2=1;for(k=0; k=N; k+)if(k!=i)(sum1*=xx-x2k;sum2*=x2i-x2k;t=t*sum1/sum2;z+=y2i*t;return z; int main()print

13、f(the answer is %fn,f1(0.596);printf(the answer is %fn,f1(0.99);printf(the answer is %fn,f2(1.8);printf(the answer is %fn,f2(6.15);Exp6#include #include #define Num 10float sNum,xNum,yNum,faiNumNum,afaNum;float beidaNum,aNum,xfaiNum,ydNum,max,pcpfh;float inside(float bNum,float cNum)(int p;float j=0

14、;for(p=1; pNum; p+)j+=cp*bp;return j; int main()int i,j,k,n,index,flag;char conti;conti=;printf(please input data length:n);scanf(%d,&n);printf(please input x and yn);for(i=1; i=n; i+)(printf(x%d=,i);scanf(%f,&xi);printf(y%d=,i);scanf(%f,&yi);while(conti=)(printf(please input compute timen);scanf(%d

15、,&index);pcpfh=0;afa1=0;a0=0;for(i=1; i=n; i+)afa1+=xi;a0+=yi;fai0i=1;)afa1=afa1/n;a0=a0/n;for(i=1; i=n; i+)(fai1i=xi-afa1;)a1=inside(fai1,y)/inside(fai1,fai1);for(k=1; kindex; k+)(for(i=1; i=n; i+)xfaii=xi*faiki;afak+1=inside(faik,xfai)/inside(faik,faik);beidak=inside(faik,faik)/inside(faik-1,faik-

16、1);for(j=1; j=n; j+)faik+1j=(xj-afak+1)*faikj-beidak*faik-1j;ak+1=inside(faik+1,y)/inside(faik+1,faik+1);printf( the result of %d computingn,index);for(i=1; i=index; i+)printf(a%d=%fn,i,ai);for(i=1; i=index; i+)printf(afa%d=%fn,i,afai);for(i=1; i=index; i+)printf(beida%d=%fn,i,beidai);for(i=1; i=n;

17、i+)for(k=0; k=index; k+)si+=ak*faiki;ydi=float(fabs(yi-si);pcpfh+=ydi*ydi;si=0;max=0;for(i=1; imax)max=ydi;flag=i;printf(當(dāng) x=%f 時，偏差最大=%f,偏差平方和為 %fn,xflag,max,pcpfh);printf(continue with space,quit with else button);conti=getchar();conti=getchar();Exp 7#include #include float C(float x)(return(exp(x)/(4+x*x);/return (log(1+x)/(1+x*x);int main()(float a=0,b=1;float h=b-a,M1,M2,s,x;int i;M1=h/2*(C(a)+C(b);for(i=0;i10;i+)printf(h=%f,T=%fn,h,M1);s=0;x=a+hwhile(xb)(s+=C(x);x+=h;M2=M1/2+h/2*s;h/=2;M1=M2;printf(h=%f,T=%fn,h,M1);return 0;Exp 8#include float C(float x,float y)return(x*x-y*y);/ret