CN114113780A - 一种线性系统频域特性测试方法及系统 - Google Patents

一种线性系统频域特性测试方法及系统 Download PDF

Info

Publication number
CN114113780A
CN114113780A CN202111480953.9A CN202111480953A CN114113780A CN 114113780 A CN114113780 A CN 114113780A CN 202111480953 A CN202111480953 A CN 202111480953A CN 114113780 A CN114113780 A CN 114113780A
Authority
CN
China
Prior art keywords
frequency domain
frequency
input signal
linear system
testing
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
Application number
CN202111480953.9A
Other languages
English (en)
Inventor
李理
贺晨
洪权
刘伟良
熊尚峰
蔡昱华
吴晋波
刘志豪
龚禹生
肖俊先
李林山
陈胜春
曾林俊
牟秀君
吴雪琴
张伦
肖纳敏
吴春
阳隋
金梦汝
向东
肖阳
罗旖旎
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
State Grid Hunan Electric Power Co Ltd
Training Center of State Grid Hunan Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd
State Grid Hunan Electric Power Co Ltd
Training Center of State Grid Hunan Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Hunan Electric Power Co Ltd, State Grid Hunan Electric Power Co Ltd, Training Center of State Grid Hunan Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202111480953.9A priority Critical patent/CN114113780A/zh
Publication of CN114113780A publication Critical patent/CN114113780A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Testing Electric Properties And Detecting Electric Faults (AREA)

Abstract

本发明涉及电力系统测试技术领域，公开了一种线性系统频域特性测试方法及系统，该方法在特定频率点对系统频域特性进行测量，通过插值得到系统在频率段的频域特性，可为设计系统控制器提供快速、准确的参考。

Description

