MATLAB符号工具箱使用简介.docVIP

MATLAB符号工具箱使用简介 一、符号对象的创建 例1 s=hello s = hello 例2： sym(‘x’) ans = x y=sym(‘2+cos(x)’) y = 2+cos(x) 例3: syms y u p=exp(-y/u) p = exp(-y/u) q=y^2+u^3+u*y q = y^2+u^3+u*y 4: syms a b c k t y real g=q+a*(2*p-1)^3-b*sin(4*y) g = y^2+u^3+u*y+a*(2*exp(-y/u)-1)^3-b*sin(4*y) findsym (g) % 找出所有的符号变量 ans = a, b, u, y findsym (g,2) % 找出距x最近的2个符号变量 ans = y, u 例5 考察符号对象和普通数据对象之间的差别: sqrt(2) ans = 1.4142 x=sqrt(sym(2)) x = 2^(1/2) double(x) ans = 1.4142 二、微积分基本运算 例1 求极限 syms x a f=(x-a)/(x^2-a^2); limit(f,x,2) % 求x-2时的极限 ans = 1/(2+a) limit(f,a,2) % 求a-2时的极限 ans = 1/(x+2) limit(f) % 默认 x-0 ans = 1/a limit(1/x) % 求x-0时极限不存在 ans = NaN limit(1/x,x,0,‘left’) % 左极限 ans = -inf limit(1/x,x,0,‘right % 右极限 ans = inf 例2 求导数 syms x a t ; f=sin(a*x)-cos(a*t) ; g1=diff(f) % 默认对 x求导 g1 = cos(a*x)*a g2=diff(f,a) % 对 a求导 g2 = cos(a*x)*x+sin(a*t)*t g3=diff(f,x,2) % 对 x求二阶导数 g4 = -sin(a*x)*a^2 例3 对符号矩阵求导, 是作用于矩阵的每个元素. syms x a t ; A=[sin(a*x),cos(t*x);a^2-x*t,t^2+t*x] A = [ sin(a*x), cos(t*x)] [ a^2-x*t, t^2+t*x] dy=diff(A) % 默认对 x求导 dy = [ cos(a*x)*a, -sin(t*x)*t] [ -t, t] dy=diff(A, t) % 对 t求导 dy = [ 0, -sin(t*x)*x] [ -x, 2*t+x] 例4 求不定积分 syms x a t ; f=sin(a*x); g=int(f) % 默认对 x求积分 g = -1/a*cos(a*x) g=int(f,a) % 对 a求积分 g = -1/x*cos(a*x) 例5 求定积分 syms x a t ; f=sin(a*x);g=int(f,0,pi)