MATLAB的符号计算培训讲义.pptVIP

solve(s1,s2,…,sn,v1,v2,…,vn)：求解符号表达式s1,s2,…,sn组成的代数方程组，求解变量分别v1,v2,…,vn。 说明： 若不指明符号表达式 s1,s2,...,sn的值，系统默认为0。例如给出一个表达式x^2-3*x-8,则系统将按x^2-3*x-8=0进行运算. 例、解代数方程：a*x2-b*x-6=0. syms a b x solve(a*x^2-b*x-6) ans = [ 1/2/a*(b+(b^2+24*a)^(1/2))] [ 1/2/a*(b-(b^2+24*a)^(1/2))] 即该方程有两个根: x1=1/2/a*(b+(b^2+24*a)^(1/2))； x2=1/2/a*(b-(b^2+24*a)^(1/2)) 例. f = ax2+bx+c 求解 f=a*x^2+b*x+c; solve(f) —— 对缺省变量x求解 ans =[1/2/a*(-b+(b^2-4*a*c)^(1/2))] [1/2/a*(-b-(b^2-4*a*c)^(1/2))] solve(f , b ) —— 对指定变量b求解 ans = -(a*x^2+c)/x 计算机 格式 一般格式 例. 求符号方程tan(2*x)=sin(x)的解 f2=solve(tan(2*x)=sin(x)) f2 = %matlab4.2的解 [ 0] [acos(1/2+1/2*3^(1/2))] [acos(1/2 -1/2*3^(1/2))] f3=solve(tan(2*x)=sin(x)) f3= matlab6.5的解 [ 0] [ pi] [ atan(1/2*(-2*3^(1/2))^(1/2),1/2+1/2*3^(1/2))] [ atan(-1/2*(-2*3^(1/2))^(1/2),1/2+1/2*3^(1/2))] [ atan(1/2*2^(1/2)*3^(1/4)/(1/2-1/2*3^(1/2)))+pi] [ -atan(1/2*2^(1/2)*3^(1/4)/(1/2-1/2*3^(1/2)))-pi] numeric(f3) ans = 0 3.1416 0 + 0.8314i 0 - 0.8314i 1.9455 -1.9455 numeric(f2) ans = 0 0 + 0.8314i 1.9455 matlab4.2版与6.5版 的对比 使用solve函数求解一般的符号代数方程组 [x,y] = solve(x^2 + x*y + y = 3,x^2 - 4*x + 3 = 0) x = [ 1] [ 3] y = [ 1] [ -3/2] 例. 解方程组 x+y+z=1 x-y+z=2 2x-y-z=1 f=solve(x+y+z=1,x-y+z=2,2*x-y-z=1) f = x: [1x1 sym] y: [1x1 sym] z: [1x1 sym] [x,y,z]=solve(x+y+z=1,x-y+z=2,2*x-y-z=1) x = 2/3; y =-1/2; z =5/6 2、符号常微分方程求解 在MATLAB中，用大写字母D表示导数。例如，Dy表示 y，D2y表示y，Dy(0)=5表示y(0)=5。D3y+D2y+Dy-x+5=0 表示微分方程y‘’‘+y’‘+y’-x+5=0。符号常微分方程求解可以 通过函数dsolve来实现，其调用格式为： dsolve(e,c,v) 该函数求解常微分方程 e 在初值条件 c 下的特解。参数 v 描述 方程中的自变量，省略时按缺省原则处理，若没有给出初值 条件 c