实验数学一： MATLAB的基本操作.ppt

ans = 4*x+11+2*x^2 f-g ans = 4*x-1-2*x^2 f*g ans = (4*x+5)*(2*x^2+6) f/g ans = (4*x+5)/(2*x^2+6) f^g ans = (4*x+5)^(2*x^2+6) （1）MATLAB提供了simple和 simplify命令，用于对表达式的化简。 例10 化简函数 解 输入命令 ： f=sym(cos(x)^2-sin(x)^2); f=simple(f) 5、几种初等数学符号表达式的运算 f = cos(2*x) 例11 化简 解 输入命令 ： syms x; f=(3*x^4+20*x^3+29*x^2-20*x-32)/(3*x^4-4*x^3-35*x^2+4*x+32); simplify(f) ans= (x+4)/(x-4) 这里注意: simplify是运用多种恒等式转换对符号表达式进行综合化简， simple是找出符号表达式的代数上的最简形式，多次使用可找到最少字母的简化式。 如化简 输入命令 ： syms x f=(1/x^3+6/x^2+12/x+8)^(1/3); f1=simple(f) f1 = (2*x+1)/x 若再输入命令 ： f2=simple(f1) f2 = 2+1/x (2) MATLAB提供了factor命令，用于对表达式进行分解。 例12 将 分解因式。 解 输入命令 ： syms x; f=x^12-1; factor(f) ans = (x-1)*(1+x^2+x)*(1+x)*(1-x+x^2)*(1+x^2)*(x^4-x^2+1) 说明：factor(f)中的f可以是正整数，但在MATLAB6.5 不超过2^32. 例如 factor(12345) ans = 3 5 823 所以有 12345=3 5 823 factor(1234567890) ans = 2 3 3 5 3607 3803 所以有 1234567890= 2 3 3 5 3607 3803 (3) MATLAB提供了collect命令，用于合并符号表达式 的同类项。 例13 合并 同类项。 解 输入命令 ： syms x; f= (x-1)^2*(x-3)*(x-5)*(x-7); collect(f) ans = -105+x^5-17*x^4+102*x^3-262*x^2+281*x (4) MATLAB提供了expand命令，用于对符号表达 式进行展开。 例14 将 和 进行展开。 解 输入命令 ： syms x y; v=[(x+1)^3 sin(x+y)]; expand(v) ans= [x^3+3*x^2+3*x+1,sin(x)*cos(y)+cos(x)*sin(y)] (5) MATLAB提供了numden命令，用于对符号表达 式的通分。 格式：[N,D]= numden(S) %将符号表达式S的每一 元素进行通分，其中N为通分后的分子表达式，D为通分后的分母表达式。 例15 将式子 通分。 解 输入命令 ： syms x y [N,D]=numden(x/y+y/x) N = x^2+y^2 D = x*y 所以 = (6) MATLAB提供了compose命令，用于对符号表 达式的复合运算。 例16 已知 , ,求f(g(y)). 解 输入命令 ： syms x y; f=1/x^3; g=tan(y); compose(f,g) ans = 1/tan(y)^3 例17已知 ，求f(g(y)). ans = (2+sin(y)^2)^3 解 输入命令 ： syms x y f=x^3; g=2+(sin(y))^2; compose(f,g) (7) MATLAB提供了finverse命令，用于求符号 表达式的反函数。 例18 求 的反函数。 ans = sin(1/x) 解 输入命令 ： f=sym(1/sin(x)); finverse(f) 例19 求 的反函