一种线性系统频域特性测试方法及系统
技术领域
本发明涉及电力系统测试技术领域，尤其涉及一种线性系统频域特性测试方法及系统。
背景技术
线性系统的信号具有时域和频域的特性，通过对时域、频域特性进行分析，可以获得线性系统的本质特征。系统控制器的设计常以系统频域特性为主，时域特性为辅，主要利用系统波特图、奈奎斯特图指导控制器的设计。
系统的频域特性可通过动态信号分析仪进行测取，其原理为向系统施加一定频率的扰动信号，测取系统的输出信号，根据输入和输出信号辨识系统的频域特性。以电力系统稳定器(PSS)试验为例，试验中通过输入含有多个频段信息的白噪声信号对电力系统进行扰动，根据白噪声输入信号及系统反馈输出信号进行频域辨识，获得电力系统在该频段的频域特性。采用白噪声信号进行频域特性测量的结果较为精确，但是计算量大，耗时较长，对动态信号分析仪的功能要求较高，因此可以针对具体系统要求对频域特性测试方法进行一定的优化。
发明内容
本发明提供了一种线性系统频域特性测试方法及系统，以解决现有技术中存在的问题。
为了实现上述目的，本发明通过如下的技术方案来实现：
第一方面，本发明提供一种线性系统频域特性测试方法，包括：
S1：构造特定离散频率的正弦波扰动信号作为进行线性系统频域特性辨识的输入信号；
S2：施加输入信号对线性系统进行扰动，测取与输入信号对应的系统输出信号；
S3：对输入信号及输出信号在特定频率点进行频域特性辨识；
S4：计算线性系统在特定频率点的幅频特性、相频特性。
可选地，所述输入信号满足如下公式：
Figure BDA0003395241270000011
式中，input表示时域输入信号，fn为离散的频率点，Ts为采样周期，n为离散的时间序列，π为常数。
可选地，所述方法还包括构造连续的低通滤波环节的传递函数，所述测取与输入信号对应的系统输出信号，包括：
根据所述传递函数测取与输入信号对应的系统输出信号。
可选地，所述输出信号满足如下关系式：
Figure BDA0003395241270000012
式中，output表示输出信号，h表示传递函数的时域表达。
第二方面，本申请实施例提供一种线性系统频域特性测试系统，包括：包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述第一方面所述方法的步骤。
有益效果：
本发明提供的线性系统频域特性测试方法，在特定频率点对系统频域特性进行测量，通过插值得到系统在频率段的频域特性，可为设计系统控制器提供快速、准确的参考。
附图说明
图1为本发明优选实施例的一种线性系统频域特性测试方法的流程图。
具体实施方式
下面对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
除非另作定义，本发明中使用的技术术语或者科学术语应当为本发明所属领域内具有一般技能的人士所理解的通常意义。本发明中使用的“第一”、“第二”以及类似的词语并不表示任何顺序、数量或者重要性，而只是用来区分不同的组成部分。同样，“一个”或者“一”等类似词语也不表示数量限制，而是表示存在至少一个。“连接”或者“相连”等类似的词语并非限定于物理的或者机械的连接，而是可以包括电性的连接，不管是直接的还是间接的。“上”、“下”、“左”、“右”等仅用于表示相对位置关系，当被描述对象的绝对位置改变后，则该相对位置关系也相应地改变。
应理解，由于线性系统的频域特性一般是连续的，整个系统的频域特性可由特定点的频域特性插值代替。比如PSS试验中对0.1～2Hz之间的频段的测量，可由0.1～2Hz之间的频段的以0.1Hz为步长的20个频率点的频域特性代替。基于以上方法，仍以PSS试验为例，可以对频域特性测试过程进行简化，如构造20个特定频率的正弦波信号作为输入扰动源对系统进行扰动，然后辨识这20个频率点的输入、输出频域特性，即可完成对电力系统在以上频段的频域特性辨识。因此为了更加简便的进行系统频域特性的测量，可基于离散频率正弦波扰动信号对线性系统频域特性测试，为评判系统频域性能提供快速、准确的参考。基于此，本申请实施例提供一种线性系统频域特性测试方法。
请参见图1，本申请实施例提供一种线性系统频域特性测试方法，包括：
S1：构造特定离散频率的正弦波扰动信号作为进行线性系统频域特性辨识的输入信号。
在该步骤中，系统在特定频率点附近含有丰富的信息，设定输入扰动信号为特定离散频率正弦波的叠加，叠加的信号作为白噪声扰动信号的简化，用作系统扰动的输入。
S2：施加输入信号对线性系统进行扰动，测取与输入信号对应的系统输出信号。
在该步骤中，施加输入扰动信号后，同时测取系统输出信号，输入、输出的测量点在时域上应能一一对应，方便后期进行快速傅里叶变换处理。
S3：对输入信号及输出信号在特定频率点进行频域特性辨识。
在该步骤中，利用快速傅里叶变换对输入、输出信号进行频域分析，得到特定频率点的频域特性，包含幅频特性及相频特性。
S4：计算线性系统在特定频率点的幅频特性、相频特性。
在该步骤中，基于输入、输出信号的频域特性，对特定频率点对应进行频域幅值相除、频域相位相减，获得系统在特定频率的频域特性。
上述的线性系统频域特性测试方法，在特定频率点对系统频域特性进行测量，通过插值得到系统在频率段的频域特性，可为设计系统控制器提供快速、准确的参考。
可选地，所述输入信号满足如下公式：
Figure BDA0003395241270000031
式中，input表示时域输入信号，fn为离散的频率点，Ts为采样周期，n为离散的时间序列，π为常数。
可选地，所述方法还包括构造连续的低通滤波环节的传递函数，所述测取与输入信号对应的系统输出信号，包括：
根据所述传递函数测取与输入信号对应的系统输出信号。
可选地，所述输出信号满足如下关系式：
Figure BDA0003395241270000032
式中，output表示输出信号，h表示传递函数的时域表达。
在一示例中，根据以上方法对某一线性系统的频域特性进行仿真测取，时域输入信号input如式(1)，其中fn为离散的频率点，Ts为采样周期，n为离散的时间序列。
Figure BDA0003395241270000033
构造连续的低通滤波环节的传递函数为H＝1/(1+Tr*s)，其中Tr＝0.02，H的时域表达式h如式(2)。
Figure BDA0003395241270000041
离散化后h[n]如式(3)，其中Ts为采样周期，n为离散的时间序列。
h[n]＝50e-50nTs (3)
根据卷积计算时域输出信号如式(4)。
Figure BDA0003395241270000042
分别对输入信号input及输出信号output进行快速傅里叶分析(fft)。分别在0.1～2Hz之间的以0.1Hz为步长的20个频率点，分别取输入、输出信号的频率特性，计算低通滤波环节H在该点的频率特性。即在特定的频率下，计算幅值特性如式(5)及相位特性如式(6)，其中amp表示幅值特性函数，deg表示相位特性函数，k为特定频率点，Input为输入信号频域特性，Output为输出信号频域特性。
amp(H[k])＝amp(Output[k])/amp(Input[k]) (5)
deg(H[k])＝deg(Output[k])-deg(Input[k]) (6)
本申请实施例还提供一种线性系统频域特性测试系统，包括：包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述方法的步骤。
该线性系统频域特性测试系统可以实现上述线性系统频域特性测试方法的各个实施例，且能达到相同的有益效果，此处，不做赘述。
以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。

