利用中心差分格式数值求解导数.pdfVIP

利用中心差分格式数值求解导数 目录 一、问题描述 . 错误 ! 未定义书签。 二、格式离散 . 错误 ! 未定义书签。 二阶导数中心差格式离散 错误 ! 未定义书签。 追赶法求解线性方程组简述 错误 ! 未定义书签。 计算流程图 错误 ! 未定义书签。 三、程序中主要符号和数组意义 错误 ! 未定义书签。 四、计算结果与讨论 错误 ! 未定义书签。 五、源程序 . 错误 ! 未定义书签。 一、问题描述 2 2 利用中心差分格式近似导数 d y / dx ，数值求解 2 d y sin 2x 0 x 1 2 dx y / x 0 0, y / x 1 1 步长分别取 x 0.05, 0.01, 0.001, 0.0001 二、格式离散 将 x 轴上 [0,1] 之间的线段按上述步长，等步长的离散为 n 个小段，包括端点， 共 n+1 个网格节点，示意图如下： 线段上边的数字表示 x 轴上的坐标值，线段下边的数字表示节点编号，从 0 到 n 编号。 二阶导数中心差格式离散 2 d y yi 1 2yi yi 1 sin(2x) 2 2 整理为线性方程形式 dx x 2 yi 1 2 yi yi 1 x sin(2x) 其中， x 为空间离散步长； i=1,2 ，……， n-1 包括边界条件的线性方程组如下： y0 0 边界条件 2 y0 2y1 y2 x sin(2 * 1* x) 2 yi 1 2 yi yi 1 x sin(2 * i * x) 2 yn 2 2yn 1 yn x sin(2 * (n 1) * x) yn 0 边界条件 改写成矩阵形式： Ay f 1 0 y0 f 0 1 2 1 y1 f 1 1 2 1 其中， y y ，f