Matlab：符号计算.PPT

数学软件 Matlab 主要内容 符号运算 Matlab 符号运算 符号运算举例 内容提要 符号对象 符号对象的建立 符号对象的建立 符号表达式 基本符号运算 符号矩阵 内容提要 symvar subs subs 举例 vpa 内容提要 因式分解 函数展开 合并同类项 函数简化 函数简化 函数简化 horner 多项式 计算极限 计算导数 计算积分 符号级数求和 代数方程求解 微分方程求解 其它运算 上机作业 上机作业 上机作业 上机要求 常微分方程求解 y=dsolve(eq1,eq2,...,cond1,cond2,... ,v) 这里 eq1、eq2、. . . 为微分方程  cond1、cond2、...为初值条件，v 为自变量 y 为输出的解。 略（详见 “实验四：求微分方程的解”） 反函数 finverse(f,v) % 求 f 关于指定变量 v 的反函数 finverse(f) % 求 f 关于默认变量的反函数 syms x t; f=x^2+2*t; g1=finverse(f,x) g2=finverse(f,t) 例：计算函数 的反函数 1、指出下面的 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]) 试试 eval(M2) 和 eval(M3) 的输出结果分别是什么？ （自学 eval 的用法） a1=1e10; b1=1e-10; c1=(a1+b1-a1)/b1; a2=sym(a1); b2=sym(b1); c2=(a2+b2-a2)/b2; 2、下面语句计算出来的 c1，c2 相等吗，为什么？上机验证 注：class(x) % 查看变量 x 的数据类型 3、 简化表达式 4、设 A 是一个符号矩阵 (定义如下)，试指出 symvar(A,1) 的输出结果，并由此能得出什么结论？ syms a b t u v x y; A=[a+b*x, sin(t)+u; x*exp(-t), log(y)+v] 5、试指出下面两条命令的结果是否相同，如果不同，哪个是正确的？为什么？ factor(sym(12345678901234567890)) factor(sym(12345678901234567890)) 6、教材第 54 页 第 12 题 * * —— Matlab 符号运算 Matlab 符号运算介绍 符号对象与基本符号运算 symvar、 subs 和 vpa 常见的符号计算 计算以推理方式进行，不受计算误差累积所带来的困扰 符号计算指令的调用比较简单，与教科书上的公式相近 符号计算可以给出完全正确的封闭解，或任意精度的数值解（封闭解不存在时 ） 符号计算所需的运行时间相对较长 符号运算的特点 Matlab 符号运算是通过符号数学工具箱（Symbolic Math Toolbox）来实现的。 Matlab 的符号数学工具箱可以完成几乎所有得符号运算功能，如：符号表达式的运算，符号矩阵的运算，符号微积分，符号作图，符号代数方程求解，符号微分方程求解等。 此外，该工具箱还支持可变精度运算，即支持以指定的精度返回结果。 Matlab 符号运算 求一元二次方程 ax2 + bx + c = 0 的根 solve(a*x^2+b*x+c=0) 求的根 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 符号运算介绍 符号对象与基本符号运算 symvar、 subs 和 vpa 常见的符号计算 在进行符号运算时，必须先定义基本的符号对象，可以是 符号变量、符号表达式等 符号对象是一种数据结构 符号表达式：含有符号对象的表达式称 符号矩阵/数组：元素为符号表达式的矩阵/数组 Matlab 符号对象 sym 用来建立单个符号对象，一般调用格式为： 符号对象的定义/声明：sym、syms 例： a=sym(a) 符号变量 = sym(x) 参数 x 可以是一个数或数值矩阵，也可以是字符串 a 是符号变量 b 是符号常量 b=sym(1/3) C 是符号矩阵 C=sym([1 ab; c d]) syms 符号变量1 符号变量2 ... 符号变量n 例： syms a b c; a=sym(a); b=sym(b); c=sym(c); sy