kbhmnAAA北航数值分析第三次大作业.docVIP

kbhmnAAA北航数值分析第三次大作业 kbhmnAAA北航数值分析第三次大作业 PAGE kbhmnAAA北航数值分析第三次大作业 数值分析第三次大作业 算法的设计方案 1、求解非线性方程组 将题目中给出的当作已知量代入题目给定的非线性方程组，求出与相对应的数组te[i][j]，ue[i][j]，此处采用的是牛顿法解非线性方程组，其算法如书上91页所示。 分片二次代数插值 对所求出的数组te[i][j]，ue[i][j]，通过分片二次代数插值运算，得到与数组te[11][21]，ue[11][21]对应的数组ze[11][21]，从而得到二元函数z=，此处采用如书上101页例2中所示的分片二次代数插值。 曲面插值 利用x[11]，y[21]，ze[11][21]建立二维函数表，进行曲面插值计算，逐步提高k值，计算其精度，看其是否满足要求，以此来确定循环结束的时刻，并得到曲面拟合的系数矩阵C[r][s]，此处的算法如书142页所示，只需将所需矩阵给出，然后按公式进行计算即可。 比较 观察和逼近的效果。观察逼近效果只需要利用新给的点列重复上面（1）和（2）的过程，得到与新的插值节点对应的，再与对应的比较即可，这里求解可以直接使用（3）中的C[r][s]和k。 几点说明 分片二次插值的结果x[i]，y[j]，ze[i][j]输出到一个文件中，方便结果的复制与粘贴。 曲面插值的结果输出