Matlab第四讲：符号计算概述.ppt

数学软件 Matlab Matlab 符号运算 本讲主要内容 符号运算 Matlab 符号运算 符号运算举例 符号对象 符号对象的建立 符号对象的建立 符号表达式 基本符号运算 符号对象的基本运算 findsym findsym 举例 subs subs 举例 符号矩阵 六类常见符号运算 因式分解 函数展开 合并同类项 函数简化 函数简化 函数简化举例 分式通分 计算极限 计算导数 计算积分 符号求和 其它运算 符号求解 符号求解 dsolve 求解析解 dsolve 的使用 dsolve 举例 dsolve 举例 dsolve 举例 上机作业 solve 也可以用来解方程组 solve( f1 , f2 , ... , fN , v1 , v2 , ... , vN) 求解由 f1 , f2 , ... , fN 确定的方程组关于 v1 , v2 , ... , vN 的解 例：解方程组 [x,y,z]=solve(x+2*y-z=27,x+z=3, ... x^2+3*y^2=28,x,y,z) 输出变量的顺序要书写正确！ solve 在得不到解析解时，会给出数值解 dsolve 的使用 y=dsolve(eq1,eq2, ... ,cond1,cond2, ... ,v) 其中 y 为输出， eq1、eq2、...为微分方程，cond1、cond2、...为初值条件，v 为自变量。 例 1：求微分方程 的通解，并验证。 y=dsolve(Dy+2*x*y=x*exp(-x^2),x) syms x; diff(y)+2*x*y - x*exp(-x^2) 几点说明 如果省略初值条件，则表示求通解； 如果省略自变量，则默认自变量为 t dsolve(Dy=2*x,x); ％ dy/dx = 2x dsolve(Dy=2*x); ％ dy/dt = 2x 若找不到解析解，则返回其积分形式。 微分方程中用 D 表示对 自变量 的导数，如： Dy y； D2y y； D3y y 例 2：求微分方程 在初值条件 下的特解，并画出解函数的图形。 y=dsolve(x*Dy+y-exp(x)=0,y(1)=2*exp(1),x) ezplot(y); 例3：求微分方程组 在初值条件 下的特解，并画出解函数的图形。 [x,y]=dsolve(Dx+5*x+y=exp(t),Dy-x-3*y=0, ... x(0)=1, y(0)=0, t) ezplot(x,y,[0,1.3]); 注：解微分方程组时，如果所给的输出个数与方程个数相同，则方程组的解按词典顺序输出；如果只给一个输出，则输出的是一个包含解的结构（structure）类型的数据。 例： [x,y]=dsolve(Dx+5*x=0,Dy-3*y=0, ... x(0)=1, y(0)=1,t) r = dsolve(Dx+5*x=0,Dy-3*y=0, ... x(0)=1, y(0)=1,t) 这里返回的 r 是一个 结构类型 的数据 r.x %查看解函数 x(t) r.y %查看解函数 y(t) 只有很少一部分微分方程（组）能求出解析解。大部分微分方程（组）只能利用数值方法求数值解。 dsolve的输出个数只能为一个 或 与方程个数相等 * * Matlab 符号运算介绍 符号对象与基本符号运算 findsym 和 subs 几类常见的符号计算 计算以推理方式进行，不受计算误差累积所带来的困扰 符号计算指令的调用比较简单，与教科书上的公式相近 符号计算可以给出完全正确的封闭解，或任意精度的数值解（封闭解不存在时 ） 符号计算所需的运行时间相对较长 符号运算的特点 Matlab 符号运算是通过符号数学工具箱（Symbolic Math Toolbox）来实现的。该工具箱是建立在功能强大的 Maple 软件的基础上的。 Matlab 的符号数学工具箱可以完成几乎所有得符号运算功能，如：符号表达式的运算，符号矩阵的运算，符号微积分，符号作图，符号代数方程求解，符号微分方程求解等。此外，该工具箱还支持可变精度运算，即支持以指定的精度返回结果。 Matlab 符号运算介绍 求一元二次方程 ax2 + bx + c = 0 的根