第八讲、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 Dy a*y,y 0 b ans b*exp a*t dsolve Df f+sin t ans -1/2*cos t -1/2*sin t +exp t *C1 dsolve Df f+sin