意大利数学家Z.高津托
意大利伟大数学家Sire Zepartzatt Gozinto的生卒年代是一个谜[1],但是他发明的 “高筋图” 在 制造资源管理、物料清单(BOM)管理、智能阅读、科学文献影响因子计算 等方面具有重要应用。
高津托图
下图是一个制造业物料需求高津托图,节点FP1、FP2分别表示最终产品的需求量,边上的数值表示组装部件所需要的上游零部件的数量,物料清单(BOM)系统需要知道所有零部件的总需求。图中:
Primary Demand(主需求) -- 市场对零部件的需求数量
Secondary Demand(次需求) -- 因产品组装产生的对零部件的需求
Total Demand(总需求)-- 以上两个需求之和
Product No. (产品(拓扑次序)编号)-- 根据组装约束对零部件产品进行拓扑排序的次序数
数学模型
设图中的零部件类型数为n,装配关系(边)数为m
设pd[i]为节点i的主需求(常量)
sd[i]为节点i的次需求(决策变量)
td[i]为节点i的总需求(被动变量)
pd[i]为节点i的产品拓扑次序编号(决策变量)
根据装配逻辑,对任何边k,如果边k的起始节点为a[k],终止节点为b[k],权值为c[k],则:
sd[i]=sum{k=1,...,m;a[k]==i}(c[k]td[b[k]]) | i=1,...,n
td[i]=sd[i]+pd[i]|i=1,...,n
把零部件从装配上游到下游排序:
pn[b[k]] >= pn[a[k]] + 1 | k=1,...,m
pn[i]>=1|i=1,...,n
pn[i]<=n|i=1,...,n
+Leapms模型:
min sum{i=1,...,n}pn[i]
subject to sd[i]=sum{k=1,...,m;a[k]==i}(c[k]td[b[k]]) | i=1,...,n
td[i]=sd[i]+pd[i]|i=1,...,n pn[b[k]] >= pn[a[k]] + 1 | k=1,...,m
pn[i]>=1|i=1,...,n
pn[i]<=n|i=1,...,n where
m,n are numbers
e,pd are sets
a[k],b[k],c[k] are numbers | k=1,...,m
sd[i],td[i] are variables of nonnegative numbers|i=1,...,n
pn[i] is a variable of nonnegative number|i=1,...,n data_relation
m=_$(e)/3
n=_$(pd)
a[k]=e[3k-2]|k=1,...,m
b[k]=e[3k-1]|k=1,...,m
c[k]=e[3k] |k=1,...,m
data
pd={150 50 20 230 0 0 0 0}
e={
3 1 1
4 1 2
4 2 3
4 3 3
4 5 2
5 2 4
6 3 4
6 4 5
7 4 3
7 5 1
8 5 2
}
求解:
+Leapms>load
Current directory is "ROOT".
.........
gozinto.leap
.........
please input the filename:gozinto
================================================================
1: min sum{i=1,...,n}pn[i]
2: subject to
3:
4: sd[i]=sum{k=1,...,m;a[k]==i}(c[k]td[b[k]]) | i=1,...,n
5: td[i]=sd[i]+pd[i]|i=1,...,n
6:
7: pn[b[k]] >= pn[a[k]] + 1 | k=1,...,m
8: pn[i]>=1|i=1,...,n
9: pn[i]<=n|i=1,...,n
10:
11: where
12: m,n are numbers
13: e,pd are sets
14: a[k],b[k],c[k] are numbers | k=1,...,m
15: sd[i],td[i] are variables of nonnegative numbers|i=1,...,n
16: pn[i] is a variable of nonnegative number|i=1,...,n
17:
18: data_relation
19: m=_$(e)/3
20: n=_$(pd)
21: a[k]=e[3k-2]|k=1,...,m
22: b[k]=e[3k-1]|k=1,...,m
23: c[k]=e[3k] |k=1,...,m
24: data
25: pd={150 50 20 230 0 0 0 0}
26: e={
27: 3 1 1
28: 4 1 2
29: 4 2 3
30: 4 3 3
31: 4 5 2
32: 5 2 4
33: 6 3 4
34: 6 4 5
35: 7 4 3
36: 7 5 1
37: 8 5 2
38: }
================================================================
>>end of the file.
Parsing model:
1D
2R
3V
4O
5C
6S
7End.
..................................
number of variables=24
number of constraints=43
..................................
+Leapms>solve
The LP is solved to optimal.
找到线性规划最优解.非零变量值和最优目标值如下:
.........
pn1*=4
pn2*=4
pn3*=3
pn4*=2
pn5*=3
pn6*=1
pn7*=1
pn8*=1
sd3*=150
sd4*=1360
sd5*=200
sd6*=8630
sd7*=4970
sd8*=400
td1*=150
td2*=50
td3*=170
td4*=1590
td5*=200
td6*=8630
td7*=4970
td8*=400
.........
Objective*=19
.........
+Leapms>
结果
参考文献
[1] Rousseau, R. . (1987). The gozinto theorem: using citations to determine influences on a scientific publication. Scientometrics, 11(3-4), 217-229.