Matlab课件第三讲符号计算.ppt

第三讲 MATLAB符号计算 3.1 符号计算基础 符号类：sym类，符号类的实例即符号对象 符号对象：符号变量、符号常量、符号表达式 1.建立符号变量和符号常数 1）用函数sym（‘’）表示。 2）用命令syms表示多个变量。 2.建立符号表达式：包括代数式、符号方程、符号矩阵等。 3.基本符号运算 与数值运算的区别 ※ 数值运算中必须先对变量赋值，然后才能参与运算。 　 ※ 符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达。 1）符号表达式的四则运算：同数字运算 2）符号矩阵建立 数值矩阵A=[1,2;3,4] A=[a,b;c,d] —— 不识别 用matlab函数sym创建矩阵 命令格式：A=sym([ ]) 3）因式分解与展开 若S为符号表达式或符号矩阵。 factor(S) 对S分解因式 expand(S) 对S进行展开 collect(S) 对S合并同类项 collect(S,v) 对S按变量v合并同类项 simple(S)调用MATLAB的其他函数对S进行综合化简，并显示化简过程 simplify(S) 应用函数规则对S化简（ex01.m） 3.2 符号函数运算 对符号函数进行运算，首先必须： （1）定义符号变量 （2）定义符号函数 注意：定义符号函数之前必须先定义符号变量 1.符号表达式中变量确定 findsym(S,n)——查找函数f中的符号变量，n是查找变量的个数 2.函数求极限 3.3 求解符号方程和符号微分方程 1. 符号代数方程（组）的求解solve 1）定义符号变量 2）定义符号方程组 格式1：[x1,x2,…,xn]=solve(eq1,eq2,…,eqn,x1,x2,…,xn) 格式2： V= solve(eq1,eq2,…,eqn,x1,x2,…,xn) eq1,eq2,…,eqn为符号方程， x1,x2,…,xn为符号变量 2.符号微分方程的求解dsolve 格式1： [y1,y2,…,yn]=dsolve(‘方程1, 方程2’,’条件1, 条件2’,’自变量’) 格式2： S=dsolve(‘方程1, 方程2’,’条件1, 条件2’,’自变量’) 注：当自变量和初条件省略时，默认自变量为t，并含有积分常数 微分符以大写字母D表示，高阶为D2，D3 3.符号变量与数值变量的转换 符号变量的运算结果仍是符号变量，需数值转换 1）vpa(f,d):将符号函数进行数值转换，有效位为d 2) double(f)将不含自由变量的符号函数f作数值变换，有效位固定为5位 注意：符号函数中的自由变量指只含阿拉伯数值字符的符号变量 求下列函数的根 （ex07.m） 3.4 符号级数运算 1）符号的级数求和(ex08.m) symsum(s,v,a,b) v是求和变量，a和b分别为开始项和结束项 2）函数的泰勒级数 taylor(f,v,n,a) 将函数f按变量v展开为泰勒级数，展开到第n项（也就是n-1次幂） * * 3.1 符号计算基础 3.2 符号函数运算 3.3 符号方程和符号微分方程 3.4 符号级数运算 如：syms x y z; f=2*x+x^2*x-5*x+x^3 f=2*x/(5*x) f=(x+y)*(x-y) f =-3*x+2*x^3 f =2/5 f =(x+y)*(x-y) 表达式 输入命令格式 备注 limit (f,x,a) 若a=0，且是对x求极限， 可简写为limit(f) limit (f,x,a,’left’) 趋于a的左极限 limit (f,x,a,’right’) 趋于a的右极限 例：求 ex02.m 3.符号函数微分与积分 1）导数：diff(f,x,n)，n阶导数，缺省值是1 2）不定积分：int(f,x) 3）定积分：int(f,x,a,b) 求 的导数和[0,1]区间上的积分 ex03.m 求解方程组 ex04.m （1）求解微分方程(ex05.m) （2）求解微分方程组( ex06.m) 需要注意，不是任何微分方程都能进行符号微分方程的求解。无解析解的微分方程要用数值求解方法进行