发明名称 一种以机器特征属性为区间数的作业车间瓶颈识别方法
摘要 本发明提出了一种以机器特征属性为区间数的作业车间瓶颈识别方法,首先采用优化方法获得调度问题的调度优化方案,其次建立区间属性值的多属性瓶颈识别模型,将瓶颈识别问题转化为根据多个评价属性综合评价候选机器的多属性决策问题,然后采用改进的三参数区间TOPSIS进行瓶颈的最终判定。本发明采用改进的三参数区间TOPSIS瓶颈识别方法增强了区间的信息表示能力,提高了区间的信息精度,改进的三参数相离度公式解决了属性区间分布概率不相同的权重确定问题,能有效地确定机器属性的权重;区间多属性瓶颈识别方法的贴近度C<sub>i</sub>充分挖掘设备属性蕴涵的信息,相比确定性多属性瓶颈识别方法,具有更高的瓶颈识别质量。
申请公布号 CN102789598B 申请公布日期 2015.07.01
申请号 CN201210232167.1 申请日期 2012.07.06
申请人 西北工业大学 发明人 王军强;陈剑;王烁;郭银州;崔福东;张映锋;孙树栋
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 西北工业大学专利中心 61204 代理人 陈星
主权项 一种以机器特征属性为区间数的作业车间瓶颈识别方法,其特征在于:包括以下步骤:步骤1:建立生产系统及加工任务的仿真模型,并在仿真模型中设置随机扰动;采用优化算法优化生产投料次序,得到G个优化调度方案{S<sub>1</sub>,S<sub>2</sub>,…,S<sub>G</sub>};步骤2:建立所有m台加工机器的机器集M={M<sub>1</sub>,M<sub>2</sub>,...,M<sub>m</sub>},建立机器的n个评价属性的评价属性集X={X<sub>1</sub>,X<sub>2</sub>,...,X<sub>n</sub>};步骤3:根据步骤1中的G个优化调度方案,计算每个优化调度方案下每台机器的评价属性值,得到每台机器下每个评价属性的属性值区间,其中第i台机器的第j个评价属性的属性值区间为<img file="FDA00001857433600011.GIF" wi="281" he="64" /><img file="FDA00001857433600012.GIF" wi="46" he="64" />和<img file="FDA00001857433600013.GIF" wi="51" he="64" />分别表示根据G个优化调度方案计算得到的第i台机器第j个评价属性中G个属性值的最小值和最大值;建立用于瓶颈识别的区间多属性决策矩阵<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>D</mi><mo>~</mo></mover><mo>=</mo><msub><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>~</mo></mover><mi>ij</mi></msub><mo>)</mo></mrow><mrow><mi>m</mi><mo>&times;</mo><mi>n</mi></mrow></msub><mo>=</mo><msub><mrow><mo>{</mo><mo>[</mo><msubsup><mi>x</mi><mi>ij</mi><mi>L</mi></msubsup><mo>,</mo><msubsup><mi>x</mi><mi>ij</mi><mi>U</mi></msubsup><mo>]</mo><mo>}</mo></mrow><mrow><mi>m</mi><mo>&times;</mo><mi>n</mi></mrow></msub><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>m</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA00001857433600014.GIF" wi="1137" he="67" /></maths>步骤4:采用改进的三参数区间TOPSIS进行瓶颈的最终判定:步骤4.1:采用向量标准化方法将区间多属性决策矩阵<img file="FDA00001857433600015.GIF" wi="531" he="68" />转化为规范化决策矩阵<img file="FDA00001857433600016.GIF" wi="527" he="67" />其中:当评价属性x<sub>ij</sub>类型为效益型属性时:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>r</mi><mi>ij</mi><mi>L</mi></msubsup><mo>=</mo><msubsup><mi>x</mi><mi>ij</mi><mi>L</mi></msubsup><mo>/</mo><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>(</mo><msubsup><mi>x</mi><mi>ij</mi><mi>U</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mtd></mtr><mtr><mtd><msubsup><mi>r</mi><mi>ij</mi><mi>M</mi></msubsup><mo>=</mo><msubsup><mi>x</mi><mi>ij</mi><mi>M</mi></msubsup><mo>/</mo><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>(</mo><msubsup><mi>x</mi><mi>ij</mi><mi>M</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mtd></mtr><mtr><mtd><msubsup><mi>r</mi><mi>ij</mi><mi>U</mi></msubsup><mo>=</mo><msubsup><mi>x</mi><mi>ij</mi><mi>U</mi></msubsup><mo>/</mo><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>(</mo><msubsup><mi>x</mi><mi>ij</mi><mi>L</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00001857433600017.GIF" wi="453" he="473" /></maths>当评价属性x<sub>ij</sub>类型为成本型属性时:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>r</mi><mi>ij</mi><mi>L</mi></msubsup><mo>=</mo><mrow><mo>(</mo><mn>1</mn><mo>/</mo><msubsup><mi>x</mi><mi>ij</mi><mi>U</mi></msubsup><mo>)</mo></mrow><mo>/</mo><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>(</mo><mn>1</mn><mo>/</mo><msubsup><mi>x</mi><mi>ij</mi><mi>L</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mtd></mtr><mtr><mtd><msubsup><mi>r</mi><mi>ij</mi><mi>M</mi></msubsup><mo>=</mo><mrow><mo>(</mo><msubsup><mrow><mn>1</mn><mo>/</mo><mi>x</mi></mrow><mi>ij</mi><mi>M</mi></msubsup><mo>)</mo></mrow><mo>/</mo><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>(</mo><mn>1</mn><mo>/</mo><msubsup><mi>x</mi><mi>ij</mi><mi>M</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mtd></mtr><mtr><mtd><msubsup><mi>r</mi><mi>ij</mi><mi>U</mi></msubsup><mo>=</mo><mrow><mo>(</mo><msubsup><mrow><mn>1</mn><mo>/</mo><mi>x</mi></mrow><mi>ij</mi><mi>L</mi></msubsup><mo>)</mo></mrow><mo>/</mo><msqrt><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mrow><mo>(</mo><msubsup><mrow><mn>1</mn><mo>/</mo><mi>x</mi></mrow><mi>ij</mi><mi>U</mi></msubsup><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00001857433600018.GIF" wi="559" he="473" /></maths><img file="FDA00001857433600019.GIF" wi="56" he="63" />为区间<img file="FDA000018574336000110.GIF" wi="151" he="63" />的区间中心,表示区间<img file="FDA000018574336000111.GIF" wi="150" he="63" />中评价属性x<sub>ij</sub>的数学期望;步骤4.2:将评价属性权重向量W=(ω<sub>1</sub>,ω<sub>2</sub>,...,ω<sub>j</sub>,...,ω<sub>n</sub>)纳入规范化决策矩阵<img file="FDA00001857433600021.GIF" wi="209" he="67" />中,得到加权规范化决策矩阵<img file="FDA00001857433600022.GIF" wi="241" he="67" /><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mover><mi>Y</mi><mo>~</mo></mover><mo>=</mo><msub><mrow><mo>(</mo><msub><mover><mi>y</mi><mo>~</mo></mover><mi>ij</mi></msub><mo>)</mo></mrow><mrow><mi>m</mi><mo>&times;</mo><mi>n</mi></mrow></msub><mo>=</mo><msub><mrow><mo>(</mo><msub><mi>&omega;</mi><mi>j</mi></msub><msub><mover><mi>r</mi><mo>~</mo></mover><mi>ij</mi></msub><mo>)</mo></mrow><mrow><mi>m</mi><mo>&times;</mo><mi>n</mi></mrow></msub></mrow>]]></math><img file="FDA00001857433600023.GIF" wi="448" he="68" /></maths>所述评价属性权重向量满足单位化约束条件ω<sub>j</sub>≥0,j=1,2,...,n,且<img file="FDA00001857433600024.GIF" wi="219" he="120" />步骤4.3:确定评价属性的正、负理想解<img file="FDA00001857433600025.GIF" wi="188" he="51" />第j个评价属性的正理想解<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mover><mi>y</mi><mo>~</mo></mover><mi>j</mi><mo>+</mo></msubsup><mo>=</mo><mo>[</mo><msubsup><mi>y</mi><mi>j</mi><mrow><mo>+</mo><mi>L</mi></mrow></msubsup><mo>,</mo><msubsup><mi>y</mi><mi>j</mi><mrow><mo>+</mo><mi>U</mi></mrow></msubsup><mo>]</mo><mo>=</mo><mo>[</mo><munder><mi>max</mi><mi>i</mi></munder><mrow><mo>(</mo><msubsup><mi>y</mi><mi>ij</mi><mi>L</mi></msubsup><mo>)</mo></mrow><mo>,</mo><munder><mi>max</mi><mi>i</mi></munder><mrow><mo>(</mo><msubsup><mi>y</mi><mi>ij</mi><mi>U</mi></msubsup><mo>)</mo></mrow><mo>]</mo><mo>;</mo></mrow>]]></math><img file="FDA00001857433600026.GIF" wi="775" he="72" /></maths>第j个评价属性的负理想解<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mover><mi>y</mi><mo>~</mo></mover><mi>j</mi><mo>-</mo></msubsup><mo>=</mo><mo>[</mo><msubsup><mi>y</mi><mi>j</mi><mrow><mo>-</mo><mi>L</mi></mrow></msubsup><mo>,</mo><msubsup><mi>y</mi><mi>j</mi><mrow><mo>-</mo><mi>U</mi></mrow></msubsup><mo>]</mo><mo>=</mo><mo>[</mo><munder><mi>min</mi><mi>i</mi></munder><mrow><mo>(</mo><msubsup><mi>y</mi><mi>ij</mi><mi>L</mi></msubsup><mo>)</mo></mrow><mo>,</mo><munder><mi>min</mi><mi>i</mi></munder><mrow><mo>(</mo><msubsup><mi>y</mi><mi>ij</mi><mi>U</mi></msubsup><mo>)</mo></mrow><mo>]</mo><mo>;</mo></mrow>]]></math><img file="FDA00001857433600027.GIF" wi="758" he="72" /></maths>区间中心的正理想解<img file="FDA00001857433600028.GIF" wi="347" he="73" />区间中心的负理想解<img file="FDA00001857433600029.GIF" wi="337" he="72" />步骤4.4:确定每台机器分别到正理想解和负理想解的距离:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msubsup><mi>D</mi><mi>i</mi><mo>+</mo></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>|</mo><mo>|</mo><msub><mover><mi>y</mi><mo>~</mo></mover><mi>ij</mi></msub><mo>-</mo><msubsup><mover><mi>y</mi><mo>~</mo></mover><mi>j</mi><mo>+</mo></msubsup><mo>|</mo><mo>|</mo><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>[</mo><mo>|</mo><msubsup><mi>y</mi><mi>ij</mi><mi>L</mi></msubsup><mo>-</mo><msubsup><mi>y</mi><mi>j</mi><mrow><mo>+</mo><mi>L</mi></mrow></msubsup><mo>|</mo><mo>+</mo><mo>|</mo><msubsup><mi>y</mi><mi>ij</mi><mi>M</mi></msubsup><mo>-</mo><msubsup><mi>y</mi><mi>j</mi><mrow><mo>+</mo><mi>M</mi></mrow></msubsup><mo>|</mo><mo>+</mo><mo>|</mo><msubsup><mi>y</mi><mi>ij</mi><mi>U</mi></msubsup><mo>-</mo><msubsup><mi>y</mi><mi>j</mi><mrow><mo>+</mo><mi>U</mi></mrow></msubsup><mo>|</mo><mo>]</mo></mrow>]]></math><img file="FDA000018574336000210.GIF" wi="1180" he="120" /></maths><maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msubsup><mi>D</mi><mi>i</mi><mo>-</mo></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>|</mo><mo>|</mo><msub><mover><mi>y</mi><mo>~</mo></mover><mi>ij</mi></msub><mo>-</mo><msubsup><mover><mi>y</mi><mo>~</mo></mover><mi>j</mi><mo>-</mo></msubsup><mo>|</mo><mo>|</mo><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>[</mo><mo>|</mo><msubsup><mi>y</mi><mi>ij</mi><mi>L</mi></msubsup><mo>-</mo><msubsup><mi>y</mi><mi>j</mi><mrow><mo>-</mo><mi>L</mi></mrow></msubsup><mo>|</mo><mo>+</mo><mo>|</mo><msubsup><mi>y</mi><mi>ij</mi><mi>M</mi></msubsup><mo>-</mo><msubsup><mi>y</mi><mi>j</mi><mrow><mo>-</mo><mi>M</mi></mrow></msubsup><mo>|</mo><mo>+</mo><mo>|</mo><msubsup><mi>y</mi><mi>ij</mi><mi>U</mi></msubsup><mo>-</mo><msubsup><mi>y</mi><mi>j</mi><mrow><mo>-</mo><mi>U</mi></mrow></msubsup><mo>|</mo><mo>]</mo></mrow>]]></math><img file="FDA000018574336000211.GIF" wi="1177" he="120" /></maths>并确定每台机器对理想解的贴近度C<sub>i</sub>:<img file="FDA000018574336000212.GIF" wi="412" he="60" />C<sub>i</sub>为区间多属性瓶颈识别评价值,C<sub>i</sub>值最大的机器为系统的瓶颈机器。
地址 710072 陕西省西安市友谊西路127号