仿真工具教程_02MATLAB符号计算精要.ppt

第2章 MATLAB符号计算 2.1 符号对象和符号表达式 2.2 符号数字及表达式的操作 2.3 符号微积分 2.4 符号微分方程的求解 2.5 符号变换和符号卷积 2.6 符号代数方程的求解 什么是符号计算？ 符号计算是指： 解算数学表达式、方程不是在离散化的数值点上进行，而是凭借一系列恒等式，数学定理，通过推理和演绎，获得解析结果。 符号计算的特点： 1）符号表达式计算前必须先定义符号变量。 2）符号计算是精确计算。 3）符号计算的计算速度较慢。 4）符号计算的运算符和函数与数值计算中的运算符和函数几乎完全相同 2.1 符号对象和符号表达式 符号对象的创建 （1）sym(‘变量’,参数) 功能：把变量定义为符号对象。 参数有三种选择： ’positive’ 表示为“正”符号变量， ’real’ 表示为“实”符号变量， ’unreal’ 表示为“非实”符号变量。 （2）syms函数 格式：syms arg1 arg2 …参数 功能：创建多个符号变量。 例如： syms r_l real syms a b c d positive 一、符号数字与数值数字 【例2.1-1】符号（类）数字与数值（类）数字之间的差异。 a=pi+sqrt(5) sa=sym(pi+sqrt(5)) Ca=class(a) Csa=class(sa) vpa(sa-a) 二、符号变量 1、符号表达式允许使用自由变量。确定自 由变量的原则： 1）小写字母i和j不能作为自由变量。 2）符号表达式中如果有多个字符变量，则按照以下顺序选择自由变量： 首先选择x作为自由变量；如果没有x，则选择在字母顺序中最接近x的字符变量；如果与x相同距离，则在x后面的优先。 【例2.1-2】 用符号计算研究方程 的解。 （1）不指定变量（默认） syms u v w z Eq=u*z^2+v*z+w; result_1=solve(Eq) findsym(Eq,1) 2）指定变量为Z result_2=solve(Eq,z) 2、确定自由变量的指令 findsym (11a版为symvar)的格式为： findsym(EXPR,n) 功能：确定EXPR中的自由变量。 其中EXPR可以是符号表达式或符号矩 阵；n为按顺序得出符号变量的个数， 当n省略时，则不按顺序给出EXPR中所 有的符号变量。 【例】findsym确定自由变量是对整个矩阵进行的。 syms a b t u v x y A=[a+b*x,sin(t)+u;x*exp(-t),log(y)+v] findsym(A,1) 2.1.1符号计算中的算符 与数字计算的符号相同 2.1.2符号计算中的函数指令 见表2.1-1，P45 注意与数值计算的函数和指令的异同 2.1.3 符号对象的识别 例2.1-5 数据对象及其识别指令的使用。 （1） clear a=1;b=2;c=3;d=4; Mn=[a,b;c,d] ; Mc=[a,b;c,d]; Ms=sym(Mc) ; 例：应用digits和vpa函数设置运算精度。 a=sym(2*sqrt(5)+pi) %创建符号对象 digits %显示默认的有效位数 vpa(a) %用默认的位数计算并显示 vpa(a,20) %按指定的精度计算并显示 digits(15) %改变默认的有效位数 vpa(a) %按digits指定的精度计算并显示 2.2.2符号表达式的基本操作 指令：simple(EXPR) 【例2.2-2】简化。 syms x f=(1/x^3+6/x^2+12/x+8)^(1/3); g1=simple(f) g2=simple(g1) 2.2.3表达式中的置换操作 1、通用置换指令 subs (S) 功能：用MATLAB工作空间中的变量替换S符号 表达式中的所有变量。 subs (S,NEW) 功能：用变量NEW替换符号表达式S中的自由 变量。 subs (S,OLD,NEW) 功能：用变量NEW替换符号表达式S中的变量 OLD。 【例2.2-4】演示subs的置换规则。 （1）syms a x f=a*sin(x)+5 （2）f1=subs(f