东北大学 MATLAB符号运算2014.pptVIP

syms x y a s=log(2*x/y); simplify(s) ans = log(2)+log(x/y) s=(-a^2+1)/(1-a); simplify(s) ans = a+1 s=sym([2*cos(x)^2-sin(x)^2,sqrt(16)]); simplify(s) ans = [ 3*cos(x)^2-1, 4] 符号表达式的化简 s1=simple(s) 从各种方法中自动选择最简格式 ，并显示化简过程 [s1,how]=simple(s) 化简并返回实际采用的化简方法 其中，s为原始表达式，s1为化简后表达式，how为采用的化简方法。 simple( ) 用几种不同的算术简化规则对符号表达式进行简化，使其用最少的字符来表示。 虽然并非表达式中的字符越少，表达式就越简单，但采用这个标准往往能够得到满意的结果，尤其是对于包含三角函数的表达式。 例：sym x simple(cos(x)^2+sin(x)^2) 从结果看出，simple比较这些不同函数的结果，最终把最少字符作为标准。 符号表达式的化简 syms x s=cos(x)+sqrt(-sin(x)^2) s = cos(x)+(-sin(x)^2)^(1/2) simple(s) simplify: cos(x)+(-sin(x)^2)^(1/2) radsimp: cos(x)+i*sin(x) combine(trig): cos(x)+1/2*(-2+2*cos(2*x))^(1/2) factor: cos(x)+(-sin(x)^2)^(1/2) expand: cos(x)+(-sin(x)^2)^(1/2) combine: cos(x)+1/2*(-2+2*cos(2*x))^(1/2) convert(exp): 1/2*exp(i*x)+1/2/exp(i*x)+1/2*((exp(i*x)-1/exp(i*x))^2)^(1/2) convert(sincos): cos(x)+(-sin(x)^2)^(1/2) convert(tan): (1-tan(1/2*x)^2)/(1+tan(1/2*x)^2)+(-4*tan(1/2*x)^2/(1+tan(1/2*x)^2)^2)^(1/2) collect(x): cos(x)+(-sin(x)^2)^(1/2) mwcos2sin: cos(x)+(-sin(x)^2)^(1/2) ans = cos(x)+i*sin(x) 例：化简多项式 syms s; P=(s+3)^2*(s^2+3*s+2)*(s^3+12*s^2+48*s+64) ； P = (s+3)^2*(s^2+3*s+2)*(s^3+12*s^2+48*s+64) simple(P) % 一系列化简尝试，得出计算机认为的最简形式 ans = (s+3)^2*(s+2)*(s+1)*(s+4)^3 [a,m]=simple(P) % 返回化简方法为因式分解方法，用 factor( ) 函数将得同样结果 a = (s+3)^2*(s+2)*(s+1)*(s+4)^3 m = factor expand(P) ans = s^7+21*s^6+185*s^5+883*s^4+2454*s^3+3944*s^2+3360*s+1152 变量替换 其中，f为原表达式，用x*替换x得出新的。 例：求其 Taylor 幂级数展开 syms a b c d t; % 假设这些变量均为符号变量 f=cos(a*t+b)+sin(c*t)*sin(d*t); % 定义给定函数 f(t) f1=subs(f,{a,b,c,d,t},{0.5*pi,pi,0.25*pi,0.125*pi,4}) f1 = -1.0000 符号型数值可采用变精度函数求值 vpa(A)：显示表达式A在当前精度下的值 vap(A,n)：显示符号表达式A在精度n下的值 vpa(pi) ans = 3.1415926535897932384626433832795 vpa(pi,60) ans = 3.14159265358979323846264338327950288419716939937510582097494 符号型变量数据类型 串演算函数 eval函数 利用字符串建立运算式后，再用eval命令执行它，可以使程序设计更加灵活。 功能是模拟字符串运算式，表达式一