第章Matlab符号运算V.ppt

数学软件 Matlab 本讲主要内容 符号运算 Matlab 符号运算 符号运算举例 符号对象 符号对象的建立 符号对象的建立 符号表达式 基本符号运算 符号对象的基本运算 findsym findsym 举例 subs subs 举例 常见符号运算 因式分解 函数展开 合并同类项 函数简化 函数简化举例 分式通分 计算极限 计算导数 计算积分 代数方程求解 代数方程求解 微分方程求解 微分方程求解 微分方程求解 微分方程求解 其它运算 上机作业 网络资源： The MathWorks - Matlab主页 水木社区上的MathTools版，有关Matlab的问题可以在上面提问 Matlab 中文论坛 例3：求微分方程组 在初值条件 下的特解，并画出解函数的图形。 [x,y]=dsolve(Dx+5*x+y=exp(t),Dy-x-3*y=0, ... x(0)=1, y(0)=0, t) 注：解微分方程组时，如果所给的输出个数与方程个数相同，则方程组的解按词典顺序输出；如果只给一个输出，则输出的是一个包含解的结构（structure）类型的数据。 反函数 finverse(f,v)：求 f 关于指定变量 v 的反函数 finverse(f )：求 f 关于默认变量的反函数 syms x t; f=x^2+2*t; g1=finverse(f,x) g2=finverse(f,t) 例：计算函数 的反函数 a1=1e10; b1=1e-10; c1=(a1+b1-a1)/b1; a2=sym(a1); b2=sym(b1); c2=(a2+b2-a2)/b2; 1、下面语句计算出来的 c1，c2 相等吗，为什么？上机验证。 2、设 A 是一个符号矩阵 (定义如下)，试指出 findsym(A,1) 的输出结果，并由此能得出什么结论？ syms a b t u v x y; A=[a+b*x,sin(t)+u; x*exp(-t),log(y)+v] 3、下面两个命令的结果是否相同，若不同，哪个是正确的？为什么？ factor(sym(12345678901234567890)) factor(sym(12345678901234567890)) * * —— Matlab 符号运算 Matlab 符号运算介绍 符号对象与基本符号运算 findsym 和 subs 常见符号运算 计算以推理方式进行，不受计算误差累积所带来的困扰 符号计算指令的调用比较简单，与教科书上的公式相近 符号计算可以给出完全正确的封闭解，或任意精度的数值解（封闭解不存在时 ） 符号计算所需的运行时间相对较长 符号运算的特点 Matlab 符号运算是通过符号数学工具箱（Symbolic Math Toolbox）来实现的。 Matlab 的符号数学工具箱可以完成几乎所有得符号运算功能，如：符号表达式的运算，符号矩阵的运算，符号微积分，符号作图，符号代数方程求解，符号微分方程求解等。 此外，该工具箱还支持可变精度运算，即支持以指定的精度返回结果。 Matlab 符号运算 求一元二次方程 ax2 + bx + c = 0 的根 solve(a*x^2+b*x+c=0) 求的根 f (x) = (cos x)2 的一次导数 x=sym(x); diff(cos(x)^2) 计算 f (x) = x2 在区间 [a, b] 上的定积分 syms a b x; int(x^2,a,b) 符号对象：在进行符号运算时，必须先定义基本的符号对象，可以是符号常量、符号变量、符号表达式等 符号对象是一种数据结构 符号表达式：含有符号对象的表达式称 符号矩阵/数组：元素为符号表达式的矩阵/数组 sym 函数用来建立单个符号变量，一般调用格式为： 符号对象的定义/声明：sym、syms 例： a=sym(a) 符号变量 = sym(x) 参数 x 可以是一个数或数值矩阵，也可以是字符串 a 是符号变量 b 是符号常量 b=sym(1/3) C 是符号矩阵 C=sym([1 ab; c d]) syms 符号变量1 符号变量2 ... 符号变量n 例： syms a b c a=sym(a); b=sym(b); c=sym(c); 符号对象的定义/声明：sym、syms syms 命令用来建立多个符号变量，一般调用格式为： 例： 建立符号表达式通常有以下 2 种方法：