MATLAB与Mathematica
文章平均质量分 94
数学计算语言
陆嵩
中国科学院数学与系统科学研究院 计算数学与科学工程计算研究所 科学与工程计算国家重点实验室
State Key Laboratory of Scientific and Engineering Computing,
Institute of Computational Mathematics and Scientific/Engineering Computing,
Academy of Mathematics and Systems Science,
Chinese Academy of Sciences
Email: lusong@lsec.cc.ac.cn
展开
-
MATLAB 工具箱傻瓜式求解 NS(Navier Stoke)方程
可惜的是,MATLAB 的自带的工具箱无法对它进行很好地求解,主要原因在于非线性项(对流项)无法处理,如果把它和右端载荷项放在一块处理,显然 MATLAB 工具箱是没办法支持这件事情的。打开工具箱,选择圆柱绕流算例,就是这个例子,点击 Run,它会自动地给你演示如何做。这个工具够傻瓜够强大,能求解很多 CFD 的问题,能和 FEniCS 和 Openfoam 联动,求足够完整的文档,以及完备的社区。安装就不用说,一搜就有了。需要指明的是,这份代码的运行需要在工具箱打开,变量载入空间的时候才能 work。原创 2022-10-18 13:05:46 · 5235 阅读 · 3 评论 -
核酸检测多少人为一组混检合适?
但是对于被检测者来说,不管是混检还是单检,做一次核酸真的就是要付出做一次核酸的成本,排一次队,刮一次鼻子,遭一次罪。有人说,可以再采用的时候,把每个人的样本分两份,当检测出患病的时候,再把备份拿出来单独测,这样就不需要喊人来车第二次了。考虑到患病率极低,这种分两份的人工成本是极高的,已经远远超过了把一起混检的人喊来再测一次的成本,更远远超过了不采用最优参数。其实很简单,2个管子,6 个人,每个管子 3 个人,2 个患病的人进入同一个管子的概率不是 1/2,而是。是一件非常麻烦的事情,不是初等数学可以做的。原创 2022-10-14 10:13:32 · 15598 阅读 · 22 评论 -
MATLAB Tensor 和 N-Way 工具箱张量补全功能评测
既然现有的方法这么垃圾,我们就基于比较快的分解方法,做一个改进。一个小改进就是在做 cp_als 最小二乘交替步的时候,直接在已知数值的位置上,每步迭代用已知值去替代它。Nway 工具箱实现的非常 stupid,使得速度非常慢,为什么 Nway 实现得这么 stupid,举个例子,它算一堆矩阵 KR 除某个之外的乘积的时候,一个一个算的,一点都不明智。ALS 方法就是先固定 B 和 C 以求解 A,然后再固定 A 和 C 求解 B,接着固定 A 和 B 以求解 C,并继续重复整个过程,直到收敛标准被满足。原创 2022-10-11 13:16:58 · 1906 阅读 · 2 评论 -
MATLAB 张量工具箱的使用
中的所有列都拿出来,做 Kronecker 内积之后求和,求和的系数来自于张量$ \mathcal{G}$ 的。运行的结果上可以看出,原来缺失的位置补全的效果可能不是非常地好。MATLAB 的 Tensor Toolbox 是一个强大的张量工具,它提供了一些强大的张量函数。是完全的,没被噪声污染的数据,我们称之为精确解,用上述的张量形式存储。给出了补全矩阵和原来的矩阵的一个误差表示,它不仅仅考虑缺失的部分。表示的是 CP 分解求和符号的上标,也是因子的个数。的 分量,并且可以执行许多操作,例如。......原创 2022-08-04 08:12:29 · 7652 阅读 · 2 评论 -
零基础使用 MATLAB 求解偏微分方程(建议收藏)
零基础使用 MATLAB 求解偏微分方程(建议收藏)文章目录零基础使用 MATLAB 求解偏微分方程(建议收藏)偏微分开源工具介绍PDE 工具箱函数汇总介绍0 基础:GUI 界面操作示例问题工具箱求解导出为代码形式代码导出相关数据0.1 基础:编程调用 PDE 工具箱偏微分开源工具介绍百分之九十以上的重要的工程和数学科学研究,和偏微分方程都脱不开关系。在所有的偏微分方程中,百分之九十九都是没有解析解的。没有解析解怎么办,我们只能通过有限元或者有限差分等方法,求解偏微分方程数值解。如果您有一些代码基础,原创 2021-10-02 01:30:07 · 45148 阅读 · 4 评论 -
MATLAB 提取图片曲线
文章目录我用 MATLAB 提取图片曲线行文动机图像的读入与裁剪颜色拾取颜色转换与色差计算分离曲线二值化,提取数据数据点分类与排序后话我用 MATLAB 提取图片曲线给你一张图片,如何提取里面曲线的数据,从而利用这些数据进行图像重绘、加工处理、测距、拟合得到函数表达式等操作呢?行文动机前段时间,有个朋友问了我一个问题,大概意思就是要给图像的流线测距离,在我的印象里面,MATLAB 是似乎没有这种直接的功能的。那么换个角度来理解一下这个问题,如果给你一张图像,如何提取里面点的数据?其实,有了曲原创 2021-08-16 23:09:19 · 11571 阅读 · 20 评论 -
常见的使用水平集函数的零水平集表示的曲面汇总(3D)(MATLAB 画三维图)
常见的使用水平集函数的零水平集表示的曲面汇总(3D)文章目录常见的使用水平集函数的零水平集表示的曲面汇总(3D)Donut甜甜圈Atom原子Popcorn爆米花Tooth牙齿Heart爱心Apple苹果Orthocircle正交圆Swiss cheese block瑞士奶酪块Heart心形Flower花代码感谢周同学和张师兄的公式提供。Donut甜甜圈水平集函数:ϕ(x,y,z)=(x2+y2−0.6)2+z2−0.09\phi(x, y, z)=\left(\sqrt{x^{2}+y^{2原创 2021-07-25 21:27:02 · 2763 阅读 · 8 评论 -
【小波分析】七、小波分析与非线性逼近(下)
【小波分析】七、小波分析与非线性逼近(下)文章目录【小波分析】七、小波分析与非线性逼近(下)洛伦兹空间小波系数刻画函数空间向 LpL_pLp 扩展小波系数刻画 LpL^pLp小波系数刻画 Besov 空间非线性小波逼近n 项逼近引入n 项逼近的 Jackson 不等式洛伦兹范数下的 Jackson 不等式Besov 范数下的 Jackson 不等式n 项小波逼近的 Bernstein 不等式一般的 n 项逼近推广到多变量情形n 项逼近和插值空间n 项逼近和 k 泛函有界区域上的 n 项逼近n 项逼近 T原创 2021-05-29 22:22:07 · 706 阅读 · 2 评论 -
MATLAB 求函数极值的内置函数一览表(实则优化算法函数汇总)
MATLAB 求函数极值的内置函数一览表收集了几乎所有的 MATLAB 内置的优化函数,可收藏,需要时查阅。文章目录MATLAB 求函数极值的内置函数一览表简介概览求函数极值(或最大值最小值)求零点(解方程)最小二乘问题求极值fminbnd:单变量fmincon:约束、非线性、多变量fminunc:无约束、多变量fminsearch:无约束、多变量、无导数linprog:线性规划quadprog:二次规划fminimax:minmax 问题fgoalattain:目标达到fseminf:半无限约束、原创 2021-05-22 17:38:09 · 5824 阅读 · 4 评论 -
【小波分析】三、正交小波的构造
【小波分析】三、正交小波的构造文章目录【小波分析】三、正交小波的构造内容回顾傅里叶变换常用性质正交小波典例: 哈尔小波香农小波shannon 采样地理考虑 B 为 PI考虑 B 为 PI 的 2 的整数次幂倍逼近关系再议正交小波基的转移正交小波基的构造正交多分辨分析内容回顾回顾一下之前的内容。前面我们从线性代数谈起,引入了"变换即找基"的观念。接着我们谈了一下小波分析的历史,从傅里叶得到法国科学院的大奖却没拿到奖金一脸委屈到哈尔基,再到加窗傅里叶变换,再到小波分析的萌芽。再之后,我们就开始正式地介绍原创 2021-05-17 00:32:55 · 4098 阅读 · 0 评论 -
[有限元方法阶段汇总篇] 有限元入门简单 1D 示例程序(Helmholtz 方程)
[有限元方法阶段汇总篇] 有限元入门简单 1D 示例程序(Helmholtz 方程)文章目录[有限元方法阶段汇总篇] 有限元入门简单 1D 示例程序(Helmholtz 方程)前言一些链接动机问题描述关于特征值问题亥姆霍兹方程程序与结果问题一问题二写在后面的话前言一些链接之前写过三篇基础的有限元基础入门级别的程序介绍,引起比较大的反响,很多人在 CSDN 私信和添加我微信讨论了相关问题,有限元入门见如下链接。有限元方法入门:有限元方法简单的一维算例有限元方法入门:有限元方法简单的二维算例(三角形原创 2020-12-26 02:05:39 · 4165 阅读 · 8 评论 -
[MATLAB 在科学计算中的应用] 使用MATLAB 进行非线性拟合
[MATLAB 在科学计算中的应用] 使用MATLAB 进行非线性拟合文章目录[MATLAB 在科学计算中的应用] 使用MATLAB 进行非线性拟合前言引述MATLAB 曲线拟合函数简述一二维数据非线性拟合一维数据拟合例子二维数据拟合例子高维数据非线性拟合lsqcurvefitnlinfit 函数数据拟合工具箱的使用写在后面的话前言引述所谓曲线拟合,就是给定数据和数据可能满足的带参函数表达,希望能找到一组最好的参数,是的这个函数能够“最大限度”地刻画给定的数据。在一维(自变量维数)的情况下,换句话说,原创 2020-12-22 01:07:39 · 6349 阅读 · 0 评论 -
数据科学和机器学习中的优化理论与算法(上)
数据科学和机器学习中的优化理论与算法(上)数据科学和机器学习当前越来越热,其中涉及的优化知识颇多。很多人在做机器学习或者数据科学时,对其中和优化相关的数学基础,包括随机梯度下降、ADMM、KKT 条件,拉格朗日乘数法、对偶问题等,不太了解,一知半解地用,用着用着就出错了。本文希望从基础知识的角度,尽可能全地对最简单的优化理论和算法做一个小结。内容涵盖以下几个方面:优化简介、无约束优化、线搜索方法、信赖域方法、共轭梯度方法、拟牛顿方法、最小二乘问题、非线性方程、约束优化理论、非线性约束优化算法、二次规划原创 2020-12-21 15:28:46 · 1532 阅读 · 2 评论 -
数学分析中的点荟萃(持续补充修缮)
数学分析中的点荟萃(持续补充)文章目录数学分析中的点荟萃(持续补充)导数篇常用求导公式导数技巧微分中值定理基本认知积分篇积分技巧常用积分表处理原则积分表格林公式和高斯公式格林公式高斯公式积分等式积分中值定理积分第一中值定理积分第二中值定理广义积分柯西判别法无穷限积分瑕积分极限和连续篇极限技巧基本认知二阶曲面方程函数篇泰勒展开定义常见函数泰勒级数一些重要的特定函数符号函数黎曼函数狄利克雷函数分段三角函数分段指数函数其他函数双曲函数二元函数Γ\GammaΓ和β\betaβ函数常用等价无穷小一般stirling原创 2020-10-04 20:20:27 · 775 阅读 · 0 评论 -
关于两个周期函数的和的周期性的讨论
关于两个周期函数的和的周期性的讨论因为排版和敲数学公式的局限性,很多地方写得并不是非常严格,或者有些跳跃,望海涵。初衷想这个问题的初衷是在给同学们习题课的时候(华东师大版的数学分析),里面有一道题,如下:求下列函数的周期:cosx2+2sinx3\cos \frac{x}{2}+2 \sin \frac{x}{3}cos2x+2sin3x。这道题本身比较简单,显然 12π12\pi12π 是它的一个周期,如果这里的周期理解为基本周期(最小正周期)的话(有同学发问了),我们还得 chec原创 2020-10-03 23:51:50 · 22041 阅读 · 9 评论 -
汉诺塔游戏玩法介绍(攻略和编程实现)
汉诺塔游戏玩法介绍(攻略和编程实现)汉诺塔游戏介绍汉诺塔是根据一个传说形成的数学问题。描述如下:有三根杆子 A、B、C 。 A 杆上 N 个 穿孔圆盘,盘的尺寸由下到上依次变小。要求按下列规则将所有圆盘移至 C 杆:每次只能移动一个圆盘大盘不能叠在小盘上面可将圆盘临时置于 B 杆,也可将从 A 杆移出的圆盘重新移回 A 杆,但都必须遵循上述两条规则。问:如何移?最少要移动多少次?最...原创 2020-03-04 03:08:52 · 9049 阅读 · 0 评论 -
交替方向乘子法(ADMM)的数学基础
交替方向乘子法(ADMM)网上的一些资料根本就没有把ADMM的来龙去脉说清楚,发现只是一个地方简单写了一下流程,别的地方就各种抄,共轭函数,对偶梯度上升什么的,都没讲清楚,给跪了。下面我来讲讲在机器学习中用得很多的ADMM方法到底是何方神圣。共轭函数给定函数f:Rn→Rf: \mathbb{R}^{n} \rightarrow \mathbb{R}f:Rn→R,那么函数f∗(y)=max...原创 2019-12-31 19:23:18 · 4661 阅读 · 6 评论 -
函数型数据主成分分析(FPCA)
本文主要介绍了以下几个方面的内容:简单介绍了经典的主成分分析方法,包括其数学推导,算法步骤,和几个实际算例;简单介绍了其它的数据降维方法,譬如局部线性嵌入以及它的简单算例;更近一步,我们介绍了函数型主成分分析方法(FPCA),包括其基本思想、数学推导、算法描述等,最为重要的是,我们将该方法和本领域进行结合,有了一些新的思考,感谢"数据科学与矩阵优化"课程给带来的灵感。文章目录前言主成分分析(...原创 2019-11-24 19:57:27 · 18939 阅读 · 21 评论 -
zzu数学 实验一泰勒级数
zzu数学 实验一泰勒级数泰勒逼近(*任意函数的泰勒逼近*)n = 10;(*定义逼近程度*)f = Tan[x];(*任给一个函数*)taylor[f_, x0_, n_] := Sum[(D[f, {x, k}]; D[f, {x, k}] /. x -> x0)/k!*(x - x0)^k, {k, 0, n}];(*定义泰勒展开函数*)figures = Table[taylor[f原创 2016-08-24 21:45:54 · 1078 阅读 · 0 评论 -
zzu数学 实验零入门测试
zzu数学 实验零入门测试入门练习代码版本一Animate[Clear[d, a, t];(*清除变量*) a = Graphics[Line[{{-1, 2}, {0, 0}, {2, 4}}]];(*绘制两根初始线段*) l = {-1, 2}/n; r = {2, 4}/n;(*将两线段n等分*) t = Table[{ i*l, (n - i + 1)*r}, {i, 1, n原创 2016-08-24 21:50:10 · 506 阅读 · 0 评论 -
zzu数学 实验二圆周率pi的计算
zzu数学 实验二圆周率pi的计算示例代码:Select[Eigenvalues[RandomReal[1, {5, 5}]], Abs[#] < 1 &](*挑选特征值小于1的矩阵特征值*)app[n_] := 4 N[Length[Select[RandomReal[{-1, 1}, {n, 2}], Dot[#, #] <= 1 &]]/n](*dot这里表示的就是点的长度其实因为l原创 2016-08-24 22:02:13 · 1035 阅读 · 0 评论 -
zzu数学 实验三最佳分数近似国歌
zzu数学 实验三最佳分数近似国歌国歌编曲15 _ = 16 _ = 2^(2/12)7 _ = 2^(4/12)do = 2^(5/12)re = 2^(7^12)mi = 2^(9/12)sol_ = 1la_ = 2^(2/12)si_ = 2^(4/12)do = 2^(5/12)re = 2^(7^12)mi = 2^(9/12)sol = 1la = 2^(2/1原创 2016-08-24 22:08:08 · 583 阅读 · 0 评论 -
zzu数学 实验四数列之3n+1问题
zzu数学 实验四数列之3n+1问题源代码示例series[a_] := Module[{a1 = a, an = {a}}, Do[a1 = If[Mod[a1, 2] == 1, 3 a1 + 1, a1/2]; AppendTo[an, a1]; If[Length[an] >= 3 && an[[-1]] == 1 && an[[-2]] == 2, Break[]原创 2016-08-24 22:33:59 · 1263 阅读 · 0 评论 -
zzu数学 实验五素数问题
zzu数学 实验五素数问题素数问题DivPrime[n_]:=Block[{primeset={2,3}},For[k=4,k<=n,k++,Block[{i,flag=True},For[i=1,primeset[[i]]<=Sqrt[k]&&i<=Length[primeset],i++,If[Mod[k,primeset[[i]]]==0,flag=False;Break[]]];I原创 2016-08-24 22:39:53 · 607 阅读 · 0 评论 -
zzu数学 实验六骰子问题
zzu数学 实验六骰子问题骰子实验程序prob[m_, n_] := Tuples[Range[n], m];A = prob[3, 6];For[i = 1, i <= Length[A], i++, A[[i]] = Sort[A[[i]]]];Length[Union[A]]Tally[A]prob[m_, n_] := Tuples[Range[n], m];A = prob[3,原创 2016-08-24 22:41:44 · 850 阅读 · 0 评论 -
zzu数学 实验七几何变换
zzu数学 实验七几何变换在如photoshop等的各大图形处理软件中,我们能看到各种图形变换,斜切,对称,翻转等,其妙无穷,我们却对其原理一无所知。其实,并不难,今天我以mathematica为例。做两个简单的示例。示例1 ** Linear Transformation **Clear[f]; t = Pi/6;a1 = Cos[t]; b1 = -Sin[t]; a2 = Si原创 2016-11-06 10:07:11 · 1164 阅读 · 0 评论 -
zzu数学 实验八物理现象之模拟电场线
zzu数学 实验八物理现象之模拟电场线 大家都学过高中物理,看到两个点电荷之间的电场线,等势线等,有没想过这些图是怎么画出来的?如果是随手画出来的,为什么会这么精准?今天就带大家走进,电场线的绘制。提供两块代码作为对比参考。代码一 ** Physics ** ** wave **x0 = 0.0; v0 = 9原创 2016-11-06 10:22:41 · 2300 阅读 · 1 评论 -
zzu数学 实验九迭代一方程求解
方程求解的方法有很多啦,最常用,使用最广泛的就是迭代法,原理好理解,效果也不差。一般方程求解g[x_] := x^3 - 12 x + 2;f = x - g[x]/D[g[x], x];x0 = -3;series = {x0};For[i = 1, i <= 5, i++, x0 = f /. x -> x0; AppendTo[series, x0]]series = N[ser原创 2016-11-06 10:33:11 · 814 阅读 · 0 评论 -
zzu数学 实验十一最速降线
zzu数学 实验十一最速降线 意大利科学家伽利略在1630年提出一个分析学的基本问题──“一个质点在重力作用下,从一个给定点到不在它垂直下方的另一点,如果不计摩擦力,问沿着什么曲线滑下所需时间最短”。他说这曲线是圆,可是这是一个错误的答案。 瑞士数学家约翰.伯努利在1696年再提出这个最速降线的问题(problem of brachistochrone),征求解答。次年已有多位数学家得到正确答案原创 2016-11-06 10:55:20 · 4597 阅读 · 2 评论 -
Mathematica画图函数命令大全
二维作图 Plot[f,{x,xmin,xmax}] 一维函数f[x]在区间[xmin,xmax]上的函数曲线 Plot[,f2.{f1.},{x,xmin,xmax}] 在一张图上画几条曲线 ListPlot[{y1,y2,..}] 绘出由离散点对(n,yn)组成的图 ListPlot[{{x1,y1},{x2,y2},..}] 绘出由离散点对(xn,yn)组成的图 Plarametri原创 2016-11-06 11:05:31 · 36397 阅读 · 0 评论 -
Mathematica精华
Clear[d, a, n, t];n = 50;a = Graphics[Line[{{n - 1, 0}, {0, 0}, {0, n - 1}}]];d = Array[f, n - 1];t = Table[{{0, i}, {n - 1 - i, 0}}, {i, 1, n - 1}];For[i = 1, i <= n - 1, i++, d[[i]] = Graphics[L原创 2016-11-06 11:07:41 · 812 阅读 · 0 评论 -
Mathematica做拟合
打开mathematica键入如下内容:advol={{0.49,100},{0.69,150},{0.88,200},{1.06,250}}; line = Fit[advol, line = Fit[advol, line = Fit[advol, line = Fit[advol, line = Fit[advol, line = Fit[advol, {1,x,x^2},x] 按 Shi原创 2016-11-06 11:10:07 · 14412 阅读 · 0 评论 -
函数列的一致收敛性
设fn(x)=(1+xn)n,n=1,2,3…f_n(x)=(1+\frac{x}{n})^n,n=1,2,3…研究{fn(x)f_n(x)}在[0,σ](σ>0)[0,\sigma](\sigma>0)上的一致收敛性。Clear[x,n];f[x_,n_]:= (1+x/n)^n;Animate[Plot[{f[x,n],E^x},{x,-3,3},PlotRange->{-3,10},Pl原创 2017-02-07 22:17:38 · 7090 阅读 · 1 评论 -
从《货殖列传》中学习炒股
从《货殖列传》中可以看出很多经济学的基本思想,是一篇作为学金融学经济的人不可不读的千古好文。其中有一句话是这样说的:”是故其智不足与权变,勇不足以决断,仁不能以取予,强不能有所守,虽欲学吾术,终不告之矣。”这句话,告诉我们股票投资最重要的四个品质:智、勇、仁和强。 智指的是权变,就是权衡变数,要求我们能够判断何时为“利好出尽”何时为“利空出尽”。勇指的是决断,这就不必多说了。仁就是买卖,要懂得买卖原创 2017-02-22 23:22:11 · 882 阅读 · 0 评论 -
matlab入门练习之杂例二
留数计算函数residue的使用clcclearb=[3 5 2 9]a=[4 6 0 -5][r,p,k]=residue(b,a)逻辑运算符的使用a=[1,2,-3,0,0]b=[0,1,0,2,0]xor(a,b)any(a)%一真即真,相当于或运算all(a)isequal(a,b)ismember(a,b)matlab中的循环及其控制clear; clc ;c原创 2017-01-25 23:56:03 · 793 阅读 · 0 评论 -
小观matlab插值函数
插值与拟合是数值分析中的重要内容,我们来看看matlab中的插值函数。空间插值[x,y,z,v]=flow(8); % flow是MATLAB自带的测试函数slice(x,y,z,v,[3,5],2,[-2,3]) % 画切片图title('插值前')[xi,yi,zi]=meshgrid(0.1:0.25:10,-3:0.25:3,-3:0.25:3原创 2017-01-25 23:22:49 · 862 阅读 · 0 评论 -
matlab入门练习之杂例一
提供几个入门杂例,供第一次摸matlab的人琢磨。matlab之美在于其代码的简洁与其在各门各业中的使用。累和乘积与上下三角clcclearA=magic(9)cumsum(A)prod(A,2)%按行求乘积triu(A)%上三角 tril(ones(4,4),-1)A=[2 1 -1;2 1 0;1 -1 1]B=[1 -1 3;4 3 2]X=B/AX*AA*B'因式分原创 2017-01-25 23:05:11 · 1271 阅读 · 0 评论 -
粒子群优化算法与遗传脚本
粒子群算法脚本:PSO.m%% PSO%% 清空环境clcclearclose allwarning off%% 参数初始化%粒子群算法中的两个参数c1 = 1.49445;c2 = 1.49445;maxgen=200; % 进化次数 sizepop=200; %种群规模%粒子更新速度Vmax=1;Vmin=-1;%种群popmax=50;popmin=-原创 2017-01-25 22:18:36 · 1529 阅读 · 7 评论 -
数学建模模拟退火法MATLAB程序参考模板
数学建模模拟退火法MATLAB程序参考模板想说的其他话都备注在程序中了,这里不多说。主程序% 模 拟 退 火 算 法 ( Simulated Annealing Algorithm ) MATLAB 程 序%模拟火算法(MATLAB 实现)clear ;% 程 序 参 数 设 定Coord = ... % 城 市 的 坐 标 Coordinates[ 0.6683 0.6195 0.4原创 2017-01-25 21:04:23 · 4761 阅读 · 0 评论 -
DIJKSTRA程序算法示例
DIJKSTRA程序算法示例 Dijkstra算法是典型最短路算法,用于计算一个节点到其他所有节点的最短路径。主要特点是以起始点为中心向外层层扩展,直到扩展到终点为止。Dijkstra算法能得出最短路径的最优解,但由于它遍历计算的节点很多,所以效率低。 Dijkstra算法是很有代表性的最短路算法,在很多专业课程中都作为基本内容有详细的介绍,如数据结构,图论,运筹学等等。Dijkstra原创 2017-01-22 13:59:50 · 2128 阅读 · 1 评论