2016年11月_心态与习惯

原创 latex 引用公式

latex 运用引用功能，给公式添加标签 \label，就不用因为公式的编号变动而麻烦地更改大量数字了。1. 一般情况下，可以用 \ref 命令例如，下面一个公式：\ begin{equation}I_{t}=I_{0}+\sum_{i=1}^{t}(Q_{i}-d_{i})\label{con:inventoryflow}\ end{equation}Equation

2016-11-24 03:16:04 168385 8

原创根据一组数据拟合出最佳分布函数, arena, Python, Matlab, R语言

不论是 matlab 还是 R 语言，根据一组数据拟合分布函数时，只能拟合出特定分布的参数值，并不能给出最适合这组数据的分布函数。例如，matlab 有 normfit 函数来计算一组数据拟合成正态分布时的均值和方差，wblfit 函数来给出拟合成威布尔分布的参数值，还有 gamfit, binofit 函数等。R 语言中也有 fitdistr函数。它们都不能自动给出最佳拟合分布，个人觉得可...

2016-11-22 02:38:53 32803 23

原创收敛速度：线性收敛，超线性收敛，r 阶收敛

最优化理论中，评价一个算法的收敛速度有两个衡量尺度，Q-收敛与 R-收敛，我们一般用到的是 Q-收敛，它包括：线性收敛，超线性收敛，r 阶收敛。设相邻两个迭代点：x_(k+1), x_{k}，最优值点 x*，若存在实数 q>0，满足：1. 若 0线性收敛2. 若 q=0，则表示算法超线性收敛若并且 r>1, q>=0，则称算法 r 阶收敛。

2016-11-21 20:28:30 44688 16

原创最优化理论中的方法概述

最优化方法首先可以分为：随机搜索方法与确定型优化方法。1. 随机搜索方法：元启发式方法，例如遗传算法，人工神经网络，模拟退火算法等。2. 确定型方法：模式搜索方法，梯度搜索方法。2.1 模式搜索方法：根据函数值的变化规律寻找最优值，例如单纯形法2.2 梯度搜索方法：根据函数的梯度求解，可以分为线搜索方法与信赖域搜索方法2.2.1 线搜索方法：利用黑塞矩阵，梯度，

2016-11-21 19:26:32 1384

原创 Matlab代码，一维搜索用进退法确定搜索区间

在使用二维优化算法时，经常需要调用一维搜索算法求步长，而一维搜索算法需要确定搜索区间。进退法不仅可以直接当作一维搜素算法使用，也可以用来确定一维搜索区间。进退法确定搜索区间的原理是：假设一维目标函数为一个单峰函数，初始点 x0, 步长 h, 迭代点 x1= x0+h。 lamda初始值设为1.(1) 若 f(x0)> f(x1)，则 x2=x1+lamda*h. 找到一个 lamda

2016-11-21 03:22:53 18415 12

原创 static dynamic 情形下的（R, T）策略及 Matlab 代码

基于论文： Tarim S A, Kingsman B G. The stochastic dynamic production/inventory lot-sizing problem with service-level constraints[J]. International Journal of Production Economics, 2004, 88(1): 105-119.

2016-11-19 04:19:32 749

原创 static dynamic策略下的随机批量生产问题

基于论文：Bookbinder J H, Tan J Y. Strategies for the probabilistic lot-sizing problem with service-level constraints[J]. Management Science, 1988, 34(9): 1096-1108.随机批量生产问题中："static dynamic" 策略是指生产周

2016-11-19 03:28:17 776

原创最优化中的牛顿法，二阶收敛性

最近偶尔翻阅一本写的不错的最优化理论教材，该书讲得很详细，很透彻。我对非线性规划理论又有了全新的认识，发现牛顿法可以说是是无约束优化中最重要的方法，其他方法：LM方法，高斯牛顿，拟牛顿法，共轭梯度法可以说是对牛顿法的扩展。准备闲来无事时将牛顿法的原理以及求解过程用数例好好再过一遍。适用对象：二阶可微函数牛顿法的几何意义本质：在原函数的某一点处用一个二次函数近似

2016-11-17 02:05:35 33635 4

原创 static uncertainty策略下的随机批量生产问题（stochastic lot sizing）的建模与求解,Matlab

基于论文：Bookbinder J H, Tan J Y. Strategies for the probabilistic lot-sizing problem with service-level constraints[J]. Management Science, 1988, 34(9): 1096-1108.对于单产品批量生产问题，问题模型为：文中提出了三种订货策略。

2016-11-15 19:39:11 1469

原创证明牛顿法在极小点附近收敛

