计算机科学 泰勒级数,一阶常微分方程泰勒级数解法的计算机实现.pdf

一阶常微分方程泰勒级数解法的计算机实现

第2l卷 第 西 安 冶 金 建 筑 学 院 学 报 V01．2lNo．4

1989年 12月 J．Xian Inst．ofMetal1．＆ Cons．Eng． Dec． 19B9

一 阶常微分方程泰勒级数解法的计算机实现

陈 泽 宏

摘 要

术文在微分运算蒋确解算法杓基硎上提出≈ 算公式和算法，盐一所常微舒方程的泰 勒级数

解法得 在计算讥上实现，为求解一 常微箭方程提供了一种新的比较可靠的缸值方法。文 中还

给出了几个例子，把泰勒圾敛解法与龙格一胄塔法的解题络果进行了比较。

关键词：人工智能}计算方法；微分方程；泰勒缎数 ／微分运 捧精确解；

Turbo—prolog语 言

前 言

在很多工程领域中的物理定律是借助微分方程表述的，而能用手算方法满意处理的微分

方程数量是相当有限的。因此，计算机解微分方程就成为一个十分重要 昀课题 。

目前，在计算机上求一阶常微分方程 =，(，)满足初值条件 ( )= 。的辩 值解，通

常采用龙格～库塔等方法。从理论上讲，把 ()在各结点以泰勒级数 展 开，即司求‘得它

在下一个结点的数值解 。但是，在泰勒级数展开式

= ()= ( 。) ，(。。)(— )+ ( 一 。) +

+ ( 一 。) +… (1)

中除 0o)和 (o)是 已知的外，要求 出各阶导数 ， …并且要计算出它们在点 。处的

值。这是实现泰勒级效解法的关键。尽管可以用差分代替微分，但是产生的氓差将壤其失去

实用性，因此很少器用。由于这个问题的存在，使泰勒级数解法至今未能实用化。

本文收到 日期l1989—08—2a

西 安 冶 盎 建 筑 学 院 学 报 第2l卷

1 公式推导

本文作者应用人工智能中的符号微分技术提出的徽分运算精 确 解法 ，可以计算出一

已知函数的各阶导数的精确表达式和在给定点处的精确值，是在计算机上实现泰勒级数解法

的基础。但是，在一阶常微分方程 ，0， 中含有 自变量 和因变量

． 且其中的函数关

糸 0)未知a一般地说 ’它的各阶导徽 也宙有蚵 ． 因此，必须推导出它 的各阶导数

表达式，下面是推导过程。

令 =F(#， )，由复合 函数微分法可得到

c，， 塞+墓 (2)

又 因为

毖 、，() (3)

故

(4)

为简洁起见，记，=，0， )，上式成为。

嚣，+ (5)

g(， )的一阶偏导数为

a塑 堑世+ f

葶 0 8 ‘a 0 a 。 (6)

()+等，+器 (7)

同样，令 ，_ (， )，可得到