八讲 MATLAB符号运算.pptVIP

第八讲 MATLAB符号运算 * 第八讲 MATLAB符号运算 8.1 设置符号变量 8.2 微分积分运算 8.3 方程求解 8.4 线性代数运算 8.5 其它命令 8.1 设置符号变量 1.设置符号变量（sym, syms） S = sym(A) x = sym(x) x = sym(x,real) x = sym(x,unreal) syms arg1 arg2 ... syms arg1 arg2 ... real syms arg1 arg2 ... unreal 8.2 微分积分运算 2.1微分 diff(S,v) （s对变量v求偏导） diff(S,n) （对s求n次导） diff(S,v,n)（s对变量 v 求n 次偏导） 8.2 微分积分运算 2.2积分 R = int(S) R = int(S,v) R = int(S,a,b) R = int(S,v,a,b) 例如： int(-2*x/(1+x^2)^2) ans= 1/(1+x^2) 8.2 微分积分运算 int(exp(-x^2)) ans = 1/2*pi^(1/2)*erf(x) int(exp(-x^2),0,inf) ans = 1/2*pi^(1/2) 8.2 微分积分运算 int(x/(1+z^2),z) ans= atan(z)*x int(x*log(1+x),0,1) ans= 1/4 int(2*x, sin(t), 1) ans= 1-sin(t)^2 8.2 微分积分运算 int([exp(t),exp(alpha*t)]) ans= [exp(t), 1/alpha*exp(alpha*t)] 2.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) 8.2 微分积分运算 例如： syms x a t h; limit(sin(x)/x) = 1 limit(1/x,x,0,right) = inf limit(1/x,x,0,left) = -inf limit((sin(x+h)-sin(x))/h,h,0) = cos(x) 8.2 微分积分运算 v = [(1 + a/x)^x, exp(-x)]; limit(v,x,inf,left) = [exp(a), 0] 2.4.级数和 symsum(s,t,a,b) s中t从a到b求和。 symsum(x^2,x,1,n) ans = 1/3*(n+1)^3-1/2*(n+1)^2+1/6*n+1/6 8.3 方程求解 3.解方程 3.1代数方程（solve） g = solve(eq) g = solve(eq,var) g = solve(eq1,eq2,...,eqn) g = solve(eq1,eq2,...,eqn,var1,var2,...,varn) 8.3 方程求解 如： 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 solve(x + y = 1,x - 11*y = 5) y = -1/3, x = 4/3 8.3 方程求解 A = solve(a*u^2 + v^2, u - v = 1, a^2 - 5*a + 6) A = a: [1x4 sym] u: [1x4 sym] v: [1x4 sym] 这里A.a = [ 2, 2, 3, 3] 8.3 方程求解 A.u = [ 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 = [ -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)] 8.3 方程求解 3.2微分方程（dsolve） r = dsolve(eq1,eq2,..., cond1,cond2,..., v) r = dsolve(eq1,eq2,...,cond1,c