课堂授课专题：MATLAB语言基础.pptVIP

Page * MATLAB的符号计算功能 符号变量与符号表达式 符号变量 定义：syms a b c sym(‘a’) sym(‘b’) 符号表达式 定义： syms x y z; f1=x*y/z; f1=sym(‘x*y/z’) Page * MATLAB的符号计算功能 符号矩阵的创建：用Matlab函数sym创建矩阵 命令格式：A=sym([ ]) A = sym([a , 2*b ; 3*a , 0]) A = [ a, 2*b] [3*a, 0] 符号矩阵的修改和转换 符号矩阵的修改：用A1=subs(A, new, old)来修改 符号矩阵与数值矩阵的转换 将数值矩阵转化为符号矩阵，函数调用格式：sym(A) 将符号矩阵转化为数值矩阵，函数调用格式：double(A) A=[1/3,2.5;1/0.7,2/5] A = 0.3333 2.5000 1.4286 0.4000 例：A2=subs(A1, c, b) A2 =[ a, 2*c] [3*a, 4*c] sym(A) ans = [ 1/3, 5/2] [10/7, 2/5] double(A) ans = 0.3333 2.5000 1.4286 0.4000 Page * MATLAB的符号计算功能 符号矩阵的运算 7.0以后版本对符号表达式直接用+ - * / ^进行相应运算 符号运算函数 inv 逆矩阵 transpose 矩阵的转置 simple 符号矩阵简化 simplify 符号简化运算 f=sym([a,b;c,d]) f = [ a, b] [ c, d] g=sym([2*a,c;b,3*d]) g = [ 2*a, c] [ b, 3*d] f*g ans = [ 2*a^2+b^2, a*c+3*b*d] [ 2*a*c+b*d, c^2+3*d^2] f.*g ans = [ 2*a^2, c*b] [ c*b, 3*d^2] f=sym(cos(x)+sin(x)); g=sym(cos(x)-sin(x)); h=f*g h = (cos(x)+sin(x))*(cos(x)-sin(x)) simplify(h) ans = 2*cos(x)^2-1 simple(h) ans = cos(2*x) Page * MATLAB的符号计算功能 符号微积分 diff(f) 对缺省变量x求微分 diff(f,v) 对指定变量v求微分 diff(f,v,n) 对指定变量v求n阶微分 int(f) 对f表达式的缺省变量x求积分 int(f,v) 对f表达式的v变量求积分 int(f,v,a,b) 对f表达式的v变量在(a,b) 区间求定积分 F=int(int(x*exp(-x*y),x),y) F = 1/y*exp(-x*y) F=diff(sym(x^3+3*x^2*y+5*x*y^2), x,2) F = 6*x+6*y Page * MATLAB的符号计算功能 符号代数方程求解 solve 解一般的线性方程、非线性方程、代数方程和代数方程组 当方程组不存在符号解且无其他自由参数，则给出数值解 例1：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 例2：符号方程tan(2*x)=sin(x)求解 f1=solve(tan(2*x)=sin(x)) f1 = pi 0 atan(1/2*(-2*3^(1/2))^(1/2),1/2+1/2*3^(1/2)) atan(-1/2*(-2*3^(1