北航数值分析大作业.docxVIP

《数值分析》大作业（3）一、算法设计方案1．算法的总体设计根据题目中给定的关于z、t、u 的数表及方程组，可以得到关于z、x、y 的数表，以此数表作为基础，可以得到Z(x，y)的表达式。要在区域D???x,y?0?x?0.8,0.5?y?1.5?上作二元拟合函数krsp?x,y???crsxyr,s?010 20，并使p?x,y?在某一最小的k 值达到如下精度，2?????f?x,y??p?x,y???10?7，其中x?0.08i，y ?0.5?0.05j。ijijiji?0 j?0拟合节点的xi，yj 值已知，要求出相应的zij?g?tij,uij??f?xi,yj?。先解非线性方程组由一组?xi,yj?可得到一组?tij,uij?，再根据所给的数表作二元分片插值得到函数z?g?t,u?，最后将?tij,uij?代入函数z?g?t,u?得到zij?g?tij,uij??f?xi,yj?。用到的算法有：求解非线性方程组，二元二次分片插值，二元函数拟合。2．算法的具体实现（1）解非线性方程组题中的非线性方程组采用Newton 法是收敛的，且有较快的收敛速度。在x，y 已知的情况下，方程组是关于t，u，v，w的非线性方程组。求解时近似解向量x?k?满足条件x?k??x?k?1?x?k????10?12。即对向量形式为?f1(x1,x2,…,xn)?0F(x)?0，一般形式为…,?的非线性方程组，??f2(x1,x2,xn)0?????fn(x1,x2,…,xn)?0令(k)(k)(k)(k)(k)h?(h1,h2,?,hn),hj?0,j?1,2,?n1 1 1 1 f(x(k)?h(k)e)?f(x(k))h(k)1f(x(k)?h(k)e)?f(x(k))h1 n 1 1 ?(k)nJ(x(k),h(k))???f(x(k)?h(k)e)?f(x(k))f(x(k)?h(k)e)?f(x(k))n1 1 1 h(k)1nn 1 n h?(k)n其中ej是第j 个n 维基本单位向量。迭代公式x(k?1) ?x(k)?J(x(k),h(k))?1F(x(k)),k?0,1,……（2）二元二次分片插值根据数表作二元二次分片插值得到函数z?g?t,u?。分片插值能保证收敛性，且有较高的精度。根据给定的数表，可将整个插值区域分成16 个小的区域。分片插值最终要求zij ?g?tij,uij??f?xi,yj?，故先判断?tij,uij?所在的区域，再作此区域的插值，计算zij。（3）二元函数拟合拟合节点为(x,y,z),i?0,1?10,j?0,1?20。基函数为乘积型的基函ijijrs数{?(x)?(y)(r?0...k;s ?0...k)}，其中?(x)?xr,?(y)?ys 。拟合函数为rskkp(x,y)???crs?r (x)?s(y)，其中C?crs为要求的系数矩阵。实际上就是先s?0r?0求解系数C的矩阵：C?(BTB)?1BTUG(GTG)?1其中：kk?1x0?x0??1y0?y0?B??1x1?x1?,G??1y1?y1 ?,U?(f(x,y))????????????ijn?m????1x?xk1y?yk?nn??nn?令A?(BTB)?1BTU，DT ?G(GTG)?1，作变换得BTBA?BTU ，GTGD?GT ，可通过列主元素Gauss 法解得矩阵A和D，再算得C?ADT 。计算时的k 值需要程序自动确定，要求最小的k 值使精度达到：10 20iiii?=??(f(x,y)?p(x,y))2 ?10?7。i=0 j?0二、计算结果1．二元二次分片插值得到数表：?i?0,1,?,10;j?0,1,?,20?。xi，yj ，f?xi,yj?，2．二元拟合选择过程的k，?值。3 ．达到精度要求的k ，?的值及crs?r?0,1,?,k;s?0,1,?,k?。p?x,y?中的系数4．数表：*，y* ，f?x*,y*?，p?x*,y*?，?i?1,2,?,8;j?1,2,?,5?。xijijij附录全部源程序#includeiostream.h#includestdio.h#includemath.hdoublefabsmax(doublex[4]);intfabsmax_flag(double x[4],intk);voidgauss(double a_0[4][4],double b[4],double deltx[4]);voidfunction(double t_1[11][21],double u_1[11][21],double x_0[4]);voidchazhi(doublet_1[11][21],doubleu_1[11][21],doublez[11][21],double t_0