MATLAB8.X程序的设计及典型应用第五章 符号计算.pptVIP

MATLAB实现符号多项式展开的指令为expand()，调用格式为： ?expand（S） ?对符号表达式S进行多项式、三角函数、指数函数和对数函数等展开。 【例5-15】展开表达式 ，其中n为任意给定的整数.编写文件名为exm5_15的脚本文件： clear,clc,syms x y n=input(Please input n=? ); g=(x-y)^n;g=expand(g) 在指令窗中运行exm5_15.m后，以n=4为例，结果为： Please input n=? 4 g = x^4 - 4*x^3*y + 6*x^2*y^2 - 4*x*y^3 + y^4 4符号表达式的展开 【例5-16】展开矩阵 中的各表达式. 编写文件名为exm5_16的脚本文件： syms alpha beta f=[sin(alpha+beta),sin(alpha-beta);sin(2*alpha),tan(2*beta)] clear,clc expand(f) 在指令窗中运行exm5_16.m，结果为： f =[ sin(alpha + beta), sin(alpha - beta)] ans = [cos(alpha)*sin(beta)+cos(beta)*sin(alpha),cos(beta)*sin(alpha)- cos(alpha)*sin(beta)] [ 2*cos(alpha)*sin(alpha), (2*tan(beta))/(tan(beta)^2-1)] 该题也证明了函数运算中的正弦函数、正切函数的倍角公式以及正弦函数和余弦函数的两角和、差公式。 符号表达式可以用许多等价的形式来表示。为了让用户得到最简洁的表达形式，MATLAB提供了两个指令进行符号化简：simplify()和simple()，其调用格式相同： ?simplify(S)或simple(S) ?对符号表达式S化简，输出最简化形式。二者区别在于对于较为复杂的表达式化简，不带输出宗量的simple() 往往给出符号表达式的化简过程，而simplify()只给出化简后的最简形式。 5 符号表达式的化简 【例5-17】化简指令simplify()和simple()使用实例：化简表达式 。 syms x; f=2*x/(x^2-x+1)+3/(x+1); simple(f) simplify: (5*x^2 - x + 3)/(x^3 + 1) radsimp: 3/(x + 1) + (2*x)/(x^2 - x + 1) …… collect(x): (5*x^2 - x + 3)/(x^3 + 1) ans = (5*x^2 - x + 3)/(x^3 + 1) simplify(f) ans = (5*x^2 - x + 3)/(x^3 + 1) 【例5-18】化简 编写文件名为exm5_18的脚本文件： clear,clc,syms x f=(1/x^3+6/x^2+12/x+8)^(1/3); f1=simplify(f) f2=simplify(f1) g1=simple(f) g2=simple(g1) 在指令窗中运行exm5_18.m后，结果为： f1 = ((2*x+1)^3/x^3)^(1/3) f2 = ((2*x+1)^3/x^3)^(1/3) g1 = (2*x+1)/x g2 = 2+1/x simplify()利用各种类型代数恒等式，包括求和、积分和分数幂、三角、指数和常用对数、Bessel函数、超几何函数和 函数，来简化表达式；simple()使用几种不同的简化工具，然后选择在结果表达式中含有最少字符的那种形式作为输出。建议用户在化简过程中多用几次化简指令。 【说明】 【例5-19】已知矩阵A= ,计算A的平方，以及A中各元素的平方，并化简。 编写文件名为exm5_19的脚本文件： clear,syms x ;A=[cos(x) sin(x);sin(x) cos(x)]; B=simple(A^2) %计算矩阵A的2次幂 C=simple(A.^2) %计算矩阵A中各元素的2次幂 在指令窗中运行exm5_19.m后，结果为： B = [ 1, sin(2*x)] [ sin(2*x), 1] C = [ cos(x)^2, sin(x)^2] [ sin(x)^2, cos(x)^2] A= 5.3.4符号表达式的置换 MATLAB符号数学工具箱中提供了用