![地下水数值模拟基础](https://wfqqreader-1252317822.image.myqcloud.com/cover/292/40936292/b_40936292.jpg)
上QQ阅读APP看本书,新人免费读10天
设备和账号都新为新人
一、差分的概念
任意足够光滑的函数f(x)沿x的正向和负向分别用Taylor级数展开,有
![img](https://epubservercos.yuewen.com/424EE0/21277064901840706/epubprivate/OEBPS/Images/txt002_1.jpg?sign=1739373045-FZp5FP4fPpFJ0Jvm5Qt4B0Isxy262JtR-0-e59ef2ed6e2abfd6327ebd7645d0964d)
于是有
![img](https://epubservercos.yuewen.com/424EE0/21277064901840706/epubprivate/OEBPS/Images/txt002_2.jpg?sign=1739373045-YSvfavZDovM8gEewqVLO9Hy4QYCirQKk-0-bd825e3fb7fd4e612e17a69e48404617)
或
![img](https://epubservercos.yuewen.com/424EE0/21277064901840706/epubprivate/OEBPS/Images/txt002_3.jpg?sign=1739373045-AHKwoFVim50uOjHgvImcW4LhYZ9qXGxI-0-16904f664c2a61d5fe4d0e420ea85b87)
由此可给出一阶导数的近似表达式
![img](https://epubservercos.yuewen.com/424EE0/21277064901840706/epubprivate/OEBPS/Images/txt002_4.jpg?sign=1739373045-KQdlJH1tzkjkuSjZlyn6TD0MyraIfVGC-0-2d69736c6f551b0d08d568228f312987)
近似表达式把级数截断了,由此产生的误差称为截断误差E,可以用被截断级数的第一项也是最大的一项来表示
![img](https://epubservercos.yuewen.com/424EE0/21277064901840706/epubprivate/OEBPS/Images/txt002_5.jpg?sign=1739373045-jz1IKb85F2UPzaG2nDMSL7qrfDcnPDR1-0-15aaab75a382f911d36c0758a9c6ffa6)
此误差O(Δx)与Δx同阶。对于足够小的Δx来说,误差O(Δx)的绝对值将小于cΔx(式中c为任意常数)。所以为了保证采用差分近似导数时的误差足够小,必须采用足够小的Δx。
由式(2-1a)和式(2-1b)相加,可得f(x)一阶导数的“中心差”近似表达式
![img](https://epubservercos.yuewen.com/424EE0/21277064901840706/epubprivate/OEBPS/Images/txt002_6.jpg?sign=1739373045-fG44Zd84SFTTW8XUw2uHXL7wBYhAzSrr-0-05ad95c8e723fb60afb9eddf45f4f8a1)
被截去的第一项为所以式(2-3)的截断误差为O[(Δx)2]。对于足够小的Δx来说,它比前差和后差有更高的近似,所以数值近似的精度不仅和Δx的大小有关,还和导数取何种差分形式有关。
式(2-1a)和式(2-1b)相减,可得f(x)二阶导数的近似表达式
![img](https://epubservercos.yuewen.com/424EE0/21277064901840706/epubprivate/OEBPS/Images/txt002_8.jpg?sign=1739373045-anfYgaLk1sT3VBzy6MUE59hDOmA1GIFa-0-fb03369de2eae3144a6928de4539a125)
被截去的第一项为截断误差为O[(Δx)2]。