CN113624269A - 一种基于谐波的频率响应测量系统及方法 - Google Patents
一种基于谐波的频率响应测量系统及方法 Download PDFInfo
- Publication number
- CN113624269A CN113624269A CN202110865727.6A CN202110865727A CN113624269A CN 113624269 A CN113624269 A CN 113624269A CN 202110865727 A CN202110865727 A CN 202110865727A CN 113624269 A CN113624269 A CN 113624269A
- Authority
- CN
- China
- Prior art keywords
- frequency
- output
- square wave
- analog
- frequency response
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 50
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims description 36
- 230000003750 conditioning effect Effects 0.000 claims description 33
- 238000012545 processing Methods 0.000 claims description 18
- 238000003775 Density Functional Theory Methods 0.000 claims description 11
- 238000000691 measurement method Methods 0.000 claims description 11
- 238000001228 spectrum Methods 0.000 claims description 9
- 238000001914 filtration Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 7
- 230000003321 amplification Effects 0.000 claims description 3
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D18/00—Testing or calibrating apparatus or arrangements provided for in groups G01D1/00 - G01D15/00
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
本发明属于线性系统频率响应测量技术领域,公开了一种基于谐波的频率响应测量系统和方法。利用占空比为50%的方波作为测量信号,把测量系统中的模数转换器过采样率设置为半整数,通过数学运算一次得到待测线性系统在多个频点上的频率响应。本发明提出的测量系统和方法,具有测量效率高、测量成本低的优点。
Description
技术领域
本发明属于线性系统频率响应测量技术领域,尤其涉及一种基于谐波的频率响应测量系统及方法。
背景技术
线性系统广泛存在。测量线性系统频率响应很重要。
测量线性系统频率响应,常用的方法有点频法、扫频法等。点频法和扫频法都是使用正弦信号作为测量信号输入待测线性系统。施加一次正弦信号进行一次测量,只能得到一个频点上的频率响应。而且,点频法和扫频法都需要正弦波产生电路,而正弦波产生电路比起方波产生电路,具有电路更复杂、成本更高等缺点。
所以,点频法和扫频法存在着测量效率相对较低、测量成本相对较高的问题。它们不能应用在对测量效率或测量成本有高要求的场景。
发明内容
本发明目的在于提供一种基于谐波的频率响应测量系统及方法,以解决点频法和扫频法存在的测量效率相对较低、测量成本相对较高的技术问题。
为解决上述技术问题,本发明的一种基于谐波的频率响应测量系统及方法的具体技术方案如下:
一种基于谐波的频率响应测量系统,包括:
一个频率可控方波生成器,输入是计算控制单元输出的控制信号Ctrl0,输出是方波信号x(t);
一个待测线性系统,输入是频率可控方波生成器输出的方波x(t),输出m(t)接入信号调理电路;
一个信号调理电路,输入是待测线性系统的输出m(t)和计算控制单元输出的控制信号Ctrl2,输出y(t)接入模数转换器;
一个模数转换器,输入是信号调理电路的输出y(t)和计算控制单元输出的控制信号Ctrl1,输出码字y(n)接入计算处理单元;
一个计算处理单元,输入是模数转换器输出的码字y(n),输出控制信号Ctrl0、Ctrl1和Ctrl2分别接入频率可控方波生成器、模数转换器和信号调理电路,并输出测量结果。
本发明还公开了一种基于谐波的频率响应测量方法,包括若干个测量轮次,每个测量轮次包括以下四个步骤:
第一步,根据目标测量频率,计算控制单元设定方波生成器生成的方波频率f0和方波幅度等参数,设定信号调理电路的低通截止频率fC和模数转换器过采样率OSR,启动方波生成器,其中过采样率OSR是半整数;
第二步,计算处理单元调整信号调理电路的增益AV,保证模数转换器的输入处在接近满量程、未饱和的状态;
第三步,计算处理单元对模数转换器输出的结果y(n)进行K点DFT计算,得到在频率点f0、3f0、5f0、……Mf0处信号y(t)的频谱数据Y(f0)、Y(3f0)、Y(5f0)、……Y(Mf0),其中M是奇数且满足M<OSR;
第四步,计算处理单元用Y(f0)、Y(3f0)、Y(5f0)、……Y(Mf0)分别除以信号调理电路的增益AV(f0)、AV(3f0)、AV(5f0)、……AV(Mf0),再分别除以信号x(t)的K点DFT频谱数据X(f0)、X(3f0)、X(5f0)、……X(Mf0),分别得到待测线性系统在频率点f0、3f0、5f0、……Mf0处的频率响应H(f0)、H(3f0)、H(5f0)、……H(Mf0)。
进一步地,所述频率可控方波生成器,生成的方波幅度、频率f0和相位等参数精确定义,方波频率f0可调,由计算控制单元输出的控制信号Ctrl0决定。
进一步地,所述信号调理电路,对输入信号m(t)进行线性放大和低通滤波,它的增益AV可调、低通滤波的截止频率fC可调,由计算控制单元输出的控制信号Ctrl2决定。
进一步地,所述模数转换器,其采样频率fs由计算控制单元输出的控制信号Ctrl1决定,与方波频率f0之间满足fs=K×f0,其中K是奇数。
进一步地,所述计算处理单元,可以进行数学运算,可以输出控制信号,可以读取外部数字码字输入,它输出控制信号Ctrl0接入频率可控方波生成器,控制方波频率等参数;输出控制信号Ctrl1接入模数转换器,控制模数转换器采样频率;输出控制信号Ctrl2接入信号调理电路,控制信号调理电路的增益Av和低通截止频率fC;并输出测量结果。
进一步地,所述M<K/2。
本发明的一种基于谐波的频率响应测量系统及方法具有以下优点:通过一次测量可以得到待测线性系统多个频点上的频率响应;使用方波生成器而不是正弦波产生电路,简化了电路,降低了成本。
附图说明
图1是本发明的一种基于谐波的频率响应测量系统基本框架示意图。
图2是本发明一种基于谐波的频率响应测量方法主要步骤示意图。
图3是本发明具体实施例基本框架示意图。
具体实施方式
为了更好地了解本发明的目的、结构及功能,下面结合附图,对本发明一种基于谐波的频率响应测量系统及方法做进一步详细的描述。
如图1所示,一种基于谐波的频率响应测量系统包括:
一个频率可控方波生成器,输入是计算控制单元输出的控制信号Ctrl0,输出是方波信号x(t);
一个待测线性系统,输入是频率可控方波生成器输出的方波x(t),输出m(t)接入信号调理电路;
一个信号调理电路,输入是待测线性系统的输出m(t)和计算控制单元输出的控制信号Ctrl2,输出y(t)接入模数转换器;
一个模数转换器,输入是信号调理电路的输出y(t)和计算控制单元输出的控制信号Ctrl1,输出码字y(n)接入计算处理单元;
一个计算处理单元,输入是模数转换器输出的码字y(n),输出控制信号Ctrl0、Ctrl1和Ctrl2分别接入频率可控方波生成器、模数转换器和信号调理电路,并输出测量结果。
如图2所示,本发明的一种基于谐波的频率响应测量方法可以包括若干个测量轮次,每个测量轮次主要包括以下四个步骤:
第一步,根据目标测量频率,计算控制单元设定方波生成器生成的方波频率f0和方波幅度等参数,设定信号调理电路的低通截止频率fC和模数转换器过采样率OSR,启动方波生成器,其中过采样率OSR是半整数;
第二步,计算处理单元调整信号调理电路的增益AV,保证模数转换器的输入处在接近满量程、未饱和的状态;
第三步,计算处理单元对模数转换器输出的结果y(n)进行K点DFT计算,得到在频率点f0、3f0、5f0、……Mf0处信号y(t)的频谱数据Y(f0)、Y(3f0)、Y(5f0)、……Y(Mf0),其中M是奇数且满足M<OSR;
第四步,计算处理单元用Y(f0)、Y(3f0)、Y(5f0)、……Y(Mf0)分别除以信号调理电路的增益AV(f0)、AV(3f0)、AV(5f0)、……AV(Mf0),再分别除以信号x(t)的K点DFT频谱数据X(f0)、X(3f0)、X(5f0)、……X(Mf0),分别得到待测线性系统在频率点f0、3f0、5f0、……Mf0处的频率响应H(f0)、H(3f0)、H(5f0)、……H(Mf0)。
优选地,所述频率可控方波生成器,生成的方波幅度、频率f0和相位等参数精确定义,方波频率f0可调,由计算控制单元输出的控制信号Ctrl0决定。
优选地,所述信号调理电路,对输入信号m(t)进行线性放大和低通滤波。它的增益AV可调、低通滤波的截止频率fC可调,由计算控制单元输出的控制信号Ctrl2决定。
优选地,所述模数转换器,其采样频率fs由计算控制单元输出的控制信号Ctrl1决定,与方波频率f0之间满足fs=K×f0,其中K是奇数。
优选地,所述计算处理单元,可以进行数学运算,可以输出控制信号,可以读取外部数字码字输入。它输出控制信号Ctrl0接入频率可控方波生成器,控制方波频率等参数;输出控制信号Ctrl1接入模数转换器,控制模数转换器采样频率;输出控制信号Ctrl2接入信号调理电路,控制信号调理电路的增益AV和低通截止频率fC;并输出测量结果。
优选地,上述测量方法中的参数M是奇数,且满足M<K/2。
上述测量系统和方法具有如下优点:通过一次测量可以得到待测线性系统多个频点上的频率响应;使用方波生成器而不是正弦波产生电路,简化了电路,降低了成本。
如果需要测量某线性系统许多频点上的频率响应,可以重复上述步骤,每次设定不同的方波频率f0和低通截止频率fC等参数。
上述测量系统和方法能一次测出待测线性系统多个频点上频率响应的原理是,基于线性电路的叠加定理和占空比为50%的方波的频谱特点,利用方波的基波和低次谐波来测量系统频率响应,通过对模数转换器过采样率的设置来消除可能的高次谐波频谱混叠。
记待测线性系统频率响应为H(f),记信号x(t)、y(t)在频率点jf0处的傅里叶系数分别为X(jf0)、Y(jf0),记信号调理电路在频率点jf0处的增益为AV(jf0),其中j=1,2,……M。X(jf0)、Y(jf0)、AV(jf0)是复数,包含幅度和相位。有
Y(jf0)=AV(jf0)X(jf0)H(jf0)
其中,j=1,2,......M (公式1)
记y(n)作K点DFT计算之后的结果是YF(jf0),j=1,2,......M。
于是有
……
由于我们限定测量信号x(t)的占空比为50%,信号x(t)的频谱具有偶次谐波能量为0的特点,即
X(mf0)=0,m=2,4,…∞ (公式5)
联立公式1和公式5,我们可以得到信号y(t)的偶数次谐波能量也为0,即
Y(mf0)=0,m=2,4,...∞ (公式6)
由于设定K和M都是奇数,所以y(t)中的K-1、K+1、3K-1、3K+1等偶次谐波能量为0。于是有
类似地,有
……
根据公式1可得
由于设定测量信号x(t)的占空比为50%,根据信号与系统知识,我们可以得到
由于低通滤波,|AV((2·i·K-1)f0)|<<|AV(f0)|,且公式11中的比值是一个很小的数,于是公式10中的比值接近于0。于是公式7中y(t)的2K-1、2K+1以及更高次的奇数次谐波可以忽略。类似的,公式8和公式9中的奇数次谐波也可以忽略。
所以,频谱混叠带来的影响可以忽略不记。
于是有
……
于是利用以下公式15即可求得待测线性系统的频率响应。
其中,j=1,2,......M (公式15)
本实施例中,我们测试一个线性系统在1kHz-10kHz内的频率响应。实施例的系统结构框图如图3所示。频率可控方波生成器采用555定时器结合电阻阵列来实现;信号调理电路采用程控放大器和程控滤波器,计算控制单元采用stm32单片机。
使用本发明提出的测量系统和方法,具体测量过程如下。
第一轮:
第一步,stm32单片机设定方波生成器生成的方波频率为1kHz,设定程控滤波器低通截止频率为33kHz,设定模数转换器过采样率为49.5,启动方波生成器;
第二步,stm32单片机调整程控放大器的增益,保证模数转换器处在接近满量程、未饱和的状态,记录此时程控放大器和程控滤波器的级联增益AV;
第三步,stm32单片机对数字信号y(n)做99点DFT,得到y(t)在1kHz、3kHz、5kHz、7kHz、9kHz这5个频率点上的频谱数据Y1、Y3、Y5、Y7、Y9;
第四步,stm32单片机用Y1、Y3、Y5、Y7、Y9分别除以程控放大器和程控滤波器在1kHz、3kHz、5kHz、7kHz、9kHz这5个频率点上的级联增益AV1、AV3、AV5、AV7、AV9,再分别除以对x(t)作99点DFT得到的在1kHz、3kHz、5kHz、7kHz、9kHz这5个频率点上x(t)的频谱数据X1、X3、X5、X7、X9,分别得到待测线性系统在频率点1kHz、3kHz、5kHz、7kHz、9kHz这5个频率点处的频率响应H1、H3、H5、H7、H9。
第二轮:
第一步,stm32单片机设定方波生成器生成的方波频率为2kHz,设定程控滤波器低通截止频率为66kHz,设定模数转换器过采样率为49.5,启动方波生成器;
第二步,stm32单片机调整程控放大器的增益,保证模数转换器处在接近满量程、未饱和的状态,记录此时程控放大器和程控滤波器的级联增益AV;
第三步,stm32单片机对数字信号y(n)做99点DFT,得到y(t)在2kHz、6kHz、10kHz这3个频率点上的频谱数据Y2、Y6、Y10;
第四步,stm32单片机用Y2、Y6、Y10分别除以程控放大器和程控滤波器在2kHz、6kHz、10kHz这3个频率点上的级联增益AV2、AV6、AV10,再分别除以对x(t)作99点DFT得到的在2kHz、6kHz、10kHz这3个频率点上x(t)的频谱数据X2、X6、X10,分别得到待测线性系统在频率点2kHz、6kHz、10kHz这3个频率点处的频率响应H2、H6、H10。
上述两轮测量,得到了待测线性系统在1kHz、2kHz、3kHz、5kHz、6kHz、7kHz、9kHz、10kHz这8个频率点上的频率响应。
使用本发明所示测量系统和方法,只需2次测量即可获得8个频率点上的频率响应。按照传统测量方法,需要测8次。本例中,本发明提出的测量系统和方法效率是传统点频法、扫频法的4倍。同时,本实施例中用555定时器结合电阻阵列取代正弦波生成电路,具有成本更低的优点。
由以上实施例可以看出,本发明提出的测量系统和方法,具有测量效率高、测量成本低的优势。
可以理解,本发明是通过一些实施例进行描述的,本领域技术人员知悉的,在不脱离本发明的精神和范围的情况下,可以对这些特征和实施例进行各种改变或等效替换。另外,在本发明的教导下,可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此,本发明不受此处所公开的具体实施例的限制,所有落入本申请的权利要求范围内的实施例都属于本发明所保护的范围内。
Claims (7)
1.一种基于谐波的频率响应测量系统,其特征在于,包括:
一个频率可控方波生成器,输入是计算控制单元输出的控制信号Ctrl0,输出是方波信号x(t);
一个待测线性系统,输入是频率可控方波生成器输出的方波x(t),输出m(t)接入信号调理电路;
一个信号调理电路,输入是待测线性系统的输出m(t)和计算控制单元输出的控制信号Ctrl2,输出y(t)接入模数转换器;
一个模数转换器,输入是信号调理电路的输出y(t)和计算控制单元输出的控制信号Ctrl1,输出码字y(n)接入计算处理单元;
一个计算处理单元,输入是模数转换器输出的码字y(n),输出控制信号Ctrl0、Ctrl1和Ctrl2分别接入频率可控方波生成器、模数转换器和信号调理电路,并输出测量结果。
2.一种利用如权利要求1所述的基于谐波的频率响应测量系统进行频率响应测量的方法,其特征在于,包括若干个测量轮次,每个测量轮次包括以下四个步骤:
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110865727.6A CN113624269A (zh) | 2021-07-29 | 2021-07-29 | 一种基于谐波的频率响应测量系统及方法 |
PCT/CN2022/083899 WO2023005248A1 (zh) | 2021-07-29 | 2022-03-30 | 一种基于谐波的频率响应测量系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110865727.6A CN113624269A (zh) | 2021-07-29 | 2021-07-29 | 一种基于谐波的频率响应测量系统及方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113624269A true CN113624269A (zh) | 2021-11-09 |
Family
ID=78381667
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110865727.6A Pending CN113624269A (zh) | 2021-07-29 | 2021-07-29 | 一种基于谐波的频率响应测量系统及方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN113624269A (zh) |
WO (1) | WO2023005248A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023005248A1 (zh) * | 2021-07-29 | 2023-02-02 | 浙江大学 | 一种基于谐波的频率响应测量系统及方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101034130A (zh) * | 2007-01-26 | 2007-09-12 | 上海欣泰通信技术有限公司 | 采用组合方波扫频测试频域特性的方法和装置 |
CN105606900A (zh) * | 2016-03-18 | 2016-05-25 | 华南理工大学 | 一种基于方波信号的单相谐波阻抗测量方法 |
CN105974343A (zh) * | 2016-06-20 | 2016-09-28 | 吉林大学 | 具有增益自动调节功能的地面磁共振信号检测装置及检测方法 |
US20180180652A1 (en) * | 2016-12-22 | 2018-06-28 | Texas Instruments Incorporated | Square-wave-based impedance analysis |
CN110297126A (zh) * | 2018-03-21 | 2019-10-01 | 中国科学院金属研究所 | 仪器放大器频响特性的频谱倍频分析测试法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3599546B2 (ja) * | 1997-01-20 | 2004-12-08 | アキュート株式会社 | 波形検出装置およびその装置を利用した状態監視システム |
CN203772956U (zh) * | 2014-05-21 | 2014-08-13 | 南京信息工程大学 | 一种基于正交调制原理的频率特性测试仪 |
CN113624269A (zh) * | 2021-07-29 | 2021-11-09 | 浙江大学 | 一种基于谐波的频率响应测量系统及方法 |
-
2021
- 2021-07-29 CN CN202110865727.6A patent/CN113624269A/zh active Pending
-
2022
- 2022-03-30 WO PCT/CN2022/083899 patent/WO2023005248A1/zh unknown
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101034130A (zh) * | 2007-01-26 | 2007-09-12 | 上海欣泰通信技术有限公司 | 采用组合方波扫频测试频域特性的方法和装置 |
CN105606900A (zh) * | 2016-03-18 | 2016-05-25 | 华南理工大学 | 一种基于方波信号的单相谐波阻抗测量方法 |
CN105974343A (zh) * | 2016-06-20 | 2016-09-28 | 吉林大学 | 具有增益自动调节功能的地面磁共振信号检测装置及检测方法 |
US20180180652A1 (en) * | 2016-12-22 | 2018-06-28 | Texas Instruments Incorporated | Square-wave-based impedance analysis |
CN110297126A (zh) * | 2018-03-21 | 2019-10-01 | 中国科学院金属研究所 | 仪器放大器频响特性的频谱倍频分析测试法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023005248A1 (zh) * | 2021-07-29 | 2023-02-02 | 浙江大学 | 一种基于谐波的频率响应测量系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2023005248A1 (zh) | 2023-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JPH0750136B2 (ja) | 周波数測定方法 | |
US20210194464A1 (en) | Fri sparse sampling kernel function construction method and circuit | |
WO2023005248A1 (zh) | 一种基于谐波的频率响应测量系统及方法 | |
US3833797A (en) | Statistical noise processor | |
Van Der Ouderaa et al. | Some formulas and applications of nonuniform sampling of bandwidth-limited signals | |
US5537344A (en) | High-speed processing apparatus and method, signal analyzing system, and measurement apparatus and method | |
CN114019235B (zh) | 一种基于同步离散傅里叶变换的频率响应测量方法 | |
Tavares | A comparison of integration and low-pass filtering | |
US4296374A (en) | Wideband digital spectrometer | |
CN114252899B (zh) | 一种核信号的级联冲激卷积成形方法和装置 | |
JP2867769B2 (ja) | 音響測定方法およびその装置 | |
CN107124678A (zh) | 一种音频谐波失真的测量系统 | |
CN111487476B (zh) | 一种准峰值检波方法和准峰值检波器 | |
Polge et al. | Evaluating analog-to-digital converters | |
CN104267259B (zh) | 电力谐波测量方法和系统 | |
CN112146746B (zh) | 声级计量方法和声级计 | |
CN117330939B (zh) | 基于ate的开关电容滤波器调测方法、系统及存储介质 | |
JP3664812B2 (ja) | 歪み率測定方法およびその装置 | |
Kvedaras et al. | Method of dynamic parameters measurement of high-speed D/A converters | |
Bradford | On errors in published formulas for the noise-response integrals I n | |
Tiwari et al. | FPGA implementation of wavelet filters for power system harmonics estimation | |
Huang et al. | Multi-Channel FRI Sampling System Based on Non-Ideal Filters | |
RU164241U1 (ru) | Устройство управления коэффициентом усиления передающего устройства рлс на основе восстановления параметров усиливаемого сигнала | |
RU2251791C2 (ru) | Способ фильтрации измерительных сигналов | |
HOLT et al. | Resampling the output of switched networks |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211109 |