一种基于距离约束副本交换的蛋白质结构预测方法
【专利摘要】一种基于距离约束副本交换的蛋白质结构预测方法，首先在每个温度层通过对查询序列进行随机的折叠和变换生成一个初始种群；在种群更新中，以Rosetta Score3为优化目标函数，基于Anfinsen提出的蛋白质天然态结构时自由能最低的结构，在每一个温度层中依次将种群中的每个个体作为目标个体，然后随机选取和目标个体不一样的两个个体进行变异、交叉产生变异个体，再将另一个个体随机选取一段与变异个体进行交换产生测试个体，再将测试个体与目标个体进行能量值比较，对于能量升高的测试个体引入距离谱的知识，将相邻的温度层的对应个体进行副本交换。本发明构象空间采样能力较好、预测精度高。
【专利说明】
一种基于距离约束副本交换的蛋白质结构预测方法
技术领域
[0001] 本发明涉及生物信息学、计算机应用领域，尤其涉及的是一种基于距离约束副本 交换的蛋白质结构预测方法。
【背景技术】
[0002] 蛋白质分子在生物细胞化学反应过程中起着至关重要的作用。它们的结构模型和 生物活性状态对我们理解和治愈多种疾病有重要的意义。蛋白质只有折叠成特定的三维结 构才能产生其特有的生物学功能。因此，要了解蛋白质的功能，就必须获得其三维空间结 构。
[0003] 蛋白质三级结构预测是生物信息学的一个重要任务。蛋白质构象优化问题现在面 临最大的挑战是对极其复杂的蛋白质能量函数曲面进行搜索。蛋白质能量模型考虑了分子 体系成键作用以及范德华力、静电、氢键、疏水等非成键作用，致使其形成的能量曲面极其 粗糙，构象对应局部极小解数目随序列长度的增加呈指数增长。而蛋白质构象预测算法能 够找到蛋白质稳定结构的机理是，大量的蛋白质亚稳定结构构成了低能量区域，所以能否 找到蛋白质全局最稳定结构的关键是算法能够找到大量的蛋白质亚稳定结构，即增加算法 的种群多样性。因此，针对更加精确的蛋白质力场模型，选取有效的构象空间优化算法，使 新的蛋白质结构预测算法更具有普遍性和高效性成为生物信息学中蛋白质结构预测的焦 点问题。
[0004] 目前，蛋白质结构预测方法大致可以分为两类，基于模板的方法和不基于模板的 方法。其中，不基于模板的从头预测(Ab-inito)方法应用最为广泛。它适用于同源性小于 25%的大多数蛋白质，仅从序列产生全新结构，对蛋白质分子设计及蛋白质折叠的研究等 具有重要意义。当前有以下几种比较成功的从头预测方法:张阳与Jeffrey Skolnick合作 的TASSER(Threading/Assembly/Refinement)方法、David Baker及团队设计的Rosetta方 法、Shehu等设计的FeLTr方法等。但是到目前还没有一种十分完善的方法来预测蛋白质的 三维结构，即使获得了很好的预测结果，但也只是针对某些蛋白质而言的，目前主要的技术 瓶颈在于两个方面，第一方面在于采样方法，现有技术对构象空间采样能力不强，另一方面 在于构象更新方法，现有技术对构象的更新精度仍然不足。
[0005] 因此，现有的构象空间优化方法存在采样效率、复杂度及预测精度方面存在不足， 需要改进。
【发明内容】
[0006] 为了克服现有群体构象空间优化方法的构象空间搜索维数较高、收敛速度较慢、 预测精度较低的不足，本发明提出一种构象空间采样能力较好、预测精度高的基于距离约 束副本交换的蛋白质结构预测方法。
[0007] 本发明解决其技术问题所采用的技术方案是：
[0008] -种基于距离约束副本交换的蛋白质结构预测方法，所述优化方法包括以下步 骤：
[0009] 1)给定输入序列信息；
[0010] 2)根据序列信息从QUARK服务器上获取距离谱文件DP，rpk为距离谱中记录下的残 基对，Dk为该残基对之间的距离，其中k e (1，N)，N为距离谱中残基对数量；
[0011] 3)初始化：设置群体规模popSize，交叉概率CR，8个温度层T，能量函数选用 Rosetta Score3,首先在每个温度层通过对查询序列随机的折叠和变换，生成一个规模为 popSize的初始种群，初始群体为Ρτ= {x11 i e 1}，计算能量函数值f^x1)，i e I，并设 ?/P =min/C〇，其中i为种群个体编号，I为种群个体编号集合，1 = {1，2,. . .，popSize}，/；： 为对应温度层能量函数值中的最小值，T为温度层；
[0012] 4)开始迭代，依次对每个温度层种群中的每个个体进行以下操作：
[0013] 4.1)设i = l，其中斤{1，2,3，."，陬};令卩响(^ = ?1，其中1为序号，？恤_表示目标 个体；
[0014] 4.2)随机生成正整数randl，rand2，rand3e {1，2,3,......NP}，且randl辛rand2 辛rand3辛i ;再生成4个随机整数randrangel，randrange2，randrange3，randrange4;其中 randrangel辛randrange2，randrange3辛randrange4G {1，2，···，Length}，Length为序列长 度；
[0015] 4.3)针对个体?進(11做变异操作；
[0016] 4·3· 1)令a=min(randrangel，randrange2)，b=max(randrangel，randrange2)，k e [a,b] ；^c = min(randrange3,randrange4) ,d=max(randrange3,randrange4) ,pe [c, d];其中min表示取两个数的最小值，max表示取两个数的最大值；
[0017] 4.3.2)用卩-2上位置&至丨」位置13的片段的氨基酸所对应的二面角口11；[、口8；[、011^&替 换Prandl的相同位置所对应的二面角口11；^8；[、011^3;再使用?-3上位置0到位置(1的片段的 氨基酸所对应的二面角口11；[48；[、01]1683替换？1^(]1上相同位置所对应的二面角口11；[、口8；[、 omega，再将所得Prandl进行片段组装得到变异个体Pmutant;
[0018] 4.4)针对变异个体Pmutant做交叉操作；
[0019] 4·4· 1)生成随机数rand4，rand5,其中rand4e (〇, 1)，rand5e (1，Length);
[0020] 4.4.2)若随机数rand4〈 = CR，则使用目标个体ptarget的片段rand5替换变异个体 P?tant中对应的片段生成测试个体Ptraii，否则Ptraii直接等于变异个体P? tant;
[0021 ] 4.5)针对每个温度层中的目标个体Ptarg(3t和测试个体Ptrall进行选择操作；
[0022] 4 · 5 · 1)计算Ptarget和Ptraii的能量:E(Ptarget)和E(P trai 1);
[0023] 4.5.2)若E(Ptarget)>E(Ptrail)则用P trail替换Ptarget，并且进入4.6);否则继续执行 4.5.3)；
[0024] 4.5.3)计算目标个体Ptarget和测试个体Ptrail在残基对rPi的距离相对于距离谱中 Di 之间的距离差▽ targc'tD,·和▽丨.railD,·:;.
[0025] 4.5.4)判断▽ targetD,和▽ trailD,的值，若大于6，则令其等于6，从而去除一些距离 差值较大的异常数值；
[0026] 4. 5. 5)计算距离谱中所有位置的目标个体和测试个体的距离差值总和
[0027] 4 · 5 · 6)若'_ > 响，则进入5);
[0028] 4.5 · 7)若vA?a <vAa,get，则产生一个(0，1)之间的随机数 rand6，若rand6 小于0 · 2， 贝1J用Ptraii替换Ptarget，否则进行5);
[0029] 5)在每个温度层都进行变异、交叉和选择更新完种群后，将相邻的温度层的对应 个体进行副本交换：
[0030] 5.1)首先选择种群中第i个个体，ie I(i = l，2,…，popSize);
[0031] 5.2)选择两个相邻的温度层Tj和Tj+1，j = 1，2，…，7;
[0032] 5.3)产生一个随机数randxe (〇，1)，和判别数judge
[0033]
[0034] 其中/(x()(i = l，2, · · ·，p〇pSiZe)为种群个体的能量值，K为玻尔兹曼常数， 1. .1 1'心=1，2，~，7)为第」个温度层的开尔文温度，」为索引编号^为自然常数；
[0035] 5.4)若randx< judge，将相邻温度层中的么和七.个体进行交换；
[0036] 6)判断是否满足终止条件，如不满足则转至4);若满足则输出结果。
[0037] 进一步，所述步骤3)中，设置最大迭代次数iteration;所述步骤6)中，对每个温度 层种群中的每个个体都执行完步骤4)_5)以后，迭代次数G = G+1，终止条件为迭代次数G达 到预设最大迭代次数iteration。
[0038]本发明的技术构思为:在基本DE算法的框架下，首先在每个温度层通过对查询序 列进行随机的折叠和变换生成一个初始种群;在种群更新中，以Rosetta Sc〇re3为优化目 标函数，基于Anf insen提出的蛋白质天然态结构时自由能最低的结构，在每一个温度层中 依次将种群中的每个个体作为目标个体，然后随机选取和目标个体不一样的两个个体进行 变异、交叉产生变异个体，再将另一个个体随机选取一段与变异个体进行交换产生测试个 体，再将测试个体与目标个体进行能量值比较，对于能量升高的测试个体引入距离谱的知 识，使得种群中保留能量虽高但结构更为合理的构象。在更新完种群后，将相邻两个温度层 之间种群个体进行副本交换来增加种群的多样性，通过不断更新种群以及副本交换获得一 系列亚稳态构象。本发明提供了一种预测精度较高、复杂度较低的构象空间优化方法。
【附图说明】：
[0039]图1是测试序列2頂U在种群更新过程中RMSD和能量值的关系示意图。
[0040] 图2是测试序列2IMU在本算法预测所得诱导构象与基本差分进化预测所得构象的 采样概率分布图。
【具体实施方式】：
[0041] 下面结合附图对本发明作进一步描述。
[0042] 参照图1和图2,一种基于距离约束副本交换的蛋白质结构预测方法，所述构象空 间优化方法包括以下步骤：
[0043] 1)给定输入序列信息；
[0044] 2)根据序列信息从QUARK服务器上获取距离谱文件DP，rpk为距离谱中记录下的残 基对，Dk为该残基对之间的距离，其中k e (1，N)，N为距离谱中残基对数量；
[0045] 3)初始化：设置群体规模popSize，交叉概率CR，8个温度层T，能量函数选用 Rosetta Score3,首先在每个温度层通过对查询序列随机的折叠和变换，生成一个规模为 popSize的初始种群，初始群体为Ρτ= {x11 i e 1}，计算能量函数值f^x1)，i e I，并设 ^ ，其中i为种群个体编号，〗为种群个体编号集合^ = …，popSize}，< 为对应温度层能量函数值中的最小值，T为温度层；
[0046] 4)开始迭代，依次对每个温度层种群中的每个个体进行以下操作：
[0047] 4.1)设1 = 1，其中斤{1，2,3，."，陬};令?恤_ = ?1，其中1为序号，卩*_表示目标 个体；
[0048] 4.2)随机生成正整数randl，rand2，rand3e {1，2,3,......NP}，且randl辛rand2 辛rand3辛i ;再生成4个随机整数randrangel，randrange2，randrange3，randrange4;其中 randrangel辛randrange2，randrange3辛randrange4G {1，2，···，Length}，Length为序列长 度；
[0049] 4.3)针对个体?^"(11做变异操作；
[0050] 4·3· 1)令a=min(randrangel，randrange2)，b=max(randrangel，randrange2)，k e [a,b] ；^c = min(randrange3,randrange4) ,d=max(randrange3,randrange4) ,pe [c, d];其中min表示取两个数的最小值，max表示取两个数的最大值；
[0051 ] 4.3.2)用Prand2上位置a到位置b的片段的氨基酸所对应的二面角ph i、p s i、ome ga替 换Prandl的相同位置所对应的二面角口11；^8；[、011^3;再使用?-3上位置0到位置(1的片段的 氨基酸所对应的二面角口11；[48；[、01]1683替换？1^(]1上相同位置所对应的二面角口11；[、口8；[、 omega，再将所得Prandl进行片段组装得到变异个体Pmutant;
[0052] 4.4)针对变异个体Pmutant做交叉操作；
[0053] 4·4· 1)生成随机数rand4，rand5,其中rand4e (〇, 1)，rand5e (1，Length);
[0054] 4.4.2)若随机数rand4〈 = CR，则使用目标个体Ptarget的片段rand5替换变异个体 Pmutant中对应的片段生成测试个体Ptraii，否则Ptraii直接等于变异个体P mutant;;
[0055] 4.5)针对每个温度层中的目标个体?__和测试个体?^11进行选择操作；
[0056] 4 · 5 · 1)计算Ptarget和Ptraii的能量:E(Ptarget)和E(P traii);
[0057] 4.5.2)若E(Ptarget)>E(Ptrail)则用P trail替换Ptarget，并且进入4.6);否则继续执行 4.5.3)；
[0058] 4.5.3)计算目标个体Ptarget和测试个体Ptrail在残基对r Pi的距离相对于距离谱中 Di之间的距离差▽ targetD,.和▽ trailD,.;
[0059] 4.5.4)判断▽ targetD,:和▽ trailD,·的值，若大于6，则令其等于6，从而去除一些距离 差值较大的异常数值；
[0060] 4. 5. 5)计算距离谱中所有位置的目标个体和测试个体的距离差值总和
[0061 ] 4 · 5 · 6)若▽lu > ▽&辆，则进入5);
[0062] 4.5.7)若ν%?Γ< _￥，则产生一个(0，1)之间的随机数rand6，若rand6小于0.2， 贝1J用Ptraii替换Ptarget，否则进行5);
[0063] 5)在每个温度层都进行变异、交叉和选择更新完种群后，将相邻的温度层的对应 个体进行副本交换：
[0064] 5.1)首先选择种群中第i个个体，ie I(i = l，2,…，popSize);
[0065] 5 · 2)选择两个相邻的温度层Tj和Tj+1，j = 1，2，…，7;
[0066] 5.3)产生一个随机数randxe (〇, 1)，和判别数judge
[0067]
[0068] 其中/(<)(i = l，2, · · ·，p〇pSiZe)为种群个体$的能量值，K为玻尔兹曼常 数，L(j = l，2,…，7)为第j个温度层的开尔文温度，j为索引编号，e为自然常数；
[0069] 5.4)若randx< judge，将相邻温度层中的'和·个体进行交换；
[0070] 6)判断是否满足终止条件到达迭代次数iteration，如不满足则转至4);若满足则 输出结果。
[0071] 进一步，所述步骤3)中，设置最大迭代次数iteration;所述步骤6)中，对每个温度 层种群中的每个个体都执行完步骤4)_5)以后，迭代次数G = G+1，终止条件为迭代次数G达 到预设最大迭代次数iteration。
[0072] 本实施例以序列长度为46的蛋白质2MU为实施例，一种基于距离约束副本交换的 蛋白质结构预测方法，其中包含以下步骤：
[0073] 1)输入蛋白质序列2頂U;
[0074] 2)根据序列信息从QUARK服务器上获取距离谱文件DP，rpk为距离谱中记录下的残 基对，Dk为该残基对之间的距离，其中ke (1，N)，N为距离谱中残基对数量，测试序列2MU的 距离谱的N=38;
[0075] 3)初始化:设置群体规模popSize = 50，交叉概率CR = 0.5，8个温度层T，每个温度 层的温度参数分别为200，239，286，342，409，489，585，700，单位为开尔文，最大迭代次数 iteration = 10000，能量函数选用Rosetta Score3，首先在每个温度层通过对查询序列随 机的折叠和变换，生成一个规模为popSize的初始种群，初始群体为Ρτ= {χ11 i e 1}，计算能 量函数值fVU1)，i e I，并设./丨=mip /YV)，其中i为种群个体编号，I为种群个体编号集合， 1 = {1，2,. . .，p〇pSiZe}，./；*为对应温度层能量函数值中的最小值，T为温度层；
[0076] 4)开始迭代，依次对每个温度层种群中的每个个体进行以下操作：
[0077] 4.1)设1 = 1，其中1￡{1，2,3，."，陬};令?恤_ =卩1，其中1为序号上虹_表示目标 个体；
[0078] 4.2)随机生成正整数randl，rand2，rand3e {1，2,3,......NP}，且randl辛rand2 辛rand3辛i ;再生成4个随机整数randrangel，randrange2，randrange3，randrange4;其中 randrangel辛randrange2，randrange3辛randrange4G {1，2，···，Length}，Length为序列长 度；
[0079] 4.3)针对个体?進(11做变异操作；
[0080] 4·3· 1)令a=min(randrangel，randrange2)，b=max(randrangel，randrange2)，k e [a,b] ；^c = min(randrange3,randrange4) ,d=max(randrange3,randrange4) ,pe [c, d];其中min表示取两个数的最小值，max表示取两个数的最大值；
[0081 ] 4.3.2)用Prand2上位置a到位置b的片段的氨基酸所对应的二面角phi、psi、omega替 换Prandl的相同位置所对应的二面角口11；^8；[、011^3;再使用?-3上位置0到位置(1的片段的 氨基酸所对应的二面角口11；[48；[、01]1683替换？1^(]1上相同位置所对应的二面角口11；[、口8；[、 omega，再将所得Prandl进行片段组装得到变异个体Pmutant;
[0082] 4.4)针对变异个体Pmutant做交叉操作；
[0083] 4·4· 1)生成随机数rand4，rand5,其中rand4e (〇, 1)，rand5e (1，Length);
[0084] 4.4.2)若随机数rand4〈 = CR，则使用目标个体Ptarget的片段rand5替换变异个体 P?tant中对应的片段生成测试个体Ptraii，否则Ptraii直接等于变异个体P? tant;
[0085] 4.5)针对每个温度层中的目标个体?__和测试个体?^11进行选择操作；
[0086] 4 · 5 · 1)计算Ptarget和Ptraii的能量:E(Ptarget)和E(P traii);
[0087] 4.5.2)若E(Ptarget)>E(Ptrail)则用P trail替换Ptarget，并且进入4.6);否则继续执行 4.5.3)；
[0088] 4.5.3)计算目标个体Ptarget和测试个体Ptrail在残基对r Pi的距离相对于距离谱中 Di之间的距离差▽ targeiD,.fllV trailDy
[0089] 4.5.4)判断▽ targets和▽ trailD,的值，若大于6，则令其等于6，从而去除一些距离 差值较大的异常数值；
[0090] 4. 5. 5)计算距离谱中所有位置的目标个体和测试个体的距离差值总和
[0091] 4 · 5 · 6)若'_ > 响，则进入5);
[0092] 4.5.7)若，则产生一个(0，1)之间的随机数rand6，若rand6小于0.2， 贝1J用Ptraii替换Ptarget，否则进行5 );
[0093] 5)在每个温度层都进行变异、交叉和选择更新完种群后，将相邻的温度层的对应 个体进行副本交换：
[0094] 5.1)首先选择种群中第i个个体，ie I(i = l，2,…，popSize);
[0095] 5.2)选择两个相邻的温度层Tj和Tj+1，j = 1，2，…，7;
[0096] 5.3)产生一个随机数randxe (〇，1)，和判别数judge
[0097] (1)
[0098] 其中/C4,)(i = l，2, · · ·，p〇pSiZe)为种群个体4的能量值，K为玻尔兹曼常数， 1'心=1，2，~，7)为第」个温度层的开尔文温度，」为索引编号^为自然常数；
[0099] 5 · 4)若randx< judge，将相邻温度层中的'和4;个体进行交换；
[0100] 6)对每个温度层种群中的每个个体都执行完步骤4)-5)以后，迭代次数G = G+1，终 止条件为迭代次数G达到预设最大迭代次数iteration，判断是否满足终止条件，如不满足 则转至4);若满足则输出结果。
[0101]以上阐述的是本发明给出的一个实施例表现出来的优良效果，显然本发明不仅适 合上述实施例，在不偏离本发明基本精神及不超出本发明实质内容所涉及内容的前提下可 对其做种种变化加以实施。
【主权项】
1. 一种基于距离约束副本交换的蛋白质结构预测方法，其特征在于:所述预测方法包 括W下步骤： 1) 给定输入序列信息； 2) 根据序列信息从卵A服服务器上获取距离谱文件DP，rpk为距离谱中记录下的残基对， Dk为该残基对之间的距离，其中ke(l，N)，N为距离谱中残基对数量； 3) 初始化：设置群体规模popSize，交叉概率CR，8个溫度层T，能量函数选用Rosetta Score3,首先在每个溫度层通过对查询序列随机的折叠和变换，生成一个规模为popSize的 初始种群，初始群体为Ρτ= {χ? I i e 1}，计算能量函数值?·τ(χ?)，i e I，并设.片=min/CO，其 中i为种群个体编号，I为种群个体编号集合，1 = {1，2,. . .，popSize}，若为对应溫度层能量 函数值中的最小值，T为溫度层； 4) 开始迭代，依次对每个溫度层种群中的每个个体进行W下操作： 4.1) 设1 = 1，其中1￡{1，2,3，。'，肥};令口*3巧日* = 口1，其中1为序号，口*3巧日*表示目标个体； 4.2) 随机生成正整数randl ,rand2 ,rand3 e {1,2,3 ,......NP}，且randl声rand2声 rand3 声 i ;再生成4个随机整数ran化ange 1, randrange2 , ran化ange3 , randrange4;其中 r曰ndr曰ngel声r曰ndr曰nge2，r曰ndr曰nge3声r曰ndr曰nge4E (1,2,·.. .Length} .Length为序列长 度； 4.3 )针对个体Prandl做变异操作； 4.3.1) 令a = min(randrangel ,randrange2)，b = max(randrangel ,randrange2)，kE [曰，b]；令c=min(r曰ndr曰nge3,r曰ndr曰nge4)，d=m曰x(r曰ndr曰nge3,r曰ndr曰nge4)，pE k，d]; 其中min表示取两个数的最小值，max表示取两个数的最大值； 4.3.2) 用Prand2上位置a到位置b的片段的氨基酸所对应的二面角phi、psi、omega替换 Prandl的相同位置所对应的二面角地1、931、01116旨日;再使用?望(13上位置(3到位置(1的片段的氨 基酸所对应的二面角曲1、口3；[、0111日旨日替换?望(11上相同位置所对应的二面角地：[、93；[、01116邑日， 再将所得Prandl进行片段组装得到变异个体Pmutant; 4.4)针对变异个体Pmutant做交叉操作； 4.4.1) 生成随机数rand4，rands，其中rand4 e (0,1 )，rands e (1, Length); 4.4.2) 右随机数rand4〈 = CR,则使用目柄;个体Ptarget的片段rands替换变异个体Pmutant 中对应的片段生成测试个体Ptrail，否则Ptrail直接等于变异个体Pmutant; 4.5 )针对每个溫度层中的目标个体Ptarget和测试个体Ptrail进行选择操作； 4.5.1 )计算Ptarget和Ptrail的能里：E ( Ptarget)和E ( Ptrail); 4.5.2) 若E(Ptarget)〉E(Ptrail)则用Ptrail替换Ptarget，并且进入4.6);否则继续执行 4.5.3); 4.5.3) 计算目标个体Ptarget和测试个体Ptrail在残基对巧1的距离相对于距离谱中化之间 的距离差▽化巧et化和▽trailDi; 4.5.4) 判断Vtarget化和▽化ailDi的值，若大于6，则令其等于6，从而去除一些距离差 值较大的异常数值； 4.5.5) 计算距离谱中所有位置的目标个体和测试个体的距离差值总和4.5.6) 右〉7013,861:,则进入5); 4.5.7) 若▽Dtrail<VDtarget，则产生一个(0,1 )之间的随机数 Γ3Π (16，若 Γ3Π (16小于0.2， 贝ij用Ptrail替换Ptarget,否则进化5 ); 5) 在每个溫度层都进行变异、交叉和选择更新完种群后，将相邻的溫度层的对应个体 进行副本交换： 5.1) 首先选择种群中第i个个体，ie I(i = l，2,…，popSize); 5.2) 选择两个相邻的溫度层Tj和Tj+i，j = 1，2，…，7; 5.3) 产生一个随机数randx e (0，1)，和判别数化dge(1) 其中，/(.却）α = l，2，...，popSize)为种群个体馬的能量值，K为玻尔兹曼常数，TJ (j = l，2,…，7)为第j个溫度层的开尔文溫度，j为索引编号，e为自然常数； 5.4) 若randxC ^dge，将相邻溫度层中的持和為，个体进行交换； 6) 判断是否满足终止条件，如不满足则转至4);若满足则输出结果。2.如权利要求1所述的一种基于距离约束副本交换的蛋白质结构预测方法，其特征在 于:所述步骤3)中，设置最大迭代次数iteration;所述步骤6)中，对每个溫度层种群中的每 个个体都执行完步骤4)-5)?后，迭代次数G = G+1，终止条件为迭代次数G达到预设最大迭 代次数 iteration。
【文档编号】G06F19/16GK105975806SQ201610268630
【公开日】2016年9月28日
【申请日】2016年4月26日
【发明人】张贵军, 俞旭锋, 周晓根, 郝小虎, 王柳静, 张丹
【申请人】浙江工业大学