Claims (5)

1.一种线性系统频域特性测试方法，其特征在于，包括：
S1：构造特定离散频率的正弦波扰动信号作为进行线性系统频域特性辨识的输入信号；
S2：施加输入信号对线性系统进行扰动，测取与输入信号对应的系统输出信号；
S3：对输入信号及输出信号在特定频率点进行频域特性辨识；
S4：计算线性系统在特定频率点的幅频特性、相频特性。
2.根据权利要求1所述的线性系统频域特性测试方法，其特征在于，所述输入信号满足如下公式：
Figure FDA0003395241260000011
式中，input表示时域输入信号，fn为离散的频率点，Ts为采样周期，n为离散的时间序列，π为常数。
3.根据权利要求1所述的线性系统频域特性测试方法，其特征在于，所述方法还包括构造连续的低通滤波环节的传递函数，所述测取与输入信号对应的系统输出信号，包括：
根据所述传递函数测取与输入信号对应的系统输出信号。
4.根据权利要求3所述的线性系统频域特性测试方法，其特征在于，所述输出信号满足如下关系式：
Figure FDA0003395241260000012
式中，output表示输出信号，h表示传递函数的时域表达。
5.一种线性系统频域特性测试系统，包括：包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述权利要求1至4中任一所述方法的步骤。
CN202111480953.9A 2021-12-06 2021-12-06 一种线性系统频域特性测试方法及系统 Pending CN114113780A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111480953.9A CN114113780A (zh) 2021-12-06 2021-12-06 一种线性系统频域特性测试方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111480953.9A CN114113780A (zh) 2021-12-06 2021-12-06 一种线性系统频域特性测试方法及系统

Publications (1)

Publication Number Publication Date
CN114113780A true CN114113780A (zh) 2022-03-01

Family

ID=80367198

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111480953.9A Pending CN114113780A (zh) 2021-12-06 2021-12-06 一种线性系统频域特性测试方法及系统

Country Status (1)

