02Matlab符号运算.ppt

数学实验 邓化宇 上海电力学院数理系dhy0826@126.com MATLAB符号运算 Matlab 符号运算介绍 Matlab 符号运算特点 Matlab 符号运算举例 符号对象与符号表达式 符号对象的建立 符号对象的建立 符号表达式的建立 符号表达式的建立： 符号对象的基本运算 符号表达式的替换 用给定的数据替换符号表达式中的指定的符号变量 subs 举例 计算函数值 符号表达式的数值化 对符号求值的命令为vpa，即Variable precision arithmetic 上机作业 因式分解 函数展开 合并同类项 函数简化 函数简化 函数简化举例 分式通分 计算极限 计算极限 计算导数 计算导数 计算积分 计算积分 符号求和 Taylor级数展开 Taylor级数展开 代数方程和微分方程求解 其它运算 上机作业 上机作业 1. 2. 3. 4. 5. 上机作业提高练习 1.一元函数或二元函数最值问题：驻点 2.利用函数的幂级数展开式来求pi的近似值。 数值积分 (1)自适应Simpson法 函数quad()采用递推自适应Simpson法计算积分，精度较高，较常使用。 q = quad(fun,a,b) q = quad(fun,a,b,tol) 功能：求一元函数fun的积分ｑ，积分上限和下限分别为a和b，绝对误差为tol，默认时取值10^-6 ，函数fun是待求解的对象，可以是字符串、内联函数、Ｍ函数文件名的函数句柄。 数值积分 数值积分 数值积分 数值积分 (2)梯形法（被积函数由一个表格定义） trapz函数采用梯形法求取数值积分，适用于由表格形式定义的函数关系的求定积分问题，求值速度快，但精度差。 q = trapz(Y) q = trapz(X,Y) 其中向量X,Y定义函数关系Y=f(X)。 【例】求exp(-x2)在[0,1]上的积分。 上机作业 取一定数量的点，利用梯形数值积分公式近似计算圆周率pi 卫星的轨道长度 我国第一颗人造地球卫星近地点 距地球表面为h=439km，，远地点距地球表面H=2384km，地球半径为R=6371km，求该卫星的轨道长度。 卫星的轨道长度 卫星的轨道长度 例：试对正弦函数y=sinx进行泰勒幂级数展开，观察不同阶次下的近似效果 x0=-2*pi:0.01:2*pi;y0=sin(x0); syms x;y=sin(x);plot(x0,y0); axis([-2*pi 2*pi -1.5 1.5]);hold on for n=8:2:16 p=taylor(y,x,n), y1=subs(p,x,x0); line(x0,y1) end 代数方程求解（见后） solve(f,v)：求方程关于指定自变量的解，f 可以是用字符串表示的方程、符号表达式或符号方程； solve 也可解方程组(包含非线性)； 得不到解析解时，给出数值解。 微分方程求解（见后） dsolve和ode45 反函数 finverse(f,v)：求 f 关于指定变量 v 的反函数 finverse(f)：求 f 关于默认变量的反函数 syms x t; f=x^2+2*t; g1=finverse(f,x) g2=finverse(f,t) 例：计算函数 的反函数 简化表达式 试指出下面两条命令的结果是否相同，如果不同，哪个是正确的？为什么？ factor(sym(12345678901234567890)) factor(12345678901234567890) 【例】求exp(-x2)在[0,1]上的积分。 ①f=inline(exp(-x.^2));％内联函数 q1=quad(f,0,1) 【例】求exp(-x2)在[0,1]上的积分。 ②编写被积函数表达式，函数名为quad1.m function f = quad1(x) f=exp(-x.^2); q2=quad(quad1,0,1) 【例】求exp(-x2)在[0,1]上的积分。 ③ q3=quad(exp(-x.^2),0,1) X=0:0.01:1; Y=exp(-X.^2); q=trapz(X,Y) 问题 * * 符号对象的建立 1 表达式简化 2 符号微积分 3 级数求和 4 卫星的轨道长度 5 Matlab 符号运算是通过符号数学工具箱（Symbolic Math Toolbox）来实现的。Matlab 符号数学工具箱是建立在功能强大的 Maple 软件的基础上的，当 Matlab 进行符号运算时，它