江苏师范大学数值分析实验报告结果

实验一:拉格朗日插值法#include void main()static float Lx10,Ly10;intn,i,j;floatx,y=0,p;printf(enter n=);输入点的个数scanf(%d,&n);printf(enter xin);输入点对应的x的值for(i=1;i=n;i+)scanf(%f,&Lxi);printf(enter yin);输入点对应的y的制for(i=1;i=n;i+)scanf(%f,&Lyi);printf(enter x=);输入验证的x的值scanf(%f,&x);for(i=1;i=n;i+)p=1;for(j=1;j=n;j+)if(i!=j)p=p*(x-Lxj)/(Lxi-Lxj);y=y+p*Lyi;printf(y=%fn,y);/入验证的x对应的y的值n=bxi019E0.401 0.Benter0 0.1entepX=0.15 D:VC6.0C+ +6.DM crosoft Vhual Stu d ioMyProjectstS H喑蛆日聲1 启entep enter 0 0.1 enter 0,39894 0.39695 0.39142 0.38138 0.36812 9.35206 entei* x=0Bi 15y=0.394473Press anj key to continueX=0.31* DtVC6 .JC - + 6.DMi crosoft /isusl Stu d io MyProject5fi.T IS BBn=bxi0.195 0.3 B.4B1 0.5yiB-39B94 0.39695 8.39142 0.3S1.38 0.36812 S_352(3t &ntep x=0-31y=0.380219Ppess mny ke to continue实验二：最小二乘法#include stdio.hfloatgs(float a2020,float b20,int n )inti,j,k,l;float s;k=1;while(k!=n+1)if(akk!=0)for(i=k+1;i=n+1;i+)aik=aik/akk;bi=bi-aik*bk;for(j=k+1;j=1;k-)s=0;for(l=k+1;l=n+1;l+)s=s+akl*bl;bk=(bk-s)/akk;return 0;int main()float a2020=0.0;定义 a 矩阵float c2020;定义 c 矩阵float ct2020;定义 ct 矩阵 float刈20;/定义数组用于存放x的数据 float y20;/定义数组用于存放y的数据 float b20=0.0;定义 b 矩阵inti,j,k,m,n;printf（输入所求函数的最高次数n:n）;输入n（求线性的函数输入1。）scanf(%d,&n);printf(输入测试数据的组数m:n);输入测试数据的组数 scanf(%d,&m);printf(输入x的测试数据d个:n,m);输入x的测试数据m个 for(i=1;i=m;i+) scanf(%f,&xi);printf(输入y的测试数据d个:n,m);输入y的测试数据m个 for(i=1;i=m;i+)scanf(%f,&yi);for(i=1;i=m;i+)/c 矩阵第一列赋值为 1ci1=1.0;/求 Cfor(j=2;j=n+1;j+)for(i=1;i=m;i+)cij=xi*cij-1;/输出 Cprintf(C 矩阵如下：n);for(i=1;i=m;i+)for(j=1;j=n+1;j+)printf(%f ,cij);if(j=n+1) printf(n);求c的转置矩阵CT for(i=1;i=m;i+)for(j=1;j=n+1;j+)ctji=cij;/输出 CTprintf(CT 矩阵如下：n);for(i=1;i=n+1;i+)for(j=1;j=m;j+)printf(%f ,ctij);if(j=m)printf(n);/求 afor(i=1;i=n+1;i+)for(j=1;j=n+1;j+)for(k=1;k=m;k+) aij+=ctik*ckj;/输出 aprintf(a 矩阵如下：n);for(i=1;i=n+1;i+)for(j=1;j=n+1;j+)printf(%f ,aij);if(j=n+1)printf(n);/求 bfor(i=1;i=n+1;i+)for(k=1;k=m;k+)bi+=ctik*yk;/输出 bprintf(b 矩阵如下： n);for(i=1;i=n+1;i+)printf(%f ,bi);gs(a,b, n);调用高斯函数求方程组的解printf(nn);pri ntf(输出求得的函数的系数为:n); for(i=1;i=n+1;i+)输出求得的函数的系数 pri ntf(a%d=%f ,i-1,bi); printf(nn);return 0;i0123456Ti (C)19.125.030.136.040.045.150.0Ri()76.3077.8079.2580.8082.3583.9085.10实验三：数值积分：定步长复合梯形#include stdio.hfloat f(float x)return 4/(1+x*x);int main()inti,n;floata,b,h,xi;float T=0.0; printf(input a,b,n:n); scanf(%f%f%d,&a,&b,&n); T=T+f(a)+f(b);h=(b-a)/n;xi=a;for(i=1;i=0.000001)n=n*2;Tn=T2n; T2n=Tn/2.0;for(k=1;k=n;k+) T2n=T2n+f(a+(2*k-1)*(b-a)/(2*n)*(b-a)/(2*n);printf(T%d=%lfn,2*n,T2n);printf(n满足条件的结果为：lfn,T2n);return 0;实验四非线性方程求解#include stdio.h #include math.h double f(double x) returnpow(x,6)-x-1;double f_(double x)return 6*pow(x,5)-1; int main()double x0,x,eps;inti,k,N;printf(请输入初始x0的值：n);scanf(%lf,&x0);printf(迭代次数：n);scanf(%d,&N);printf(请输入误差e： n); scanf(%lf,&eps);printf(ktxkn);printf(0t%.8lfn,x0);x=x0-f(x0)/f_(x0);if(!(f_(x0)=0)for(k=1;k=N;k+)printf(%dt%.8lfn,k,x); if(fabs(x-x0)N)i=2;elsei=1;if(i=0) printf(nx=%.8lfn,x); printf(f(x)=%.30lfn,f(x); printf(k=%dn,k);else if(i=1)printf(给定的x0导致计算中断!n);else if(i=2)printf(迭代N次后仍不满足精度要求!n);return 0;严舞囂#inc-ude =stdioh=doub-e f(doub-e xi、doub-eyi)宀return yrint main()宀int Fl- doub-e a、b、yssssspcccc二aaaa二nnnn2hhhh*-input a、b、n、y9rrr%_f=20ar%_f=20br%d=20nr%_f=20yordoub-e hu?a- doub-e x-ii3 doub-eyliyg pr5tf(=xit 埒甘姦HyM=r pr5tf(=%_ft%_fn=、xi、yirintdoub-e xilyiljfor(lilxun=+)宀xiluxi+hj yiluyi+h*f(xi、yirx-iixiljy-iiyEprintfrfxtfxrrXLyir returns一耳巧忙穴、bFri 1*011SIH10 1 劈is講%昌1 000060E 9SB830sIBB 蛊 0 肿 i00g ss 宦国蛊B 500038 黒SSSB 誇00當 08ASS0 P290S0 国 6561SI0 .sss 4出s 4-JIXI2百P s ss 国3”M6P曲 no cQnrtia匚 e2.改进欧拉法#include stdio.hdouble f(double xi,doubleyi)return -yi;int main()int n;double a,b,y0;printf(input a,b,n,y0:n);scanf(%lf,&a);scanf(%lf,&b);scanf(%d,&n);scanf(%lf,&y0);double h=(b-a)/n;double xi=a;doubleyi=y0;printf(xit 改进的欧拉格式 yin);printf(%lft%lfn,xi,yi);int i;double xi1,yi1,yp,yc;for(i=1;i=n;i+)xi1=xi+h;yp=yi+h*f(xi,yi);yc=yi+h*f(xi1,yp);yi1=(1.0/2.0)*(yp+yc);xi=xi1;yi=yi1; printf(%lft%lfn,xi,yi);return 0;实验六高斯消元法#include stdio.h#include math.hint main()double A1010,b10;double m10;doubleeps=0.000001;inti,j,n,k,ik;printf(输入矩阵的阶数n： n);scanf(%d,&n);printf(输入矩阵 A:n);for(i=1;i=n;i+)for(j=1;j=n;j+) scanf(%lf,&Aij);printf(输入矩阵 b:n);for(i=1;i=n;i+)scanf(%lf,&bi);double max;double temp;for(k=1;kn;k+)max=0;for(i=k;i=n;i+)if(max-fabs(Aik)0) max=Aik;ik=i;if(k!=ik)for(j=k;j=n;j+)temp=Akj;Akj=Aikj;Aikj=temp;temp=bk;bk=bik; bik=temp;for(i=k+1;i=n;i+) mi=Aik/Akk; for(j=k;j=n;j+)Aij=Aij-Akj*mi;for(j=k+1;j=n;j+) bj=bj-bk*mj;printf(A上三角矩阵：n);for(i=1;i=n;i+) for(j=1;j=n;j+)printf(%12lf,Aij);if(j=n) printf(n);printf(b 矩阵:n);for(i=1;i=1;i-)for(j=n;ji;j-)bi=(bi-Aij*bj);bi=bi/Aii;printf(n求得得结果为：n);for(i=1;i=n;i+)printf(x%d=%lft,i,bi);return 0;ii4_0湘0胴x2=-7.0(dW(ll!l(!lm3Press smy key to continue1.000000B上倔7O.0O600B6.096096 -1.SaGQBG -0.596006输八无阻的阴数n：3 16 2 13 1 输入矩阵h：2 ? 4A丄二角矩阵:3.0OS0OQ0.0OS0O6b矩阵:2 .nnpmnn求得得结果対:xi =1 .MUidMUtdI D:VC6.0C +6.0M; crosoft Visual Stu d iaMyP rectsVraWf jTGDebijgXMWf JTD.ee 1=1 丨回