优·04_MATLAB符号运算2011.pptVIP

* MATLAB和R软件 * Matlab 符号作图 空间曲线绘图： ezplot3 缺省的绘图区间为 [0, 2?] ezplot3(x,y,z,[a,b]) ezplot3(x,y,z) * MATLAB和R软件 * Matlab 符号作图 空间曲面绘图： ezmesh、ezsurf ezmesh(z(x,y),[a,b,c,d]) ezmesh(z(x,y),[a,b]) ezmesh(z(x,y)) ezmesh(x(s,t),y(s,t),z(s,t),[a,b,c,d]) ezmesh(x(s,t),y(s,t),z(s,t),[a,b]) ezmesh(x(s,t),y(s,t),z(s,t)) ezsurf 的用法 与 ezmesh 相同 * MATLAB和R软件 * 上机作业 a1=1e10; b1=1e-10; c1=(a1+b1-a1)/b1; a2=sym(a1); b2=sym(b1); c2=(a2+b2-a2)/b2; 指出下面的 M1，M2，M3 分别是什么，并上机验证。 a=1; b=2; c=3;d=4; M1=[a,b;c,d]; M2=[a,b;c,d]; M3=sym([a,b;c,d]); 下面语句计算出来的 c1，c2 相等吗，为什么？上机验证。 补充：class(x) 查看指定变量 x 的类型 * MATLAB和R软件 * 上机作业 简化表达式 设 A 是一个符号矩阵 (定义如下)，试指出 findsym(A,1) 的输出结果，并由此能得出什么结论？ syms a b t u v x y; A=[a+b*x,sin(t)+u; x*exp(-t),log(y)+v] 试指出下面两条命令的结果是否相同，如果不同，哪个是正确的？为什么？ factor(sym(12345678901234567890)) factor(12345678901234567890) * Matlab 演示 * Matlab 演示 * Matlab 演示 * Matlab 演示 MATLAB和R软件 MATLAB符号运算 * MATLAB和R软件 * Matlab 符号运算介绍 Matlab 符号运算是通过符号数学工具箱（Symbolic Math Toolbox）来实现的。Matlab 符号数学工具箱是建立在功能强大的 Maple 软件的基础上的，当 Matlab 进行符号运算时，它就请求 Maple 软件去计算并将结果返回给 Matlab。 Matlab 的符号数学工具箱可以完成几乎所有的符号运算功能。主要包括：符号表达式的运算，符号表达式的复合、化简，符号矩阵的运算，符号微积分、符号作图，符号代数方程求解，符号微分方程求解等。此外，该工具箱还支持可变精度运算，即支持以指定的精度返回结果。 * MATLAB和R软件 * Matlab 符号运算特点 计算以推理方式进行，因此不受计算误差累积所带来的困扰。 符号计算指令的调用比较简单，与数学教科书上的公式相近。 符号计算可以给出完全正确的封闭解，或任意精度的数值解（封闭解不存在时）。 符号计算所需的运行时间相对较长。 * MATLAB和R软件 * Matlab 符号运算举例 求一元二次方程 ax2 + bx + c = 0 的根 solve(a*x^2+b*x+c) 求函数 f (x) = (cos x)2 的一次导数 x=sym(x); diff(cos(x)^2) 计算 f (x) = x2 在区间 [a, b] 上的定积分 syms a b x; int(x^2,a,b) * MATLAB和R软件 * 在进行符号运算时，必须先定义基本的符号对象，可以是符号常量、符号变量、符号表达式等。符号对象是一种数据结构。 符号对象与符号表达式 含有符号对象的表达式称为符号表达式，Matlab 在内部把符号表达式表示成字符串，以与数字变量或运算相区别。 符号矩阵/数组：元素为符号表达式的矩阵/数组。 * MATLAB和R软件 * sym 函数用来建立单个符号变量，一般调用格式为： 符号对象的建立：sym 和 syms 符号对象的建立 例： a=sym(a) 符号变量 = sym(A) 参数 A 可以是一个数或数值矩阵，也可以是字符串 a 是符号变量 b 是符号常量 b=sym(1/3) c 是符号矩阵 c=sym([1 ab; c d]) * MATLAB和R软件 * 符号对象的建立：sym 和 syms 符号对象的建立 syms 命令用来建立多个符号变量，一般调用格式为： syms 符号变量1 符号变量2 ... 符号变量n 例： sym