第8次课 符号计算.pptVIP

符号计算 【学习目标】 ● 掌握符号对象的定义方法以及符号表达式的运算法则。 ● 掌握微积分的符号计算方法。 ● 掌握级数求和的方法以及将函数展开为泰勒级数的方法。 ● 掌握代数方程和微分方程符号求解的方法。 建立符号对象 sym函数用来建立单个符号量 调用格式为 符号量名?=?sym(符号字符串) 注：符号字符串可以是常量、变量、函数表达式。 函数syms一次可以定义多个符号变量。 调用格式为 syms 符号变量名1 符号变量名2…符号变量名n 注： 变量间用空格而不要用逗号分隔。 例如，用syms函数定义4个符号变量a、b、c、d，命令如下： syms a b c d 建立符号表达式 （1）用sym函数建立符号表达式。例如： U=sym(3*x^2-5*y+2*x*y+6) U = 3*x^2-5*y+2*x*y+6 M=sym([a,b;c,d]) M = [ a, b] [ c, d] （2）使用已经定义的符号变量组成符号表达式。例如： syms x y; V=3*x^2-5*y+2*x*y+6 V = 3*x^2 + 2*y*x - 5*y + 6 1．符号表达式的四则运算 用+、?、*、/、^ 运算符实现，运算结果依是符号表达式。 f= sym(2*x^2+3*x-5) g= sym(x^2-x+7) f+g ans= 3*x^2 + 2*x + 2 2．符号表达式的提取分子和分母运算 如果符号表达式是一个有理分式或可以展开为有理分式，可利用numden函数来提取符号表达式中的分子或分母。 调用格式为 [n,d]=numden(s) 分子和分母存放在n与d中 如果符号表达式是一个符号数组（矩阵），将获得两个新数组（矩阵）n和d，其中n是分子数组（矩阵）， d是分母数组（矩阵）。例如： h=sym([3/2,(2*x+1)/3;a/x+a/y,3*x+4]) ; [n,d]=numden(h) n = [ 3, 2*x + 1] [a*x + a*y, 3*x + 4] d = [ 2, 3] [ x*y, 1] 3．符号表达式的因式分解与展开 ● factor(s)：对符号表达式s分解因式。 ● expand(s)：对符号表达式s进行展开。 ● collect(s,v)：对符号表达式s按变量v合并同类项。 syms x y; s1=x^3-y^3; factor(s1) %对s分解因式 ans = (x - y)*(x^2 + x*y + y^2) s2=(-7*x^2-8*y^2)*(-x^2+3*y^2); expand(s2) %对s展开 ans = 7*x^4 - 13*x^2*y^2 - 24*y^4 s3=(x+y)*(x^2+y^2+1) collect(s3,y) %按y合并同类项 ans = y^3 + x*y^2 + (x^2 + 1)*y + x*(x^2 + 1) 4．符号表达式的化简 ● simplify(s)：应用MuPAD简化规则对s进行化简。 ● simple(s)：调用MATLAB的其他函数对表达式进行综合化简，并显示化简过程。 s= sym((x^2+5*x+6)/(x+2)); simplify(s) ans= x + 3 s= sym(cos(3*acos(x))); simple(s) … ans = 4*x^3 - 3*x 5．符号表达式与数值表达式之间的转换 函数sym可以将数值表达式变换成它的符号表达式。 函数eval可以将符号表达式变换成数值表达式。 6．符号多项式与多项式系数向量之间的转换 函数sym2poly可以将符号多项式转换为多项式系数向量。 函数poly2sym可以将多项式系数向量转换为符号多项式。 符号矩阵 函数作用于符号矩阵时，是分别作用于矩阵的每一个元素。 求函数极限 求函数极限的函数是limit。 对于极限值为“没有定义”的极限，给出的结果为NaN； 极限值为无穷大时，给出的结果为inf。 ● limit(f,x,a)：求符号函数f(x)的极限值，即计算当变量x趋近于常数a时，f(x)函数的极限值。 ● limit(f,x,a,‘right’)或limit(f,x,a,‘left’)：求符号函数f的极限值。right表示变量x从右边趋近于a，left表示变量x从左边趋近于a。 极限1： f=sym((sin(x+h)-sin(x))/h); limit(f,sym(h),0) ans = cos(x) 极