計算实习报告3.docVIP

数 值 分 析（B） 大 作 业（三） 姓名： 学号： 电话： 1、算法设计： ①整体思路 这一次的作业和前面两次相比，感觉要难很多，最明显的就是书上没有完整的算法实现流程。再加上插值和逼近这一段在课上也没听太懂什么意思，所以刚开始思考这个程序时感觉很吃力。不过经过后来好好看书和向同学请教，渐渐对插值和分片元二次拟合有点懂了，也大概知道了解决本次问题分别要完成的子问题。 要完成本次任务，主体上有三部分问题要解决： （1）解非线性方程组。将D内的)当作已知量代入题目给定的非线性方程组，得到与相对应的t[i][j],u[i][j]。 （2）分片双二次插值。采用分片双二次插值，得到与t[11][21],u[11][21]]对应的z[11][21]，也即建立与的对应关系，得到二元函数。 （3）曲面拟合。利用x[i],y[j],z[11][21]建立二维函数表，再根据精度的要求选择适当k值，并得到曲面拟合的系数矩阵C[k][k]。 （4）观察逼近效果。观察逼近效果只需要重复上面（1）和（2）的过程，得到与新的插值节点对应的，再与对应的比较即可。 程序示意图如下所示： ②各部分问题的算法实现 （1）解非线性方程组 非线性方程组常用数值解法有简单迭代法和牛顿法。牛顿法收敛快，一般都能达到平方收敛，因此这里选择New