对流迎风matlab,16.2.1.1 迎风格式

16.2.1.1  迎风格式

迎风格式的形式如下所示：

其中 为时间步长，h为空间步长。

以 为例，当用迎风格式求解对流方程时，在计算求解区间的左端点处的下一个时刻的函数值时，要用到左端点的左边一个节点的值，因此必须向左延拓一个节点，才能计算下一个时刻的左端点的函数值，如此得出，M个时间步的迎风格式，应向左延拓M个节点函数值，这点在程序中可以看出。

在MATLAB中编程实现的迎风格式的函数为：peHypbYF

功能：用迎风格式解对流方程

调用格式：u = peHypbYF (a,dt,n,minx,maxx,M)

其中，a：方程中的常数；

dt：时间步长；

n：空间节点个数；

minx：求解区间的左端；

maxx：求解区间的右端；

M：时间步的个数；

u：求解区间上的数值解。