第4章 符号运算功能.pptVIP

例 先看数值积分，例 符号表达式与多项式的转换 sym2poly 多项式的符号表达式转换为按降幂排列的行向量，该函数只能对含有一个变量的符号表达式进行转换 poly2sym 按降幂排列把行向量转换为符号表达式 例 设有连接点O(0,0)和A(1,1)的一段向上凸的曲线弧OA，对于曲线弧上的任一点P(x,y)，曲线弧OP与线段OP所围图形的面积为x2，求曲线弧的方程。 函数taylor不支持字符定义方式 sym2poly函数不支持字符定义格式 4.3 符号方程的求解 4.3.1 代数方程的求解 solve函数求解涉及符号的一般代数方程（组）。调用格式为： solve(f)：f为符号方程，该函数可以求f关于系统默 认变量为自变量的符号方程的解。 solve(f,v)：求关于指定变量v的解。 [a1,a2,…,an]=solve(f1,f2,…,fn) [a1,a2,…,an]=solve(f1,f2,…,fn,’v1,v2,…,vn’) 如果表达式不是一个方程式(不含等号)，则在求解之前自动将表达式置成0。 变量为x 变量为b 解方程组 默认变量为x,y,z 4.3.2 常微分方程的求解 求解常微分方程（组）的函数：dsolve，它不但可以求解无附加条件的常微分方程（组），而且可以求解有附加条件的常微分方程（组）。 MATLAB中用D表示对变量求导数，Dy表示对y求一阶导数， Dny表示对y求n阶导数。 1．求解无附加条件的常微分方程 用dsolve函数求常微分方程通解的调用格式为： (1) dsolve(’f’) (2) dsolve(’f’,’v’) 说明：f为待求解的常微分方程， 第一种调用格式视变量t为自变量进行求解； 第二种调用格式中v为指定变量 C1是积分常数 dsolve函数只支持字符定义方式，而不支持符号定义方式。 dsolve(D2y+2*x=2*y, x) ans = exp(2^(1/2)*x)*C2+exp(-2^(1/2)*x)*C1+x t是自变量 x是自变量 练习 2．求解有初始条件的常微分方程 求解有初始条件的常微分方程时，dsolve函数的调用格式为： dsolve(’f’,’c1,c1,…,cn’,’v’) dsolve(’f’,’c1’,’c2’,…,’cn’,’v’) 其中，f为常微分方程；c1，c2，…，cn 为初始条件；v为指定变量。 * * 第4章 MATLAB的符号运算功能 符号计算对未赋值的符号对象（常数、变量、表达式）进行运算和处理。 数值运算中必须先对变量赋值，然后才能参与运算。符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达。 符号运算不产生数值运算的运算误差，可以在运算最后将数字代入结果，避免了中间运算的误差，能够以指定的精度返回结果。 数值运算与符号运算的区别 4.1.1 创建符号常量 sym 创建单个符号变量， syms 多个符号变量。 4.1 符号对象的创建 创建符号表达式 在创建符号表达式之前，必须把符号表达式所包含的全部符号变量都创建完毕。 4.1.2 创建符号方程 字符串数组 符号方程 4.1.3 创建符号矩阵 创建符号矩阵有以下几种方法： 用sym 直接创建符号矩阵； 用类似创建普通数值矩阵的方法创建符号矩阵； 由数值矩阵转换为符号矩阵。 subs 通用置换, 求符号表达式的值 syms a b c x y; m=sym([1/10 x y; x+y tan(a/b) 2*c ; a b*c 3^2]); m m = [ 1/10, x, y] [ x+y, tan(a/b), 2*c] [ a, b*c, 3^2] mn=subs(m,{a,b,c,x,y},{2,2,2,2,2}) mn = 0.1000 2.0000 2.0000 4.0000 1.5574 4.0000 2.0000 4.0000 9.0000 1．用 sym 直接创建符号矩阵 矩阵行与行之间用“；”隔开，各矩阵元素之间用“，”或“空格”隔开。 可变精度的算术运算 符号运算是对符号表达式的运算，因为不需要进行数值运算，所以无舍入误差。求符号运算结果用函数 eval求其数值，仅在结果转换时会引入一次性的舍入误差。 数值计算精度受到每次计算结果所保留的位数限制。如果保留位数是 16 位，则第 17 位以后的数据将被舍去，所以任何数值运算都会引入舍入误差。重复多次的数值运算还会造成累计误差。 原理上，符号运算可以实现任何数位的运算，但当保留位数增加时，每次计算就需要增