计算方法

课程简介

张瑞
中国科学技术大学数学科学学院
http://staff.ustc.edu.cn/~rui

课程简介

ieee-float

如何求解数学模型问题？

  • 最理想的结果，可以找到数学模型问题的解析解
  • 在大部分情形下，求解数学模型的解析解很困难，或者解析解很复杂，不利于计算、分析
  • 随着计算机的发展，可以用计算机来求解

计算方法是一种研究并求解数学问题的数值近似解的方法。

具体来说：

  1. 简单或复杂的数学问题。如$\sqrt x$, $\ln x$, $Ax=b$, $\displaystyle\int_a^b f(x)dx$, 微分方程(组) $...$ 等等问题
  2. 构造数值计算算法，将问题变为$+$, $-$, $\times$, $\div$的操作
  3. 编程让计算机计算，并得到结果

课程特点

计算方法是连接数学模型问题和计算机的桥梁，它具有

  • 理论性: 需要对模型问题有些简单的认知（如：解的存在唯一性）；需要对算法的一些特性进行分析
  • 实践性: 需要编程实现算法，

随着计算机技术的发展，计算方法已经广泛应用到包括物理、化学、生物学、经济学等各行各业。

程序作业

本课程有编程的要求：

  1. 编程语言不限。可以使用如：C, Python, Fortran, R, Matlab, Mathematica等， 但不能使用库函数完成主要功能。
  2. 至少使用双精度计算，给出小数点后12位以上的数字。
  3. 每次作业有截止日期，不许拖延。
  4. 提交网址： http://num.hpc.ac.cn

课程内容

本课程只是简单地介绍一些基础的数值分析算法，包括

  • 数值分析: 高等数学问题的数值求解。 如定积分的Newton-Leiberniz公式
    \[\int_a^b f(x)dx=F(b)-F(a) \]
    其中的原函数$F(x)$是理论上存在，但实际中很难找到，甚至不能表达。如$\int \sin\frac1x dx$
  • 数值代数： 线性代数问题的数值求解。如线性方程组的Gram法则
    \[Ax=b\,\quad \Rightarrow \quad x_i=\frac{D_i}{D} \]
    算法可以用计算机编程实现。但计算量太大。当$n=20$时，需要用10亿/秒的计算机算30年。
  • 常微分方程数值解: 只能得到很小部分常微分方程的解析解。对于一般的微分方程，只能寻求数值解。

课程基础

本课程需要掌握一门编程语言，一元微积分和线性代数的相关知识。

  • 连续函数的介值定理
  • (Rolle) 函数$f(x)$在$[a,b]$上连续，在$(a,b)$内可导，且$f(a)=f(b)$， 则存在点$\xi\in(a,b)$满足
    \[f'(\xi)=0 \]
  • (Taylor公式) 设$f(x)$在区间$I$内具有$n+1$阶导数，$x$, $x_0$是$I$内任意两点，存在$x$与$x_0$之间的数$\xi$，满足
    \[\begin{aligned} f(x)=f(x_0)&+f'(x_0)(x-x_0)+\dfrac{f''(x_0)}{2!}(x-x_0)^2+\cdots \\ &+\dfrac{f^{(n)}(x_0)}{n!}(x-x_0)^n+R_n(x) \end{aligned} \]
    其中
    \[R_n(x)=\frac{f^{(n+1)}(\xi)}{(n+1)!}(x-x_0)^{n+1} \]
    $R_n(x)$称为Lagrange余项。
  • (积分中值定理) $f(x)$在$[a,b]$上连续，则存在$\xi\in(a,b)$，使得
    \[\int_a^b f(x)dx=f(\xi)(b-a) \]
  • Gram法则，初等矩阵，初等变换等

目录

本节读完

例 1.

1. 

~~ snippet
% 实际问题的模型抽象
\begin{tikzpicture}[scale=2]
  \node (u) at (0,0)[shape=rectangle, draw] {实际问题};
  \node (y) at (3,0)[shape=rectangle, draw] {数学模型};
  \node (z) at (3,-1.5)[shape=rectangle, draw] {数学结论};
  \node (x) at (0,-1.5)[shape=rectangle, draw] {解释实际问题};

  \draw[->] (u) to node[auto, blue] {简化抽象} (y);
  \draw[->] (y) -- node[auto, blue] {求解} (z);
  \draw[->] (z) -- node[auto, blue] {结果分析} (x);
  \draw[->] (x) -- node[auto, blue] {检验} (u);
\end{tikzpicture}
~~
  1. need to install dvisvgm 2.8 to convert pdf to svg
  2. install 
    ####### compile dvisvgm 2.8
sudo apt install libkpathsea-dev pkg-config libfreetype6-dev libbrotli-dev libwoff-dev
make
  3. install 
    ###### xelatex metapost
sudo apt install texlive-xetex texlive-metapost
####### free chinese fonts
sudo apt install fontconfig fonts-arphic-ukai latex-cjk-chinese-arphic-gbsn00lp fonts-arphic-gbsn00lp latex-cjk-chinese-arphic-gkai00mp fonts-wqy-microhei fonts-wqy-zenhei fonts-noto