对外经济贸易大学MATLAB基础与应用课件-第1章计算机数学语言概述.pptxVIP

对外经济贸易大学信息学院 * 对外经济贸易大学信息学院 * 第 1 章计算机数学语言概述 对外经济贸易大学《MATLAB基础与应用》 * * 本章主要内容 数学问题计算机求解概述 计算机数学语言概述 关于本书及相关内容 本章要点简介 网络资源简介 * * 1.1 数学问题计算机求解概述 1.1.1为什么要学习计算机数学语言？ 1.1.2 数学问题的解析解与数值解 1.1.3数学运算问题软件包发展概述 * * 1.1.1 为什么要学习计算机数学语言？ 数学问题求解 手工推导 借助计算机 用数值分析技术，从底层编写起 应用现成软件进行计算机求解 解析解与数值解 * * 先考虑下面一些例子 【例1-1】高等数学问题：已知函数，如何 求导及高阶导数？ 思路：① 由分式求导公式，得出 ② 逐次求导则可以得出 问题：求导过程很繁杂，容易出错 * * 计算机求解结果 不是最简 * * 基于计算机的化简结果 靠手工推导的方法难以准确得出 手工无从推导，计算机能，1秒 * * 【例1-2】代数方程求根 古典方法 一、二、三、四阶：直接方法 五阶或以上 Abel定理，认为无解 现代数值方法 林士谔-Bairstrow 算法，又称为劈因子法 具体实例 * * 双精度变量的数值方法 代入方程后的误差 精确解 * * 【例1-3】 矩阵行列式求解问题 代数余子式 1 个 n 阶行列式可以表示成 n 个 n-1 阶行列式的和，… 可以将高阶矩阵行列式转换成1阶矩阵行列式 结论：任意矩阵行列式解析解存在 问题：忽略了可计算性 n=20, 银河机，3000年 * * 实例：Hilbert 矩阵，n=20 传统数值分析结论：矩阵奇异 双精度级别下的数值解 -1.1004e-195 * * 该矩阵行列式的精确结果 近似值 计算时间 0.2 秒 * * 【例1-4】非线性常微分方程组 没有解析解 传统数值方法 计算步长选择 计算时间 变步长 刚性方程专用方法 Van der Pol方程 * * 传统数学书上没有的微分方程解法 延迟微分方程例子 分数阶微分方程 解决方法：计算机数学语言＋算法 * * 【例1-5】 最优化问题 求解方法：一般线性规划 若找不到全局最优解，可以试遗传算法 若要求 为整数 －－ 整数规划 * * 其他课程的关系 其他数学分支 积分变换、复变函数、偏微分方程、数据插值与拟合、概率论与数理统计、数值分析 * * 1.1.2 数学问题的解析解与数值解 数学家和其他科学技术工作者的区别 数学家：理论严格证明、存在性 工程技术人员：如何直接得出解 解析解不能使用的场合 不存在 数学家解决方法，引入符号erf(a) 工程技术人员更感兴趣积分的值 数值解 * * 解析解不能使用的场合 解析解不存在：无理数，无限不循环小数 p 数学家：尽量精确地取值，小日本60亿位 工程技术人员：足够精确即可 祖充之 3.1415926，阿基米德的~3.1418 解析解存在但不实用或求解不可能 高阶矩阵行列式 * * 数值解应用场合 在力学领域，常用有限元法求解偏微分方程； 在航空、航天与自动控制领域，经常用到数值线性代数与常微分方程的数值解法等解决实际问题； 工程与非工程系统的计算机仿真中，核心问题的求解也需要用到各种差分方程、常微分方程的数值解法； 在高科技的数字信号处理领域，离散的快速Fourier 变换 (FFT) 已经成为其不可或缺的工具。 …… …… * * 1.1.3 数学运算问题软件包发展概述 享有国际声望的软件包 线性代数LINPACK 矩阵特征值计算 LINPACK NAG (Oxford: Numerical Algorithm Group) Press W H, Flannery B P, Teukolsky S A, and Vitterling W T. Numerical recipes, the art of scientific computing. Cambridge: Cambridge University Press, 1986 * * 软件包作用 从历史发展角度，起了不可替代的作用 对计算机数学语言的强有力支持 但不能过多依赖 使用烦琐 应该在计算机数学语言的意义下利用之 * * 举例：求取矩阵特征值 EISPACK 软件包解法 计算机数学语言解法：eig(A) * * 考虑一个实际编程例子 如何编写一个能求出两个矩阵相乘的计算机通用子程