CN116595730A - 一种电磁暂态仿真误差评价处理方法及装置 - Google Patents
一种电磁暂态仿真误差评价处理方法及装置 Download PDFInfo
- Publication number
- CN116595730A CN116595730A CN202310489287.8A CN202310489287A CN116595730A CN 116595730 A CN116595730 A CN 116595730A CN 202310489287 A CN202310489287 A CN 202310489287A CN 116595730 A CN116595730 A CN 116595730A
- Authority
- CN
- China
- Prior art keywords
- calculating
- difference
- normal distribution
- error index
- distribution function
- 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.)
- Pending
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 59
- 238000004088 simulation Methods 0.000 title claims abstract description 53
- 230000001052 transient effect Effects 0.000 title claims abstract description 48
- 238000003672 processing method Methods 0.000 title claims abstract description 17
- 238000005315 distribution function Methods 0.000 claims abstract description 66
- 238000000034 method Methods 0.000 claims abstract description 52
- 238000012545 processing Methods 0.000 claims abstract description 31
- 230000009466 transformation Effects 0.000 claims abstract description 18
- 238000004590 computer program Methods 0.000 claims description 23
- 238000003860 storage Methods 0.000 claims description 8
- 238000010586 diagram Methods 0.000 description 11
- 230000006870 function Effects 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000011158 quantitative evaluation Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000006872 improvement Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013507 mapping 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
- 238000010206 sensitivity analysis Methods 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Biology (AREA)
- Computer Hardware Design (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Computation (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种电磁暂态仿真误差评价处理方法及装置,涉及电力系统仿真技术领域。所述方法包括:获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果。所述装置执行上述方法。本发明实施例提供的方法及装置,能够提高方法的通用性。
Description
技术领域
本发明涉及电力系统仿真技术领域,具体涉及一种电磁暂态仿真误差评价处理方法及装置。
背景技术
电磁暂态仿真在时域下求解系统电路,常用于电力系统电磁暂态现象研究、控制保护算法研究等。随着电磁暂态仿真逐渐成为评估电力系统动态行为的重要手段,需要首先保证仿真的有效性,如何评价仿真误差,成为保证仿真准确性的前提条件。
仿真误差评价方法可分为定性评价与定量评价。定性评价指由相关专家与技术人员根据经验采用目测比对的方式作出主观判断,这种方法要求判断者对模型原理、测试条件和数据可信度理解全面,判断结果带有明显的主观偏向,缺乏稳定性。定量评价指根据仿真的数据信息,借助计算机定量计算误差评价指标,评价结果具有客观性与可信性,能够为模型优化与改进提供定量化的建议。
定量评价主要使用数值分析方法与特征选择验证法(Feature SelectiveValidation,FSV)。数值分析方法主要包括均方误差分析、灵敏度分析法、相似性系数法、假设检验法等,这类方法体现的误差信息相对单一,缺少对模型进行改进的指导信息,无法对误差细节进行溯源分析。为了使误差评价更加全面,FSV方法应运而生。FSV方法根据用户关心的数据特征进行差异提取,能够全面反映数据在整体及细节上的差异,同时结合专家目测评价过程,将定量的误差评价结果定性表述,符合实际工程使用需要,但是,用户无法直观方便的选择关心的数据,因此上述方法的通用性欠佳。
发明内容
针对现有技术中的问题,本发明实施例提供一种电磁暂态仿真误差评价处理方法及装置,能够至少部分地解决现有技术中存在的问题。
一方面,本发明提出一种电磁暂态仿真误差评价处理方法,包括:
获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;
根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;
根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果。
其中,所述根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,包括:
根据所述正态分布函数和所述幅值差异量计算第一总体误差指标;
根据所述正态分布函数和所述特征差异量计算第二总体误差指标;
根据所述正态分布函数和所述全局差异量计算第三总体误差指标;
将所述第一总体误差指标、所述第二总体误差指标和所述第三总体误差指标确定为所述总体误差指标。
其中,所述根据所述正态分布函数和所述幅值差异量计算第一总体误差指标,包括:
根据如下公式计算所述第一总体误差指标:
其中,ADMtot为第一总体误差指标、ADM(i)为幅值差异量、f(i)为正态分布函数、i为第i个频点、N为频点个数。
其中,所述根据所述正态分布函数和所述特征差异量计算第二总体误差指标,包括:
根据如下公式计算所述第二总体误差指标:
其中,FDMtot为第二总体误差指标、FDM(i)为特征差异量、f(i)为正态分布函数、i为第i个频点、N为频点个数。
其中,所述根据所述正态分布函数和所述全局差异量计算第三总体误差指标,包括:
根据如下公式计算所述第三总体误差指标:
其中,GDMtot为第三总体误差指标、GDM(i)为全局差异量、f(i)为正态分布函数、i为第i个频点、N为频点个数。
其中,所述根据所述总体误差指标确定误差评价结果,包括:
根据所述总体误差指标的数值所在的预设数值区间确定误差评价结果。
一方面,本发明提出一种电磁暂态仿真误差评价处理装置,包括:
获取单元,用于获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;
计算单元,用于根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;
确定单元,用于根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果。
再一方面,本发明实施例提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如下方法:
获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;
根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;
根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果。
本发明实施例提供一种计算机可读存储介质,包括:
所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现如下方法:
获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;
根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;
根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果。
本发明实施例还提供一种计算机程序产品,所述计算机程序产品包括计算机程序,所述计算机程序被处理器执行时实现如下方法:
获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;
根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;
根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果。
本发明实施例提供的电磁暂态仿真误差评价处理方法及装置,获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果,不但能够全面反映数据在整体及细节上的差异,方便进行误差细节溯源,而且使得用户能够直观方便的选择关心的数据,提高方法的通用性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1是本发明一实施例提供的电磁暂态仿真误差评价处理方法的流程示意图。
图2是本发明实施例提供的FSV方法基本原理图。
图3是本发明另一实施例提供的电磁暂态仿真误差评价处理方法的流程示意图。
图4是本发明一实施例提供的电磁暂态仿真误差评价处理装置的结构示意图。
图5为本发明实施例提供的计算机设备实体结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互任意组合。
图1是本发明一实施例提供的电磁暂态仿真误差评价处理方法的流程示意图,如图1所示,本发明实施例提供的电磁暂态仿真误差评价处理方法,包括:
步骤S1:获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量。
步骤S2:根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量。
步骤S3:根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果。
在上述步骤S1中,装置获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量。装置可以是执行该方法的计算机设备等,例如为服务器。本申请技术方案中对数据的获取、存储、使用、处理等均符合国家法律法规的相关规定。
如图2所示,获取电磁暂态仿真的待分析数据,具体包括:选取待分析的两组数据,其中一组作为待评估数据,另一组作为可信参考数据。
如图3所示,对上述两组数据分别进行傅里叶变换分解,得到直流分量DC、低频分量Lo和高频分量Hi。三组分量的分界点采用IEEE标准中的定义方式,频率最低的4个频点定义为直流分量DC,除去直流分量后,将剩余数据总体能量40%作为高低频率的分界点。通过公式说明如下:
ip=i40%+5
式中,S为除直流分量外,剩余频点的频谱值之和;N为频点个数;F(i)为第i个频点的频谱值;i40%为除直流分量外,从第5个频点累加到S的40%所对应的频点;ip为高低频率分量的分界点,除直流分量外,频率在ip以下的定义为低频分量,频率在ip以上的定义为高频分量。
在上述步骤S2中,装置根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量。
说明如下:
幅值差异量ADM反映数据在总体数值上的差异性,通过直流分量和低频分量计算得出,计算公式为:
ODM(i)通过如下公式表示:
其中,ADM(i)为第i个频点的幅值差异量,ODM(i)为第i个频点的直流偏移差异量,未提及的内容为本领域常规内容。
特征差异量FDM反映数据细节特征的差异性,通过低频分量和高频分量的导数计算得出,计算公式为:
FDM(i)=2(|FDM1(i)+FDM2(i)+FDM3(i)|)
其中,未提及的内容为本领域常规内容。
计算全局差异量GDM,具体方法如下:
GDM为全局差异量,是平衡整体差异与细节差异之后的误差评价结果,计算公式为:
kA×kF=1
式中,kA表示ADM所占权重,kF表示FDM所占权重,本实施例中二者都可以取为1。
在上述步骤S3中,装置根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果。
所述根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,包括:
根据所述正态分布函数和所述幅值差异量计算第一总体误差指标;
根据所述正态分布函数和所述特征差异量计算第二总体误差指标;
根据所述正态分布函数和所述全局差异量计算第三总体误差指标;
将所述第一总体误差指标、所述第二总体误差指标和所述第三总体误差指标确定为所述总体误差指标。
所述根据所述正态分布函数和所述幅值差异量计算第一总体误差指标,包括:
根据如下公式计算所述第一总体误差指标:
其中,ADMtot为第一总体误差指标、ADM(i)为幅值差异量、f(i)为正态分布函数、i为第i个频点、N为频点个数。
所述根据所述正态分布函数和所述特征差异量计算第二总体误差指标,包括:
根据如下公式计算所述第二总体误差指标:
其中,FDMtot为第二总体误差指标、FDM(i)为特征差异量、f(i)为正态分布函数、i为第i个频点、N为频点个数。
所述根据所述正态分布函数和所述全局差异量计算第三总体误差指标,包括:
根据如下公式计算所述第三总体误差指标:
其中,GDMtot为第三总体误差指标、GDM(i)为全局差异量、f(i)为正态分布函数、i为第i个频点、N为频点个数。
正态分布是一种的连续概率分布,其函数曲线为钟形曲线,若随机变量服从一个均值参数为μ、方差参数为σ2的概率分布,且其概率密度函数满足公式如下:
根据不同研究关注频点的不同,基于正态分布函数对所有频点误差加权求平均值,对于重点关注频点误差施加较大系数,对于边缘频点误差施加较小系数,从而得到数据总体误差水平。
式中,f(i)是均值为μ,方差可选为1的正态分布概率密度函数,μ的取值与待分析的关注频点相同,保证在关注频点处误差权重较大,放大关注频点及周围频点的误差特征。
所述根据所述总体误差指标确定误差评价结果,包括:
根据所述总体误差指标的数值所在的预设数值区间确定误差评价结果。
按照表1的映射关系将定量结果使用定性描述的方式表现,直观反映数据误差的判断结果。
表1
本发明实施例提供的电磁暂态仿真误差评价处理方法,获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果,不但能够全面反映数据在整体及细节上的差异,方便进行误差细节溯源,而且使得用户能够直观方便的选择关心的数据,提高方法的通用性。
进一步地,所述根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,包括:
根据所述正态分布函数和所述幅值差异量计算第一总体误差指标;可参照上述说明,不再赘述。
根据所述正态分布函数和所述特征差异量计算第二总体误差指标;可参照上述说明,不再赘述。
根据所述正态分布函数和所述全局差异量计算第三总体误差指标;可参照上述说明,不再赘述。
将所述第一总体误差指标、所述第二总体误差指标和所述第三总体误差指标确定为所述总体误差指标。可参照上述说明,不再赘述。
进一步地,所述根据所述正态分布函数和所述幅值差异量计算第一总体误差指标,包括:
根据如下公式计算所述第一总体误差指标:
其中,ADMtot为第一总体误差指标、ADM(i)为幅值差异量、f(i)为正态分布函数、i为第i个频点、N为频点个数。可参照上述说明,不再赘述。
进一步地,所述根据所述正态分布函数和所述特征差异量计算第二总体误差指标,包括:
根据如下公式计算所述第二总体误差指标:
其中,FDMtot为第二总体误差指标、FDM(i)为特征差异量、f(i)为正态分布函数、i为第i个频点、N为频点个数。可参照上述说明,不再赘述。
进一步地,所述根据所述正态分布函数和所述全局差异量计算第三总体误差指标,包括:
根据如下公式计算所述第三总体误差指标:
其中,GDMtot为第三总体误差指标、GDM(i)为全局差异量、f(i)为正态分布函数、i为第i个频点、N为频点个数。可参照上述说明,不再赘述。
进一步地,所述根据所述总体误差指标确定误差评价结果,包括:
根据所述总体误差指标的数值所在的预设数值区间确定误差评价结果。可参照上述说明,不再赘述。
图4是本发明一实施例提供的电磁暂态仿真误差评价处理装置的结构示意图,如图4所示,本发明实施例提供的电磁暂态仿真误差评价处理装置,包括获取单元401、计算单元402和确定单元403,其中:
获取单元401用于获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;计算单元402用于根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;确定单元403用于根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果。
具体的,装置中的获取单元401用于获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;计算单元402用于根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;确定单元403用于根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果。
本发明实施例提供的电磁暂态仿真误差评价处理装置,获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果,不但能够全面反映数据在整体及细节上的差异,方便进行误差细节溯源,而且使得用户能够直观方便的选择关心的数据,提高方法的通用性。
进一步地,所述确定单元403具体用于:
根据所述正态分布函数和所述幅值差异量计算第一总体误差指标;
根据所述正态分布函数和所述特征差异量计算第二总体误差指标;
根据所述正态分布函数和所述全局差异量计算第三总体误差指标;
将所述第一总体误差指标、所述第二总体误差指标和所述第三总体误差指标确定为所述总体误差指标。
进一步地,所述确定单元403还具体用于:
根据如下公式计算所述第一总体误差指标:
其中,ADMtot为第一总体误差指标、ADM(i)为幅值差异量、f(i)为正态分布函数、i为第i个频点、N为频点个数。
进一步地,所述确定单元403还具体用于:
根据如下公式计算所述第二总体误差指标:
其中,FDMtot为第二总体误差指标、FDM(i)为特征差异量、f(i)为正态分布函数、i为第i个频点、N为频点个数。
进一步地,所述确定单元403还具体用于:
根据如下公式计算所述第三总体误差指标:
其中,GDMtot为第三总体误差指标、GDM(i)为全局差异量、f(i)为正态分布函数、i为第i个频点、N为频点个数。
进一步地,所述确定单元403具体用于:
根据所述总体误差指标的数值所在的预设数值区间确定误差评价结果。
本发明实施例提供电磁暂态仿真误差评价处理装置的实施例具体可以用于执行上述各方法实施例的处理流程,其功能在此不再赘述,可以参照上述方法实施例的详细描述。
图5为本发明实施例提供的计算机设备实体结构示意图,如图5所示,所述计算机设备包括:存储器501、处理器502及存储在存储器501上并可在处理器502上运行的计算机程序,所述处理器502执行所述计算机程序时实现如下方法:
获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;
根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;
根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果。
本实施例公开一种计算机程序产品,所述计算机程序产品包括计算机程序,所述计算机程序被处理器执行时实现如下方法:
获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;
根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;
根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果。
本实施例提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现如下方法:
获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;
根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;
根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果。
本发明实施例与现有技术中的技术方案相比,获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果,不但能够全面反映数据在整体及细节上的差异,方便进行误差细节溯源,而且使得用户能够直观方便的选择关心的数据,提高方法的通用性。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在本说明书的描述中,参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种电磁暂态仿真误差评价处理方法,其特征在于,包括:
获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;
根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;
根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果。
2.根据权利要求1所述的电磁暂态仿真误差评价处理方法,其特征在于,所述根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,包括:
根据所述正态分布函数和所述幅值差异量计算第一总体误差指标;
根据所述正态分布函数和所述特征差异量计算第二总体误差指标;
根据所述正态分布函数和所述全局差异量计算第三总体误差指标;
将所述第一总体误差指标、所述第二总体误差指标和所述第三总体误差指标确定为所述总体误差指标。
3.根据权利要求2所述的电磁暂态仿真误差评价处理方法,其特征在于,所述根据所述正态分布函数和所述幅值差异量计算第一总体误差指标,包括:
根据如下公式计算所述第一总体误差指标:
其中,ADMtot为第一总体误差指标、ADM(i)为幅值差异量、f(i)为正态分布函数、i为第i个频点、N为频点个数。
4.根据权利要求2所述的电磁暂态仿真误差评价处理方法,其特征在于,所述根据所述正态分布函数和所述特征差异量计算第二总体误差指标,包括:
根据如下公式计算所述第二总体误差指标:
其中,FDMtot为第二总体误差指标、FDM(i)为特征差异量、f(i)为正态分布函数、i为第i个频点、N为频点个数。
5.根据权利要求2所述的电磁暂态仿真误差评价处理方法,其特征在于,所述根据所述正态分布函数和所述全局差异量计算第三总体误差指标,包括:
根据如下公式计算所述第三总体误差指标:
其中,GDMtot为第三总体误差指标、GDM(i)为全局差异量、f(i)为正态分布函数、i为第i个频点、N为频点个数。
6.根据权利要求1至5任一所述的电磁暂态仿真误差评价处理方法,其特征在于,所述根据所述总体误差指标确定误差评价结果,包括:
根据所述总体误差指标的数值所在的预设数值区间确定误差评价结果。
7.一种电磁暂态仿真误差评价处理装置,其特征在于,包括:
获取单元,用于获取电磁暂态仿真的待分析数据,对待分析数据进行变换处理,得到直流分量、低频分量和高频分量;
计算单元,用于根据所述直流分量、所述低频分量和所述高频分量计算每个频点的幅值差异量和特征差异量,根据所述幅值差异量和所述特征差异量计算每个频点的全局差异量;
确定单元,用于根据正态分布函数、所述幅值差异量、所述特征差异量和所述全局差异量计算总体误差指标,并根据所述总体误差指标确定误差评价结果。
8.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至6任一所述方法。
9.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现权利要求1至6任一所述方法。
10.一种计算机程序产品,其特征在于,所述计算机程序产品包括计算机程序,所述计算机程序被处理器执行时实现权利要求1至6任一所述方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310489287.8A CN116595730A (zh) | 2023-05-04 | 2023-05-04 | 一种电磁暂态仿真误差评价处理方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310489287.8A CN116595730A (zh) | 2023-05-04 | 2023-05-04 | 一种电磁暂态仿真误差评价处理方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116595730A true CN116595730A (zh) | 2023-08-15 |
Family
ID=87603763
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310489287.8A Pending CN116595730A (zh) | 2023-05-04 | 2023-05-04 | 一种电磁暂态仿真误差评价处理方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116595730A (zh) |
-
2023
- 2023-05-04 CN CN202310489287.8A patent/CN116595730A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Peters et al. | Performance of the trim and fill method in the presence of publication bias and between‐study heterogeneity | |
KR102024953B1 (ko) | 동적 특이치 바이어스 감소 시스템 및 방법 | |
WO2008088652A2 (en) | Method and system for generating a predictive analysis of the performance of peer reviews | |
CN111143981B (zh) | 虚拟试验模型验证系统及方法 | |
CN105374367A (zh) | 异常帧检测方法和装置 | |
CN113360983B (zh) | 一种边坡可靠度分析与风险评估方法 | |
CN114089255B (zh) | 一种电容式电压互感器稳定性评估方法 | |
CN110825627B (zh) | 基于网格区域密度的适应性随机测试用例生成方法 | |
CN110567576B (zh) | 一种变电站厂界噪声超标原因的确定方法及装置 | |
CN116595730A (zh) | 一种电磁暂态仿真误差评价处理方法及装置 | |
JP2012230072A (ja) | 破壊検査システム、破壊検査方法、データ処理プログラムおよびプログラム記録媒体 | |
CN110098610B (zh) | 故障扰动下电力系统振荡主导模式的实时辨识方法及系统 | |
CN116933034A (zh) | 生态系统碳储量变化的确定方法、系统、电子设备及介质 | |
CN111914424A (zh) | 一种基于短期测风资料的设计风速取值方法及系统 | |
CN111382908A (zh) | 一种考虑大震时间相关性的地震随机事件集模拟方法 | |
CN115759263A (zh) | 一种基于因果推断的策略效果评估方法及装置 | |
CN105719181A (zh) | 风险等级评估方法及装置 | |
CN110717244B (zh) | 基于平均偏离度算法的数据信任度分析计算机仿真方法 | |
US8589138B2 (en) | On-the-fly improvement of certainty of statistical estimates in statistical design, with corresponding visual feedback | |
US8027764B2 (en) | Method and system for automatic test data generation for lookup tables | |
CN113552444A (zh) | 一种泄漏电流特征阈值在线整定方法及装置 | |
CN106772306B (zh) | 一种目标物的检测方法及服务器 | |
CN112653129B (zh) | 暂态功角稳定裕度估算方法、装置及系统 | |
CN111142134A (zh) | 一种坐标时间序列处理方法及装置 | |
CN113643754B (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 |