CN111538314B - 一种交流伺服***控制性能分析方法 - Google Patents
一种交流伺服***控制性能分析方法 Download PDFInfo
- Publication number
- CN111538314B CN111538314B CN202010336155.8A CN202010336155A CN111538314B CN 111538314 B CN111538314 B CN 111538314B CN 202010336155 A CN202010336155 A CN 202010336155A CN 111538314 B CN111538314 B CN 111538314B
- Authority
- CN
- China
- Prior art keywords
- value
- evaluation index
- servo system
- formula
- torque
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 27
- 238000011156 evaluation Methods 0.000 claims abstract description 206
- 230000004044 response Effects 0.000 claims abstract description 35
- 230000006870 function Effects 0.000 claims description 77
- 239000013598 vector Substances 0.000 claims description 64
- 238000000034 method Methods 0.000 claims description 48
- 239000011159 matrix material Substances 0.000 claims description 36
- 101100177269 Arabidopsis thaliana HCAR gene Proteins 0.000 claims description 15
- 238000010845 search algorithm Methods 0.000 claims description 13
- 230000001052 transient effect Effects 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000010606 normalization Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 238000005520 cutting process Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 description 9
- 238000004590 computer program Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000010219 correlation analysis Methods 0.000 description 4
- 239000000047 product Substances 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000306 component Substances 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
- G05B23/0243—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults model based detection method, e.g. first-principles knowledge model
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/20—Pc systems
- G05B2219/24—Pc safety
- G05B2219/24065—Real time diagnostics
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种用于交流伺服***性能分析方法,通过步骤:采集预设台数交流伺服***的控制性能的评价指标数据集,所述评价指标包括动态性能评价指标、稳态性能评价指标以及频域响应评价指标;根据预设的最优值取值条件,确定所述评价指标数据集中的最优值;获取所述评价指标的权重值;根据所述权重值,以及每台所述交流伺服***的评价指标数据与所述最优值之间的关联系数,分别计算所述与所述之间的关联度;根据所述关联度的大小,获取所述交流伺服***的评价结果;其中,所述关联度的值越大对应的所述交流伺服***的综合性能评价越高。
Description
技术领域
本发明涉及交流伺服***技术领域,尤其涉及一种用于交流伺服***控制性能分析的优化评价方法。
背景技术
交流伺服***作为工业生产机械的核心部件,是决定工业生产效率及产品品质的重要因素。随着科技的进步,交流伺服***逐渐向轻量化、高速化、精确化和智能化发展。高性能的交流伺服***在工业自动化领域的使用逐渐增多,为了将交流伺服***的控制性能调试到最优状态,就必须采取一种用于交流伺服***控制性能分析方法。但是,交流伺服***控制性能评价指标复杂、层次多样、关联性强、具有较强的非线性特征,并且在不同工作情况下对交流伺服***的控制性能要求也不同。目前,传统的评价方法中评价指标权重主观性较强,评价过程中定性成分较多,定量数据较少,需要一种科学、客观的方法对其进行改进。
发明内容
本发明的主要目的在于提供一种科学、客观的交流伺服***控制性能分析方法,解决现有技术中交流伺服***控制性能评价方法中存在的问题。
为实现上述目的,本发明提供的交流伺服***控制性能分析方法,包括步骤:
S1,采集预设台数Z台交流伺服***的控制性能的评价指标数据,得到所述Z台交流伺服***的评价指标数据集其中,i表示第i台交流伺服***,表示第i台交流伺服***的评价指标数据,3≤j≤M,ri,j分别表示第i台交流伺服***的j个评价指标的值,M表示所述评价指标的个数,所述评价指标包括动态性能评价指标(A)、稳态性能评价指标(B)以及频域响应评价指标(C),其中,所述动态性能评价指标选(A)选自以下指标中的一个或者多个:超调量、调整时间、正反转速差率以及转速调整率,所述稳态性能评价指标(B)选自以下指标中的一个或者多个:转矩平均误差、转矩波动系数以及转速波动系数,所述频域响应评价指标(C)选自以下指标中的一个或者多个:截止频域、幅值裕度以及相角裕度;
S3,获取所述评价指标的权重值WS(1×j)=[w1,w2,w3...,wj],其中,所述权重值WS(1×j)表示j个所述评价指标的权重值的集合,所述wj表示第j个所述评价指标的权重值;
优选地,所述步骤S1,包括以下步骤中的一个或者多个:
对被测交流伺服***突加20N·m的负载,记录所述被测交流伺服***的阶跃响应曲线,获取所述阶跃响应曲线中的所述超调量和所述调整时间;
对被测交流伺服***在额定电压220V下空载运行,将转速设定为1000r/min,测量所述被测交流伺服***的正反两个方向的转速平均值nccw和ncw,根据公式计算正反转速差率Kn,其中,Kn为正反转速差率,ncw为交流伺服***顺时针旋转时的转速平均值,单位为r/min,nccw为交流伺服***逆时针旋转时的转速平均值,单位为r/min;
对被测交流伺服***的转速设定为1000r/min,负载设定为10N·m,每隔10s将负载增大2N·m直到负载为30N·m,记录所述被测交流伺服***的转速变化值与额定转速nN的百分比,按公式计算得到所述被测交流伺服***转速调整率Δn,其中i表示改变负载的次数,ni表示所述被测交流伺服***的实际转速,单位为r/min,nN表示所述被测交流伺服***的额定转速,单位为r/min;
对被测交流伺服***施加20N·m的恒定负载稳定运行,记录所述被测交流伺服***实际产生的电磁转矩曲线,按照公式计算转矩平均误差其中为转矩平均误差,Ti为实际转矩大小,T*为设定转矩大小,N为采样次数;
对被测交流伺服***施加20N·m的恒定负载稳定运行,记录所述被测交流伺服***瞬时转矩的最大值Tmax,最小转矩值Tmin,按照公式进行转矩波动系数的计算KfT,其中KfT为转矩波动系数,Tmax为瞬态转矩的最大值,单位为N·m,Tmin为瞬态转矩的最小值,单位为N·m;
对被测交流伺服***施加20N·m的恒定负载稳定运行,记录所述被测交流伺服***瞬时转速的最大值nmax,最小转速值nmin,按照公式进行转速波动系数Kfn的计算,其中Kfn为转速波动系数,nmax为瞬态转速的最大值,单位为r/min,nmin为瞬态转速的最小值,单位为r/min;
优选地,所述步骤S4,包括:
按照公式Δij=|ri,j-r0,j|对归一化后的所述ri,j的值计算绝对差值矩阵Δ(N×M),其中,Δij为所述绝对差值矩阵第i行第j列元素的绝对值差值;
其中,z为所述ri,j的原始实际值,所述公式f1和f2的输出值为归一化后的所述ri,j的值;u,v分别为对应的所述评价指标的下界和上界,z∈[u,v]。
优选地,所述j的所述动态性能(A)的具体评价指标包括所述超调量、所述调整时间、所述正反转速差率以及所述转速调整率,所述稳态性能(B)的具体评价指标包括所述转矩平均误差、所述转矩波动系数以及所述转速波动系数,所述频域响应(C)的具体评价指标包括所述截止频域、所述幅值裕度以及所述相角裕度;
所述获取所述评价指标的权重值WS(1×j)=[w1,w2,w3...,wj]的步骤,包括:
获取所述动态性能评价指标(A)、稳态性能评价指标(B)、频域响应评价指标(C)中的评价指标权重值WA(1×4)、WB(1×3)以及WC(1×3);
优选地,所述获取动态性能评价指标(A)、稳态性能评价指标(B)、频域响应评价指标(C)中的评价指标权重值WA(1×4)、WB(1×3)以及WC(1×3)的步骤,包括:
根据模糊层次分析理论中五标度法则构造模糊判断矩阵MA(4×4)、MB(3×3)以及MC(3×3);
根据公式构建一致性指标CIF,其中,n为模糊判断矩阵MA(4×4)、MB(3×3)以及MC(3×3)的维度;wp满足公式中的约束条件;a=1,为两评价指标间重要程度差异的度量单位;bpq为模糊判断矩阵中的第p行第q列元素;ypq=b1q-bpq,为模糊判断矩阵的第一行第q列元素与第p行第q列元素之差,为ypq的平均值;
将所述构建一致性指标CIF的公式作为目标函数,用和声搜索算法求出一组(w1,w2,...wn,bpq)使得目标函数值小于0.1;
根据K组解向量将和声记忆库HM设置为一个K×H的二维矩阵,其中,H为(w1,w2,...wn,bpq)的中求解值的个数,xkh表示第k组解向量(w1,w2,...wn,bpq)的中对应的第h个求解值;
根据公式xkh=La+(Lb-La)×rand()分别初始化所述HM中的K组解向量,并将所述K组解向量分别代入所述CIF目标函数,得到K个目标函数值f1,f2,f3,...fk;其中,La=0.1为xkh变量的下界,Lb=0.9为变量xkh的上界,rand()为产生[0,1]之间随机数的函数;
获取预设的和声记忆库候选概率HMCR,调频宽度BW,基音调整概率PAR和即兴创作的次数N,判断所述rand()是否小于所述HMCR;
在所述rand()小于所述HMCR时,在所述HM内随机选择一组所述解向量作为新的解向量判断所述rand()是否小于所述PAR,若是,则对所述新的解向量以预设概率0.5进行扰动,并将扰动后的代入所述CIF目标函数得到的目标函数值fnew,若否,直接将所述新的解向量代入所述CIF目标函数得到的目标函数值fnew;
判断所述fnew是否小于所述fmin,若是,则将所述替换所述目标函数值f1,f2,f3,...fk中最大的目标函数值fmax对应的解向量,并产生新的解向量进行下次迭代,若否,则直接产生新的解向量进行下次迭代;
优选地,所述预设的和声记忆库候选概率HMCR=0.95,调频宽度BW=0.001,基音调整概率PAR=0.9,即兴创作的次数N=400-600,K=100。
在本发明的技术方案中,本发明提出一种用于交流伺服***性能分析方法可以从动态性能、稳态性能和频域响应方面对交流伺服***控制性能进行分析,所述评价指标包括动态性能评价指标(A)、稳态性能评价指标(B)以及频域响应评价指标(C),其中,所述动态性能(A)评价指标选自以下指标中的一个或者多个:超调量、调整时间、正反转速差率以及转速调整率,所述稳态性能评价指标(B)选自以下指标中的一个或者多个:转矩平均误差、转矩波动系数以及转速波动系数,所述频域响应评价指标(C)选自以下指标中的一个或者多个:截止频域、幅值裕度以及相角裕度。具体的,通过步骤S1,采集预设台数Z台交流伺服***的控制性能的评价指标数据,得到Z台交流伺服***的评价指标数据集步骤S2,根据预设的最优值取值条件,确定所述评价指标数据集中的最优值步骤S3,获取所述评价指标的权重值WS(1×j)=[w1,w2,w3...,wj];步骤S4,根据所述权重值,以及每台所述交流伺服***的评价指标数据与所述最优值之间的关联系数,分别计算所述与所述之间的关联度步骤S5,根据所述关联度的大小,获取所述交流伺服***的评价结果;其中,所述关联度的值越大对应的所述交流伺服***的综合性能评价越高。本发明中的用于交流伺服***性能分析方法结合具体的控制性能评价指标加权值以及每台所述交流伺服***的评价指标数据与所述最优值之间的关联度,通过关联度的大小来获取所述交流伺服***的评价结果,其中,所述关联度的值越大对应的所述交流伺服***的综合性能评价越高,为交流伺服***的方案设计、产品选型提供了参考方向。
附图说明
图1为本发明交流伺服***控制性能分析方法第一实施例的流程示意图;
图2为本发明交流伺服***控制性能分析方法的步骤S3的子流程示意图;
图3为本发明交流伺服***控制性能分析方法步骤S3的和声搜索算法的流程示意图;
图4为本发明交流伺服***控制性能分析方法的步骤S4的子流程示意图。
本发明目的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
另外,本发明各个实施例之间的技术方案可以相互结合,但是必须是以本领域普通技术人员能够实现为基础,当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在,也不在本发明要求的保护范围之内。
请参阅图1,为实现上述目的,本发明的第一实施例提供一种交流伺服***控制性能分析方法,本发明提供的交流伺服***控制性能分析方法,应用于交流伺服***转矩变化响应性能的评估,所述交流伺服***控制性能分析方法包括如下步骤:
所述评价指标包括动态性能评价指标(A)、稳态性能评价指标(B)以及频域响应评价指标(C),其中,所述动态性能评价指标(A)选自以下指标中的一个或者多个:超调量、调整时间、正反转速差率以及转速调整率,所述稳态性能评价指标(B)选自以下指标中的一个或者多个:转矩平均误差、转矩波动系数以及转速波动系数,所述频域响应评价指标(C)选自以下指标中的一个或者多个:截止频域、幅值裕度以及相角裕度。在本实施例中,以10个所述评价指标为例。
具体的,以上10个所述评价指标:超调量、调整时间、正反转速差率、转速调整率、转矩平均误差、转矩波动系数、转速波动系数、截止频域、幅值裕度和相角裕度的具体获取步骤如下:
对被测交流伺服***突加20N·m的负载,记录所述被测交流伺服***的阶跃响应曲线,获取所述阶跃响应曲线中的所述超调量和/或所述调整时间;
对被测交流伺服***在额定电压220V下空载运行,将转速设定为1000r/min,测量所述被测交流伺服***的正反两个方向的转速平均值nccw和ncw,根据公式计算正反转速差率Kn,其中,Kn为正反转速差率,ncw为交流伺服***顺时针旋转时的转速平均值,单位为r/min,nccw为交流伺服***逆时针旋转时的转速平均值,单位为r/min;
对被测交流伺服***的转速设定为1000r/min,负载设定为10N·m,每隔10s将负载增大2N·m直到负载为30N·m,记录所述被测交流伺服***的转速变化值与额定转速nN的百分比,按公式计算得到所述被测交流伺服***转速调整率Δn,其中i表示改变负载的次数,ni表示所述被测交流伺服***的实际转速,单位为r/min,nN表示所述被测交流伺服***的额定转速,单位为r/min;
对被测交流伺服***施加20N·m的恒定负载稳定运行,记录所述被测交流伺服***实际产生的电磁转矩曲线,按照公式计算转矩平均误差其中为转矩平均误差,Ti为实际转矩大小,T*为设定转矩大小,N为采样次数;
对被测交流伺服***施加20N·m的恒定负载稳定运行,记录所述被测交流伺服***瞬时转矩的最大值Tmax,最小转矩值Tmin,按照公式进行转矩波动系数的计算KfT,其中KfT为转矩波动系数,Tmax为瞬态转矩的最大值,单位为N·m,Tmin为瞬态转矩的最小值,单位为N·m;
对被测交流伺服***施加20N·m的恒定负载稳定运行,记录所述被测交流伺服***瞬时转速的最大值nmax,最小转速值nmin,按照公式进行转速波动系数Kfn的计算,其中Kfn为转速波动系数,nmax为瞬态转速的最大值,单位为r/min,nmin为瞬态转速的最小值,单位为r/min;
本领域技术人员应当知道,在JB/T 10184-200:《交流伺服驱动器通用技术条件》以及GB/T 16439-2009:《交流伺服***通用技术条件》中可以知道上述相关性能评价指标的相关内容,在此不再赘述。
在一较佳的实施例中,所述预设的最优值取值条件包括:将所述评价指标数据集中的所述动态性能评价指标(A)以及所述稳态性能评价指标(B)中的最小值作为最优值;将所述频域响应评价指标(C)中的最大值作为最优值。
通常来说,可以将控制性能评价指标分为正向指标和负向指标,正向指标中其数值越大对***的越有利,则将评价指标数据集中正向指标的最大值作为最优值,例如截止频域、相角裕度和幅值裕度;负向指标中其数值越小对***越有利,则将评价指标数据集中负向指标的最小值作为最优值,例如超调量、调整时间、正反转速差率、转速调整率、转矩平均误差、转矩波动系数、转速波动系数。
可以理解的是,在其他实施方式中,为了减小误差,可以去掉评价指标数据集中的一个或者几个最大值,以及一个或者几个最小值,然后根据预设的最优值取值条件进行选取最优值;以防止数据误测偏差带来的影响。当然,也可以在测试时采取多组测试取平均值的方式得到每台交流伺服***的评价指标数据,这样可以不用删除最大值和最小值的方式。
S3,获取所述评价指标的权重值WS(1×j)=[w1,w2,w3...,wj],其中,所述权重值WS(1×j)表示j个所述评价指标的权重值的集合,所述wj表示第j个所述评价指标的权重值。
请一并结合图2,在一具体示例中,步骤S3可以包括:
步骤S31,获取动态性能评价指标(A)、稳态性能评价指标(B)、频域响应评价指标(C)中的评价指标权重值WA(1×4)、WB(1×3)以及WC(1×3);
其中,所述获取动态性能(A)、稳态性能(B)、频域响应(C)中的评价指标权重值WA(1×4)、WB(1×3)以及WC(1×3)的步骤,包括:
根据模糊层次分析理论中五标度法则构造模糊判断矩阵MA(4×4)、MB(3×3)以及MC(3×3);
根据公式构建一致性指标CIF,其中,n为模糊判断矩阵MA(4×4)、MB(3×3)以及MC(3×3)的维度;wp满足公式中的约束条件;a=1,为两评价指标间重要程度差异的度量单位;bpq为模糊判断矩阵中的第p行第q列元素;ypq=b1q-bpq,为模糊判断矩阵的第一行第q列元素与第p行第q列元素之差,为ypq的平均值;
将所述构建一致性指标CIF的公式作为目标函数,用和声搜索算法求出一组(w1,w2,...wn,bpq)使得目标函数值小于0.1;
根据K组解向量将和声记忆库HM设置为一个K×H的二维矩阵,,其中,H为(w1,w2,...wn,bpq)的中求解值的个数,xkh表示第k组解向量(w1,w2,...wn,bpq)的中对应的第h个求解值;
根据公式xkh=La+(Lb-La)×rand()分别初始化所述HM中的K组解向量,并将所述K组解向量分别代入所述CIF目标函数,得到K个目标函数值f1,f2,f3,...fk;其中,La=0.1为xkh变量的下界,Lb=0.9为变量xkh的上界,rand()为产生[0,1]之间随机数的函数;
获取预设的和声记忆库候选概率HMCR,调频宽度BW,基音调整概率PAR和即兴创作的次数N,判断所述rand()是否小于所述HMCR;
在所述rand()小于所述HMCR时,在所述HM内随机选择一组所述解向量作为新的解向量判断所述rand()是否小于所述PAR,若是,则对所述新的解向量以预设概率0.5进行扰动,并将扰动后的代入所述CIF目标函数得到的目标函数值fnew,若否,直接将所述新的解向量代入所述CIF目标函数得到的目标函数值fnew;
判断所述fnew是否小于所述fmin,若是,则将所述替换所述目标函数值f1,f2,f3,...fk中最大的目标函数值fmax对应的解向量,并产生新的解向量进行下次迭代,若否,则直接产生新的解向量进行下次迭代;
具体的,以稳态性能评价指标(B)的评价指标权重值WB(1×3)的计算为例,来说明:
根据模糊层次分析理论中五标度法则构造模糊判断矩阵MB(3×3),将转矩平均误差(B1)、转矩波动系数(B2)、转速波动系数(B3)评价指标权重值分别用w1,w2,w3表示,并根据B1、B2、B3之间的联系,以及表1中的五标度规则构造模糊判断矩阵
MB(3×3),
M<sub>B</sub> | B1 | B2 | B3 |
B1 | 0.5 | 0.6 | 0.7 |
B2 | 0.4 | 0.5 | b<sub>23</sub> |
B3 | 0.3 | b<sub>32</sub> | 0.5 |
表1 0.1~0.9五标度
根据模糊判断矩阵MB(3×3)和
公式构建一致性指标CIF,其中,n=3为模糊判断矩阵MB(3×3)的维度;wp满足公式中的约束条件;a=1,为两评价指标间重要程度差异的度量单位;bpq为模糊判断矩阵中的第p行第q列元素;ypq=b1q-bpq,为模糊判断矩阵的第一行第q列元素与第p行第q列元素之差,为ypq的平均值。
将所述构建一致性指标CIF的公式作为目标函数,用和声搜索算法求出一组(w1,w2,...wn,bpq)使得目标函数值小于0.1;
根据100组解向量设置和声记忆库HM为一个100×4的二维矩阵,其中,(w1,w2,...wn,bpq)的中求解值的个数4,xkh表示第k组解向量(w1,w2,...wn,bpq)的中对应的第h个求解值;
根据公式xkh=La+(Lb-La)×rand()分别初始化所述HM中的100组解向量,并将所述100组解向量分别代入所述CIF目标函数,得到100个目标函数值f1,f2,f3,...fk;其中,La=0.1为xkh变量的下界,Lb=0.9为变量xkh的上界,rand()为产生[0,1]之间随机数的函数;
请一并结合图3,为本实施例中一具体示例中的和声搜索算法的流程图,包括步骤:
步骤S100,获取预设的和声记忆库候选概率HMCR=0.95,调频宽度BW=0.001,基音调整概率PAR=0.9,即兴创作的次数N=500;
步骤S101,判断所述rand()是否小于所述HMCR;
步骤S104,判断所述fnew是否小于所述fmin;若是,进入步骤S1041,若否,返回步骤S101直接产生新的解向量进行下次迭代;
步骤105,判断当前迭代次数是否达到即兴创作的次数N=500,若是,将所述目标函数值f1,f2,f3,...fk中最小的目标函数值fmin对应的解向量表示为所述若否,返回步骤S101直接产生新的解向量进行下次迭代。
进一步地,请一并结合图4,所述步骤S4,包括:
S42,按照公式Δij=|ri,j-r0,j|对归一化后的所述ri,j的值计算绝对差值矩阵Δ(N×M),其中,Δij为所述绝对差值矩阵第i行第j列元素的绝对值差值;
本发明一较佳示例中采用灰色关联度计算方式,灰色关联分析在灰色理论中是一个很重要的部分,国内外很多研究人员对构造灰色关联模型进行了长时间的探索,并取得了很多创新性的成就。根据白色***和黑色***,邓聚龙教授提出了一种解决信息不完备***的数学方法,即灰色关联分析法[31]。按照经典控制理论,颜色代表的是对于一个***已知信息的多少,白色就代表信息充足,比如传统的液压控制***,输入与输出之间的关系是能够确定的。然而,对于一个黑色***代表信息量完全缺失,对***的内部结构与元素间关系是完全不确定的***,通常称为这个叫做黑箱或者黑盒[54-57]。灰色***就是介于这两者之间,表示该***的信息量有部分缺失,对该***只有部分了解。机器人交流伺服电机的综合性能受到各种因素影响,因素包含的方面较多,而且不同因素间的耦合性较强。对于这些因素的评估通常是以统计学和运筹学为基础的方法进行描述,并不适用于复杂多变的场合。灰色关联分析法是一种多因素统计的分析方法,常用与多属性决策领域。它以各因素的样本数据为依据,用灰色关联度来描述各因素间的关系强弱、大小和次序。
在其他实施方式中,所述关联度的计算方式可以是其他实施方式,以上灰色关联度计算方式为本发明一较佳示例。还有基于Apriori算法、Hash算法等其他方式的关联度算法,在此不再赘述。
其中,z为所述ri,j的原始实际值,所述公式f1和f2的输出值为归一化后的所述ri,j的值;u,v分别为对应的所述评价指标的下界和上界,z∈[u,v]。
具体的,当控制性能评价指标数据越大对交流伺服***影响越有利时,将这种控制性能评价指标称为正向指标采用公式f1作为该指标的隶属度函数,如截止频率等性能指标。当控制性能评价指标数据越大对交流伺服***影响越有害时,将这种控制性能评价指标数据称为负向指标采用公式f2作为该指标的隶属度函数,如超调量等。然后,将控制性能评价指标数据通过对应的隶属度函数进行映射,将各个控制性能评价指标数据映射到[0,1]之间,统一评价指标的量纲,以便后续关联度计算。
具体的,通常可以定综合性能评估级别,例如设定机器人交流伺服电机综合性能级别为好、较好、中、较差、差五个级别,并设定综合性能级别对应的灰色关联度分段值为0.8,0.7,0.6,0.45。综合性能级别与灰色关联度之间的对应关系如表2所示。
表2综合性能分级
将机器人交流伺服电机灰色关联度与表2中的分段值进行比较,对机器人交流伺服电机的综合性能进行分级,确定机器人交流伺服电机的综合性能评级。
具体实例:
在一具体计算实例中,交流伺服***控制性能指标权重为:
w<sub>1</sub> | w<sub>2</sub> | w<sub>3</sub> | w<sub>4</sub> | w<sub>5</sub> |
0.046 | 0.046 | 0.051 | 0.086 | 0.046 |
w<sub>6</sub> | w<sub>7</sub> | w<sub>8</sub> | w<sub>9</sub> | w<sub>10</sub> |
0.109 | 0.072 | 0.128 | 0.221 | 0.192 |
选取I1,I2,…,I100中最大的交流伺服***控制性能评价指标数据,按照归一化处理后,具体见隶属度函数f1和f2,得到控制性能评价指标隶属度;控制性能评价指标的加权值为隶属度乘以对应的权重值;分段值为对比参照的理论区间。
可以看到选取I1,I2,…,I100中最大的交流伺服***控制性能较为优秀。具体表现为在稳态性能评价指标加权值上远远超过分段值。但是在超调量、转速调整率、相角裕度等评价指标加权值上低于分段值。反映出交流伺服***控制性能中稳态性能较优秀,动态性能与频域响应较差。
在本发明的技术方案中,本发明提出一种用于交流伺服***性能分析方法可以从动态性能、稳态性能和频域响应方面对交流伺服***控制性能进行分析,所述评价指标包括动态性能评价指标(A)、稳态性能评价指标(B)以及频域响应评价指标(C),其中,所述动态性能(A)评价指标选自以下指标中的一个或者多个:超调量、调整时间、正反转速差率以及转速调整率,所述稳态性能评价指标(B)选自以下指标中的一个或者多个:转矩平均误差、转矩波动系数以及转速波动系数,所述频域响应评价指标(C)选自以下指标中的一个或者多个:截止频域、幅值裕度以及相角裕度。具体的,通过步骤S1,采集预设台数Z台交流伺服***的控制性能的评价指标数据,得到Z台交流伺服***的评价指标数据集步骤S2,根据预设的最优值取值条件,确定所述评价指标数据集中的最优值步骤S3,S3,获取所述评价指标的权重值WS(1×j)=[w1,w2,w3...,wj];步骤S4,根据所述权重值,以及每台所述交流伺服***的评价指标数据与所述最优值之间的关联系数,分别计算所述与所述之间的关联度步骤S5,根据所述关联度的大小,获取所述交流伺服***的评价结果;其中,所述关联度的值越大对应的所述交流伺服***的综合性能评价越高。本发明中的用于交流伺服***性能分析方法结合具体的控制性能评价指标加权值以及每台所述交流伺服***的评价指标数据与所述最优值之间的关联度,通过关联度的大小来获取所述交流伺服***的评价结果,其中,所述关联度的值越大对应的所述交流伺服***的综合性能评价越高,为交流伺服***的方案设计、产品选型提供了参考方向。
针对机器人交流伺服电机综合性能评估中指标数据信息量缺失导致其综合性能难以定量评估的问题,本章采用灰色关联分析法对机器交流伺服电机的综合性能评估展开研究。首先,建立机器人交流伺服电机综合性能评估指标数据的参考序列和比较序列。其次,结合参考序列与比较序列建立关联系数矩阵。然后,根据指标权重与关联系数矩阵计算机器人交流伺服电机的关联度大小,并根据关联度对机器人交流伺服电机综合性能进行定量评估,确定机器人交流伺服电机的综合性能评级。最后,结合机器人交流伺服电机综合性能评估指标数据集对所提出的评估方法进行仿真实验,并将基于灰色关联度的机器人交流伺服电机的综合性能评估结果与第三章中的基于模糊层次分析法的机器人交流伺服电机综合性能评估结果进行对比,实验结果表明,两种方法的评估结果一致,能够有效的对机器人交流伺服电机的综合性能进行评估。
本发明还提供一种交流伺服***控制性能分析装置,所述交流伺服***控制性能分析装置包括存储在所述存储器上并可在所述处理器上运行的操作***以及计算机程序;所述计算机程序被所述处理器执行时实现如上中任一项所述的交流伺服***控制性能分析方法的步骤。
可以理解的是,由于本实施例中的包含交流伺服***控制性能分析方法的计算机程序至少包括上述交流伺服***控制性能分析方法实施例的全部技术方案,同时至少具有以上实施例的全部技术效果,此处不再一一赘述。
此外,为实现上述目的,本发明还提供一种计算机可读存储介质,所述计算机可读存储介质上存储有计算机程序,所述计算机程序被处理器执行时实现如上述任一实施例交流伺服***控制性能分析方法的的步骤。
由于本实施例计算机可读存储介质的技术方案至少包括上述交流伺服***控制性能分析方法实施例的全部技术方案,因此至少具有以上实施例的全部技术效果,此处不再一一赘述。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在如上所述的一个计算机可读存储介质(如ROM/RAM、磁碟、光盘)中,包括若干指令用以使得一台终端设备进入本发明各个实施例所述的方法。
需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者***不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者***所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括该要素的过程、方法、物品或者***中还存在另外的相同要素。
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
以上仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (10)
1.一种交流伺服***控制性能分析方法,其特征在于,包括步骤:
S1,采集预设台数Z台交流伺服***的控制性能的评价指标数据,得到所述Z台交流伺服***的评价指标数据集其中,i表示第i台交流伺服***,表示第i台交流伺服***的评价指标数据,3≤j≤M,ri,j分别表示第i台交流伺服***的j个评价指标的值,M表示所述评价指标的个数,所述评价指标包括动态性能评价指标(A)、稳态性能评价指标(B)以及频域响应评价指标(C),其中,所述动态性能评价指标选(A)选自以下指标中的一个或者多个:超调量、调整时间、正反转速差率以及转速调整率,所述稳态性能评价指标(B)选自以下指标中的一个或者多个:转矩平均误差、转矩波动系数以及转速波动系数,所述频域响应评价指标(C)选自以下指标中的一个或者多个:截止频域、幅值裕度以及相角裕度;
S3,获取所述评价指标的权重值WS(1×j)=[w1,w2,w3...,wj],其中,所述权重值WS(1×j)表示j个所述评价指标的权重值的集合,所述wj表示第j个所述评价指标的权重值;
2.根据权利要求1中所述的交流伺服***控制性能分析方法,其特征在于,所述步骤S1,包括以下步骤中的一个或者多个:
对被测交流伺服***突加20N·m的负载,记录所述被测交流伺服***的阶跃响应曲线,获取所述阶跃响应曲线中的所述超调量和所述调整时间;
对被测交流伺服***在额定电压220V下空载运行,将转速设定为1000r/min,测量所述被测交流伺服***的正反两个方向的转速平均值nccw和ncw,根据公式计算正反转速差率Kn,其中,Kn为正反转速差率,ncw为交流伺服***顺时针旋转时的转速平均值,单位为r/min,nccw为交流伺服***逆时针旋转时的转速平均值,单位为r/min;
对被测交流伺服***的转速设定为1000r/min,负载设定为10N·m,每隔10s将负载增大2N·m直到负载为30N·m,记录所述被测交流伺服***的转速变化值与额定转速nN的百分比,按公式计算得到所述被测交流伺服***转速调整率Δn,其中i表示改变负载的次数,ni表示所述被测交流伺服***的实际转速,单位为r/min,nN表示所述被测交流伺服***的额定转速,单位为r/min;
对被测交流伺服***施加20N·m的恒定负载稳定运行,记录所述被测交流伺服***实际产生的电磁转矩曲线,按照公式计算转矩平均误差其中为转矩平均误差,Ti为实际转矩大小,T*为设定转矩大小,N为采样次数;
对被测交流伺服***施加20N·m的恒定负载稳定运行,记录所述被测交流伺服***瞬时转矩的最大值Tmax,最小转矩值Tmin,按照公式进行转矩波动系数的计算KfT,其中KfT为转矩波动系数,Tmax为瞬态转矩的最大值,单位为N·m,Tmin为瞬态转矩的最小值,单位为N·m;
6.根据权利要求1中所述的交流伺服***控制性能分析方法,其特征在于,所述j的所述动态性能(A)的具体评价指标包括所述超调量、所述调整时间、所述正反转速差率以及所述转速调整率,所述稳态性能(B)的具体评价指标包括所述转矩平均误差、所述转矩波动系数以及所述转速波动系数,所述频域响应(C)的具体评价指标包括所述截止频域、所述幅值裕度以及所述相角裕度;
所述获取所述评价指标的权重值WS(1×j)=[w1,w2,w3...,wj]的步骤,包括:
获取所述动态性能评价指标(A)、稳态性能评价指标(B)、频域响应评价指标(C)中的评价指标权重值WA(1×4)、WB(1×3)以及WC(1×3);
7.根据权利要求6中所述的交流伺服***控制性能分析方法,其特征在于,所述获取动态性能评价指标(A)、稳态性能评价指标(B)、频域响应评价指标(C)中的评价指标权重值WA(1×4)、WB(1×3)以及WC(1×3)的步骤,包括:
根据模糊层次分析理论中五标度法则构造模糊判断矩阵MA(4×4)、MB(3×3)以及MC(3×3);
根据公式构建一致性指标CIF,其中,n为模糊判断矩阵MA(4×4)、MB(3×3)以及MC(3×3)的维度;wp满足公式中的约束条件;a=1,为两评价指标间重要程度差异的度量单位;bpq为模糊判断矩阵中的第p行第q列元素;ypq=b1q-bpq,为模糊判断矩阵的第一行第q列元素与第p行第q列元素之差,为ypq的平均值;
将所述构建一致性指标CIF的公式作为目标函数,用和声搜索算法求出一组(w1,w2,...wn,bpq)使得目标函数值小于0.1;
根据K组解向量将和声记忆库HM设置为一个K×H的二维矩阵,其中,H为(w1,w2,...wn,bpq)的中求解值的个数,xkh表示第k组解向量(w1,w2,...wn,bpq)的中对应的第h个求解值;
根据公式xkh=La+(Lb-La)×rand()分别初始化所述HM中的K组解向量,并将所述K组解向量分别代入所述CIF目标函数,得到K个目标函数值f1,f2,f3,...fk;其中,La=0.1为xkh变量的下界,Lb=0.9为变量xkh的上界,rand()为产生[0,1]之间随机数的函数;
获取预设的和声记忆库候选概率HMCR,调频宽度BW,基音调整概率PAR和即兴创作的次数N,判断所述rand()是否小于所述HMCR;
在所述rand()小于所述HMCR时,在所述HM内随机选择一组所述解向量作为新的解向量判断所述rand()是否小于所述PAR,若是,则对所述新的解向量以预设概率0.5进行扰动,并将扰动后的代入所述CIF目标函数得到的目标函数值fnew,若否,直接将所述新的解向量代入所述CIF 目标函数得到的目标函数值fnew;
判断所述fnew是否小于所述fmin,若是,则将所述替换所述目标函数值f1,f2,f3,...fk中最大的目标函数值fmax对应的解向量,并产生新的解向量进行下次迭代,若否,则直接产生新的解向量进行下次迭代;
10.根据权利要求9中所述的交流伺服***控制性能分析方法,其特征在于,所述预设的和声记忆库候选概率HMCR=0.95,调频宽度BW=0.001,基音调整概率PAR=0.9,即兴创作的次数N=400-600,K=100。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010336155.8A CN111538314B (zh) | 2020-04-26 | 2020-04-26 | 一种交流伺服***控制性能分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010336155.8A CN111538314B (zh) | 2020-04-26 | 2020-04-26 | 一种交流伺服***控制性能分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111538314A CN111538314A (zh) | 2020-08-14 |
CN111538314B true CN111538314B (zh) | 2021-06-18 |
Family
ID=71973235
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010336155.8A Active CN111538314B (zh) | 2020-04-26 | 2020-04-26 | 一种交流伺服***控制性能分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111538314B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2022065759A (ja) * | 2020-10-16 | 2022-04-28 | セイコーエプソン株式会社 | ロボットのパラメーターセットを調整する方法、プログラム、および情報処理装置 |
CN112711237B (zh) * | 2020-12-29 | 2022-09-27 | 华润电力技术研究院有限公司 | 一种火电机组自动控制品质在线评估方法和*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101989827A (zh) * | 2010-11-18 | 2011-03-23 | 东南大学 | 基于惯量辨识的交流伺服***速度环控制参数自整定方法 |
CN102497156A (zh) * | 2011-12-27 | 2012-06-13 | 东南大学 | 永磁同步电机速度环的神经网络自校正控制方法 |
CN103338003A (zh) * | 2013-06-28 | 2013-10-02 | 西安交通大学 | 一种电机负载转矩及惯量同时在线辨识的方法 |
CN104749956A (zh) * | 2015-03-24 | 2015-07-01 | 北京理工大学 | 基于和声搜索算法的工业机器人结构优化方法 |
CN107390512A (zh) * | 2017-07-27 | 2017-11-24 | 中国船舶重工集团公司第七二四研究所 | 一种冗余驱动伺服***的抗干扰无间断容错控制方法 |
-
2020
- 2020-04-26 CN CN202010336155.8A patent/CN111538314B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101989827A (zh) * | 2010-11-18 | 2011-03-23 | 东南大学 | 基于惯量辨识的交流伺服***速度环控制参数自整定方法 |
CN102497156A (zh) * | 2011-12-27 | 2012-06-13 | 东南大学 | 永磁同步电机速度环的神经网络自校正控制方法 |
CN103338003A (zh) * | 2013-06-28 | 2013-10-02 | 西安交通大学 | 一种电机负载转矩及惯量同时在线辨识的方法 |
CN104749956A (zh) * | 2015-03-24 | 2015-07-01 | 北京理工大学 | 基于和声搜索算法的工业机器人结构优化方法 |
CN107390512A (zh) * | 2017-07-27 | 2017-11-24 | 中国船舶重工集团公司第七二四研究所 | 一种冗余驱动伺服***的抗干扰无间断容错控制方法 |
Non-Patent Citations (1)
Title |
---|
并行混沌与和声搜索的多目标混合优化算法;袁小芳等;《湖南大学学报(自然科学版)》;20201231;第45卷(第4期);第96页至第103页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111538314A (zh) | 2020-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111538314B (zh) | 一种交流伺服***控制性能分析方法 | |
CN112785119B (zh) | 基于聚类和层次分析算法的配网电压越限原因分析方法 | |
CN110365053B (zh) | 基于延迟优化策略的短期风电功率预测方法 | |
CN113240261A (zh) | 一种区域电能质量监测分析*** | |
CN111832839B (zh) | 基于充分增量学习的能耗预测方法 | |
CN112732786A (zh) | 金融数据处理方法、装置、设备及存储介质 | |
Kumar et al. | Comparative analysis of SOM neural network with K-means clustering algorithm | |
CN112712240A (zh) | 一种台区线损成因分析方法及装置 | |
CN111598313A (zh) | 一种风功率预测方法 | |
CN115496249A (zh) | 一种基于聚类算法的工业可调负荷潜力分析方法及*** | |
Liu et al. | Evaluation on regional science and technology resources allocation in China based on the zero sum gains data envelopment analysis | |
CN110796303A (zh) | 一种基于ewt和odbscan的短期电力负荷预测方法 | |
CN113591322A (zh) | 一种基于极端梯度提升决策树的低压台区线损率预测方法 | |
CN107808245A (zh) | 基于改进决策树方法的管网调度*** | |
CN109580918B (zh) | 一种用于预测石脑油的分子组成的方法 | |
CN115936303A (zh) | 基于机器学习模型进行暂态电压安全分析方法 | |
CN111539275B (zh) | 一种基于负载特性可视化的电气负载分类方法和*** | |
CN114004525A (zh) | 基于模糊综合评价法面向用户侧的电能替代效益评价方法 | |
CN113011779A (zh) | 一种基于模糊综合评价的能源消费价格补偿方法及装置 | |
Klein-Helmkamp et al. | Cluster-Based Analysis of Feasible Operating Regions for Aggregated Energy Storage Systems | |
CN114970719B (zh) | 一种基于改进svr模型的物联运行指标预测方法 | |
Zheng et al. | Construction of Diagnosis and Evaluation Model of “Total Factor” Resource Input | |
CN115833135B (zh) | 一种机会约束交流最优潮流自适应求解方法、装置及介质 | |
CN117421882A (zh) | 一种基于量测信息的综合负荷模型在线建模方法 | |
CN117893092A (zh) | 一种基于改进熵权法的用户能效分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |