数值分析第三次大作业

A2：李（一）开始（二）计算数xi0.08i,yj0.05j

（三）（xiyj）,（i0,12,,10j0,12,20）代入题目中的非线性方程组，使用牛顿迭代法，求解对应（tu）；（四）t,zf(xi,yj)11*21个数值解；ji（五）从k=0开始，逐渐增大k的值，并使用最小二乘法曲面拟合法z=f(x,y)进行拟合，得到每次的k,。当107时结束计算，否则置k+1新计算，输出拟合结果；ji（六）计算

),

x

的值并输出结果，括f与p的差值，以观察p(x,y)x*0.1i,y*0.50.2

f(x,y)的效果。其 （七）结束#include"stdafx.h"//********************调用函数doublenewton_steffenson(doublexi,doubleyj,doubleeps,int doublef0(doublet,doubleu,doublev,doublew,doublex,doubley); doublef1(doublet,doubleu,doublev,doublew,doublex,doubley);//非线性方程doublef2(doublet,doubleu,doublev,doublew,doublex,doubley);//非线性方程等doublef3(doublet,doubleu,doublev,doublew,doublex,doubley);//非线性方程等voidlua(doublelua_mat[][21],intlc);LU分解子函doublelx(doublexl,inti,intk);//xdoublely(doubleyl,inti,intk);//y基函doublezuidazhi(doublex1,doublex2,doublex3,doublex4);//求最大值函数voidgta_mat(doubleinv_matrix[][21],intr,intc);voidsolve_c(doubleinv_matrix[][21],intr,intc);//求取intarray_M[21];doublea_matrix[21][21]={0};doubledoubledoublearray_t[6]={0,0.2,0.4,0.6,0.8,1.0};doublearray_z[6][6]=doublearray_x[11];doublearray_y[21];doublevoidmain(){inti,j,k,r,s;doubleepsilon=1e-12;intM=1000;for(i=0;i<11;i++)//求解插值数表{printf("x[%d]=%lf,y[%d]=%lf,}//P(x,y)doublematrix_b[21][21]={0};doublematrix_g[21][21]={0};doublesum=0;doubledoublematrix_temp[21][21]={0};{for(j=0;j<=k;j++)//Bfor(j=0;j<=k;j++)//G求//****************************求σK******************************doublesigma=0;{} {printf("\nk的值是：%d\nσ的值是%17.12e\n\n",k,sigma);}}printf("****C\n");{}doublex[8],y[5];doublefxyxin[8][5];doublepxyxin[8][5];doubledelta[8][5];{delta[i][j]=pxyxin[i][j]-}{{ }}}//求解非线性方doublenewton_steffenson(doublexi,doubleyj,doubleeps,int{intintdouble doubledoublematrix_j[21][21]={0};doublehu=0.4,ht=0.2;doublearray_f[4];doubledelta_x[4];doubleb[4];double{//求F（x）//求向量h{}{}//求矩阵inte[4]={0};{{

{}}else{

{}}else{

{}}else{{}}} //对矩阵J进行选主元的LU分for(k1=0;k1<=n-{double}{doublesum=0;}{doublesum=0;}//ΔX求解结doubledeltasum,xsum;

