第六章 边值问题差分法简介.docVIP

第六章 边值问题差分法简介 我们考虑如下简单边值问题。 其中，为已知函数，（为了解方程收敛）。 为已知常数。 我们要求函数，实质是关心两个数域之间的对应关系： 所以我们知道一系列对应的就足够了， 所以这类问题解决思路是将方程离散化，转化为解一个现行方程组。 首先，用分点 将区间n等分，h为步长，称为节点。 其次在[a，b]内每个内部节点 上用数值微分公式。 替代原方程中的二阶导数得在节点满足的关系是式： 注意：上式是关于未知数为的线性方程。 我们有n-1个这样的方程，组成方程组。 未知数，共n+1个。 再加上边值条件，n+1个未知数，n+1个方程，方程组封闭。 解此线性方程组，得到问题的数值解。 方程组为： 矩阵形式： 由于 所以，对角占优矩阵，用追赶法求解。 解决实际问题， 首先，用分点 将区间n等分； 再直接计算得到线性方程组出发，解之。 例：一维无限深势阱 边值问题为： 令，则有： 设a=1 故区间为[-1，1] 离散化 N=10，步长 即 常数 * 非零解= ** 解出近似本征值代入*，解出yi