数值分析大作业三课案.docVIP

《数值分析B》 计算实习 题目三 学号：ZY1207227 姓名：丁少行 2012-11-25 目录 一．算法设计方案 1 1.1 根据x、y求解u、t 1 1.2 通过已知的二维数表，进行分片二次插值 1 1.3 进行曲面拟合和求系数矩阵 2 1.4 计算k值及相应的系数矩阵 2 二．源程序及注释 3 三、运行结果 24 3.1 数表 24 3.2 k的选择过程 32 3.3 系数 33 3.4 数表 33 一．算法设计方案 算法共分为四步，首先是求解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 PiecewiseSecondOrderInterpolation(double x[], double y[], double z[], int n, int m, double u, double v); //拟牛顿法中方程组表达式 void NewtonMethodfunc(double x[], double y[], double b[]); //曲面拟合 void SurfaceFitting(double x[], double y[], double z[], int n, int m, double a[], int ff, int q, double dt); int main() { int it = 0; int i = 0, j = 0; double b[4]; double z[11*21]; double x[11], y[21]; double t[11][21], u[11][21]; double tuvw[4] = {1.0, 1.0, 1.0, 1.0}; //迭代初始值 f