MATLAB语言及应用-第六章.pptVIP

第6章 MATLAB符号计算 6.1 符号计算基础 6.2 符号导数及其应用 6.3 符号积分 6.4 级数 6.5 代数方程的符号求解 6.6 常微分方程的符号求解 6.1 符号计算基础 6.2 符号导数及其应用 6.3 符号积分 6.4 级数 6.5 代数方程的符号求解 6.6 常微分方程的符号求解 6.1 符号计算基础 6.1.1 符号对象 1. 建立符号变量和符号常数 (1)sym函数 sym函数用来建立单个符号量，例如，a=sym(a)建立符号变量a，此后，用户可以在表达式中使用变量a进行各种运算。 例6.1考察符号变量和数值变量的差别。 在 MATLAB命令窗口，输入命令： a=sym(a);b=sym(b);c=sym(c);d=sym(d); %定义 4个符号变量 w=10;x=5;y=-8;z=11; %定义4个数值变量 A=[a,b;c,d] %建立符号矩阵A B=[w,x;y,z] %建立数值矩阵B det(A) %计算符号矩阵A的行列式 det(B) %计算数值矩阵B的行列式 (2)syms函数 syms函数的一般调用格式为： syms var1 var2 … varn 函数定义符号变量var1,var2,…,varn等。用这种格式定义符号变量时不要在变量名上加字符分界符()，变量间用空格而不要用逗号分隔。 2. 建立符号表达式 符号表达式可以采用单引号或sym/syms函数定义 例6.2用两种方法建立符号表达式。 在MATLAB窗口，输入命令： U=sym(3*x^2+5*y+2*x*y+6) %定义符号表达式U syms x y; %建立符号变量x、y V=3*x^2+5*y+2*x*y+6 %定义符号表达式V 2*U-V+6 %求符号表达式的值 例6.3计算3阶范得蒙矩阵行列式的值。设A是一个由符号变量a,b,c确定的范得蒙矩阵。 命令如下： syms a b c; U=[a,b,c]; A=[[1,1,1];U;U.^2] %建立范得蒙符号矩阵 det(A) %计算A的行列式值 例6.4建立x,y的一般二元函数。 在MATLAB命令窗口，输入命令： syms x y; f=sym(f(x,y)) 6.1.2 基本的符号运算 1. 符号表达式运算 (1)符号表达式的四则运算 例6.5符号表达式的四则运算示例。 在 MATLAB命令窗口，输入命令： syms x y z; f=2*x+x^2*x-5*x+x^3 %结果为最简形式 f=2*x/(5*x) %符号表达式的结果为最简形式 f=(x+y)*(x-y) %结果不是x^2-y^2，而是(x+y)*(x-y) (2)因式分解与展开 factor(S) 对S分解因式,S是符号表达式或符号矩阵。 expand(S) 对S进行展开,S是符号表达式或符号矩阵。 collect(S) 对S合并同类项,S是符号表达式或符号矩阵。 collect(S,v) 对S按变量v合并同类项，S是符号表达式或符号矩阵。 例6.6 对符号矩阵A的每个元素分解因式。 命令如下： syms a b x y; A=[2*a^2*b^3*x^2-4*a*b^4*x^3+10*a*b^6*x^4,3*x*y-5*x^2;4,a^3-b^3]; factor(A) %对A的每个元素分解因式 例6.7 计算表达式S的值。 命令如下： syms x y; s=(-7*x^2-8*y^2)*(-x^2+3*y^2); expand(s) %对s展开 collect(s,x) %对s按变量x合并同类项(无同类项) factor(ans) % 对ans分解因式 练习： 1 展开符号表达式 f1=(x+1)7和 f2=cos(x+y)。 expand(S) 2 展开符号表达式(x+3)(x+t)(y-1)和(x-2)2(x+1)-x2。 3 因式分解符号表达式 x5-y5。 factor(S) 4 对符号表达式 (x-ex)(x+y)进行同类项合并。collect(S) 5 按照不同方式合并表达式 (x2-aey)(