- 博客(371)
- 资源 (3)
- 收藏
- 关注
RSS订阅原创 一点浩然气，千里快哉风
最后几个月，感觉效率特别低，可能心理上有些倦怠的缘故，英国整体上也是一个效率低的国家，这边的导师和合作者做事也普遍比较慢。每次一想到国内那个比较内卷的环境，就想把论文赶快写完投出去，赶快把数量刷上去，然而现实的英国，什么都很慢，心理难免有些焦虑。商科领域的顶级期刊基本全是美国的，新加坡或香港大部分教师毕业于美国常青藤，在文章发表上普遍强于英国这边的教师。英国这边的大学老师是铁饭碗，工资收入也不高，在科研上没有太大压力或动力。总体上，来英国的访问学者很少是国内名牌985高校的教师，普通211和一本的居多。
2023-05-23 03:31:17
16
原创 苹果电脑 Macbook 使用 cplex
如果是 intel 芯片，选择后缀是 osx 但没有 atm 的文件。如果是 M 系列芯片，要选择后缀带 osx-atm 的文件。
2023-04-22 19:28:56
78
原创 幂函数与指数函数的近似
幂函数1xα可以近似为指数函数eαx，甚至可以进一步近似为1αx。在一本书中指数平滑方法的介绍中见到了这个近似，总结一下。
2023-04-04 18:30:35
521
原创 Introduction to linear optimization exercise 第二章课后题答案
Introduction to linear optimization》这本书很不错，但是课后题没答案，准备通过网络搜集以及自己的练习，整理一下。
2023-03-26 20:36:44
206
原创 拟凸函数，拟凹函数，单峰函数
为单变量函数时，单峰函数（只有一个局部极小值或局部极大值的函数）要么是拟凸函数，要么是拟凹函数，而多变量函数时则不一定。拟凸（quasi-convex）函数很早就听说过，但是标准定义一直不太了解，现在总结一下。几何意义是函数任意两点连线上的点，在该函数上的值小于这两点对应函数值的最大值。是拟凸函数：若对于其定义域内的任意两个点。上图就是一个拟凸函数。一个定义在凸集上的实数函数。
2023-02-24 05:26:37
312
原创 全期望值定理与全方差定理
全期望值定理（law of total expectation）比较熟悉，竟然还有个全方差定理（law of total variance），关于条件期望与条件方差的，总结一下。
2023-02-17 07:26:06
477
原创 学术生涯中曾经的小目标
utd 期刊上的一些文章，有的没啥理论难度，我看也是在灌水，但大部分还是挺有深度的，确实比欧洲的一些期刊文章质量高些。工作后，想着反正没有非升即走压力，就把论文慢慢投，EJOR, Omega, IJPE, IJPR, COR, JORS, AOR, CIE 这些期刊一个一个投，拒了就换下一个。刚开始准备读博时，特别想发表一篇 sci。后来忍了一年多，实在坚持不下去了，因为 sci 迟迟发不出来，等到博士快毕业时才发出来。或许未来能够突破下，发一个弱的 utd，但是 utd 中的强顶刊，现在还不敢奢望。
2023-02-08 01:12:44
258
原创 《Introduction to linear optimization》 看了一遍了
从 2019 年 4 月份买到这本书，到把这本书完全看一遍，都过了 3 年多了。准备顺着第一遍的笔记慢慢看第二遍，再做做课后题。
2023-02-01 00:38:08
304
原创 一个变量在一个范围内取值，线性化
偶尔碰到一个线性化的问题：一个变量在一个区域内为一个值，在其他区域内为另一个值，即。引入一个 0-1 变量。
2022-12-27 03:53:07
522
9
原创 关于矩阵求导的几个性质
在高级的理论教材中，经常能见到针对一组变量的求导，往往涉及到矩阵。并不能快速手动推导，需要熟悉记住一些结论，于是决定写篇博客总计一下。这些性质都很容易通过对矩阵展开运算验证。对于 m×1m\times 1m×1 列向量 y\mathbf{y}y， m×nm\times nm×n 矩阵 A\mathbf{A}A， n×1n\times 1n×1 列向量 xxx，若它们的关系如下：y=Ax\mathbf{y=Ax}y=Ax则∂y∂x=A\frac{\partial \mathbf{y}}{\parti
2022-12-07 20:42:58
285
原创 凸包(convex hull)，凸包络面(convex envelope), 凸低估计量(convex underestimator), 图上方(epigraph)，
凸分析中经常见到这些概念，目前这方面的中文资料似乎不太多，决定写篇博客总结一下。凸包在文献中比较常见些，也称作凸包络面 convex envelope。凸包一般针对某个集合（函数也可以有凸包，但我看到一些文献将函数的凸包称作凸包络面 convex envelope）。凸包的定义为：对于某个有限集合 {v1,v2,…,vn}\{v_1, v_2, \dots, v_n\}{v1,v2,…,vn}，它的凸包为conv{v1,v2,…,vn}={θ1v1+θ2v2+⋯+θnvn∣θi≥0,∑nθi=1,∀i
2022-11-17 00:47:51
793
2
原创 关于组合数（二项系数）的一个递推公式
对于组合数，一些英文资料里也称为二项系数（二项分布的系数），有一个递推公式。这个非常容易证明，将表达式展开为阶乘计算即可。最近听到一位教授提起，发下自己忘了，于是写篇博客记录下。
2022-11-16 06:03:07
402
原创 将一个整数划分为几部分，考虑顺序，java 代码
最近又无意中遇到了一个整数划分问题，将一个整数划分为几个数的和，考虑顺序。网上搜到一些方案是 python 的，没有给出考虑顺序的结果，自己适当修改了下，用 java 实现了。
2022-10-25 05:27:47
1161
原创 1024 程序员节
今年在英国访学，最近感觉做科研好难啊，困难大于预期。但困难总可以一步一步简化，然后慢慢解决的。突然意识到今天是10月24日，程序员节。我比较喜欢编程，每年这一天都想写点什么。只要自己一直在奋进，相信未来一定会更好！
2022-10-24 23:11:59
128
原创 一个凸函数关于一个随机变量的期望仍然是凸函数
在证明 base stock 策略时用到了这个性质，设有一个关于 $x$ 的凸函数 $f(x)$，另有一个随机变量 $\xi$，该随机变量的概率密度函数为 $\phi(\xi)$。则 $f(x)$ 对这个随机变量一个非负函数的期望得到的另外一个函数，即 $g(x)=\int f(x-\xi)\phi(\xi)d\xi$，也为凸函数。
2022-10-03 23:31:27
570
原创 到底什么是拓扑空间，拓扑
拓扑可以理解为一个给定集合内元素所组成的，另外一个满足上面三个条件的集合。拓扑空间有时候也直接称为拓扑，一个集合。，下面的六个新集合中，最后两个并不属于。下面，用维基百科的一个图形例子来说明。
2022-09-08 00:48:08
692
3
原创 推土距离， Wasserstein distance
在鲁棒优化中见到一个名词：Wasserstein distance, 又称作 earth mover's distance，堆土距离。用来测量两个随机分布的距离。之所以叫堆土距离，是因为它的物理含义是：将一个随机分布的直方图，变化成另一个随机分布的直方图，直方图变动的最小距离和。直方图中的每个小柱子可以视作一个土堆。...
2022-05-26 18:50:54
855
1
原创 克隆巴赫系数 Cronbach‘s alpha 及 R, Python 实现
克隆巴赫系数 Cronbach’s alpha 是度量问卷信度（Reliability）的一个指标，一般大于 0.7 表示问卷的信度可以接受，大于 0.8 表示问卷的信度良好，大于 0.9 表示问卷的信度优秀。信度即测量的一致性，表示若问卷中的问题重复很多次问同一个人，这个人的回答总是一致的。效度是衡量问卷中的问题能否测量出它应该测量的事物的成都，例如网上的一些 IQ 测试问题真的能反映一个人的真实智力吗？这个就属于效度要衡量的东西。信度的标准定义为真值方差与测量值的方差之比。测量值 xxx 与真值 τ
2022-04-21 14:48:33
8935
原创 eclipse 中的 git unauthorized 2022
eclipse 不能与 github 连接，最后决定放弃用 eclipse 进行 push 或 pull ，用 github 自己的 github desktop 软件进行 git 操作比较好
2022-04-18 21:30:54
403
4
原创 随机对偶动态规划 SDDP，报童模型的一个 python 例子
将 benders 分解与 SAA 方法结合在一起，就产生了随机对偶动态规划方法。这两年在论文里看到了这个方法，一般用在运输或能源规划问题中，在库存管理问题中还没见别人用过。
2022-03-31 21:42:34
3008
4
原创 读了下神书 Introduction to Linear Optimization
这本 MIT 教授 Dimitris Bertsimas 写的神书 《Introduction to Linear Optimization》，早在 2019 年我一个博士同学就推荐过给我。我当时从淘宝买了影印版，读了几页，之后就把这本书遗忘了。后来在 2021 年看到一位 MIT 博士在公众号中推荐了这本运筹学著作，并评这本书为三大神书之一。于是又拿起了这本书，从头开始细看。的确如评价所说，“抽吸剥茧”，将线性规划中的很多原理讲解地明明白白，后悔没有早点读完。这本书比较厚，有 500 多页，目前我才看
2022-02-28 18:42:03
904
原创 两阶段与多阶段随机规划
这段时间发现一本不错的随机规划教材：《Handbooks in Operations Research and Management Science: Stochastic Programming》。这本书比之前看的一些教材更通俗易懂一些。总结一下这几天看到的知识点。1. 两阶段规划two-stage model 中，xxx 为第一阶段的决策变量，必须在不确定性发生之前作出决定，yyy 为第二阶段的决策变量，在不确定性发生之后作出决定。ξ\xiξ 为随机变量，而 www 为随机变量的一个具体表现值。
2022-01-04 23:08:20
2585
原创 总结2021
2021年很快就过去了，这一年申请到了国家自科，也录用了一篇领域内还可以的论文，生活也算有点安定。博士毕业工作已 3 年，目前的学术成果中规中矩，仍然有些遗憾与不足。目前单位的科研环境不好，但最大好处是压力不大。在新的一年里，希望自己不忘初心，继续在学术道路上探索，工作与生活都更健康，更有规律一些。...
2022-01-01 18:40:59
881
2
原创 关于条件方差的一个性质
今天看到关于条件方差的一条性质，记录下来：Var(Y)=E(Var(Y∣X))+Var(E(Y∣X))\text{Var} (Y) = \text{E}(\text{Var}(Y|X))+\text{Var}(\text{E}(Y|X)) Var(Y)=E(Var(Y∣X))+Var(E(Y∣X))而条件均值的性质为：E(Y)=E(E(Y∣X))\text{E}(Y)=\text{E}(\text{E}(Y|X))E(Y)=E(E(Y∣X))...
2021-12-29 12:21:31
2695
1
原创 python Numpy 的一维数组 1D array 既是行向量也是列向量，矩阵相乘符号@
我发现 Numpy 中的一维数组既可以是行向量也可以是列向量：import numpy as npa = np.array([1, 2, 3, 4])In [18]: a.shape(4,)IPdb [19]: a.T.shape(4,)可以发现，一维数组的转置就是它自身。若在定义数组时，里面有两个中括号，就是一个行向量了。b = np.array([[1, 2, 3 ,4]])In [24]: b.shape(1, 4)但若 numpy 的二维数组与一维数组做矩阵乘法，一维
2021-12-26 20:09:19
3156
原创 统计模型评价准则 AIC
统计模型评价时，经常见到一个准则 AIC， 全称是 Akaike information criterion，是以日本的统计学家 Akaike 命名的，它的计算公式为：AIC=2k−ln(L^)AIC=2k-\ln(\hat{L})AIC=2k−ln(L^)其中，kkk 是模型中的待估参数数量， L^\hat{L}L^ 是该模型极大似然估计的最大值。AIC 值越小，说明该统计模型损失的信息越少。...
2021-12-03 17:39:03
2131
原创 详解 Benders 分解与一个算例的 python 代码
听说过 benders 分解几年了，一直没有细看，最近论文里面也见到，还是有必要了解一下它的基本思想与用法的。
2021-11-15 01:19:35
4039
19
原创 max-min 不等式(minimax inequality)
在一篇鲁棒优化相关的论文中，看到了 minimax inequality，又称作 max-min 不等式。查了一下，原来在拉格朗日对偶中应用过，写篇博客总结一下。
2021-11-09 13:25:33
1602
原创 录用了一篇 IJPR
昨天录用了一篇 IJPR，若是博士毕业之前能发 IJPR，不知有多高兴啊。这篇文章从投到录将近一年：一审两个半月，大修，状态是：Major Revision but allow Resubmission二审一个月多点，小修，让修改一些语句三审花了四个多月，又是小修，原来这个小修又发给原审稿人了，而审稿人似乎暑假度假了，特别慢，又修改了几句话三审一个月，小修，领域编辑让补充了一些数据说明之类的话四审/五审，半个多月，小修/录用，主编那里让引用几篇期刊的新论文，领域编辑让稍微扩充一下论文结论读
2021-11-05 20:46:41
2083
12
原创 内卷的国内学术圈 1024
最近浏览了一下国内一些高校相近研究领域的青年教师，发现大家的科研成果水平越来越高了，不禁感叹学术圈越来越内卷。一些普通211高校的新入职教师，可能要有 utd 顶刊才能留下来一些普通一本学校的新入职教师，科研成果或许比10年前一些985高校教授的成果都好觉得自己论文写的太慢，但也不想那么急功近利，可是还得有好成果啊。读博期间要是能发一篇领域内中上等的期刊，不知道要多高兴。但是现在，似乎没有那么振奋人心（仅仅才过了三年），可能有 utd 才会特别高兴些。。。搞科研，最主要还得靠自己，最起码一个科
2021-10-27 23:47:01
487
2
原创 指示约束条件 (indicator constraint) 与大M (big M) 条件
数学规划软件 cplex 或 gurobi 都有专门的指示约束条件：约束一个变量与另外一些变量的关系，例如对于一个 0-1 变量 zzz，它与 xxx 有以下关系：x>0⇒z=1x=0⇒z=0\begin{aligned}x>0 \Rightarrow z=1\\x=0\Rightarrow z=0\end{aligned}x>0⇒z=1x=0⇒z=0一般情况下可以使用大 M 构造成线性规划模型：x≤Mzx≥0\begin{aligned}x &\leq Mz
2021-09-08 21:42:18
3127
2
原创 无限阶段平稳马尔科夫过程的求解算法：策略迭代与值迭代
无限阶段的马尔科夫决策过程，相关论文一般比较难，初步总结一下知识点。若无限阶段的马尔科夫过程为平稳的，那么各阶段的状态转移概率与回报函数都是相同的，这样**马尔科夫过程最终会达到一个平稳的状态**。（在库存问题中，平稳过程意味着各阶段的随机需求函数相同，收益与成本参数也相同）
2021-08-27 14:29:14
699
原创 大致看了下伍德里奇的《计量经济学导论》
我的研究方向也有不少做实证研究的，而计量经济学是实证研究的重要武器，因此，在 8 月份时看了这本经典的计量经济学教材。的确比不少国内教材好，优点是比较详细，很多知识点讲解地很清楚看到第10章之后，略微有点吃力了。因为看的速度太快，没有完全消化，很多知识点需要再看一遍这本书有很多例子和课后题，配合一些计量软件练习，应该能学地更好，stata 是做计量最好的软件，有必要学一学...
2021-08-22 16:39:08
601
Cplex, Ilog, OPL 使用手册，manual
2017-11-14
Cplex OPL 使用手册，OPL 编程
2017-11-14
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人