- 博客(371)
- 资源 (3)
- 收藏
- 关注
原创 一点浩然气,千里快哉风
最后几个月,感觉效率特别低,可能心理上有些倦怠的缘故,英国整体上也是一个效率低的国家,这边的导师和合作者做事也普遍比较慢。每次一想到国内那个比较内卷的环境,就想把论文赶快写完投出去,赶快把数量刷上去,然而现实的英国,什么都很慢,心理难免有些焦虑。商科领域的顶级期刊基本全是美国的,新加坡或香港大部分教师毕业于美国常青藤,在文章发表上普遍强于英国这边的教师。英国这边的大学老师是铁饭碗,工资收入也不高,在科研上没有太大压力或动力。总体上,来英国的访问学者很少是国内名牌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关注的人