第三讲符号计算优秀课件.pptVIP

第三讲 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 符号级数运算 3.5 积分变换 如：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 需要注意，不是任何微分方程都能进行符号微分方程的求解。无解析解的微分方程要用数值求解方法进行 ex08.m