设二阶连续可微函数 f 的一个局部极小点为 x*，此处的一阶倒数我 g(x*)，二阶导数为 G(x*)，并且G(x*)非奇异。若牛顿法的初始点充分靠近 x*，求证：\begin{equation}\lim\limits_{k\rightarrow+\infty}x^{k}=x^{*}\end{equation}证明：因为 x* 为局部极小点，所以\begin{equation...

2016-11-15 07:41:28 3726

原创最速下降法的迭代方向为什么相互垂直

证明：最速下降法的迭代公式：\begin{equation}x_{k+1}=x_{k}-\lambda\Delta f(x)\end{equation}其中，迭代步长 lambda 由一维搜索得到，lambda 满足\begin{equation}\lambda=\arg~\max\limits_{\lambda} f(x_{k}-\lambda\Delta f(x))\end{equation}...

2016-11-15 02:35:12 6241 4

原创为何最优化中的牛顿法是椭球范数下的最速下降法

梯度与黑塞矩阵分别由下列符号表示：\begin{equation}g(x)=\nabla f(x_k)\end{equation}\begin{equation}g(x)=\nabla ^2 f(x_k)\end{equation}牛顿法的迭代关系式为：\begin{equation}x_{k+1}=x_{k}-\lambda\nabla ^{2}f(x_{k})^{-1}\nab...

2016-11-14 04:20:29 3659 9

原创 matlab心形图大全，几个经典的函数图像，有趣的函数图像，matlab画图

收藏几个经典的函数，用来当作数学优化方法求解的例子。再收藏几个有趣的函数图像。1 平面心形图像经过网上搜索，平面坐标系中，下面函数最想心形图像，而且函数形式简单，不需要分段。改变右面的数值可以变化心形图像的大小。> ezplot('x^2+(y-(x^2)^(1/3))^2=9')

2016-11-11 21:25:35 160556 7

原创初始库存不为0 时的 Wager Whitin 算法 Matlab 代码

最近读几篇老论文，发现Wager-Whitin 算法在初始库存不为 0 时仍然成立。自己没有想到如何证明。原 wagner whitin的论文上又看了看，似乎证明时没有考虑到初始库存不我0时定理 1 不成立的情况。但经过编程测试，用了五六百个算例，发现仍然成立。此时，最优解有如下特点：第一个生产周期刚好将初始库存与生产量消化完，之后的阶段仍然满足“”零库存“”性质。Matla

2016-11-09 19:38:03 1802 1

原创 matlab中的 peaks 函数的含义

matlab 中的 peaks 函数是一个典型的多元函数。它本质上是一个二元高斯分布的概率密度函数，函数表达式为：\begin{equation}f(x,y)=3(1-x)^{2}e^{-x^2-(y+1)^2}-10(\frac{1}{5}x-x^3-y^5)e^{-x^2-y^2}-\frac{1}{3}e^{-(x+1)^2-y^2}\end{equation}在 matlab 命令

2016-11-07 20:30:01 86353 1

原创 matlab 中画箭头

matlab 画图时，没有直接画箭头的函数。网上搜索解决方案时，有人提到用 annotation 函数，但该函数的位置坐标并不是传统的坐标值，而是在整个图形位置的单位坐标，不好调整；另一种方法用 fill 函数填充一个三角形区域表示箭头，代码比较多并且繁琐。解决方案：有外国大神直接编写了一个 arrow 函数，坐标值是真实坐标，可以像使用 plot 画直线那样使用该函数画箭头。文件

2016-11-07 01:16:08 23149 5

原创期望库存量，期望缺货量的线性拟合，基本思想与 matlab 代码

基于论文：Rossi R, Tarim S A, Prestwich S, et al. Piecewise linear lower and upper bounds for the standard normal first order loss function[J]. Applied Mathematics and Computation, 2014, 231: 489-502.在

2016-11-07 00:56:51 3219

原创最优化中的鞍点

在学习最优化课程时，不时听到“鞍点”这个名词。老师很快提了这个词，但没有详细介绍鞍点的含义。鞍点（saddle point）的数学含义是：目标函数在此点上的梯度（一阶导数）值为 0，但从该点出发的一个方向是函数的极大值点，而在另一个方向是函数的极小值点。判断鞍点的一个充分条件是：函数在一阶导数为零处（驻点）的黑塞矩阵为不定矩阵。半正定矩阵：所有特征值为非负，或主子式...

2016-11-06 23:13:03 21682 4

原创遗传算法求解旅行商 TSP 问题，matlab代码

学习启发式算法时，旅行商问题是一个经典的例子。其中，遗传算法可以用来求解该问题。遗传算法是一种进化算法，由于其启发式算法的属性，并不能保证得到最优解。求解效果与初始种群选取，编码方法，选择方法，交叉变异规则有关。上课时，老师不知从哪里找了一个非常粗糙的程序，自己将不少错误修正，增加了一些注释方便理解，并增加了一些代码使程序更完美。该代码可以动态显示每一代的路线，非常直观！（备注：对于一般的 ...

2016-11-02 01:24:11 76436 190

天天向上的专栏