Ch5_MatLab的符号计算与数值计算.pptxVIP

MatLab的符号计算与数值计算一、MatLab下的符号计算 1.求极限 在MatLab中, 求极限的命令是 基本格式 由于Matlab偏重于数值计算, 而求极限是一个符号运算的过程, 因此在此之前要赋予命令例 求极限 在命令窗口执行命令输出结果为1.例 求极限 求该极限的命令为计算结果: 某些情况下, 我们要求函数的单侧极限. 相应的命令为:例 求极限命令计算结果若要求极限则执行命令计算结果 在高等数学中, 我们知道函数在点 处的极限不存例 求极限在, 但左右极限都存在, 且分别有在MatLab中相应的结果是:极限不存在左极限右极限例 求极限首先定义函数, 再求极限.定义函数求极限对表达式及指定变量求 阶导数.2.求导数 基本格式 功能例 求下列导数:⑴求:⑵求:⑶求:⑴定义函数, 然后求导: 得⑵本题是由参数方程确定的导数问题. 由求导公式:在MatLab下执行操作:定义符号求一阶导数求二阶导数⑶本题为多元函数求导分别求出对两个变量的偏导 由偏导的定义, 知因此, 在MatLab下执行命令3.求积分 在这里所指的积分均指的是不定积分, 即满足及定积分 基本格式 功能对表达式进行积分, 积分变量为 功能对表达式计算定积分.⑴ ⑵例 求积分⑶解 ⑴在MatLab下执行命令: 积分结果为的表达式即:⑵积分结果即:⑶ 我们知道, 在定积分中有公式是奇数.例 计算广义积分命令例 做出函数 的图形, 再计算积分:程序如下:输出计算结果:输出近似值 其中 表示级数求和中的初始项和终止项.4.级数 ⑴级数求和 基本格式例 求下列级数的和.⑴⑵命令如下:输出结果:对函数 在点 处展开到 阶.例 将 在点 处展开到9阶. ⑵函数展开成Taylor级数 基本格式 功能: 命令如下:例 对函数 在点 展开到12阶.的多项式命令如下: Taylor级数逼近计算器 在MatLab中, 有一个用来图示Taylor级数逼近情况的计算器. 在命令窗口中输入:是一个交互式的Taylor级数逼近计算器.输入函数展开的阶数例 求函数 在区间 的阶Taylor级数.5.符号方程求解 ⑴求解代数方程 基本格式: 功能: 求解符号表达式的代数方程.例 求解下列方程:⑴ ⑵ 命令如下: 在MatLab中, 以 表示求导, 而 表示对 的二阶 ⑵符号微分方程求解导数. 例如:表示微分方程说明: 格式中的 表示定解条件. 无此项时表示求方程 基本格式功能:对表达式确定的微分方程求解.的通解.例 求解下列微分方程:⑴⑵方程⑴为贝努里方程, 特解为:方程⑵为二阶常系数线性微分方程. 通解为例 求解微分方程组其中 由6.应用举例例 求二重积分围成.解 先作出区域的图形.在MatLab下求曲线的交点.得交点后分别做两个两次积分, 即得到问题之解 设进货量为 需求量为 则收益函数为:3.一煤炭部门煤的进价为每吨65元, 零售价为每吨70元,若当年卖不出去, 则第二年要削价20%, 若供应缺货, 有关部门每吨罚款10元, 已知顾客对煤炭需求量求最优策略.这里 均满足其它再由已知条件:故相应的概率密度函数为由此得到收益的期望值得积分结果为:对其求导并令其为零, 得函数惟一的驻点:相应程序如下:每年进货 万吨.即最优策略为:二、数值解法1.数值积分 数值积分方法很多, 经常使用的是梯形法, 辛普森法, 及Newton—Cotes方法. 上章介绍的积分法就是属于梯形法. 数值积分的意义: 在上章中, 我们用梯形公式计算函数在区间上的定积分,其方法是将区间做了一千等分, 在每个小区间上用梯形面积近似替代曲边梯形面积. 对于一个较为复杂的函数而言, 其计算量将是非常巨大的, 数值积分的思想是:对每一次的计算结果加以判定: 本次的计算结果是否足够精确, 若没有达到预定目标, 则进行下一次的计算. 并尽可能利用已有的结果. 设函数 积分区间为 ⑴变步长梯形法记:时的变步长梯形法图示若 则迭代结束