精讲多练Mtlab3.ppt

精讲多练MATLAB 杨琦 西安交通大学计算机教学中心 2003-05-14 第3章 MATLAB的符号计算 3.1 符号变量和符号表达式(Symbolic expression) 3.2 微积分 3.3 方程求解 3.4 符号表示式的运算 3.5 数学式的化简 3.6 sym函数 3.7 求反函数和复合函数 3.8 表达式替换 3.9 任意精度计算 3.10 符号积分变换 3.11 Maple函数的使用 3.1 符号变量和符号表达式 (Symbolic expression) 1. 建立符号变量和符号常数 (1)sym函数 sym函数用来建立单个符号量，例如，a=sym(‘a’)建立符号变量a，此后，用户可以在表达式中使用变量a进行各种运算。 (2)syms函数 syms函数的一般调用格式为： syms var1 var2 … varn 函数定义符号变量var1,var2,…,varn等。用这种格式定义符号变量时不要在变量名上加字符分界符()，变量间用空格而不要用逗号分隔。 例3.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的行列式 例3.2 计算二次方程f=ax^2+bx+c的积分、微分。 syms a b c x f=sym(a*x^2+b*x+c) df=diff(f) nf=int(f) 表3-1默认自变量 默认原则：自变量为除了i和j之外并且在字母位置上最接近x的小写字母；如果式子中没有上述字母，则x会被视为默认的自变量。 findsym (S,n)——符号表达式中变量的确定 MATLAB中的符号可以表示符号变量和符号常数。findsym可以帮助用户查找一个符号表达式中的的符号变量。该函数的调用格式为： findsym(S,n) 函数返回符号表达式S中的n个符号变量，若没有指定n，则返回S中的全部符号变量。 在求函数的极限、导数和积分时，如果用户没有明确指定自变量，MATLAB将按缺省原则确定主变量并对其进行相应微积分运算。可用findsym(S,1)查找系统的缺省变量，事实上，MATLAB按离字符x最近原则确定缺省变量。 2. 建立符号表达式 例3.3用两种方法建立符号表达式。 在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 %求符号表达式的值 例3.4计算3阶范得蒙矩阵行列式的值。设A是一个由符号变量a,b,c确定的范得蒙矩阵。 命令如下： syms a b c; U=[a,b,c]; A=[[1,1,1];U;U.^2]%建立范得蒙符号矩阵 det(A) %计算A的行列式值 例3.5建立x,y的一般二元函数。 在MATLAB命令窗口，输入命令： syms x y; f=sym(f(x,y)); 3.5 数学式的化简 1. 符号表达式运算 (1)符号表达式的四则运算 例3.6符号表达式的四则运算示例。 在 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进