Country Link
CN (1) CN114113780A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116599537A (zh) * 2023-05-18 2023-08-15 哈尔滨市科佳通用机电股份有限公司 一种铁路移频信号译码算法的单元测试方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005078559A (ja) * 2003-09-03 2005-03-24 Fuji Electric Holdings Co Ltd 特性不明システムの同定装置
CN105242111A (zh) * 2015-09-17 2016-01-13 清华大学 一种采用类脉冲激励的频响函数测量方法
CN108932215A (zh) * 2018-05-21 2018-12-04 华南理工大学 一种用于电力系统多输入多输出线性化模型辨识的低频段多正弦信号设计方法
CN109858170A (zh) * 2019-02-14 2019-06-07 苏州同元软控信息技术有限公司 一种强非线性Modelica系统模型的频率特性估算系统及方法
US20190383700A1 (en) * 2018-06-19 2019-12-19 Palo Alto Research Center Incorporated Model-based diagnosis in frequency domain
CN113612522A (zh) * 2021-08-13 2021-11-05 陕西航天技术应用研究院有限公司 一种频域偏离的处理方法及装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005078559A (ja) * 2003-09-03 2005-03-24 Fuji Electric Holdings Co Ltd 特性不明システムの同定装置
CN105242111A (zh) * 2015-09-17 2016-01-13 清华大学 一种采用类脉冲激励的频响函数测量方法
CN108932215A (zh) * 2018-05-21 2018-12-04 华南理工大学 一种用于电力系统多输入多输出线性化模型辨识的低频段多正弦信号设计方法
US20190383700A1 (en) * 2018-06-19 2019-12-19 Palo Alto Research Center Incorporated Model-based diagnosis in frequency domain
CN109858170A (zh) * 2019-02-14 2019-06-07 苏州同元软控信息技术有限公司 一种强非线性Modelica系统模型的频率特性估算系统及方法
CN113612522A (zh) * 2021-08-13 2021-11-05 陕西航天技术应用研究院有限公司 一种频域偏离的处理方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
周子强: "《机械工程测控技术》", 国防工业出版社, pages: 71 - 72 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116599537A (zh) * 2023-05-18 2023-08-15 哈尔滨市科佳通用机电股份有限公司 一种铁路移频信号译码算法的单元测试方法

Similar Documents

Publication Publication Date Title
US20120007583A1 (en) Method and device for frequency response measurement
CN103457603A (zh) 一种基于平均频谱测试adc动态参数的方法
CN102955068A (zh) 一种基于压缩采样正交匹配追踪的谐波检测方法
CN114113780A (zh) 一种线性系统频域特性测试方法及系统
KR100911685B1 (ko) 코히어런트하지 않게 샘플링된 데이타의 파워 스펙트럼을측정하기 위한 저누설 방법
Zamora et al. Multi-dimensional ringdown modal analysis by filtering
CN110020452B (zh) 电磁兼容性测试仿真分析方法及其系统
US20150276840A1 (en) Noise modulation for on-chip noise measurement
Tarasiuk Comparative study of various methods of DFT calculation in the wake of IEC Standard 61000-4-7
JP5712255B2 (ja) フーリエ解析による周波数測定方法および周波数測定装置
CN104931777B (zh) 一种基于两条dft复数谱线的信号频率测量方法
Betta et al. A multi-application FFT analyzer based on a DSP architecture
CN106124081A (zh) 永磁同步电机精确多点实时测温方法及系统
Yavorskyj et al. Discrete LS estimates of correlation function of bi-periodically correlated random signals
Goh et al. Hybrid FFT‐ADALINE algorithm with fast estimation of harmonics in power system
Yi et al. Instantaneous harmonic decomposition technique for three‐phase current based on multiple reference coordinates
Aiello et al. A PC-based instrument for harmonics and interharmonics measurement in power supply systems
Önal et al. Hilbert–Huang transform based approach for measurement of voltage flicker magnitude and frequency
Li et al. Digital prolate spheroidal window-based S-transform for time-varying harmonic analysis
CN115494303A (zh) 一种emi接收机信号转化方法、装置及存储介质
JP5251291B2 (ja) 微弱信号解析装置、微弱信号解析方法、及び微弱信号解析プログラム
CN111551785B (zh) 基于无迹卡尔曼滤波的频率与谐波检测方法
CN114487589A (zh) 电网宽频信号自适应测量方法、装置及系统
CN104849551B (zh) 一种谐相角分析方法
US8984342B2 (en) Unit testing and analysis using a stored reference signal

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