申请/专利权人：中国科学院合肥物质科学研究院

申请日：2021-04-28

公开（公告）日：2024-04-12

公开（公告）号：CN113486489B

主分类号：G06F30/20

分类号：G06F30/20;G06F17/12;G06Q50/06;G06F113/04

优先权：

专利状态码：有效-授权

法律状态：2024.04.12#授权;2021.10.26#实质审查的生效;2021.10.08#公开

摘要：本发明公开了一种适用于EAST电源装置脉冲有功功率的分解拟合方法，包括：1采集EAST电源装置有功功率信号序列；2、采用快速傅里叶变换对原信号序列进行分解，取幅值较大分量作为周期振荡分量序列的拟合值；3、除去原信号序列中的周期振荡分量序列的拟合值，从而得到待拟合的衰减分量序列并采用Prony算法进行拟合，得到待拟合的衰减分量序列的拟合值；4、由周期振荡分量序列的拟合值和待拟合的衰减分量序列的拟合值共同组成原信号序列的拟合值，从而完成有功功率的分解拟合。本发明能够建立精确EAST电源装置有功功率模型，并维持电网安全稳定运行。

主权项：1.一种适用于EAST电源装置脉冲有功功率的分解拟合方法，其特征是按如下步骤进行：步骤1、设采集EAST电源装置所得到的N个有功功率信号组成的原信号序列记为{Xk|k＝0,1,…,N-1}，其中，Xk表示第k个有功功率信号；每个有功功率信号之间的间隔时间为Δt，则采样频率为所述原信号序列{Xk|k＝0,1,…,N-1}是由周期振荡分量序列{Fk|k＝0,1,…,N-1}和衰减分量序列{Pk|k＝0,1,…,N-1}组成，其中，Fk表示第k个周期振荡分量，Pk表示第k个衰减分量；步骤2、采用快速傅里叶变换对原信号序列{Xk|k＝0,1,…,N-1}进行分解，取幅值较大的前NF个周期振荡分量作为周期振荡分量序列{Fk|k＝0,1,…,N-1}的拟合值其中，表示第k个周期振荡分量Fk的拟合值，并通过式1得到： 式1中，Ym表示经过快速傅里叶变换所得幅值较大的NF个周期震荡分量中的第m个周期震荡分量的幅值；fm表示经过快速傅里叶变换所得幅值较大的NF个周期震荡分量中的第m个周期震荡分量的频率；步骤3、除去原信号序列{Xk|k＝0,1,…,N-1}中的周期振荡分量序列{Fk|k＝0,1,…,N-1}的拟合值从而得到待拟合的衰减分量序列并采用Prony算法进行拟合；其中，表示待拟合的第k个衰减分量；步骤3.1、设待拟合的衰减分量序列所包含的衰减分量数目为步骤3.2、利用式2构造线性差分方程组，并利用最小二乘法进行求解，从而得到线性差分方程组的系数向量[ai|i＝1,…NP]T，其中，ai表示由待拟合的衰减分量组成的线性差分方程组的系数向量的第i个值： 步骤3.3、根据系数向量[ai|i＝1,…,NP]，利用式3构造特征多项式gz，并求取gz＝0时的根集合{zi|i＝1,…,NP}，其中，zi表示gz＝0时的第i个根： 式3中，z表示特征多项式gz中的未知数；表示未知数z的NP次方；步骤3.4、判断第i个根zi是否满足式4，若是，则执行步骤3.5，否则，令NP+1赋值给NP，并返回步骤3.2顺序执行； 式4中，zm表示gz＝0时的第m个根；表示第i个根zi的共轭复数；ziN-1表示第i个根zi的N-1次方；Re{ziN-1}表示ziN-1的实部；Im{ziN-1}表示ziN-1的虚部；R表示双精度形式的最大有限浮点数；步骤3.5、利用根集合{zi|i＝1,…,NP}构造范德蒙德行列式，利用式5构造待拟合的衰减分量的范德蒙德方程组，并用最小二乘法进行求解，从而得到范德蒙德行列式的系数向量{Bi|i＝1,…,NP}，其中，Bi表示范德蒙德行列式的系数向量中的第i个值： 步骤3.6、利用式6计算待拟合的衰减分量序列的拟合值其中，表示第k个衰减分量的拟合值： 式6中，zik表示第i个根zi的k次方；步骤4、由周期振荡分量序列{Fk|k＝0,1,…,N-1}的拟合值和待拟合的衰减分量序列的拟合值共同组成原信号序列{Xk|k＝0,1,…,N-1}的拟合值，从而完成有功功率的分解拟合。

全文数据：

权利要求：

百度查询： 中国科学院合肥物质科学研究院 适用于EAST电源装置脉冲有功功率的分解拟合方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD