高等应用数学问题MATLAB求解》.pps
《高等应用数学问题MATLAB求解》.pps》由会员分享,可在线阅读,更多相关《高等应用数学问题MATLAB求解》.pps(33页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第第 1 1 章章计算机数学语言概述计算机数学语言概述薛定宇、陈阳泉著薛定宇、陈阳泉著高等应用数学问题的高等应用数学问题的MATLABMATLAB求解求解,清华大学出版社,清华大学出版社20042004CAICAI课件开发:薛定宇课件开发:薛定宇5/28/20231高等应用数学问题的 MATLAB 求解东北大学信息学院 本章主要内容本章主要内容n n数学问题计算机求解概述n n计算机数学语言概述n n关于本书及相关内容n n本章要点简介n n网络资源简介5/28/20232高等应用数学问题的 MATLAB 求解1.1 1.1 数学问题计算机求解概述数学问题计算机求解概述n n1.1.1为什么要
2、学习计算机数学语言?n n1.1.2 数学问题的解析解与数值解n n1.1.3数学运算问题软件包发展概述5/28/20233高等应用数学问题的 MATLAB 求解1.1.1 为什么要学习为什么要学习计算机数学语言?计算机数学语言?n数学问题求解n n手工推导n n借助计算机n n用数值分析技术,从底层编写起n n应用现成软件进行计算机求解n n解析解与数值解5/28/20234高等应用数学问题的 MATLAB 求解先考虑下面一些例子先考虑下面一些例子【例1-1】高等数学问题:已知函数,如何求导及高阶导数?思路:由分式求导公式,得出 逐次求导则可以得出问题:求导过程很繁杂,容易出错5/28/20
3、235高等应用数学问题的 MATLAB 求解计算机求解结果计算机求解结果不是最简5/28/20236高等应用数学问题的 MATLAB 求解基于计算机的化简结果基于计算机的化简结果n n靠手工推导的方法难以准确得出靠手工推导的方法难以准确得出手工无从推导,计算机能,1秒5/28/20237高等应用数学问题的 MATLAB 求解【例1-2】代数方程求根n n古典方法n n一、二、三、四阶:直接方法一、二、三、四阶:直接方法n n五阶或以上五阶或以上 AbelAbel定理,认为无解定理,认为无解n n现代数值方法n n林士谔林士谔-BairstrowBairstrow 算法,又称为劈因子法算法,又称
4、为劈因子法n n具体实例5/28/20238高等应用数学问题的 MATLAB 求解n n双精度变量的数值方法n n代入方程后的误差n n精确解5/28/20239高等应用数学问题的 MATLAB 求解【例1-3】矩阵行列式求解问题n n代数余子式n n1 1 个个 n n 阶行列式可以表示成阶行列式可以表示成 n n 个个 n-1 n-1 阶行列式的阶行列式的和,和,n n可以将高阶矩阵行列式转换成可以将高阶矩阵行列式转换成1 1阶矩阵行列式阶矩阵行列式n n结论:任意矩阵行列式解析解存在n n问题:忽略了可计算性n nn=20,银河机,3000年5/28/202310高等应用数学问题的 MA
5、TLAB 求解n n实例:Hilbert 矩阵,n=20n n传统数值分析结论:矩阵奇异n n 双精度级别下的数值解-1.1004e-1955/28/202311高等应用数学问题的 MATLAB 求解该矩阵行列式的精确结果该矩阵行列式的精确结果n n近似值n n计算时间 0.2 秒5/28/202312高等应用数学问题的 MATLAB 求解【例1-4】非线性常微分方程组n n没有解析解n n传统数值方法n n计算步长选择计算步长选择n n计算时间计算时间n n变步长变步长n n刚性方程专用方法Van der Pol方程5/28/202313高等应用数学问题的 MATLAB 求解传统数学书上没有
6、的微分方程解法传统数学书上没有的微分方程解法n n延迟微分方程例子n n分数阶微分方程n n解决方法:计算机数学语言算法5/28/202314高等应用数学问题的 MATLAB 求解【例1-5】最优化问题n n求解方法:一般线性规划n n若找不到全局最优解,可以试遗传算法n n若要求 为整数 整数规划5/28/202315高等应用数学问题的 MATLAB 求解其他课程的关系其他课程的关系n n其他数学分支n n积分变换、复变函数、偏微分方程、数据插值积分变换、复变函数、偏微分方程、数据插值与拟合、概率论与数理统计、数值分析与拟合、概率论与数理统计、数值分析n n其他相关课程n n电路、电子技术、
7、电力电子技术、电机与拖动、电路、电子技术、电力电子技术、电机与拖动、自动控制原理自动控制原理5/28/202316高等应用数学问题的 MATLAB 求解1.1.2 数学问题的解析解与数值解数学问题的解析解与数值解n n数学家和其他科学技术工作者的区别n n数学家:理论严格证明、存在性数学家:理论严格证明、存在性n n工程技术人员:如何直接得出解工程技术人员:如何直接得出解n n解析解不能使用的场合n n不存在不存在n n数学家解决方法,引入符号数学家解决方法,引入符号erf(erf(a a)n n工程技术人员更感兴趣积分的值工程技术人员更感兴趣积分的值 数值解数值解5/28/202317高等应
8、用数学问题的 MATLAB 求解n n解析解不能使用的场合n n解析解不存在:无理数,无限不循环小数解析解不存在:无理数,无限不循环小数 p pn n数学家:尽量精确地取值,小日本数学家:尽量精确地取值,小日本6060亿位亿位n n工程技术人员:足够精确即可工程技术人员:足够精确即可n n祖冲之祖冲之 3.14159263.1415926,阿基米德的,阿基米德的3.14183.1418n n解析解存在但不实用或求解不可能n n高阶矩阵行列式高阶矩阵行列式5/28/202318高等应用数学问题的 MATLAB 求解数值解应用场合数值解应用场合n n在力学领域,常用有限元法求解偏微分方程;在力学领
9、域,常用有限元法求解偏微分方程;n n在航空、航天与自动控制领域,经常用到数值线性在航空、航天与自动控制领域,经常用到数值线性代数与常微分方程的数值解法等解决实际问题;代数与常微分方程的数值解法等解决实际问题;n n工程与非工程系统的计算机仿真中,核心问题的求工程与非工程系统的计算机仿真中,核心问题的求解也需要用到各种差分方程、常微分方程的数值解解也需要用到各种差分方程、常微分方程的数值解法;法;n n在高科技的数字信号处理领域,离散的快速在高科技的数字信号处理领域,离散的快速Fourier Fourier 变换变换 (FFT)(FFT)已经成为其不可或缺的工具。已经成为其不可或缺的工具。n
10、n 5/28/202319高等应用数学问题的 MATLAB 求解1.1.3 数学运算问题软件包发展概述数学运算问题软件包发展概述n n享有国际声望的软件包n n线性代数线性代数LINPACKLINPACKn n矩阵特征值计算矩阵特征值计算 LINPACKLINPACKn nNAG(Oxford:Numerical Algorithm Group)NAG(Oxford:Numerical Algorithm Group)n nPress W H,Flannery B P,Press W H,Flannery B P,TeukolskyTeukolsky S A,and S A,and Vitte
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 高等 应用 数学 问题 MATLAB 求解
限制150内