用MATLAB符号工具箱编程求常微分方程的通解.pdfVIP

洛阳师范学院学报 2005年第 2 期 ·81 · 用 MA TLAB 符号工具箱编程求常微分方程的通解 唐洪浪 , 桂现才 (湛江师范学院数学系 , 广东湛江 524048) 摘 　要 : 用 MA TLAB 的函数 d so lve解常微分方程时 , 只能得出解析解 , 对于解的结果只能用 隐式来表示的情况 , d so lve便无能为力. 因此 , 用 d so lve对某些齐次微分方程和全微分方程求 解时 , MA TLAB 往往得不到有用的结果. 本文利用 MA TLAB 符号工具箱 , 通过编程可求出齐 次方程的解 , 通过自动寻找积分因子 , 从而求出全微分方程的通解 , 另外指出文献 [ 3 ]中的 一个错误. 关键词 : d so lve 函数 ; 齐次方程 ; 全微分方程 ; 积分因子 ; MA TLAB ; 符号工具箱 中图分类号 : O 175. 1 文献标识码 : A 文章编号 : 1009 - 4970 (2005) 02 - 0081 - 04 　　MA TLAB 符号工具箱功能强大 , 解符号微分 LAB 的符号积分 in t函数来实现. ( ) 方程的 d so lve 函数对通常的一阶线性微分方程 , 判断 1 式是否为齐次方程以及求齐次方程 可降阶的高阶微分方程 , 二阶常系数齐次和非齐 的程序 M 文件 d_ so lve1. m 如下 次方程能给出通解和特解 , 还能解微分方程组. function C = d_ so lve1 ( P, Q ) 因此 , d so lve 函数是解微分方程的有力工具 , 一 R = ( - 1) sym div ( P, Q ) ; 　% ( 1) 式化为 ( ) 般而言 , 如果微分方程的解 , 数学上能用显式 2 式 ( ) ( ) f = sub s R , ′y ′, ′u x ′; 　%作变换 y = ux y = f x 来表示 , d so lve 函数都能给 出通解 , 但 ( ) 是 , 在教学与实践中 , 笔者发现 d so lve 对于某些 F = simp le f ; 　%化简 ( ) 较简单的齐次方程 、全微分方程与可转化为全微