一维对流方程在三种差分格式下解.pdf

一维对流方程在三种差分格式下的数值解 一、实验目的 ∆ 用数值方法计算一维对流方程在A 、B、C 三种差分格式下的解。∆x 取为0.05 ， 取值 ∆ 为0.5,1,2.并作相关讨论 + = 0, −8 ≤ ≤ 8 , 0 ( ) , 0 = 1 + , −1 ≤ ≤ 0 ( ) , 0 = 1 − , 0 ≤ ≤ 1 ( ) , 0 = 0, −1 或 1 二、实验要求 1. 学会对MS-FORTRAN 的基本操作 2. 用Fortran 编写程序计算一维对流方程在A 、B、C 三种格式下的解。 ∆ 3. 讨论各种格式下不同的 值的差分格式解的特点。 ∆ 三、实验步骤 1.写出一维对流方程的三种格式的差分方程： A 格式： +1 = + ∆ ( − ) +1 −1 2∆ B 格式： +1 = + ∆ ( − ) +1 ∆ C 格式： +1 = + ∆ ( − ) −1 ∆ 2.根据上述差分方程，写出求解的fortran 程序 （附件1,2）; 3.根据fortran 程序求出三种格式下的解; ∆ 4. 比较各种格式下不同的 的数值解，讨论它们的特点。 ∆ 四、实验结果分析 利用上述fortran 程序的计算结果，可以观察计算数值，初步判断差分格式的收敛性和稳定 性；然后利用EXCEL 或MATLAB 进行绘图，更形象地体现计算结果。以下为利用计算数据在MATLAB 所绘制的图 （为了使图更清晰，对比更明显，只画了x ∈[-3,3] 的图）： ∆ ∆ 1. = 0.5( =