二阶导数的中心差商.pdf

维普资讯 二 }敬， 中、差南 耘值简方， 1蚋9年 4月 四川太学糊 然#狮 靠住 斟崭 第 36喜第 2期 Jollr~1a[0ISichuanUniversity (Natural~icnecEdition) Vo136 No2 文章编号 ：0490-6756(1999)一02—0359一 7／州 二阶导数的中心差商 (二)z／。 向晓林 (四Jlj再 65) 二()ZZ 若 ，0)：D[ 一 R，， )E (D)．在工程技术 中，经 常需要计算 ， )的一 阶导数和二 阶 导数．，()的二阶导数 ，即海色矩 阵的近似计算 ，也为最优化计算方法 的基础．，0)在点 E intD内的一阶导数 f(z)：D[ 一 ．在数值微分中，广，， )通常通过 ，(z)向前 (后)及 中心 差商束计算 ，它们 的截断误差分别为 rJ(^)和 O(h) ，显然 中，釜商精度更高． 对于 ，(z)的二阶导数 cz)： 一 巾，一般有 下处理方式 ： (1)在假设 ，0)的一阶导数存在的情况下 ，对 re．x)的每一分量 ，使用一阶向前 (后)、中 心差商．此时要使用 f(zJ的一阶导函数 f(z)； (2)只利用 ，0)的函数值 ，对 ，()使用I旬前 (后)差商 ，得出求 ， )的海色矩阵的向前 (后)羞商公式 ，其截断误差为 a(^) ； (3)只利用 ，“)的函数值 ，对 ，cz)使用中心差 商 ，用中心差商来计算 V f(z)． 第一种方式需要利用 ，“)的导数值 ．第二种方式虽然R利用 ， )的晒数值 ，但截断误差 的阶太低 ，第三种方式虽然 已有提及 ，但未见有给出其具体计算公式及其截断误差．本文给 出 了=所导数的只利用函数值的中心差商公式及其截断误差． 引理 设 ，如)：Dc 一 在开球 D0一 {z… 0f1 d 0}c D 上三次连续可微 则对满足 PIf d，P∈ 的任意l量P，有 。+)一，)+，【)rp+击 )P斗击 耋耋塑 尸 其中 0E t0，1)， E R． 证 明 令 “)一f(z0+ tP)，则 ()：0一 0、(0 一 f(zO， (1)一 ，(z。+ P)．埘 () 使用一维函数 的泰勒展式，展开到=阶，并令 = 1．得证． 记 AE ，A= (dI_j 。．，0)：Dc 一 R． 堕 r二 !二l 则可用 (I)式计算 ， )，使 = ．(一 I，2，… ；= I．2，… ) 定理 令 ，c)：D[ 一 在 D内= 阶连续可导 ，n (I)，‘)的三阶导数在D 卜满足李普希兹条件 ．李普希照常数为 即V ，ED，，J∈f … ， ) 收稿 日期：1907一o9一I9 维普资讯 36O 四川大学学报 (自然科学版 第 36卷 ≤ … c } 、 。“ … …’ (2)^E 勾—一步 长 ．日 E D， 十 船 + heJ． 十 船 J， + 艇j一 船，· 船 船JE D； (3) = ‘)… 其中～ 用