数值分析第三次作业报告.docVIP

数值分析计算实习作业三 解题 根据题目要求，需要对进行插值，但是题目中只给出了插值区域，并未明确给出节点值，而是给出了一个四元非线性方程组，通过该方程组和给定的x，y值解出该方程组的u，v的值，再根据给定的二维数表，查出的值，所以程序第一步应该是求解非线性方程组； 根据x，y的不同选值，通过解非线性方程组得到了11*21组（t，u）的值，根据这些值并不能通过给定的二维数表得到f的值，所以还要进行插值计算。根据题目要求，采用分片二次代数插值的二元插值方法。 得到数表后，可以通过最小二乘法曲面拟合的方式对节点进行拟合。题目中要求程序自动确定最小的k值，这里需要注意的是，k值得选取是并不是任意的。当k=0时，得到的插值多项式是一个常数，一般不能满足精度要求。所以k的选值可以从1开始逐次递增，但k最大不能超过n-1(n为插值节点的个数)；在本题中，对于x，n=11;对于y，n=21,所以k的最大值为10。若k超过10，还找不到满足精度要求的拟合曲面，则拟合失败。 算法设计方案 求解非线性方程组（函数 Newton()） 求解非线性方程组的方法有简单迭代法、Newton法、离散Newton法。其中简单迭代法构造迭代公式非常困难，而Newton法对于初始向量的要求比较苛刻，所以采用离散Newton法比较可靠。其算法如下： 选取初始向量，给定精度水平ε； 循环体内： 选取向量h； 计算F(x)和