数值分析实习作业三.docVIP

数值分析大作业 第三题 学号： 姓名： 学院： 授课教师：吕淑娟 2016.11.21 题目：关于x, y, t, u, v, w的下列方程组（A.3） (A.3) 以及关于z, t, u的下列二维数表 z u t 0 0.4 0.8 1.2 1.6 2.0 0 -0.5 -0.34 0.14 0.94 2.06 3.5 0.2 -0.42 -0.5 -0.26 0.3 1.18 2.38 0.4 -0.18 -0.5 -0.5 -0.18 0.46 1.42 0.6 0.22 -0.34 -0.58 -0.5 -0.1 0.62 0.8 0.78 -0.02 -0.5 -0.66 -0.5 -0.02 1.0 1.5 0.46 -0.26 -0.66 -0.74 -0.5 确定了一个二元函数z = f(x, y)。 1、试用数值方法求出f(x, y)在区域 上的一个近似表达式 要求以最小的k值达到以下的精度 其中，。 2、计算（i = 1, 2, …，8；j = 1, 2,…，5）的值，以观察逼近的效果，其中，=0.1i, =0.5+0.2j。 说明： 1、用迭代方法求解非线性方程组时，要求近似解向量满足 2、作二元插值时，要使用分片二次代数插值。 3、要由程序自动确定最小的k值。 4、打印以下内容： （1）算法的设计方案； （2）全部源程序； （3）数表： （i = 0,1,2,…,10；j = 0,1,2,…,20）； （4）选择过程的值； （5）达到精度要求时的值以及中的系数（r = 0,1,…,k；s = 0,1,…,k）； （6）数表：（i = 1, 2, …,8；j = 1, 2,…,5）； （7）讨论：计算过程中遇到的现象。 5、采用f型输出的准确值，其余实型数采用e型输出并且至少显示12位有效数字。 算法设计方案 方程组(A.3)由四个方程组成，又为六元方程组，则给x,y赋一组值，则能解出一组（t,u,v,w）。依次代入xi=0.08i，yj=0.5+0.05j (其中i=0,1,…,10; j=0,1,…,20)， 利用二维数表进行分片二次代数插值，得到函数z=z(t,u)。对于区域D内任意x,y，代入方程组（A.3）解得(t,u)，再代入函数z=z(t,u)，此从x,y到z的映射即为二元函数 z=f(x,y)。分片二次代数插值具体算法见教材 将1中得到的(tij,uij)代入函数z=z(t,u)，得到zij=f(xi,yj) 同1、2步，代入=0.1i, =0.5+0.2j，（i = 1, 2, …，8；j= 1, 2,…，5），得到；将代入 ，得到。比较两函数值，观察逼近的效果。 结果 1、数表： （i = 0,1,2,…,10；j = 0,1,2,…,20） 2、选择过程的值 3、最终k应为5，此时crs矩阵如下 4、数表：（i = 1, 2, …,8；j = 1, 2,…,5） 遇到的问题 从插值得到的f(x,y)函数开始（包括其在内），所得到计算结果与课本答案存在1e-9左右的差异。经验证，并不是迭代初值的选择所导致的。 源程序 #include stdafx.h #include stdio.h #include math.h #include stdlib.h #define N 4 #define L 11 #define M 21 int _tmain(int argc, _TCHAR* argv[]) { void ntsl(double x, double y, double a[N]); //Newton迭代解方程组函数 void gsslm(double a[M][M], double b[M], double x[M], int n); void gssln(double a[N][N], double b[N], double x[N]); //列主元gauss迭代 void fcfz(double a[N], double f[N], double x, double y); //方程组赋值函数 void fcds(double a[N], double ff[N][N]); //方程导数矩阵函数 double amul(double a[N], double b[N]); //向量内积函数 double fpcz(double t0[3], double u0[3], double z0[3][3], double t, double u); //二次分片插值函数 void atra(double a[M][M], int n, int m, double at[M]