[理学]第三讲 MATLAB的符号运算.pptVIP

第三讲 MATLAB的符号运算 ——Matlab 不仅具有数值运算功能，还开发了在Matlab环境下实现符号计算的工具包Symbolic Math Toolbox 。符号工具箱在Matlab中引入新的数据类型——符号对象sym。 在数值计算中，包括输入、输出及中间过程，变量都是数值变量，而在符号运算中，变量都以字符形式保存和运算，即使是数字也被当作变量来处理。 符号表达式包括符号函数和符号方程，他们创建的方式是一样的，最简单易用的创建方法和MATLAB字符串变量的生成方法相同。 与数值运算的区别 ※ 数值运算中必须先对变量赋值，然后才能参与运算。 　 ※ 符号运算无须事先对独立变量赋值，运算结果以标准的符号形式表达。 3.1 符号表达式的生成 创建符号函数 f= ‘log(x)’ 创建符号方程 equation=‘a*x^2+b*x+c=0’ 创建符号微分方程 Diffeq=‘Dy-y=x’ 注:由这种方法创建的符号表达式对空格是很敏感的,因此不要在字符间乱加空格符. 符号表达式在MATLAB被看作是1ⅹ1阶的符号矩阵,因此可用sym命令来创建. f=sym(‘sin(x)’) f f=sym(‘sin(x)^2=0’) 另外一种符号函数的创建方法为syms 命令 Sym x F=sin(x)+cos(x) syms a b c %syms可同时定义多个变量 3.2 符号和数值之间的转化 有时符号运算的目的是为了得到精确的数值解，这样就需要对得到的解析解进行数值转换。在symbolic中有三种不同的算术运算： 数值类型 matlab的浮点算术运算 有理数类型 maple的精确符号运算 vpa类型 maple的任意精度算术运算 浮点算术运算 1/2+1/3 －－(定义输出格式format long) ans = 0.83333333333333 符号运算 sym(1/2)+(1/3) ans = 5/6 －－精确解 解析解(精确解)进行数值转化的两个函数. digits(n) —— 函数设置有效个数为n的近似解精度. vpa(x,n) —— 符号表达式在digits(n) 精度下的数值解. digits(25) vpa(1/2+1/3) ans = .8333333333333333333333333 numeric(s): 将不含自由变量的符号表达式转换为数值形式,其效果double(sym(s))相同.(旧版中的数值转化函数). vpa(5/6,40) ans = .8333333333333333333333333333333333333333 a=sym([1/4,exp(1);log(3),3/7]) a = [ 1/4,exp(1)] [log(3), 3/7] vpa(a,10) ans = [.2500000000, 2.718281828] [1.098612289, .4285714286] 3.2 符号矩阵的创建 数值矩阵A=[1,2;3,4] A=[a,b;c,d] —— 不识别 ?用matlab函数sym创建矩阵（symbolic的缩写) 命令格式：A=sym([ ]) ※ 符号矩阵内容同数值矩阵 ※ 需用sym指令定义 ※ 需用 标识 例如：A = sym([a , 2*b ; 3*a , 0]) A = [ a, 2*b] [3*a, 0] 这就完成了一个符号矩阵的创建。 注意：符号矩阵的每一行的两端都有方 括号，这是与 matlab数值矩阵的 一个重要区别。 ?用字符串直接创建矩阵 模仿matlab数值矩阵的创建方法 需保证同一列中各元素字符串有相 同的长度。 例：创建符号矩阵 syms a b c d e f g h A=[a b;c d],B=[e f;g h],C=A*B D=sym([1/(x+t),sin(x);1,exp(x)]) ? 符号矩阵的修改 a.直接修改 可用?、 ?键找到所要修改的矩阵，直接修改 b.指令修改 ①A(行标,列标)＝符号表达式 ②A1=subs(S, ‘new ,‘old) ? 符号矩阵与数值矩阵的转换 将数值矩阵转化为符号矩阵 函数