{printf("%d次迭代不成功，计算停止!!!\n",M);}}//ielseif(x[0]>(array_t[4]-//jelseif(x[1]>(array_u[4]-hu/2))//f(x,y)doublepxy1=0;doublepxy2=0;doublepxy=0;{{}}returnpxy;}//*********************F（X）*********************************doublef0(doublet,doubleu,doublev,doublew,doublex,doubley){doublef0;returnf0;}doublef1(doublet,doubleu,doublev,doublew,doublex,double{doublef1;returnf1;}doublef2(doublet,doubleu,doublev,doublew,doublex,double{doublereturnf2;}doublef3(doublet,doubleu,doublev,doublew,doublex,double{doublef3;returnf3;}//A进行选主元LU分voidlua(doublelua_mat[][21],int{inti,j,t,k;intn=lc+1;doubles[21]={0};doublesmax;intfor(k=0;k<=n-{for(i=k;i<=n- {doublesum=0;}for(i=k+1;i<=n- {{}} {doubletemper;{}}if(k<(n-{for(j=k+1;j<=n-{doublesum=0;}}}}//求最大doublezuidazhi(doublex1,doublex2,doublex3,double{doubleimax;returnimax;} 日基函doublelx(doublexl,inti,int{intdoublenum=1,den=1;doublevl;for(r=i-{{}}returnvl;}doublely(doubleyl,inti,int{intdoublenum=1,den=1;doublevl;for(r=i-{{}}returnvl;}//A子函voidgta_mat(doubleinv_matrix[][21],intr,int{inti,j,k,k1,k2,s;intn=c+1;doubledoublematrix[21][21]={0};doublesum=0;double{{{}{}{for(k2=0;k2<=n-{doubletemper;}{doublesum=0;}doublesum=0;}}}//C子函voidsolve_c(doubleinv_matrix[][21],intr,int{inti,j,k,k1,k2,s;intn=c+1;doublegtat_matrix[21][21]={0};doublec_matrix[21][21]={0};doublematrix[21][21]={0};doublesum=0;double{}{}{for(k2=0;k2<=n-doubletemper;}{doublesum=0;}{doublesum=0;}}}数表x0=0.000000,y0=0.500000:f(x0,y0)=4.465040184799e-001x0=0.000000,y1=0.550000:f(x0,y1)=3.246832629274e-001x0=0.000000,y2=0.600000:f(x0,y2)=2.101596866825e-001x0=0.000000,y3=0.650000:f(x0,y3)=1.030436083159e-001x0=0.000000,y4=0.700000:f(x0,y4)=3.401895562658e-003x0=0.000000,y5=0.750000:f(x0,y5)=-8.873581363801e-002x0=0.000000,y6=0.800000:f(x0,y6)=-1.733716327497e-001x0=0.000000,y7=0.850000:f(x0,y7)=-2.505346114666e-001x0=0.000000,y8=0.900000:f(x0,y8)=-3.202765063876e-001x0=0.000000,y9=0.950000:f(x0,y9)=-3.826680661097e-001x0=0.000000,y10=1.000000:f(x0,y10)=-4.377957667384e-001x0=0.000000,y11=1.050000:f(x0,y11)=-4.857589414438e-001x0=0.000000,y12=1.100000:f(x0,y12)=-5.266672548835e-001x0=0.000000,y13=1.150000:f(x0,y13)=-5.606384797965e-001x0=0.000000,y14=1.200000:f(x0,y14)=-5.877965387677e-001x0=0.000000,y15=1.250000:f(x0,y15)=-6.082697790899e-001x0=0.000000,y16=1.300000:f(x0,y16)=-6.221894528764e-001x0=0.000000,y17=1.350000:f(x0,y17)=-6.296883781856e-001x0=0.000000,y18=1.400000:f(x0,y18)=-6.308997600028e-001x0=0.000000,y19=1.450000:f(x0,y19)=-6.259561525454e-001x0=0.000000,y20=1.500000:f(x0,y20)=-6.149885466094e-001x1=0.080000,y0=0.500000:f(x1,y0)=6.380152265102e-001x1=0.080000,y1=0.550000:f(x1,y1)=5.066117551462e-001x1=0.080000,y2=0.600000:f(x1,y2)=3.821763692772e-001x1=0.080000,y3=0.650000:f(x1,y3)=2.648634911536e-001x1=0.080000,y4=0.700000:f(x1,y4)=1.547802002848e-001x1=0.080000,y5=0.750000:f(x1,y5)=5.199268349093e-002x1=0.080000,y6=0.800000:f(x1,y6)=-4.346804020491e-002x1=0.080000,y7=0.850000:f(x1,y7)=-1.316010567885e-001x1=0.080000,y8=0.900000:f(x1,y8)=-2.124310883088e-001x1=0.080000,y9=0.950000:f(x1,y9)=-2.860045510580e-001x1=0.080000,y10=1.000000:f(x1,y10)=-3.523860789794e-001x1=0.080000,y11=1.050000:f(x1,y11)=-4.116554565222e-001x1=0.080000,y12=1.100000:f(x1,y12)=-4.639049115188e-001x1=0.080000,y13=1.150000:f(x1,y13)=-5.092367247005e-001x1=0.080000,y14=1.200000:f(x1,y14)=-5.477611179623e-001x1=0.080000,y15=1.250000:f(x1,y15)=-5.795943883391e-001x1=0.080000,y16=1.300000:f(x1,y16)=-6.048572588895e-001x1=0.080000,y18=1.400000:f(x1,y18)=-6.361682484133e-001x1=0.080000,y19=1.450000:f(x1,y19)=-6.424676566900e-001x1=0.080000,y20=1.500000:f(x1,y20)=-6.426971026996e-001x2=0.160000,y0=0.500000:f(x2,y0)=8.400813957651e-001x2=0.160000,y1=0.550000:f(x2,y1)=6.997641656726e-001x2=0.160000,y2=0.600000:f(x2,y2)=5.660614423514e-001x2=0.160000,y3=0.650000:f(x2,y3)=4.391716081175e-001x2=0.160000,y4=0.700000:f(x2,y4)=3.192421380407e-001x2=0.160000,y5=0.750000:f(x2,y5)=2.063761923874e-001x2=0.160000,y6=0.800000:f(x2,y6)=1.006385238914e-001x2=0.160000,y7=0.850000:f(x2,y7)=2.060740067835e-003x2=0.160000,y8=0.900000:f(x2,y8)=-8.935402476698e-002x2=0.160000,y9=0.950000:f(x2,y9)=-1.736269688648e-001x2=0.160000,y10=1.000000:f(x2,y10)=-2.507999561599e-001x2=0.160000,y11=1.050000:f(x2,y11)=-3.209322694446e-001x2=0.160000,y12=1.100000:f(x2,y12)=-3.840977350046e-001x2=0.160000,y13=1.150000:f(x2,y13)=-4.403821754175e-001x2=0.160000,y14=1.200000:f(x2,y14)=-4.898811523126e-001x2=0.160000,y15=1.250000:f(x2,y15)=-5.326979655338e-001x2=0.160000,y16=1.300000:f(x2,y16)=-5.689418792921e-001x2=0.160000,y17=1.350000:f(x2,y17)=-5.987265495151e-001x2=0.160000,y18=1.400000:f(x2,y18)=-6.221686297503e-001x2=0.160000,y19=1.450000:f(x2,y19)=-6.393865356972e-001x2=0.160000,y20=1.500000:f(x2,y20)=-6.504993507878e-001x3=0.240000,y0=0.500000:f(x3,y0)=1.051515091801e+000x3=0.240000,y1=0.550000:f(x3,y1)=9.029274308302e-001x3=0.240000,y2=0.600000:f(x3,y2)=7.605802668593e-001x3=0.240000,y3=0.650000:f(x3,y3)=6.247151981455e-001x3=0.240000,y4=0.700000:f(x3,y4)=4.955197560009e-001x3=0.240000,y5=0.750000:f(x3,y5)=3.731340427746e-001x3=0.240000,y6=0.800000:f(x3,y6)=2.576567488723e-001x3=0.240000,y7=0.850000:f(x3,y7)=1.491505594102e-001x3=0.240000,y8=0.900000:f(x3,y8)=4.764698677337e-002x3=0.240000,y9=0.950000:f(x3,y9)=-4.684932320146e-002x3=0.240000,y10=1.000000:f(x3,y10)=-1.343567603849e-001x3=0.240000,y11=1.050000:f(x3,y11)=-2.149133449274e-001x3=0.240000,y12=1.100000:f(x3,y12)=-2.885737006348e-001x3=0.240000,y13=1.150000:f(x3,y13)=-3.554063647857e-001x3=0.240000,y14=1.200000:f(x3,y14)=-4.154913964886e-001x3=0.240000,y15=1.250000:f(x3,y15)=-4.689182499695e-001x3=0.240000,y16=1.300000:f(x3,y16)=-5.157838831247e-001x3=0.240000,y17=1.350000:f(x3,y17)=-5.561910752001e-001x3=0.240000,y18=1.400000:f(x3,y18)=-5.902469305629e-001x3=0.240000,y20=1.500000:f(x3,y20)=-6.397468392579e-001x4=0.320000,y0=0.500000:f(x4,y0)=1.271246751481e+000x4=0.320000,y1=0.550000:f(x4,y1)=1.115002018146e+000x4=0.320000,y2=0.600000:f(x4,y2)=9.646077272154e-001x4=0.320000,y3=0.650000:f(x4,y3)=8.203473694749e-001x4=0.320000,y4=0.700000:f(x4,y4)=6.824476781794e-001x4=0.320000,y5=0.750000:f(x4,y5)=5.510852085975e-001x4=0.320000,y6=0.800000:f(x4,y6)=4.263923859018e-001x4=0.320000,y7=0.850000:f(x4,y7)=3.084629956332e-001x4=0.320000,y8=0.900000:f(x4,y8)=1.973571296919e-001x4=0.320000,y9=0.950000:f(x4,y9)=9.310562085941e-002x4=0.320000,y10=1.000000:f(x4,y10)=-4.285992234034e-003x4=0.320000,y11=1.050000:f(x4,y11)=-9.483392529689e-002x4=0.320000,y12=1.100000:f(x4,y12)=-1.785729903640e-001x4=0.320000,y13=1.150000:f(x4,y13)=-2.555537790546e-001x4=0.320000,y14=1.200000:f(x4,y14)=-3.258401501575e-001x4=0.320000,y15=1.250000:f(x4,y15)=-3.895069883634e-001x4=0.320000,y16=1.300000:f(x4,y16)=-4.466382045995e-001x4=0.320000,y17=1.350000:f(x4,y17)=-4.973249517677e-001x4=0.320000,y18=1.400000:f(x4,y18)=-5.416640326994e-001x4=0.320000,y19=1.450000:f(x4,y19)=-5.797564797951e-001x4=0.320000,y20=1.500000:f(x4,y20)=-6.117062881476e-001x5=0.400000,y0=0.500000:f(x5,y0)=1.498321052481e+000x5=0.400000,y1=0.550000:f(x5,y1)=1.334998632066e+000x5=0.400000,y2=0.600000:f(x5,y2)=1.177125123739e+000x5=0.400000,y3=0.650000:f(x5,y3)=1.025024055020e+000x5=0.400000,y4=0.700000:f(x5,y4)=8.789600231743e-001x5=0.400000,y5=0.750000:f(x5,y5)=7.391451087035e-001x5=0.400000,y6=0.800000:f(x5,y6)=6.057448714871e-001x5=0.400000,y7=0.850000:f(x5,y7)=4.788838610666e-001x5=0.400000,y8=0.900000:f(x5,y8)=3.586506258818e-001x5=0.400000,y9=0.950000:f(x5,y9)=2.451022361964e-001x5=0.400000,y10=1.000000:f(x5,y10)=1.382683509285e-001x5=0.400000,y11=1.050000:f(x5,y11)=3.815486540699e-002x5=0.400000,y12=1.100000:f(x5,y12)=-5.525282116814e-002x5=0.400000,y13=1.150000:f(x5,y13)=-1.419868808137e-001x5=0.400000,y14=1.200000:f(x5,y14)=-2.220944390959e-001x5=0.400000,y15=1.250000:f(x5,y15)=-2.956352324598e-001x5=0.400000,y16=1.300000:f(x5,y16)=-3.626795115028e-001x5=0.400000,y17=1.350000:f(x5,y17)=-4.233061642240e-001x5=0.400000,y18=1.400000:f(x5,y18)=-4.776010361325e-001x5=0.400000,y19=1.450000:f(x5,y19)=-5.256554266672e-001x5=0.400000,y20=1.500000:f(x5,y20)=-5.675647436551e-001x6=0.480000,y0=0.500000:f(x6,y0)=1.731892740382e+000x6=0.480000,y1=0.550000:f(x6,y1)=1.562034577208e+000x6=0.480000,y2=0.600000:f(x6,y2)=1.397216918208e+000x6=0.480000,y3=0.650000:f(x6,y3)=1.237801006739e+000x6=0.480000,y4=0.700000:f(x6,y4)=1.084087532678e+000x6=0.480000,y5=0.750000:f(x6,y5)=9.363227723149e-001x6=0.480000,y6=0.800000:f(x6,y6)=7.947044490537e-001x6=0.480000,y7=0.850000:f(x6,y7)=6.593871980282e-001x6=0.480000,y8=0.900000:f(x6,y8)=5.304875868399e-001x6=0.480000,y9=0.950000:f(x6,y9)=4.080886854542e-001x6=0.480000,y10=1.000000:f(x6,y10)=2.922442012295e-001x6=0.480000,y11=1.050000:f(x6,y11)=1.829822068535e-001x6=0.480000,y12=1.100000:f(x6,y12)=8.030849403543e-002x6=0.480000,y13=1.150000:f(x6,y13)=-1.579041305164e-002x6=0.480000,y14=1.200000:f(x6,y14)=-1.053445516210e-001x6=0.480000,y15=1.250000:f(x6,y15)=-1.883980906096e-001x6=0.480000,y16=1.300000:f(x6,y16)=-2.650071493189e-001x6=0.480000,y17=1.350000:f(x6,y17)=-3.352378389040e-001x6=0.480000,y18=1.400000:f(x6,y18)=-3.991645038868e-001x6=0.480000,y19=1.450000:f(x6,y19)=-4.568681433016e-001x6=0.480000,y20=1.500000:f(x6,y20)=-5.084349932782e-001x7=0.560000,y0=0.500000:f(x7,y0)=1.971221786400e+000x7=0.560000,y1=0.550000:f(x7,y1)=1.795329599501e+000x7=0.560000,y2=0.600000:f(x7,y2)=1.624067113228e+000x7=0.560000,y3=0.650000:f(x7,y3)=1.457830582708e+000x7=0.560000,y4=0.700000:f(x7,y4)=1.296954649752e+000x7=0.560000,y5=0.750000:f(x7,y5)=1.141718105447e+000x7=0.560000,y6=0.800000:f(x7,y6)=9.923495333243e-001x7=0.560000,y7=0.850000:f(x7,y7)=8.490326633294e-001x7=0.560000,y8=0.900000:f(x7,y8)=7.119113522641e-001x7=0.560000,y9=0.950000:f(x7,y9)=5.810941589219e-001x7=0.560000,y10=1.000000:f(x7,y10)=4.566585132334e-001x7=0.560000,y11=1.050000:f(x7,y11)=3.386544961394e-001x7=0.560000,y12=1.100000:f(x7,y12)=2.271082557696e-001x7=0.560000,y13=1.150000:f(x7,y13)=1.220250891932e-001x7=0.560000,y14=1.200000:f(x7,y14)=2.339221963760e-002x7=0.560000,y15=1.250000:f(x7,y15)=-6.881870197104e-002x7=0.560000,y16=1.300000:f(x7,y16)=-1.546493442129e-001x7=0.560000,y17=1.350000:f(x7,y17)=-2.341526664587e-001x7=0.560000,y18=1.400000:f(x7,y18)=-3.073910919133e-001x7=0.560000,y19=1.450000:f(x7,y19)=-3.744348623481e-001x7=0.560000,y20=1.500000:f(x7,y20)=-4.353605565359e-001x8=0.640000,y0=0.500000:f(x8,y0)=2.215667863688e+000x8=0.640000,y1=0.550000:f(x8,y1)=2.034201133607e+000x8=0.640000,y2=0.600000:f(x8,y2)=1.856955143619e+000x8=0.640000,y3=0.650000:f(x8,y3)=1.684358164161e+000x8=0.640000,y4=0.700000:f(x8,y4)=1.516776352400e+000x8=0.640000,y5=0.750000:f(x8,y5)=1.354519041151e+000x8=0.640000,y6=0.800000:f(x8,y6)=1.197844086673e+000x8=0.640000,y7=0.850000:f(x8,y7)=1.046963049419e+000x8=0.640000,y8=0.900000:f(x8,y8)=9.020460838023e-001x8=0.640000,y9=0.950000:f(x8,y9)=7.632264776629e-001x8=0.640000,y10=1.000000:f(x8,y10)=6.306048219543e-001x8=0.640000,y11=1.050000:f(x8,y11)=5.042528145972e-001x8=0.640000,y12=1.100000:f(x8,y12)=3.842167155457e-001x8=0.640000,y13=1.150000:f(x8,y13)=2.705204766410e-001x8=0.640000,y14=1.200000:f(x8,y14)=1.631685723996e-001x8=0.640000,y15=1.250000:f(x8,y15)=6.214855811676e-002x8=0.640000,y16=1.300000:f(x8,y16)=-3.256661939682e-002x8=0.640000,y17=1.350000:f(x8,y17)=-1.210165348444e-001x8=0.640000,y18=1.400000:f(x8,y18)=-2.032513996228e-001x8=0.640000,y19=1.450000:f(x8,y19)=-2.793303595584e-001x8=0.640000,y20=1.500000:f(x8,y20)=-3.493199575400e-001x9=0.720000,y0=0.500000:f(x9,y0)=2.464684222660e+000x9=0.720000,y1=0.550000:f(x9,y1)=2.278058979399e+000x9=0.720000,y2=0.600000:f(x9,y2)=2.095251250840e+000x9=0.720000,y3=0.650000:f(x9,y3)=1.916718127997e+000x9=0.720000,y4=0.700000:f(x9,y4)=1.742854628776e+000x9=0.720000,y5=0.750000:f(x9,y5)=1.573998427334e+000x9=0.720000,y6=0.800000:f(x9,y6)=1.410434835231e+000x9=0.720000,y7=0.850000:f(x9,y7)=1.252401750608e+000x9=0.720000,y8=0.900000:f(x9,y8)=1.100094409628e+000x9=0.720000,y9=0.950000:f(x9,y9)=9.536698512613e-001x9=0.720000,y10=1.000000:f(x9,y10)=8.132510552489e-001x9=0.720000,y11=1.050000:f(x9,y11)=6.789307429659e-001x9=0.720000,y12=1.100000:f(x9,y12)=5.507748485043e-001x9=0.720000,y13=1.150000:f(x9,y13)=4.288256769731e-001x9=0.720000,y14=1.200000:f(x9,y14)=3.131047717398e-001x9=0.720000,y15=1.250000:f(x9,y15)=2.036155140327e-001x9=0.720000,y16=1.300000:f(x9,y16)=1.003454782409e-001x9=0.720000,y17=1.350000:f(x9,y17)=3.268565186572e-003x9=0.720000,y18=1.400000:f(x9,y18)=-8.765306591329e-002x9=0.720000,y19=1.450000:f(x9,y19)=-1.724672478188e-001x9=0.720000,y20=1.500000:f(x9,y20)=-2.512302207523e-001x10=0.800000,y0=0.500000:f(x10,y0)=2.717811109469e+000x10=0.800000,y1=0.550000:f(x10,y1)=2.526399501256e+000x10=0.800000,y2=0.600000:f(x10,y2)=2.338411386860e+000x10=0.800000,y3=0.650000:f(x10,y3)=2.154329377280e+000x10=0.800000,y4=0.700000:x10=0.800000,y5=0.750000:f(x10,y5)=1.799510579099e+000x10=0.800000,y6=0.800000:f(x10,y6)=1.629448220554e+000x10=0.800000,y7=0.850000:f(x10,y7)=1.464650043751e+000x10=0.8