实习内容(二)
要求:对比分析以下三种时间积分格式的稳定性。
1、空间采用中央差格式,时间采用欧拉格式:修改初始值或时间步长等并绘图展示试验结果,理解该格式的稳定性。(用实习一的程序)
2、空间采用中央差格式,时间采用梯形(隐式)格式:修改初始值或时间步长等并绘图展示试验结果,理解该格式的稳定性。
3、空间采用中央差格式,时间采用中央差格式:修改初始场或时间步长等并绘图展示试验结果,理解该格式的稳定性。
其中第一步计算初值可分别用下面三种格式计算,并绘图展示试验结果差异。
- 欧拉格式(A2=1,B2=0)
- 欧拉后差格式(A2=0,B2=1)
- 赫恩格式(A2=0.5,B2=0.5)
修改部分,添加地址,便于数据存储
修改时间步长,来看各自的稳定度;修改A1,B1的值,改变格式,图中所示为中央差格式
- 空间采用中央差格式,时间采用欧拉格式
修改时间步长并绘图展示试验结果,理解该格式的稳定性。
理解与分析:空间采用中央差格式,时间采用欧拉格式,仅改变时间步长,该格式下,差分解的误差不随时间增加,说明该格式下较稳定。
- 空间采用中央差格式,时间采用梯形(隐式)格式
时间步长并绘图展示试验结果,理解该格式的稳定性。
理解与分析:空间采用中央差格式,时间采用梯形格式,仅改变时间步长,该格式下,差分解的误差不随时间增加,说明该格式下较稳定。
- 空间采用中央差格式,时间采用中央差格式
时间步长并绘图展示试验结果,理解该格式的稳定性。
理解与分析:空间采用中央差格式,时间采用中央差格式,仅改变时间步长,该格式下,差分解的误差不随时间增加,说明该格式下较稳定。