实验符号运算.docVIP

实验3 符号运算 一、实验目的 符号(symbol)运算的基本功能. 二、实验内容与要求 1. 字符型变量、符号变量、符号表达式、符号方程的建立 用单引号设定字符串变量 a =u+4 %定义a为字符型变量 a = u+4 用命令sym(‘’)创建单个符号变量、符号表达式、符号方程. x= sym(m+n+i) %定义x为符号型变量 x= m+n+i y = sym（d*x^2 + x – 4） %定义y为符号表达式 y= d*x^2 + x – 4 e = sym( a*x^2+b*x+c=0) %定义e为符号方程 e= a*x^2+b*x+c=0 用命令syms创建多个符号变量、符号表达式. syms a b x y %定义a，b，x，y为符号变量，字母间必须用空格 s = a*x^4+b*cos(y)-x*y %定义s为符号表达式 s= a*x^4+b*cos(y)-x*y 注意：sym(‘’)中的单引号不要漏，syms后的符号变量之间不能用逗号，用syms不能建立符号方程. 2. 合并同类项 格式：collect(S) %按缺省变量x的次数合并系数. collect(S,v) %是对指定的变量v计算，操作同上. syms x y %定义x，y为符号变量 R1=collect((exp(x)+x)*(x+2)); 3. 复合函数计算 格式：compose(f,g,x,y)%返回复合函数f [ g (y)]，f = f (x)，g = g (y). syms x y f = 1/(1 + x^2*y); g = sin(y); C = compose(f,g,x,y) % 结果为1/(1+sin(y)^2*y) 4. 符号表达式的展开 格式：R = expand(S) %展开符号表达式S中每个因式的乘积. syms x y t E=expand((x-2)*(x-4)*(y-t)) 5. 符号因式分解 格式：factor(S) % S可以是正整数、符号表达式或符号整数. syms x y F1 = factor(x^4-y^4) %结果为(x-y)*(x+y)*(x^2+y^2) 问题：若F2=factor(sym(‘12345678901234567890’))，结果如何？ 6. 符号表达式的通分 numden(Numerator and denominator) 格式：[N, D] = numden(S) %将符号表达式S中的每一元素进行通分，其中N为分子的表达式，D为分母的表达式. syms x y [N, D] = numden(x/y + y/x) %结果为N=x^2+y^2, D= y*x 7. 符号表达式的化简 格式：R = simplify(S) %运用多种恒等式转换对符号表达式S进行综合化简. syms x a b c R = simplify(exp(c*log(sqrt(a+b)))) %结果(a+b)^(1/2*c) 8. 搜索符号表达式的最简形式 格式：r = simple(S) %运用包括simplify在内的各种指令找出符号表达式S的代数上的最简短形式，多次使用，可找到最少字母的简化式. 【例】 化简. syms x f=(1/x^3+6/x^2+12/x+8)^(1/3); f1=simple(f),f2=simple(f1) f1 = (2*x+1)/x f2 = 2+1/x 问题：分别用simple，simplify命令两次化简，试比较命令simple，simplify命令之间的区别和联系. 9. 将复杂的符号表达式显示成我们习惯的数学书写形式 格式：pretty(S) %用缺省的线型宽度79显示符号矩阵S中每一元素. y=sym(log(x)/sqrt(x)); dy = diff(y) pretty(dy) 计算结果为： dy = 1/x^(3/2)-1/2*log(x)/x^(3/2) 1 log(x) ---- - 1/2 ------ 3/2 3/2 x x 10. 函数的反函数 格式：g = finverse( f )%返回函数f的反函数，其中f为单值的一元数学函数， f = sym(1+3*x); V = finverse(f) %结果为-1/3+1/3*x 11. 符号表达式求和 格式：r = symsum(S,v,a,b) %对S中指定的符号变量v从a到b求和. syms n r = symsum(n^2