MATLAB第一章要点.ppt

Slide 41 (of 41) MATLAB 语言与科学运算 东北大学信息学院 薛定宇 xuedingyu@ 教学邮箱： matlab_world@126.com, 密码：matlab 课程内容与教材 课程内容 MATLAB语言 科学运算问题的计算机求解 课程教材 薛定宇、陈阳泉《高等应用数学问题的MATLAB 求解》，第三版，清华大学出版社，2013, (2008, 2004)， CRC Press, 2008 本书全套交互式PPT课件、例题、代码 (教学网站) MATLAB手册， 论坛：； 国家精品“资源共享课”（控制） 第1章计算机数学语言概述 第1章 计算机数学语言概述 为什么学习MATLAB语言 数学问题计算机求解概述 计算机数学语言概述 关于本书及相关内容 其他演示的内容 MATLAB语言的科学可视化 MATLAB在自动化专业课程与研究中应用 1.1 数学问题计算机求解概述 为什么要学习计算机数学语言 数学问题的解析解与数值解 数学运算问题软件包发展概述 常规计算机语言的局限性 为什么要学习科学运算？ 在理工科课程与研究中需要数学基础 很多问题直接转换成数学问题 电路、结构设计——线性代数 控制、化工、环境、物理——微分方程 管理、经济、金融、系统工程——最优化 医学、生物学、实验科学——数理统计 …… 控制科学的研究依赖于数学建模与运算 为什么要学习计算机数学语言 数学问题求解 手工推导 过程繁琐，需要强大的数学功底、经验 工作量巨大，只能求解简单问题 易出错，需要反复验证 借助计算机 用数值分析技术，从底层编写起 应用现成软件进行计算机求解 解析解与数值解 例 1.1 高等数学问题：已知函数 如何求4阶导数 例 1.2 Abel-Ruffini定理：5次或以上的多项式方程没有通用的解析解求解方法。 数值方法 林士谔-Bairstrow算法，又称为劈因子法 /BairstowsMethod.html 解如下多项式方程 双精度变量的数值方法 结果 检验 精确解 例 1.3 矩阵行列式求解问题 代数余子式 1 个 n 阶行列式可以表示成 n 个 n-1 阶行列式的代数和，… 可以将高阶矩阵行列式转换成1阶矩阵行列式 结论：任意矩阵行列式解析解存在 问题：忽略了可计算性 n=25, ，每秒二亿亿次巨型机59年 实例：Hilbert 矩阵，n=20 传统数值分析结论：矩阵奇异 双精度级别下的数值解 解析解 精确解 例 1.4 微分方程的解 Van der Pol方程，没有解析解 刚性方程 利用传统课程的微分方程求解算法不能求解 利用MATLAB语言的求解语句 其他类型的微分方程 延迟微分方程 分数阶微分方程 在现有微分方程课程中几乎未提及，但实际应用中可能用到，如何求解？ 例 1.5 线性规划问题 MATLAB代码 混合整数规划问题 整数规划的解 混合整数规划可以通过定义intlist类型的变量来求取 例 1.6 高等应用数学分支 积分变换 复变函数 偏微分方程 数据插值与拟合 概率论与数理统计 数值分析 考试之后还记得其中问题的求解方法吗？ 例 1.7 新的数学分支 模糊集合与粗糙集合 人工神经网络 其他相关课程 电路 电子技术 电力电子技术 电机与拖动 自动控制原理 1.1.2 数学问题的解析解与数值解 数学家和其他科学技术工作者的区别 数学家：理论严格证明、存在性 工程技术人员：如何直接得出解 解析解不能使用的场合 不存在，例如 数学家解决方法：引入符号erf(a) 工程技术人员解决方法：查表法，得出近似解 解析解不能使用的场合 解析解不存在：无理数，无限不循环小数 p 数学家：尽量精确地取值。例如，日本60亿位 工程技术人员：足够精确即可 祖冲之 3.1415926，阿基米德的3.1418 /Pi.html 解析解存在但不实用或求解不可能 高阶矩阵行列式 解析解与数值解数据结构：符号型、双精度 数值解应用场合 在力学领域，常用有限元法求解偏微分方程； 在航空、航天与自动控制领域，经常用到数值线性代数与常微分方程的数值解法等解决实际问题； 工程与非工程系统的计算机仿真中，核心问题的求解也需要用到各种差分方程、常微分方程的数值解法； 在高科技的数字信号处理领域，离散的快速Fourier 变换 (FFT) 已经成为其不可或缺的工具。 …… …… 1.1.3 数学运算问题软件包发展概述 享有国际声望的软件包 线性代数LINPACK 矩阵特征值计算 EISPACK NAG (Oxford: Numerical Algorithm Group) Press W H, Flannery B P, T