科学与工程计算.docVIP

（10分）考虑下面对流-扩散方程 请写出一个无条件稳定的二阶精度的差分格式。 Crank Nicholson 隐格式: 精度：O(τ2+h2) 稳定性：，其中。由，隐格式是无条件稳定 （10分）考虑下面对流方程 请写出迎风格式。 或者 （15分）讨论扩散方程 的差分格式 的稳定性。 Fourier分析方法讨论其稳定性： 差分格式改写为： ， 其中 令，代入上式有 因此，增长因子为 由von Neumann条件有 ，即，解之得，，又因，不满足，所以此差分格式不稳定。 （15分）考虑对流方程 的差分格式 试确定差分格式中的系数，使得此格式具有尽可能高的精度。 带入上式，整理得：，解之得 此时，余项为 该格式的精度为 五、（30分）考虑下面两点边值问题 （1）（8分）写出Galerkin变分形式与Ritz变分形式。 方程两边同乘以函数v(x)，并在[0，1]上积分得： 左边分部积分有： 若令v(0)=0,根据边界条件有 记, 则Galerkin变分问题是 记 Ritz变分问题是 （2）（7分）在二次函数空间中，计算变分问题的近似解。 一般二次函数为，根