课堂授课专题1：数学物理建模思想与数据可视化方法.pptVIP

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’)来修改 例：A2=subs(A, c, b) A2 =[ a, 2*c] [3*a, 0] Page * MATLAB的符号计算功能 符号矩阵与数值矩阵的转换 将数值矩阵转化为符号矩阵，函数调用格式：sym(A) 将符号矩阵转化为数值矩阵，函数调用格式：double(A) sym(A) ans = [ 1/3, 5/2] [10/7, 2/5] double(A) ans = 0.3333 2.5000 1.4286 0.4000 A=[1/3,2.5;1/0.7,2/5] A = 0.3333 2.5000 1.4286 0.4000 Page * MATLAB的符号计算功能 符号矩阵的运算 7.0以后版本对符号表达式直接用+ - * / ^进行相应运算 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] Page * MATLAB的符号计算功能 符号运算函数 inv 逆矩阵 transpose 矩阵的转置 simplify 符号简化运算 simple 符号矩阵简化 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 Page * MATLAB的符号计算功能 例2：符号方程tan(2*x)=sin(x)求解 f1=solve(tan(2*x)=sin(x)) f1 = acos(1/2 - 3^(1/2)/2) acos(3^(1/2)/2 + 1/2) 0 -acos(1/2 - 3^(1/2)/2) -acos(3^(1/2)/2 + 1/2) double(f1) ans = 1.9455 0 + 0.8314i 0 -1.9455