符号运算功能.docVIP

第5章 符号运算功能 在MATLAB中实现符号计算功能主要有以下三个主要途径 各种功能函数 Maple.m和mpa.m两个函数与Maple接口 符号函数计算器 5.1 符号表达式的生成 符号表达式包括函数和符号方程，两者的区别前者不包括等号，后者必须带等号。 创建符号函数 f=log(x) 创建符号方程 equation=a*x^2+b*x+c 创建符号微分方程 diffeq=Dy-y=x sym命令创建 f=sym(sin(x)) f=sym(sin(x)^2=0) syms命令 syms x %定义符号变量 f=sin(x)+cos(x) 5.2 符号和数值之间的转换 digits函数 digits(D) 函数设置为有效数字个数为D的近似精度 vpa函数 格式1：r=vpa(S) 符号表达式S在digits函数设置下的精度的数值解。 格式2：r=vpa(S,D) 符号表达式S在digits(D)函数设置下的精度的数值解。 S=solve(3*x^2-exp(x)) 解为： [ -2*lambertw(-1/6*3^(1/2))] [ -2*lambertw(-1,-1/6*3^(1/2))] [ -2*lambertw(1/6*3^(1/2))] 执行：vpa(S) 结果为：（数值解） [ .91000757248870906065733829575944] [ 3.7330790286328142006199540298434] [ -.45896226753694851459857243243408] 执行digits(4) vpa(S) 结果为：（数值解，精度为4位） [ .9100] [ 3.734] [ -.4590] 5.3 符号函数的运算 一、复合函数 compose(f,g)：返回当f=f(x)和g=g(y)时的符合函数f(g(y))，其中x和 y均为符号变量。 compose(f,g,z)：返回复合函数以z为自变量。 compose(f,g,x,z)：返回复合函数f(g(z))，且使x为f的独立变量。 compose(f,g,x,y,z)：返回复合函数f(g(z))，且使x为f的独立变量，y为g的独立变量。 举例： syms x y z t u ％定义符号变量 f=1/(1+x^2) ％定义符号函数 g=sin(y) ％定义符号函数 h=x^t ％定义符号函数 p=exp(-y/u) ％定义符号函数 compose(f,g) ％用sin(y)代替x，结果为1/(1+sin(y)^2) compose(f,g,t) ％用sin(y)代替x，t代替y，结果为1/(1+sin(t)^2) compose(f,g,x,z) ％返回复合函数f(g(z))，结果为sin(z)^t compose(h,g,t,z) ％结果为x^sin(z) compose(h,p,x,y,z) ％结果为exp(-z/u)^t compose(h,p,t,u,z) ％结果为x^exp(-y/z) 二、反函数 finverse(f) finverse(f,v) 举例： syms x y f=x^2+y finverse(f,y) 结果为：-x^2+y 5.4 符号矩阵的创立 1．使用sym函数直接创建 2．用创建子阵的方法创建符号矩阵 3．将数值矩阵转化为符号矩阵 5.5 符号矩阵的运算 一、基本运算 符号矩阵的四则运算（加、减、乘、除） 符号矩阵的行列式运算 符号矩阵的逆 符号矩阵的秩 符号矩阵的指数运算 二、矩阵分解 特征值函数 [x,y］=eig(b) 奇异值分解函数 svd(b) 约当标准型 三角抽取函数 diag： 对角线 tril： 上三角 triu： 下三角 三、符号矩阵的简化 因式分解 fact(S) syms x f=x^9-1 factor(f) ％结果为：(x-1)*(x^2+x+1)*(x^6+x^3+1) 符号矩阵的展开 expand(S) 举例： expand((x+1)^3) 得到： x^3+3*x^2+3*x+1 同类式合并 collect(S，v) 将符号矩阵S中的各元素的V的同幂项系数合并。 collect(S) 将符号矩阵S中的各元素的默认变量进行同类项合并。 举例： collect(