MATLAB讲稿-符号函数.pptVIP

补充章节 符号运算 MATLAB7.x的符号计算工具包为3.1版，它以Maple的内核为符号计算的引擎，依赖MAPLE已有的库函数。 MATLAB中包含两个工具包：基本Symbolic Math Toolbox和扩展Symbolic Math Toolbox。 基本Symbolic Math Toolbox收集了100多个MATLAB函数，提供了访问Maple内核的途径，它所采用的语言形式是MATLAB语言的自然扩充，它也允许访问Maple的线性代数软件包中的函数。 扩展Symbolic Math Toolbox增加了一些功能，包括访问所有的Maple非图形软件包、Maple编程特征和用户定义的过程。通过这两个工具包，用户可以书写的M文件，访问Maple函数和Maple工作空间。 为了获得Symbolic Math Toolbox的快速联机介绍，可在MATLAB命令窗口键入demos，MATLAB将显示MATLAB demos对话框，选取Symbolic Math即可。 Symbolic Math Toolbox定义了一个新的数据类型为符号对象（所有symbolic object）或sym。本身来讲，一个符号对象是一种数据结构，用来存贮代表符号的字符串。Symbolic Math Toolbox使用符号对象代表符号变量、表达式和矩阵。 §1 符号对象和符号表达式 一、创建符号变量 （1）生成符号变量 例：syms a b x X Y;k=sym(3); z=sym(‘c*sqrt(delta)+y*sin(theta)’); % 将字符串表达式变成符号变量EXPR=a*z*X+(b*x^2+k)*Y; %创建符号数学函数 说明：前两条指令定义a b x X Y是符号变量，定义后就可以参加后面的符号计算了。 例：用符号变量代表黄金分割比率： rho=sym((1+sqrt(5))/2); %下面可以对rho进行各种操作了。 f=rho^2-rho-1 simple(f) simplify(f) 例：求矩阵A=[a11,a12;a21,a22]的行列式值、逆和特征根 syms a11 a12 a13 a21 a22 a23 a31 a32 a33; A=[a11 a12 a13;a21 a22 a23; a31 a32 a33], DA=det(A),IA=inv(A),EA=eig(A) (2)实型和复型符号函数的创建 sym命令可以确定符号变量的数学属性，通过real选项实现。也就是说，下列指令包含了符号变量x和y，并有附加的数学特性为实型变量。 x=sym(x,real);y=sym(y,real); 或采用更有效的方式定义： syms x y real z=x+j*y; 特别地，这意味着表达式f=x^2+y^2严格非负，因此z是一个复变量，可以做如下操作。 conj(x),conj(z) 返回相应变量的共轭复数。 二、符号表达式中自由变量的确定： ———findsym函数syms a b x X Y;k=sym(3); z=sym(c*sqrt(delta)+y*sin(theta)); EXPR=a*z*X+(b*x^2+k)*Y;findsym(EXPR) % 除常数符号 k 外的所有独立符号变量都被列出findsym(EXPR,1) %在 EXPR 中确定一个自由符号变量 findsym(EXPR,2),findsym(EXPR,3) %在 EXPR 中确定 2 个和 3 个自由变量 §2 符号表达式和符号函数的操作 一、符号表达式的操作1、collect: 合并系数 EXPR=sym((x^2+x*exp(-t)+1)*(x+exp(-t))); expr1=collect(EXPR) % 默认合并 x 同幂项系数expr2=collect(EXPR,exp(-t)) % 合并 exp(-t) 同幂项系数 2、factor 指令的使用:分解因式和质数分解（1）除 x 外不含其它自由变量的情况syms a x;f1=x^4-5*x^3+5*x^2+5*x-6;factor(f1) （2）含其它自由变量的情况之一f2=x^2-a^2;factor(f2) （3）对正整数的质数分解factor(1025) 例：对多项式进行嵌套型分解——horner clear;syms a x;f1=x^4-5*x^3+5*x^2+5*x-6;horner(f1) 写出矩阵 各元素的分子、分母多项式，求矩阵各元素的分子、分母多项式syms s;A=[3/2,(s^2+3)/(2*s-1)