MATLAB学习课件.ppt

第五章 Matlab 符号运算 第一节 创建符号变量 一、sym函数定义符号变量 sym函数用于定义单个符号变量,使用函数sym的调用格式为： sym(x) sym(x,real) sym(x,unreal) 其中,参数real定义的变量为实型符号变量,unreal定义的是非实型符号变量。 【例5-1】 使用函数sym定义符号变量。 a=sym(a) % 定义符号变量a a = a sym(b,real) % 定义符号变量b,实型符号变量 ans = b c=sym(byebye) c = Byebye 有了符号运算后,也可将数值变量向符号变量转换,见下例： 第二节 符号表达式的建立、化简和替换 在Matlab中,符号表达式可以通过基本赋值语句进行建立, 方法可以采用单引号或sym/syms函数定义。 【例5-5】 使用单引号建立符号表达式。 y=a*x^2+b=0 % 定义符号代数方程 y = a*x^2+b=0 y=D2y-c*Dy+d=0 % 定义符号微分方程 y = D2y-c*Dy+d=0 在书写符号表达式时,需要注意以下几点： 数学符号??? 的书写形式为?pi?；虚数单位用?i?或?j?表示； 无穷大用?INF?或?inf?表示； 符号相乘必须用?*?连接； 指数运算以e为底的书写形式为exp(),在Matlab中,只有求以e为底的自然对数； 表达式需写在同一行； 与数学表达式不同,Matlab的表达式中只能用小括号。在Matlab函数嵌套调用时, 多重小括号嵌套使用,要避免出错。 【例5-7】 展开符号表达式和。 首先在命令窗口创建符号变量。 syms x y; f1=(x+1)^7; expand(f1) ans = x^7+7*x^6+21*x^5+35*x^4+35*x^3+21*x^2+7*x+1 f2=cos(x+y); f=expand(f1) f = cos(x)*cos(y)-sin(x)*sin(y) 1. 因式分解（factor） 函数factor的功能是将符号表达式进行因式分解, 其调用格式为： factor(s) 其中,s可以是正整数、符号整数、符号表达式或 符号矩阵。当s为正整数时,因式分解的结果返回 的是s的质数分解式。在这里,我们需要注意的是, 在整数数组中,有一位元素位数超过16位时,必须 用sym函数创建该元素；当s为符号表达式时,结果 返回乘积形式。 【例5-9】 因式分解整数352、整数12345678901234567890、符号表达式以及符号矩阵。 f1=factor(sym(352)) f1 = (2)^5*(11) f2=factor(sym(1234567890123456789)) f2 = (2)^8*(3)*(11)*(146137297599841) syms a b x y f3=factor(x^5-y^5) f3 = (x-y)*(x^4+x^3*y+x^2*y^2+x*y^3+y^4) f4=factor([(x-b)^2 x^2-b*x;a^2-b^2 a*x-a]) f4 = [ (x-b)^2, x*(x-b)] [ (a-b)*(a+b), a*(x-1)] 可见,factor函数不仅在多项式中有着强大的分解能力,而且在正整 数质数分解中也有着广泛的应用。 从结果可以看出,采用不同的合并条件时,同一个 符号表达式可能会得到不同的结果。因此,在实 际应用中,应根据需要选择适当的合并条件。 3. 化简（simplify和simple） 符号表达式化简的函数simplify和simple,其中函 数simplify的调用格式为： R=simplify(s) 在Matlab中,simplify函数的具体功能是根据一定 规则对符号表达式进行化简。S代表符号表达式 或符号矩阵,R是化简后的结果。 【例5-12】 试对表达式和进行化简。 syms t x real R1=simplify(csc(t)^2-cot(t)^2) R1 = 1 R2=simplify((x^5-1)/(x-1)) R2 = x^4+x^3+x^2+x+1 函数simplify是一个强大的化简工具,它可以完成对指数、 对数、三角函数等各种数学表达式的化简。而函数simple 将化简结果及所使用的方法均列出来,这些方法中就包括 simplify。Simple的调用格式为