Hybrid Discrete Teaching-learning-based Optimization Algorithm for Solving Complex Parallel Machine Scheduling Problem
-
摘要: 针对制造行业中广泛存在的一类复杂并行机调度问题, 即带到达时间、多工序、加工约束和序相关设置时间的并行机调度问题(Parallel machine scheduling problem with arrival time, multiple operations, process restraints and sequence-dependent setup times, PMSP_AMPS), 建立问题的排序模型并提出一种混合离散教与学优化算法进行求解, 优化目标为最小化最大完工时间.首先, 根据标准教与学算法(Teaching-learning-based optimization, TLBO)中两阶段个体更新公式的特点, 在保留每一阶段个体更新公式框架不变的前提下, 对公式中具体改变实数个体或向量的每个核心操作均用所设计的排列操作进行替换, 使其可直接在离散问题解空间中执行基于标准教与学算法机理的全局搜索, 从而明显提高了原算法的全局搜索效率.其次, 采用交换操作和插入操作构造了一种简洁有效地变邻域局部搜索, 对全局搜索发现的优质解区域进行细致搜索, 从而进一步增强了算法的性能.通过对不同测试问题的仿真实验和算法比较, 验证了所提算法可有效求解PMSP_AMPS.Abstract: This paper proposes a hybrid discrete teaching-learning-based optimization algorithm and set up the scheduling model for solving the parallel machine scheduling problem with arrival time, multiple operations, process restraints, and sequence-dependent setup times (PMSP_AMPS), which widely exists in the manufacturing process. The criterion is to minimize the maximum completion time. Firstly, according to the characteristics of the two stage individual updating formulas in the standard teaching-learning-based optimization (TLBO) algorithm, under the premise of keeping the framework of individual updating formulas at each stage unchanged, each core operation that specifically changes the real number or vector in the formula is replaced with proposed arrangement operation, so that it can execute the global search directly in the discrete solution space based on the TLBO mechanism, which significantly improve the TLBO's global exploration ability. Secondly, a simple but effective variable-neighborhood local search is constructed by using interchange and insert operations, which can perform a detailed search of the high-quality solution regions found by the global search, thereby further enhancing the performance of the algorithm. Simulation experiments and comparisons on different instances demonstrate that the proposed algorithm can effectively solve PMSP_AMPS.
-
Key words:
- Parallel machine scheduling problem (PMSP) /
- multiple operations /
- sequence-dependent setup times /
- arrival times /
- discrete teaching-learning-based optimization
-
随着我国综合实力的提高, 海洋在经济发展、军事安全、科学研究等领域的作用也越来越重要, 因而得到了广泛的关注[1-2]. 传感器作为获取信息的一种重要手段, 在海洋监测方面也发挥着越来越大的作用. 光纤传感器作为一种新型的无源传感器件, 具有灵敏度高、体积小、抗电磁干扰等优点[3-4], 已经在很多领域发挥作用. 近年来, 面向海洋监测应用的光纤传感器也得到了快速发展, 大量的海洋参数传感方案被提出, 目前主要涉及到海洋的温度、盐度和深度三个参数的测量[5]. 由于分布式光纤仅对温度和应变敏感, 目前光纤传感器的温度盐度深度测量以点式传感器为主.
光纤光栅(Fiber bragg grating, FBG)通常结合敏感材料实现对海洋温盐深参数的测量, 通过敏感材料实现灵敏度的放大. 中科院半导体所的Wang等[6]通过将FBG固定在毛细不锈钢管中实现了对FBG的温度增敏, 温度分辨力可以达到0.01 ℃; 通过将FBG和弹性膜片增敏结构相结合可以实现1.57 nm/MPa的压力测量灵敏度[7]; 将FBG和水凝胶相结合[8], 利用水凝胶的水分累积和扩散特性将盐度变化转换为FBG栅区的应变变化, 实现了灵敏度为9.5 nm/‰ 的盐度测量. 由于FBG的光谱形状为一个窄带宽的峰值, 光谱的大范围高分辨率解调实现较为方便, 而且可以通过波分复用技术实现多个传感器的级联. 但是基于该原理的传感器灵敏度普遍偏低.
基于该种情况, 研究者们提出了多种基于干涉原理的高灵敏海水温度盐度测量方案. 2010年, Liao等[9]通过飞秒激光刻蚀技术在光纤上形成Mach-Zehnder干涉仪, 直接将待测介质作为传感臂, 传感器的折射率灵敏度可以达到9148 nm/RIU, 可以实现nm/‰量级的盐度灵敏度. 本课题组提出了一种同一结构中两种干涉效应共存的传感方案[10], 在一个结构中就可以对温度和盐度双参数进行解耦. 之后, 为了方便传感器的布设, 降低外界拉伸对传感器的影响, 将透射式传感结构更改为反射式[11]. 这些基于干涉原理的传感器虽然具有很高的测量灵敏度, 但是传感器的复用比较困难. 基于干涉原理的光纤传感器的光谱在波长域内是准正弦分布的, 传感器的灵敏度和测量范围存在矛盾, 虽然已经通过干涉光谱解调算法解决了基于干涉原理的光纤传感器灵敏度和测量范围之间的矛盾[12-14], 但是同一系统中通过频分复用技术复用的传感器数量也是有限的. 为了解决干涉式光纤传感器的复用问题, 本文提出了一种基于调频连续波(Frequency modulated continuous wave, FMCW)技术的复用方案, 将反射端面返回的光与参考臂的光形成Mach-Zehnder干涉仪, 用于确定传感器的位置; 通过同一传感器不同反射端面间的拍频恢复传感器的光谱. 搭建了用于传感器分布测试的系统, 编写了数据处理软件用于光谱的采集与处理, 通过在系统中接入用于温度盐度测试的级联法布里−珀罗干涉仪(Fabry-Perot interferometer, FPI)探头和用于应变测试的FPI探头, 分别对应海洋环境测试中的温度、盐度、深度这三个基本参数, 并测试了实验系统中传感器的响应特性.
1. 理论介绍
图1给出了FMCW的技术原理图, FMCW技术利用可调谐激光器(Tunable laser source, TLS)发出频率随时间变化的光, 反射光因在光纤中传输表现出相对于参考光的延时特性, 通过探测器(Photoelectric detector, PD)探测到参考光和传感光形成的干涉光, 通过干涉光的频率反推得到反射端面的位置信息[15-16]. 假设入射光是调频速度为$ \gamma $ Hz/s的线性扫频光, 入射光经过耦合器分到参考臂和传感臂上, 参考臂的电场表达式可以写为:
$$ \begin{aligned} E_{\rm{ref}} = \sqrt{a}E_{0}{\rm{e}}^{{\rm{j}}(2 \pi f_{0}t+\pi \gamma t^{2})} \end{aligned} $$ (1) 其中$ f_{0} $表示扫频光的起始频率, $ E_{0} $为入射电场的幅值, $ a $表示耦合器到参考臂的分光比.
当入射光从反射面返回时, 可以表示为一个与入射光存在时间延时的扫频信号:
$$ E_{\rm{sen}} = \sqrt{(1-a)r}E_{0}{\rm{e}}^{{\rm{j}}[2\pi f_{0}(t-\tau)+\pi \gamma (t-\tau)^{2}]} $$ (2) 其中$ \tau $表示由于参考臂与传感臂臂长差造成的时延差, 具体可以计算为$ \tau = 2n\Delta l/c $. $ r $表示传感器端面的反射率. 两束光相干之后通过探测器对相干光进行探测, 探测得到的光强可以表示为:
$$ \begin{split} I=\;& (E_{\rm{ref}}+E_{\rm{sen}})(E_{\rm{ref}}+E_{\rm{sen}})^{*} = \\ &[a+r(1-a)]E_{0}^{2}\;+ \\ &2\sqrt{ar(1-a)}\cos(2\pi \gamma \tau t+2\pi f_{0}\tau-\pi\gamma\tau^{2}) \end{split} $$ (3) 可以看出, 拍频信号与参考光和传感臂间的时延存在线性对应关系, 进而可以通过该方法确定传感臂与参考臂的长度差. 图2给出了参考臂和传感臂光频率的变化情况, 两束具有时延的光形成了一个具有固定频率差的干涉信号, 这与式(3)是相同的. 当系统中的传感臂中存在多个反射面时, 会形成多个与参考臂具有不同光程差的干涉信号, 根据光程差可以确定反射面所处的位置, 通过不同反射端面与参考臂形成的干涉光谱间的拍频, 可以反推得到干涉光谱的信息, 进而可以实现单个传感器干涉光谱的还原.
图 2 参考臂和传感臂频率随时间的变化Fig. 2 Frequency changing of reference beam and sensing arm with time2. 仿真分析
本文编写了MATLAB代码对基于FMCW技术的传感器复用系统进行了仿真, 仿真中采用的系统如图3所示, 仿真中设置可调谐激光器的波长扫描范围为1530 nm ~ 1570 nm, 扫频速度大致为10 THz/s, 将99%注入到传感系统的传感臂中用于得到较强的反射光. 环形器将光注入到复用在传感臂上的传感器中并收集传感器的反射光. 用于传感器分光的耦合器分光比为95 : 5, 并在系统的4 m、5 m和6 m位置处设置三个FPI传感器, 通过平衡探测器(Balanced photoelectric detector, BPD)对相干光进行探测.
图4给出了系统的仿真光谱, 仿真光谱包含了不同位置处反射面的光谱叠加情况, 从时域光谱上很难区分不同反射面, 采用快速傅里叶变换(Fast Fourier transform, FFT)对叠加光谱进行了频谱分析, 可以得到不同反射端面的位置信息. 图5(a)给出了傅里叶分析的结果, 可以看出, 在频谱的4 m、5 m、6 m位置处出现了3个特征频率, 频谱中的1 m和2 m处的特征频率则是由于3个FPI传感器之间拍频形成的. 由于传感臂的反射光很弱, 因而传感器间拍频信号的强度也会很弱, 可以通过提高参考臂信号强度的方法进一步提高参考光与传感光拍频信号强度, 降低传感器之间光谱拍频对光谱恢复的影响. 由于每个FPI都是由相邻很近的两个反射面构成的, 为了实现对干涉光谱的恢复, 对4 m处的频谱进行了放大, 由图5(b)给出, 发现两个端面在频域是可以区分的, 因而可以实现对传感器光谱的还原.
图 5 仿真光谱的频谱特性图((a)仿真光谱频谱特性图; (b)仿真光谱频谱特性分析放大图)Fig. 5 Frequency spectrum of the simulated spectrum ((a) Frequency spectrum of the simulated spectrum; (b) Partial enlarged drawing of the frequency spectrum)之后对FPI光谱的恢复方法进行了研究, 利用矩形窗将传感器特征频谱处的复频谱信号截取出来, 补零后对其进行反傅里叶变换, 得到还原光谱, 如图6所示. 信号通过带通滤波器后会产生延时, 延时的大小与滤波器设置的参数有关, 通过将滤波后数据延时点删除可以消除滤波延时的影响. 由于恢复的传感器光谱为同一传感器两个反射面与参考臂形成的干涉光谱间的拍频, 两个信号的延时特性一致, 因而对恢复信号的影响可以忽略. 此外, 由于系统中同一传感器的延时参数是统一的, 即使带通滤波对系统响应光谱有微小影响也是可以忽略的.
图 6 还原光谱与真实光谱对比((a)周期不匹配的情况; (b)周期匹配的情况)Fig. 6 Comparison between the retrieved and real spectrum ((a) Mismatch phenomenon; (b) Match phenomenon)虽然真实光谱和还原光谱具有近似的谱形, 但是还原光谱的谷值处较为尖锐, 因为在反傅里叶变换(Inverse FFT, IFFT)后对信号进行取模运算, 导致信号没有负值部分, 这会造成信号的失真. 通过分析, 发现信号在拍频时导致了频率的减半, 具体原因可以由式(4) 给出, 可以看出两个信号在进行拍频之后会形成一个高频和低频信号的乘积, 拍频后的低频信号频率为两个信号频率差的一半. 因而可以通过倍频的方法将信号频率调整为一致, 即对信号做一个平方, 本方案中利用积化和差公式将拍频信号倍频, 使得拍频信号与真实信号频率相同.
$$ \cos \left ( a \right )+\cos \left ( b \right ) = 2\cos \left ( \frac{a+b}{2} \right) \cos \left ( \frac{a-b}{2} \right) $$ (4) 之后对系统中可以复用的传感器数量进行理论计算, 第$ { N} $个端面的反射强度为:
$$ P_{\rm{out}} = P_{\rm{in}}r_{\rm{m}}r_{ c1,N}^{2}(1-\alpha_{N})^{2}\prod\limits_{ n = 1}^{N-1} {r_{ c2,n}^{2}(1-\alpha_{n})^2} $$ (5) 其中$ P_{\rm{in}} $表示输入到参考臂中的光, $ r_{\rm{m}} $表示光纤反射端面的反射率, $ P_{\rm{out}} $是反射面反射到探测器的能量. $ r_{ c1,N} $表示第$ {N} $个耦合器第1个端口的输出能量比, $ r_{c2,n} $表示第$ {n} $个耦合器第2个端口的输出能量比, 平方表示光在耦合器中传播一个来回, $ \alpha_{N} $表示第$N $个耦合器的插入损耗. 假设本系统中采用99 : 1的耦合器将光分配给系统中的传感器, 传感器法兰间的连接没有损耗, 所有FPI传感器端面都置于盐水(折射率近似为1.33)中, 端面的反射率$ r_{\rm{m}} $大概为0.0025, 进入传感臂的光功率为10 mW, 传感器中第500个传感器的反射能量为2.2×$ 10^{-10} $ mW, 如果参考臂的输入光功率为20 μW, 那么两束光相干后的光功率约为4.2 nW. 这个光强度大于探测器的噪声等效功率, 可以通过光电探测器探测得到.
本方案中采用的方案为相干探测方案, 最大的传感长度需要综合考虑系统中光源的线宽、数据采集卡的采样率以及可调谐激光器的波长扫描速度; 根据光源线宽和相干长度之间的关系: $\Delta v = {c}/\left(2nL \right)$, 本方案中采用的可调谐激光器的线宽为60 kHz, 光源的相干长度为1.66 km. 按照本方案中设置的采样率(62.5 MHz)和波长扫描速度(80 nm/s), 为了利用参考干涉仪光谱实现等频率重采样, 参考干涉仪每个周期至少有5个采样点, 根据干涉光谱计算公式$ \Delta \lambda = \lambda^{2}/\left(2nL\right) $, 系统的最大传感距离为125 m.
3. 实验系统搭建
为了对系统的特性进行测试, 搭建了实验测试系统, 为了实现高的距离分辨力, 系统中采用的光频率扫描范围应该设置的尽可能大, 本文采用的可调谐激光器是等波长间隔扫描的, 然而干涉光谱在波长域并不是标准的正弦分布, 根据干涉仪的干涉光谱公式$ y = {\rm{cos}}(2\pi 2nl/\lambda) $, 波长位于正弦函数的分母上, 虽然波长相对于腔长较小, 得到的干涉光谱随波长是一个类正弦信号, 但是这也会导致FFT分析结果不准确, 尤其是在波长范围较大的情况下. 为了消除激光器非线性扫频的影响, 本方案中采用了附加干涉仪作为重采样的标准, 利用干涉光谱在频域是标准三角函数的特点, 三角函数在零点之间的间隔是确定的, 通过利用参考干涉仪的零点实现光谱的等频率间隔重采样.
搭建了实验系统如图7所示的带有附加干涉仪的传感系统. 采用的光源为波长范围为1480 nm ~ 1640 nm, 品牌为Santec, 型号为TSL770的可调谐激光器, 调谐速度为0 ~ 200 nm/s. 系统中使用品牌为Conquer, 型号为KG-PR-200 M的光探测器, 探测带宽为200 MHz, 用于获取参考部分的干涉信号, 所选的PD的波长探测范围为850 nm ~ 1650 nm. 通过 Thorlabs品牌的PDB570C型号的BPD将传感部分的参考臂与传感臂的干涉光进行相干探测. 本方案中采用的BPD工作波长范围可以覆盖1200 nm ~ 1700 nm, 探测带宽可以覆盖0 ~ 400 MHz.
选用Advantech公司的PCIE-1840采集卡(Data acquisition, DAQ)作为信号采集器件将探测器得到的信号读入电脑中, 该采集卡可以实现4通道16位分辨率的信号采集, 每个通道的采样率可以达到125 MHz, 在本实验中将采集卡的采样率设置为62.5 MHz, 实验中可调谐激光器的波长范围为1530 nm ~ 1570 nm, 波长调谐速度为80 nm/s, 数据每次采样时间为0.5 s. 需要用到采集卡的三个通道, 一个通道用作采集卡的触发信号, 一个通道用于参考信号的采集, 一个通道用于传感信号的采集. 在利用参考干涉仪的光谱进行重采样时, 将参考信号作为参考时钟, 信号的等波数采样可以通过硬件法或者软件法实现. 硬件法是通过将PD探测得到的参考信号作为外部时钟信号输入到采集卡中. 软件法则是将PD探测到的参考采样信号和BPD探测到的传感信号同时输入到采集卡的信号通道中. 由于每次扫描产生的参考信号不能连续稳定存在, 因而不能作为一个可靠的外部时钟. 此外, 采集卡对外部时钟的频率是有限制的, 参考时钟的频率应该为10 MHz左右才可以满足外部时钟采集需求. 故采用了软件法对光谱进行重采样.
为了实现对系统光谱的处理与单个传感器的光谱恢复, 基于LabView开发环境编写了用于数据处理的软件, 用于系统光谱的实时显示与处理. 图8给出了数据处理软件的前面板, 前面板包含多个用于设置采集系统的输入控件, 主要包括可调谐激光器的起始扫描波长、终止波长、扫描速度的输入, 采集卡的采样速率也可以通过输入控件进行设置, 信号通道下拉菜单可以对信号的输入通道进行选择. 此外, 为了提高数据的读取效率, 设置了用于调整采样段长与采样段数的输入控件. 为了直观地显示光谱信息, 图中右侧的4个波形图分别用于显示原始采样光谱、重采样后的光谱、重采样光谱的频谱特性图以及还原后传感器的干涉光谱. 理论上可以做到全部传感器光谱的显示, 由于目前复用规模较小, 为了清晰地显示单个光谱的变化情况, 显示面板中只显示了单个传感器光谱的恢复, 之后会考虑进行多个传感器光谱的同时处理, 利用还原光谱显示控件作为带通滤波的参数选择参考, 构建多个传感器光谱滤波参数数组, 实现多个传感器光谱同时恢复.
图9给出了数据采集处理软件的后面板, 后面板中对数据的主要操作包括采集卡配置、信号重采样、FFT频谱分析、IFFT光谱还原以及光谱存储几个部分, 这3个过程都涉及到大量的数据操作, 开始时采用了LabView内置的VI函数对光谱数据进行处理, 由于光谱数据量较大, 处理效率较低. 之后改用LabView内置的MATLAB脚本VI函数作为数据处理函数, 大大地提高了数据处理速度. 本实验中利用参考干涉仪作为参考时钟对传感光谱进行等频率重采样, 为了方便信号采样, 将时钟信号减去基值后通过比较器将波形从正弦波转换为方波, 当方波两个相邻采样点出现正负跳变时, 采集一个传感信号点, 通过该种方法可以实现快速的传感光谱重采样. 利用FFT对重采样的信号进行频谱特征分析, 之后通过矩形窗截取传感光谱的特征频谱, 实现传感器光谱的恢复.
4. 实验验证
在所设计的大容量传感器复用系统中进行了应变和盐度实验. 本文侧重于分布式干涉式光纤传感器的实现, 因而选取了应变模拟压力的测量效果, 海洋的压力通过增敏结构以应变的形式传递到光纤传感器结构上. 在系统中接入了3个传感器, 传感器1和传感器3为常见的基于单模−空心−单模结构的光纤FPI应变传感器, 该传感器通过在单模光纤中间熔接一段空心光纤制作; 传感器2为利用单模光纤错位熔接制作的双FPI级联的温盐传感器, 该传感器的制作方法可以参考本课题组的论文[12]. 首先将传感器接入到FMCW系统中, 经过采集软件的重采样之后, 得到的复合光谱如图10所示.
之后用FFT分析了复合光谱的频率特性, 图11展示出了光谱的频谱特性, 子图中给出了接入3个传感器的频谱情况, 对应整体频谱中蓝色椭圆圈出的部分. 可以看出, 由于本方案中采用的波长扫描范围较大, 所以光谱的频率分辨力较高, 可以区分同一传感器内的不同反射端面. 图中紫色的特征峰值是由于光纤法兰连接处的反射造成的. 可以明显地看出, 传感器3的光强远大于传感器1和传感器2的光强, 这是为了能够区分传感器位置. 传感器3接入的光为95%, 传感器3的反射光强度接近前两个传感器的20倍, 且3个传感器是等间隔分布的, 通过这种方法可以快速地找出系统中3个传感器的特征频率. 图中绿色的特征频率是由于传感器以及法兰之间的反射光拍频形成的.
图12给出了传感器的应力测试系统, 通过三维位移滑台用于固定光纤, 铁架台用于悬挂光纤, 将砝码悬挂在光纤自由端, 用于给传感结构施加定量的应力. 在弹性范围内, 光纤的应力和所施加的质量之间存在线性对应关系. 因为砝码的质量精度可以做到很高, 本文采用砝码悬挂法对光纤实现精确的应变控制.
对传感器中传感器1进行应力测试后, 从数据处理软件上得到了不同应变下的光谱, 对光谱数据进行了处理. 对获取的干涉光谱进行平滑与寻峰操作, 得到了不同质量下的干涉峰值, 对不同质量下的干涉峰值进行拟合, 得到的拟合结果如图13(a)所示, 可以看出, 随着所施加砝码质量的增加, 干涉光谱的峰值波长表现出红移响应, 传感器的灵敏度可以达到23.35 pm/g, 拟合线性度可以达到0.997. 图13(b)给出了传感器1在同一质量下连续监测38次的波长变化情况, 测量标准差(Standard deviation, SD)可以达到40.85 pm, 这可能是由于环境波动以及光源抖动等因素造成的.
图 13 应力传感器响应特性((a)不同质量下谐振波长拟合效果; (b)固定质量下传感器波长监测)Fig. 13 Responses of the strain sensor ((a) Wavelength fitting result under different weights; (b) Wavelength record under a fixed weight)之后对双FPI级联的温盐传感器的盐度(折射率)响应特性进行测试, 将传感器2放置于设计的盐度传感平台上. 通过胶头滴管向传感平台一侧滴加盐水, 通过吸水纸从另一侧吸收盐水, 将待测浓度的盐水进行3次冲刷用于减小浓度差的影响. 本方案中采用吸水纸的原因在于盐水的表面张力相对于有机溶液较大, 盐水无法可靠浸入传感结构(或者传感器的响应时间较长), 影响测试结果的准确性. 在实际应用中不需要更换液体, 只需要保证待测液体浸入到传感器即可, 可以考虑事先将传感器结构浸泡于有机溶液中. 由于本传感器结构较小, 浸泡需要的有机溶剂较少, 对实际测量产生的影响可以忽略, 测试环境由图14给出. 盐度测试采用的是海洋国家计量中心生产的中国系列标准海水.
图15给出了传感器的光谱情况, 图15(a)表示传感器的原始光谱, 这包含两个FPI传感器的混叠光谱, 通过带通滤波器将两个传感器的光谱进行恢复, 对盐度敏感的传感器光谱在图15(b)中给出, 对温度敏感的传感器光谱在图15(c)中给出. 可以看出, 使用带通滤波器可以很好地区分两个传感器光谱.
本实验中仅对盐度响应特性进行测试, 对不同盐度下的响应光谱进行了分析. 分析结果在图16中给出, 对光谱的特征峰值与盐水浓度进行了线性拟合, 光谱随着盐度的增加表现出右移响应. 传感器的盐度灵敏度可以达到242.58 pm/‰, 拟合线性度可以达到0.9996, 传感器的盐度和波长之间具有很好的响应特性. 该传感器的盐度灵敏度与文献[17]基本是一致的, 可以认为传感器的性能不受到复用系统影响. 传感器随盐度变化的波长移动量达到9.7 nm, 光谱移动量接近光谱周期的2倍, 采用光谱峰值追踪法无法进行有效的光谱处理. 本次实验中, 结合经验法对光谱进行寻峰, 之后的实验可以考虑采用干涉光谱腔长解调技术对光谱进行解调[12], 就可以实现大动态范围、高分辨力的光谱解调.
图17给出了传感器在同一盐度下光谱特征峰值连续监测的效果, 通过计算得到该传感器的波长标准差为20.68 pm, 对应盐度的标准差为0.085‰. 连续监测结果表明, 传感器在该复用系统中仍然可以实现较好的传感效果.
图 17 传感器2盐度光谱的连续监测效果Fig. 17 Continuous wavelength record of the salinity spectrum of sensor 25. 结论
本文提出了一种基于FMCW原理的干涉式光纤传感器的复用方法, 利用传感器的反射端面和参考臂形成的Mach-Zehnder干涉光谱之间的拍频实现传感器光谱还原, 搭建了用于大容量传感器复用的实验系统, 编写了用于数据采集与处理的软件, 实现了数据的实时采集以及对系统内单个传感器光谱的还原, 测试了系统中传感器的温度和盐度响应特性, 实验结果表明, 该复用系统在实现大容量传感器复用的同时不影响单个传感器的传感性能. 所提出的复用方法可以用于基于光纤传感器阵列的海洋参数立体剖面监测以及多点温度监测等应用场合.
-
图 1 工序为$ \pi $ = [1 3 2 5 4 1 3 1 3 4]时的甘特图
Fig. 1 The Gantt chart of $ \pi $ = [1 3 2 5 4 1 3 1 3 4]
图 6 各算法在不同时间参数下的均值线及95 %置信度下Tukey$'$s HSD检验的置信区间
Fig. 6 Means plot and 95 % Tukey$'$s HSD confidence intervals for the interaction between the algorithms and the different time factor $ \rho $
图 7 各算法在不同时间参数下的均值线及95 %置信度下Tukey$'$s HSD检验的置信区间
Fig. 7 Means plot and 95 % Tukey$'$s HSD confidence intervals for the interaction between the algorithms and the different time factor $ \rho $
表 1 工序加工约束、序相关设置时间和工件到达时间表($ \pi $ = [1 3 2 5 4 1 3 1 3 4])
Table 1 The schedule of process constraint, sequence setup time and arrival time ($ \pi $ = [1 3 2 5 4 1 3 1 3 4])
可执行操作的设备 序相关设置时间 首次到达设备的时间 操作1 操作2 操作3 工件1 工件2 工件3 工件4 工件5 m1 m2 m3 工件1 m1/m2 m1/m2/m3 m1/m3 — 83 38 39 47 35 34 23 工件2 m1 — — 53 — 66 45 25 38 78 98 工件3 m3 m2/m3 m2 64 57 — 72 46 52 23 32 工件4 m1/m3 m1 — 46 67 83 — 55 132 131 98 工件5 m1/m3 — — 40 66 83 77 — 114 99 112 
下载: 导出CSV
表 2 工件加工时间表($ \pi $ = [1 3 2 5 4 1 3 1 3 4]) ($t$/s)
Table 2 The processing time table ($ \pi $ = [1 3 2 5 4 1 3 1 3 4]) ($t$/s)
操作1 操作2 操作3 m1 m2 m3 m1 m2 m3 m1 m2 m3 工件1 62 44 — 78 86 26 48 — 87 工件2 41 — — — — — — — — 工件3 — — 31 — 58 65 — 42 — 工件4 32 — 31 27 — — — — — 工件5 74 — 36 — — — — — —
下载: 导出CSV
表 3 参数水平设置表
Table 3 Combinations of parameter values
参数 水平 1 2 3 4 $ popsize $ 10 20 30 40 $ TF $ 0 1 2 3 $ r_{m} $ 0.1 0.4 0.7 0.9
下载: 导出CSV
表 4 正交表和AVG统计
Table 4 Orthogonal array and AVG
参数组合 水平 $ AVG $ $ popsize $ $ TF $ $ r_{m} $ 1 1 1 1 500.65 2 1 2 2 503.90 3 1 3 3 497.30 4 1 4 4 498.55 5 2 1 2 489.35 6 2 2 3 491.20 7 2 3 4 485.50 8 2 4 1 488.35 9 3 1 3 482.55 10 3 2 4 481.80 11 3 3 1 482.25 12 3 4 2 483.25 13 4 1 4 483.35 14 4 2 1 484.00 15 4 3 2 483.60 16 4 4 3 482.20
下载: 导出CSV
表 5 各参数响应值
Table 5 Average response value and rank of each parameter
水平 $ popsize $ $ TF $ $ r_{m} $ 1 500.10 488.98 488.81 2 488.60 490.23 490.02 3 482.46 487.16 488.31 4 483.29 488.09 487.30 等级 1 3 2
下载: 导出CSV
表 6 DTLBO与PSO、DTLBO-Ⅰ、标准TLBO和CIWO的比较($ \rho =4 $)
Table 6 Comparisons of DTLBO, DTLBO-Ⅰ, PSO, CIWO, and standard TLBO ($ \rho =4 $)
测试问题 PSO DTLBO-Ⅰ 标准TLBO CIWO DTLBO BST WST AVG BST WST AVG BST WST AVG BST WST AVG BST WST AVG 10$ \times $5 229 235 230.15 228 231 229.45 228 233 230.05 227 246 229.95 227 244 229.3 20$ \times $5 560 587 576.7 546 579 560.3 558 589 574.5 503 554 526.5 506 534 521.4 30$ \times $5 682 708 694.8 638 685 662.3 675 709 697 615 665 635.5 607 638 622.4 40$ \times $5 752 784 771.4 696 748 725.7 728 787 767.75 654 716 690.05 661 716 692.8 50$ \times $5 1 132 1 171 1 157.45 1 093 1 140 1 114.3 1148 1 175 1 161.5 1 019 1 091 1 059.2 1 019 1 078 1 051.85 40$ \times $10 340 357 348.2 307 336 321.3 340 360 351.15 293 326 311.9 284 317 304.3 50$ \times $10 410 429 419.7 373 411 392.75 413 429 421.1 362 391 376.3 353 386 370.75 60$ \times $10 523 548 538.8 497 529 516 538 552 544.3 485 516 499.95 481 516 500.6 70$ \times $10 658 682 672 626 653 641.8 660 684 672.65 604 639 622 597 636 621.75 80$ \times $10 596 619 612.15 567 601 585.55 597 621 612.95 553 589 571.45 553 584 569.5 80$ \times $20 286 298 292.2 265 289 278.85 285 298 291.8 257 282 271.65 259 272 264.7 90$ \times $20 287 299 294.65 271 289 280.55 288 301 297.4 255 284 273.1 255 279 269.05 100$ \times $20 329 340 336 309 328 320.5 333 343 338 304 323 312.65 299 313 306.25 150$ \times $20 479 496 489.2 460 486 472 487 499 492.9 454 476 464.95 449 477 458.45 200$ \times $20 575 593 587.25 557 576 568.5 577 598 589.85 559 575 564.4 547 559 549.55 Average 522.53 543.06 534.71 495.53 525.4 511.32 523.66 545.2 536.19 476.27 511.53 493.97 473.13 503.27 488.84
下载: 导出CSV
表 7 HDTLBO与DPSO、DTLBO-Ⅱ、GA_DR_C和CCIWO的比较($ \rho =4 $)
Table 7 Comparisons of HDTLBO, DPSO, DTLBO-Ⅱ, GA_DR_C, and CCIWO ($ \rho =4 $)
测试问题 DPSO DTLBO-Ⅱ GA_DR_C CCIWO HDTLBO BST WST AVG BST WST AVG BST WST AVG BST WST AVG BST WST AVG 10$ \times $5 227 230 228.75 227 230 228.45 227 230 228.2 227 230 228.2 227 230 227.3 20$ \times $5 521 547 527.15 533 554 543.1 514 551 536.15 494 536 518.1 499 540 520.95 30$ \times $5 612 647 624.1 630 658 645.65 625 656 645.85 617 644 621.95 603 638 617.5 40$ \times $5 652 692 686.75 695 733 710.15 695 735 715.8 657 715 686.15 645 692 681.65 50$ \times $5 1 014 1 074 1 058.4 1 064 1 111 1 090 1 070 1107 1 094.75 1 024 1 088 1 057.05 1 006 1 055 1 035.9 40$ \times $10 303 317 306.14 303 328 316.2 305 325 316.4 292 320 307.7 290 316 303.65 50$ \times $10 365 387 370.95 370 395 387.95 374 400 388.2 362 390 373.85 355 382 369 60$ \times $10 481 503 498.85 498 519 510.9 497 584 541.15 481 515 496.6 481 501 493.05 70$ \times $10 601 639 615.4 626 656 640.15 636 716 694.55 597 639 621.7 597 635 606.75 80$ \times $10 556 590 579 571 596 584.55 561 590 582.75 548 586 568.45 554 575 564.85 80$ \times $20 274 285 277.41 269 285 278 273 285 279.9 275 280 269.3 262 279 273.7 90$ \times $20 269 290 273.33 275 288 282.1 277 292 285.8 268 285 276.25 266 282 272.58 100$ \times $20 304 323 313.7 313 327 320.85 317 328 322.85 300 323 314.6 297 323 310.65 150$ \times $20 472 489 475.34 469 486 476.45 470 488 480.55 469 486 470.7 464 483 469.89 200$ \times $20 553 579 566.5 564 579 573 562 587 578.1 557 583 570.05 549 579 565.2 Average 480.27 506.13 493.45 493.8 516.33 505.83 493.53 524.93 521.73 477.87 508 492.04 473 500.67 487.51
下载: 导出CSV
表 8 模具各工序加工约束及首次到达时间表
Table 8 The schedule of mold process constraint and arrival time
模具 可执行操作的设备 模具首次到达设备时间 操作1 操作2 操作3 m1 m2 m3 m4 m5 1 m1/m3/m5 m3/m4 m2/m3/m4 7 4 8 4 6 2 m1/m3/m6 m1 — 1 6 2 3 3 3 m3 — — 3 1 8 4 6 4 m2/m3/m5 m1/m2/m3/m4/m5 m2/m3/m4 8 1 6 8 7 5 m1/m3/m4 — — 1 2 6 7 3 6 m2/m5 m1/m4/m5 m1/m4/m5 5 8 3 1 5 7 m2/m4/m5 m2/m3/m4/m5 — 3 8 3 3 6 8 m1/m2/m3/m5 — — 2 6 8 7 7 9 m4/m5 m3 m3/m4 6 3 5 3 7 10 m2 m2/m4 m1/m2/m3 7 4 1 4 6 11 m1/m5 m1 — 6 5 1 8 4 12 m4/m5 m1/m2/m3/m4/m5 m1/m2/m4 5 3 1 8 8 13 m2/m4/m5 m3 m4/m5 7 1 6 6 6 14 m3 m1/m3/m5 2 1 4 6 3 15 m3/m4/m5 — — 8 6 5 3 2 16 m1/m3 m3 — 8 6 6 7 2 17 m3 — — 8 4 6 7 1 18 m4/m5 m2/m4 m5 4 3 5 2 6 19 m1/m2/m4/m5 — — 3 6 2 1 3 20 m1/m4 m3/m5 m2/m4/m5 4 7 5 7 7
下载: 导出CSV
表 9 各模具产品的序相关设置时间表
Table 9 The schedule of sequence setup time between each mold
模具 序相关设置时间 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 1 — 2 3 3 4 3 3 3 3 2 1 3 1 4 4 4 2 4 2 1 2 2 — 2 2 3 2 2 1 3 1 2 2 4 4 3 3 4 4 2 4 3 2 3 — 4 3 3 4 2 1 2 3 3 1 4 2 1 3 3 4 3 4 3 1 3 — 4 1 4 2 2 2 1 1 4 1 3 2 2 1 2 3 5 4 4 4 2 — 1 1 1 4 1 1 1 3 2 1 3 2 4 3 3 6 1 2 3 2 3 — 3 2 2 3 4 1 1 3 3 2 4 2 4 1 7 3 4 4 4 1 2 — 1 1 4 2 2 3 2 2 2 4 1 1 2 8 4 2 4 1 4 1 2 — 4 3 3 3 2 4 1 3 4 1 4 1 9 4 2 3 3 2 4 4 1 — 2 2 3 2 1 2 4 2 1 4 3 10 2 4 1 2 2 4 4 3 3 — 1 1 2 3 1 2 1 4 2 1 11 3 4 3 2 2 2 1 2 2 2 — 2 4 3 3 2 1 2 2 2 12 2 1 4 2 4 2 3 2 1 2 3 — 1 1 3 2 3 2 1 2 13 2 1 2 4 2 1 4 2 4 2 4 3 — 4 4 3 4 3 4 1 14 4 4 1 3 2 2 2 4 1 4 1 2 1 — 2 2 1 3 4 1 15 3 2 2 3 1 2 3 4 3 1 2 1 3 1 — 2 3 1 4 4 16 1 4 1 1 4 3 1 3 1 2 4 4 2 3 4 — 4 4 4 1 17 4 3 1 1 3 2 3 1 3 4 1 1 3 2 1 1 — 2 1 1 18 1 2 4 4 2 1 4 1 2 1 4 4 3 2 1 3 1 — 2 2 19 2 3 2 2 2 3 4 1 1 3 2 2 2 3 3 2 4 1 — 2 20 2 2 2 2 2 2 2 1 2 1 2 4 1 4 1 1 1 4 1 —
下载: 导出CSV
表 10 模具各工序的加工时间表
Table 10 The processing time table of each mold
模具 操作1 操作2 操作3 m1 m2 m3 m4 m5 m1 m2 m3 m4 m5 m1 m2 m3 m4 m5 1 27 — 18 — 29 — — 27 17 — — 23 29 19 — 2 13 — 26 — 11 13 — — — — — — — — — 3 — — 23 — — — — — — — — — — — — 4 — 10 13 — 12 18 25 11 10 26 — 21 31 27 — 5 9 — 8 16 — — — — — — — — — — — 6 — 9 — — 29 19 — — 25 10 30 — — 26 19 7 — 19 — 22 28 — 10 28 24 27 — — — — — 8 30 23 28 — 17 — — — — — — — — — — 9 — — — 31 24 — — 9 — — — — 8 28 — 10 — 25 — — — — 27 — 13 — 31 20 26 — — 11 25 — — — 24 12 — — — — — — — — — 12 — — — 23 — 17 13 16 — — 22 8 — 30 — 13 — 16 — 19 22 — — 28 — — — — — 27 23 14 28 — 16 — 20 — — 18 — — 16 — 16 — 21 15 — — 28 29 23 — — — — — — — — — — 16 27 — 10 — — — — 8 — — — — — — — 17 — — 26 — — — — — — — — — — — — 18 — — — 12 24 — 13 — 27 — — — — — 17 19 24 19 — 12 14 — — — — — — — — — — 20 17 — — 10 — — — 31 — 20 — 30 — 9 24
下载: 导出CSV
表 11 实例仿真结果
Table 11 Simulation results of the instance
运行次数 DTLBO-Ⅱ GA DPSO CCIWO HDTLBO 1 167 166 166 163 163 2 164 166 164 163 166 3 168 167 167 167 167 4 169 166 164 166 163 5 163 163 165 163 165 6 167 165 166 167 163 7 167 168 163 168 166 8 166 168 167 164 168 9 168 167 164 163 163 10 169 168 165 164 165 11 168 163 164 168 166 12 169 168 169 167 163 13 168 167 165 166 167 14 163 170 166 170 166 15 167 166 168 169 163 16 166 167 166 171 163 17 167 170 163 168 163 18 169 168 164 170 165 19 165 168 163 165 168 20 167 168 164 167 163 Average 166.85 166.95 165.15 166.45 164.8
下载: 导出CSV
-
[1] Pinedo M. Scheduling: theory, algorithms, and systems. Berlin Heidelberg: Springer, 2012. [2] Allahverdi A. Two-machine proportionate flowshop scheduling with breakdowns maximum lateness. Computers & Operations Research, 1996, 23(10): 909-916 [3] Lan X J, Su D D. Research of a mold job shop scheduling optimization based on particle swarm optimization algorithm. Applied Mechanics & Materials, 2015, 757(2): 201-207 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=10.4028/www.scientific.net/AMM.757.201 [4] Sang H Y, Duan P Y, Li J Q. An effective invasive weed optimization algorithm for scheduling semiconductor final testing problem. Swarm & Evolutionary Computation, 2018, 38: 42-53 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=f24ee1e67ed4add61940b0e69b5f0a53 [5] Hu H, Ng KKH, Qin Y C. Robust parallel machine scheduling problem with uncertainties and sequence-dependent setup time. Scientific Programming, 2016: 1-13 [6] Li C L, Cheng T C E. The parallel machine min-max weighted absolute lateness scheduling problem. Naval Research Logistics, 1994, 41(1): 33-46 doi: 10.1002/1520-6750(199402)41:1<33::AID-NAV3220410104>3.0.CO;2-S [7] Ranjbar M, Davari M, Leus R. Two branch-and-bound algorithms for the robust parallel machine scheduling problem. Computers & Operations Research, 2012, 39(7): 1652-1660 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=58dbb13e91377bee19b577c7d4561726 [8] Lee W C, Wang J Y, Lin M C. A branch-and-bound algorithm for minimizing the total weighted completion time on parallel identical machines with two competing agents. Knowledge-Based Systems, 2016, 105(C): 68-82 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=90998daca3d98ce1fbd2b119fbae1b4a [9] Wu L, Wang S. Exact and heuristic methods to solve the parallel machine scheduling problem with multi-processor tasks. International Journal of Production Economics, 2018, 201: 26-40 doi: 10.1016/j.ijpe.2018.04.013 [10] Chen Z L, Powell W B. A column generation based decomposition algorithm for a parallel machine just-in-time scheduling problem. European Journal of Operational Research, 1999, 116(1): 220-232 doi: 10.1016/S0377-2217(98)00136-2 [11] Yin Y, Ye D, Zhang G. Single machine batch scheduling to minimize the sum of total flow time and batch delivery cost with an unavailability interval. Information Sciences, 2014, 274(8): 310-322 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=057211486faaa3fa789b3f4b9a751866 [12] Li D W, Lu X W. Two-agent parallel-machine scheduling with rejection. Theoretical Computer Science, 2017, 703: 66-75 doi: 10.1016/j.tcs.2017.09.004 [13] Yin Y, Cheng S R, Cheng T C E, et al. Just-in-time scheduling with two competing agents on unrelated parallel machines. Omega, 2016, 63: 41-47 doi: 10.1016/j.omega.2015.09.010 [14] Woo Y B, Jung S, Kim B S. A rule-based genetic algorithm with an improvement heuristic for unrelated parallel machine scheduling problem with time-dependent deterioration and multiple rate-modifying activities. Computers & Industrial Engineering, 2017, 109: 179-190 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=f0a6c39ef977f0f9c653802880bf2420 [15] Avalos-Rosales O, Angel-Bello F, Alvarez A. Efficient metaheuristic algorithm and re-formulations for the unrelated parallel machine scheduling problem with sequence and machine-dependent setup times. International Journal of Advanced Manufacturing Technology, 2015, 76(9-12): 1705-1718 doi: 10.1007/s00170-014-6390-6 [16] 张嘉琦, 曹政才, 刘民.融合代理模型和差分进化算法的并行机动态调度方法.计算机集成制造系统, 2017, 23(1): 75-81 http://d.old.wanfangdata.com.cn/Periodical/jsjjczzxt201701009Zhang Jia-Qi, Cao Zheng-Cai, Liu Min. Integrated differential evolution algorithm with surrogate model for dynamic parallel machine scheduling. Computer Integrated Manufacturing Systems, 2017, 23(1): 75-81 http://d.old.wanfangdata.com.cn/Periodical/jsjjczzxt201701009 [17] Chen C L. Iterated hybrid metaheuristic algorithms for unrelated parallel machines problem with unequal ready times and sequence-dependent setup times. International Journal of Advanced Manufacturing Technology, 2012, 60(5-8): 693-705 doi: 10.1007/s00170-011-3623-9 [18] Damodaran P, Diyadawagamage D A, Ghrayeb O. A particle swarm optimization algorithm for minimizing makespan of nonidentical parallel batch processing machines. International Journal of Advanced Manufacturing Technology, 2012, 58(9-12): 1131-1140 doi: 10.1007/s00170-011-3442-z [19] Diana R O M. An immune-inspired algorithm for an unrelated parallel machines$'$ scheduling problem with sequence and machine dependent setup-times for makespan minimisation. Neurocomputing, 2015, 163(C): 94-105 [20] Sels V, Coelho J, Dias A M, et al. Hybrid tabu search and a truncated branch-and-bound for the unrelated parallel machine scheduling problem. Computers & Operations Research, 2015, 53: 107-117 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=a7d0c31a9a95e9f1ce534a14dda1cdb2 [21] Bitar A, Yugma C, Roussel R. A memetic algorithm to solve an unrelated parallel machine scheduling problem with auxiliary resources in semiconductor manufacturing. Journal of Scheduling, 2016, 19(4): 367-376 doi: 10.1007/s10951-014-0397-6 [22] Joo C M, Kim B S. Hybrid genetic algorithms with dispatching rules for unrelated parallel machine scheduling with setup time and production availability. Computers & Industrial Engineering, 2015, 85(C): 102-109 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=fcb17bbd465d8682d2400ec1e2f8025a [23] 罗家祥, 唐立新.带释放时间的并行机调度.问题的ILS & SS算法.自动化学报, 2005, 31(6): 917-924 http://www.aas.net.cn/article/id/15921Luo Jia-Xiang, Tang Li-Xin. A new ILS & SS algorithm for parallel-machine scheduling problem. Acta Automatica Sinica, 2005, 31(6): 917-924 http://www.aas.net.cn/article/id/15921 [24] Lin S W, Ying K C. ABC-based manufacturing scheduling for unrelated parallel machines with machine-dependent and job sequence-dependent setup times. Computers & Operations Research, 2014, 51(3): 172-181 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=5b16d15a9448bceeb8e453a1f542ce32 [25] Rao R V, Savsani V J. and Vakharia D P. Teaching- learning-based optimization: a novel method for constrained mechanical design optimization problems. Computer-Aided Design, 2011, 43(3): 303-315 doi: 10.1016/j.cad.2010.12.015 [26] Waghmare, G. Comments on "a note on teaching-learning-based optimization algorithm". Information Sciences, 2013, 229: 159-169 doi: 10.1016/j.ins.2012.11.009 [27] Shao S W, Pi D C, Shao Z. An extended teaching-learning based optimization algorithm for solving no-wait flow shop scheduling problem. Applied Soft Computing, 2017, 61: 193-210 doi: 10.1016/j.asoc.2017.08.020 [28] Xu Y, Wang L, Wang S Y, et al. An effective teaching--learning-based optimization algorithm for the flexible job-shop scheduling problem with fuzzy processing time. Neurocomputing, 2015, 148: 260-268 doi: 10.1016/j.neucom.2013.10.042 [29] Shao S W, Pi D C, Shao Z. A hybrid discrete optimization algorithm based on teaching-probabilistic learning mechanism for no-wait flow shop scheduling. Knowledge-Based Systems, 2016, 107: 219-234 doi: 10.1016/j.knosys.2016.06.011 [30] Li J Q, Pan Q K, Mao K. A discrete teaching-learning-based optimisation algorithm for realistic flowshop rescheduling problems. Engineering Applications of Artificial Intelligence, 2015, 37: 279-292 doi: 10.1016/j.engappai.2014.09.015 [31] 马永杰, 云文霞.遗传算法研究进展.计算机应用研究, 2012, 29(4): 1201-1206 doi: 10.3969/j.issn.1001-3695.2012.04.001Ma Yong-Jie, Yun Wen-Xia. Research progress of genetic algorithm. Application Rehash of Computers, 2012, 29(4): 1201-1206 doi: 10.3969/j.issn.1001-3695.2012.04.001 [32] Abdoun O, Abouchabaka J. A comparative study of adaptive crossover operators for genetic algorithms to resolve the traveling salesman problem. Computer Science, 2011, 31(11): 49-57 [33] 王凌, 钱斌.混合差分进化与调度算法.北京:清华大学出版社, 2012.Wang Ling, Qian Bin. Hybrid differential evolution and scheduling algorithms. Beijing: Tsinghua University Press, 2012. [34] 潘全科, 王凌, 赵保华.解决零空闲流水线调度问题的离散粒子群算法.控制与决策, 2008, 23(2): 191-194 doi: 10.3321/j.issn:1001-0920.2008.02.014Pan Quan-Ke, Wang Ling, Zhao Bao-Hua. Discrete particle swarm optimization for no-idle flow shop problem. Control and Decision, 2008, 23(2): 191-194 doi: 10.3321/j.issn:1001-0920.2008.02.014 期刊类型引用(13)
1. 胡蓉，伍星，毛剑琳，钱斌. 融入概率学习的混合差分进化算法求解绿色分布式可重入作业车间调度. 控制理论与应用. 2024(03): 512-521 . 
百度学术2. 杨媛媛，胡蓉，钱斌，张长胜，金怀平. 深度强化学习算法求解动态流水车间实时调度问题. 控制理论与应用. 2024(06): 1047-1055 . 百度学术3. 郭琴，郑巧仙. 基于历史最优与接受差解的教与学算法. 计算机工程与设计. 2024(11): 3499-3506 . 百度学术4. 王建华，杨琦，朱凯. 自适应多种群Jaya算法求解绿色并行机调度问题. 计算机集成制造系统. 2023(01): 111-120 . 百度学术5. 王娟，刘美红，祝世兴，陈文博，李遇贤，孙军锋. 指尖密封结构参数的多目标混合教与学优化方法. 机械工程学报. 2023(09): 157-170 . 百度学术6. 翟志波，贾国平，王涛，周鹏鹏，闫汝山，戴玉森. 基于拉普拉斯分布与鲍德温效应的教与学算法及其应用. 计算机集成制造系统. 2023(08): 2611-2621 . 百度学术7. 徐林浩，胡蓉，钱斌，于乃康. 绿色VRP的增强拉格朗日松弛启发式算法. 电子测量技术. 2023(19): 96-103 . 百度学术8. 陶丽文，刘建胜. 一种多班级离散型教与学算法求解路径规划问题. 南昌大学学报(工科版). 2022(03): 282-288 . 百度学术9. 宋君乐，陶翼飞，可晓东，丁小鹏. 分区编码HBGWO算法求解并行机分批调度问题. 机械设计. 2022(10): 48-59 . 百度学术10. 马晓梅，何非. 基于改进遗传算法的标签印刷生产调度技术. 计算机应用. 2021(03): 860-866 . 百度学术11. 曹劲松，熊福力. 基于离散教与学算法的分布式预制流水车间调度研究. 计算机测量与控制. 2021(12): 166-171+183 . 百度学术12. 张梓琪，钱斌，胡蓉. 混合交叉熵算法求解复杂零等待流水线调度问题. 控制理论与应用. 2021(12): 1919-1934 . 百度学术13. 李丽荣，李木子，李崔灿，王培崇. 具有动态自适应学习机制的教与学优化算法. 计算机工程与应用. 2020(19): 62-67 . 百度学术其他类型引用(17)
-
下载:
计量
- 文章访问数: 6402
- HTML全文浏览量: 1562
- PDF下载量: 200
- 被引次数: 30
