数值分析大作业三教程.docVIP

《数值分析B》 计算实习 题目三 学号：ZY1207227 姓名：丁少行 2012-11-25  目录  TOC \o 1-3 \h \z \u  HYPERLINK \l _Toc311387121 一．算法设计方案 1  HYPERLINK \l _Toc311387122 1.1 根据x、y求解u、t 1  HYPERLINK \l _Toc311387123 1.2 通过已知的二维数表，进行分片二次插值 1  HYPERLINK \l _Toc311387124 1.3 进行曲面拟合和求系数矩阵 2  HYPERLINK \l _Toc311387125 1.4 计算k值及相应的系数矩阵 2  HYPERLINK \l _Toc311387126 二．源程序及注释 3  HYPERLINK \l _Toc311387127 三、运行结果 24  HYPERLINK \l _Toc311387128 3.1 数表 24  HYPERLINK \l _Toc311387128 3.2 k的选择过程 32  HYPERLINK \l _Toc311387129 3.3 系数  PAGEREF _Toc311387129 \h 33  HYPERLINK \l _Toc311387130 3.4 数表  PAGEREF _Toc311387130 \h 33 PAGE \* MERGEFORMAT35 一．算法设计方案 算法共分为四步，首先是求解u、t，然后进行插值运算，再进行曲面拟合并求解系数矩阵C，最后计算k值及系数矩阵。 1.1 根据x、y求解u、t 将和代入非线性方程组 中，用Newton法解出和； 1.2 通过已知的二维数表，进行分片二次插值 1)选择（m，n）满足 如果或，则m=1或4； 如果或，则n=1或n=4。 2)u和t的九个点确定一个分片区域，即选为插值节点，相应的Lagrange形式的插值多项式为 其中 (k=m-1, m, m+1) (r=n-1, n, n+1) 并将和代入，便得到了数表。 1.3 进行曲面拟合和求系数矩阵 将拟合节点代入非线性方程组，解出中间变量。将对应的中间变量值代入分片插值多项式中，找出与的对应关系。即。选取拟合基函数为，先设定k=2, 将拟合节点代入拟合基函数，计算出矩阵，其中 的列数根据拟合精度要求来自动确定。 其中 , 因为为正定矩阵，所以矩阵求逆使用列主元素高斯消去法（初等行变换）将其变换为单位阵来求得。计算，与相应元素平方和即为拟和精度。 1.4 计算k值及相应的系数矩阵 根据另外给出的节点值来??察逼近效果。k从0开始逐渐增大，当，得到满足精度要求的系数和相应的k。 二．源程序及注释 #include math.h #include string #include iomanip #include fstream #include iostream using namespace std; double T_val[6]={0.0,0.2,0.4,0.6,0.8,1.0}; //初始化t值 double U_val[6]={0.0,0.4,0.8,1.2,1.6,2.0}; //初始化u值 double Z_val[6*6]={-0.50,-0.34, 0.14, 0.94, 2.06, 3.50, //对应的z的值 -0.42,-0.50,-0.26, 0.30, 1.18, 2.38, -0.18,-0.50,-0.50,-0.18, 0.46, 1.42, 0.22,-0.34,-0.58,-0.50,-0.10, 0.62, 0.78,-0.02,-0.50,-0.66,-0.50,-0.02, 1.50, 0.46,-0.26,-0.66,-0.74,-0.50}; //解线性方程组 int SolveLinearQuatation(int n, double a[], double b[]); //求解非线性方程组的一组实根 int SolveNonlinearEquation(int n, double eps, double t, double h, double x[], double r[], int k); //计算多项式的值 double Calculation(double c[], int ff, int q, double x, double y); //分片二次插值求f(x[i], y[i]) double PiecewiseSecondOrderInt