第8讲、MATLAB符号运算.docVIP

第八讲 MATLAB符号运算 8.1 设置符号变量 8.2 微分积分运算 8.3 方程求解 8.4 线性代数运算 8.5 其它命令 1、设置符号变量（sym, syms） S = sym(‘A’) x = sym(x) syms arg1 arg2 ... 2、微分积分运算 (1)微分 diff(S,v) （s对变量v求偏导） diff(S,n) （对s求n次导） diff(S,v,n)（s对变量 v 求n 次偏导） 例1： syms x a y=sin(a*x)/x y = sin(a*x)/x diff(y) ans = cos(a*x)*a/x-sin(a*x)/x^2 diff(y,2) ans = -sin(a*x)*a^2/x-2*cos(a*x)*a/x^2+2*sin(a*x)/x^3 diff(y,a) ans = cos(a*x) diff(y,a,2) ans = -sin(a*x)*x (2)积分 R = int(S) R = int(S,v) R = int(S,a,b) R = int(S,v,a,b) 例2： syms x t y=exp(t*x) y = exp(t*x) int(y) ans = 1/t*exp(t*x) int(y,t) ans = 1/x*exp(t*x) int(y,0,1) ans = (exp(t)-1)/t int(y,t,0,1) ans = (exp(x)-1)/x z=exp(-x^2) z = exp(-x^2) int(z) ans = 1/2*pi^(1/2)*erf(x) int(z,0,inf) ans = 1/2*pi^(1/2) int([x sin(x);cos(x) t*x^2]) ans = [ 1/2*x^2, -cos(x)] [ sin(x), 1/3*t*x^3] (3)极限 limit(F,x,a) limit(F,a) （x-a时F的极限） limit(F) （x-0时F的极限) limit(F,x,a,right) limit(F,x,a,left) 例3： syms x a h limit(sin(x)/x) ans = 1 limit(1/x,x,0,left) ans = -Inf limit(1/x,x,0,right) ans = Inf limit((sin(x+h)-sin(x))/h,h,0) ans = cos(x) (4)级数和 symsum(s,t,a,b) s中t从a到b求和。 例4： syms x n symsum(x^2,x,1,n) ans = 1/3*(n+1)^3-1/2*(n+1)^2+1/6*n+1/6 3、方程求解 (1)代数方程（solve） g = solve(eq) g = solve(eq,var) g = solve(eq1,eq2,...,eqn) g = solve(eq1,eq2,...,eqn,var1,var2,...,varn) 例5：solve(a*x^2 + b*x + c) ans= [ 1/2/a*(-b+(b^2-4*a*c)^(1/2)), 1/2/a*(-b-(b^2-4*a*c)^(1/2))] solve(a*x^2 + b*x + c,b) ans= -(a*x^2+c)/x S=solve(x+y=1,x-11*y=5) S = x: [1x1 sym] y: [1x1 sym] S.x ans = 4/3 S.y ans = -1/3 A=solve(a*u^2+v^2,u-v=1,a^2-5*a+6) A = a: [4x1 sym] u: [4x1 sym] v: [4x1 sym] A.a ans = 2 2 3 3 A.u ans = 1/3+1/3*i*2^(1/2) 1/3-1/3*i*2^(1/2) 1/4+1/4*i*3^(1/2) 1/4-1/4*i*3^(1/2) A.v ans = -2/3+1/3*i*2^(1/2) -2/3-1/3*i*2^(1/2) -3/4+1/4*i*3^(1/2) -3/4-1/4*i*3^(1/2) (2)微分方程（dsolve） r = dsolve(eq1,eq2,..., cond1,cond2,..., v) r = dsolve(eq1,eq2,...,cond1,cond2,...,v) 例6： dsolve(Dy=a*y) ans = C1*exp(a*t) dsolve