一种基于谐波的频率响应测量系统及方法

1.本发明属于线性系统频率响应测量技术领域，尤其涉及一种基于谐波的频率响应测量系统及方法。


背景技术：

2.线性系统广泛存在。测量线性系统频率响应很重要。
3.测量线性系统频率响应，常用的方法有点频法、扫频法等。点频法和扫频法都是使用正弦信号作为测量信号输入待测线性系统。施加一次正弦信号进行一次测量，只能得到一个频点上的频率响应。而且，点频法和扫频法都需要正弦波产生电路，而正弦波产生电路比起方波产生电路，具有电路更复杂、成本更高等缺点。
4.所以，点频法和扫频法存在着测量效率相对较低、测量成本相对较高的问题。它们不能应用在对测量效率或测量成本有高要求的场景。


技术实现要素：

5.本发明目的在于提供一种基于谐波的频率响应测量系统及方法,以解决点频法和扫频法存在的测量效率相对较低、测量成本相对较高的技术问题。
6.为解决上述技术问题，本发明的一种基于谐波的频率响应测量系统及方法的具体技术方案如下：
7.一种基于谐波的频率响应测量系统，包括：
8.一个频率可控方波生成器，输入是计算控制单元输出的控制信号ctrl0，输出是方波信号x(t)；
9.一个待测线性系统，输入是频率可控方波生成器输出的方波x(t)，输出m(t)接入信号调理电路；
10.一个信号调理电路，输入是待测线性系统的输出m(t)和计算控制单元输出的控制信号ctrl2，输出y(t)接入模数转换器；
11.一个模数转换器，输入是信号调理电路的输出y(t)和计算控制单元输出的控制信号ctrl1，输出码字y(n)接入计算处理单元；
12.一个计算处理单元，输入是模数转换器输出的码字y(n)，输出控制信号ctrl0、ctrl1和ctrl2分别接入频率可控方波生成器、模数转换器和信号调理电路，并输出测量结果。
13.本发明还公开了一种基于谐波的频率响应测量方法，包括若干个测量轮次，每个测量轮次包括以下四个步骤：
14.第一步，根据目标测量频率，计算控制单元设定方波生成器生成的方波频率f0和方波幅度等参数，设定信号调理电路的低通截止频率f
c
和模数转换器过采样率osr，启动方波生成器，其中过采样率osr是半整数；
15.第二步，计算处理单元调整信号调理电路的增益a
v
，保证模数转换器的输入处在
接近满量程、未饱和的状态；
16.第三步，计算处理单元对模数转换器输出的结果y(n)进行k点dft计算，得到在频率点f0、3f0、5f0、
……
mf0处信号y(t)的频谱数据y(f0)、y(3f0)、y(5f0)、
……
y(mf0)，其中m是奇数且满足m<osr；
17.第四步，计算处理单元用y(f0)、y(3f0)、y(5f0)、
……
y(mf0)分别除以信号调理电路的增益a
v
(f0)、a
v
(3f0)、a
v
(5f0)、
……
a
v
(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)。
18.进一步地，所述频率可控方波生成器，生成的方波幅度、频率f0和相位等参数精确定义，方波频率f0可调，由计算控制单元输出的控制信号ctrl0决定。
19.进一步地，所述信号调理电路，对输入信号m(t)进行线性放大和低通滤波，它的增益a
v
可调、低通滤波的截止频率f
c
可调，由计算控制单元输出的控制信号ctrl2决定。
20.进一步地，所述模数转换器，其采样频率f
s
由计算控制单元输出的控制信号ctrl1决定，与方波频率f0之间满足f
s
＝k
×
f0，其中k是奇数。
21.进一步地，所述计算处理单元，可以进行数学运算，可以输出控制信号，可以读取外部数字码字输入，它输出控制信号ctrl0接入频率可控方波生成器，控制方波频率等参数；输出控制信号ctrl1接入模数转换器，控制模数转换器采样频率；输出控制信号ctrl2接入信号调理电路，控制信号调理电路的增益a
v
和低通截止频率f
c
；并输出测量结果。
22.进一步地，所述m＜k/2。
23.本发明的一种基于谐波的频率响应测量系统及方法具有以下优点：通过一次测量可以得到待测线性系统多个频点上的频率响应；使用方波生成器而不是正弦波产生电路，简化了电路，降低了成本。
附图说明
24.图1是本发明的一种基于谐波的频率响应测量系统基本框架示意图。
25.图2是本发明一种基于谐波的频率响应测量方法主要步骤示意图。
26.图3是本发明具体实施例基本框架示意图。
具体实施方式
27.为了更好地了解本发明的目的、结构及功能，下面结合附图，对本发明一种基于谐波的频率响应测量系统及方法做进一步详细的描述。
28.如图1所示，一种基于谐波的频率响应测量系统包括：
29.一个频率可控方波生成器，输入是计算控制单元输出的控制信号ctrl0，输出是方波信号x(t)；
30.一个待测线性系统，输入是频率可控方波生成器输出的方波x(t)，输出m(t)接入信号调理电路；
31.一个信号调理电路，输入是待测线性系统的输出m(t)和计算控制单元输出的控制信号ctrl2，输出y(t)接入模数转换器；
32.一个模数转换器，输入是信号调理电路的输出y(t)和计算控制单元输出的控制信
号ctrl1，输出码字y(n)接入计算处理单元；
33.一个计算处理单元，输入是模数转换器输出的码字y(n)，输出控制信号ctrl0、ctrl1和ctrl2分别接入频率可控方波生成器、模数转换器和信号调理电路，并输出测量结果。
34.如图2所示，本发明的一种基于谐波的频率响应测量方法可以包括若干个测量轮次，每个测量轮次主要包括以下四个步骤：
35.第一步，根据目标测量频率，计算控制单元设定方波生成器生成的方波频率f0和方波幅度等参数，设定信号调理电路的低通截止频率f
c
和模数转换器过采样率osr，启动方波生成器，其中过采样率osr是半整数；
36.第二步，计算处理单元调整信号调理电路的增益a
v
，保证模数转换器的输入处在接近满量程、未饱和的状态；
37.第三步，计算处理单元对模数转换器输出的结果y(n)进行k点dft计算，得到在频率点f0、3f0、5f0、
……
mf0处信号y(t)的频谱数据y(f0)、y(3f0)、y(5f0)、
……
y(mf0)，其中m是奇数且满足m<osr；
38.第四步，计算处理单元用y(f0)、y(3f0)、y(5f0)、
……
y(mf0)分别除以信号调理电路的增益a
v
(f0)、a
v
(3f0)、a
v
(5f0)、
……
a
v
(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)。
39.优选地，所述频率可控方波生成器，生成的方波幅度、频率f0和相位等参数精确定义，方波频率f0可调，由计算控制单元输出的控制信号ctrl0决定。
40.优选地，所述信号调理电路，对输入信号m(t)进行线性放大和低通滤波。它的增益a
v
可调、低通滤波的截止频率f
c
可调，由计算控制单元输出的控制信号ctrl2决定。
41.优选地，所述模数转换器，其采样频率f
s
由计算控制单元输出的控制信号ctrl1决定，与方波频率f0之间满足f
s
＝k
×
f0，其中k是奇数。
42.优选地，所述计算处理单元，可以进行数学运算，可以输出控制信号，可以读取外部数字码字输入。它输出控制信号ctrl0接入频率可控方波生成器，控制方波频率等参数；输出控制信号ctrl1接入模数转换器，控制模数转换器采样频率；输出控制信号ctrl2接入信号调理电路，控制信号调理电路的增益a
v
和低通截止频率f
c
；并输出测量结果。
43.优选地，上述测量方法中的参数m是奇数，且满足m＜k/2。
44.上述测量系统和方法具有如下优点：通过一次测量可以得到待测线性系统多个频点上的频率响应；使用方波生成器而不是正弦波产生电路，简化了电路，降低了成本。
45.如果需要测量某线性系统许多频点上的频率响应，可以重复上述步骤，每次设定不同的方波频率f0和低通截止频率f
c
等参数。
46.上述测量系统和方法能一次测出待测线性系统多个频点上频率响应的原理是，基于线性电路的叠加定理和占空比为50％的方波的频谱特点，利用方波的基波和低次谐波来测量系统频率响应，通过对模数转换器过采样率的设置来消除可能的高次谐波频谱混叠。
47.记待测线性系统频率响应为h(f)，记信号x(t)、y(t)在频率点jf0处的傅里叶系数分别为x(jf0)、y(jf0)，记信号调理电路在频率点jf0处的增益为a
v
(jf0)，其中j＝1,2,
……
m。x(jf0)、y(jf0)、a
v
(jf0)是复数，包含幅度和相位。有
48.y(jf0)＝av(jf0)x(jf0)h(jf0)
49.其中，j＝1，2，......m
ꢀꢀ
(公式1)
50.记y(n)作k点dft计算之后的结果是y
f
(jf0)，j＝1，2，......m。
51.于是有
[0052][0053][0054]
……
[0055][0056]
由于我们限定测量信号x(t)的占空比为50％，信号x(t)的频谱具有偶次谐波能量为0的特点，即
[0057]
x(mf0)＝0，m＝2，4，
…
∞
ꢀꢀ
(公式5)
[0058]
联立公式1和公式5，我们可以得到信号y(t)的偶数次谐波能量也为0，即
[0059]
y(mf0)＝0，m＝2，4，...∞
ꢀꢀ
(公式6)
[0060]
由于设定k和m都是奇数，所以y(t)中的k
‑
1、k+1、3k
‑
1、3k+1等偶次谐波能量为0。于是有
[0061][0061][0062][0063]
类似地，有
[0064][0065]
……
[0066][0067]
根据公式1可得
[0068][0069]
由于设定测量信号x(t)的占空比为50％，根据信号与系统知识，我们可以得到
[0070][0071]
由于低通滤波，|a
v
((2
·
i
·
k
‑
1)f0)|＜＜|a
v
(f0)|，且公式11中的比值是一个很小的数，于是公式10中的比值接近于0。于是公式7中y(t)的2k
‑
1、2k+1以及更高次的奇数次谐波可以忽略。类似的，公式8和公式9中的奇数次谐波也可以忽略。
[0072]
所以，频谱混叠带来的影响可以忽略不记。
[0073]
于是有
[0074][0075][0076]
……
[0077][0078]
于是利用以下公式15即可求得待测线性系统的频率响应。
[0079][0080]
其中，j＝1，2，......m
ꢀꢀ
(公式15)
[0081]
本实施例中，我们测试一个线性系统在1khz
‑
10khz内的频率响应。实施例的系统结构框图如图3所示。频率可控方波生成器采用555定时器结合电阻阵列来实现；信号调理电路采用程控放大器和程控滤波器，计算控制单元采用stm32单片机。
[0082]
使用本发明提出的测量系统和方法，具体测量过程如下。
[0083]
第一轮：
[0084]
第一步，stm32单片机设定方波生成器生成的方波频率为1khz，设定程控滤波器低通截止频率为33khz，设定模数转换器过采样率为49.5，启动方波生成器；
[0085]
第二步，stm32单片机调整程控放大器的增益，保证模数转换器处在接近满量程、未饱和的状态，记录此时程控放大器和程控滤波器的级联增益a
v
；
[0086]
第三步，stm32单片机对数字信号y(n)做99点dft，得到y(t)在1khz、3khz、5khz、7khz、9khz这5个频率点上的频谱数据y1、y3、y5、y7、y9；
[0087]
第四步，stm32单片机用y1、y3、y5、y7、y9分别除以程控放大器和程控滤波器在1khz、3khz、5khz、7khz、9khz这5个频率点上的级联增益a
v1
、a
v3
、a
v5
、a
v7
、a
v9
，再分别除以对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。
[0088]
第二轮：
[0089]
第一步，stm32单片机设定方波生成器生成的方波频率为2khz，设定程控滤波器低通截止频率为66khz，设定模数转换器过采样率为49.5，启动方波生成器；
[0090]
第二步，stm32单片机调整程控放大器的增益，保证模数转换器处在接近满量程、未饱和的状态，记录此时程控放大器和程控滤波器的级联增益a
v
；
[0091]
第三步，stm32单片机对数字信号y(n)做99点dft，得到y(t)在2khz、6khz、10khz这3个频率点上的频谱数据y2、y6、y
10
；
[0092]
第四步，stm32单片机用y2、y6、y
10
分别除以程控放大器和程控滤波器在2khz、6khz、10khz这3个频率点上的级联增益a
v2
、a
v6
、a
v10
，再分别除以对x(t)作99点dft得到的在2khz、6khz、10khz这3个频率点上x(t)的频谱数据x2、x6、x
10
,分别得到待测线性系统在频率点2khz、6khz、10khz这3个频率点处的频率响应h2、h6、h
10
。
[0093]
上述两轮测量，得到了待测线性系统在1khz、2khz、3khz、5khz、6khz、7khz、9khz、10khz这8个频率点上的频率响应。
[0094]
使用本发明所示测量系统和方法，只需2次测量即可获得8个频率点上的频率响应。按照传统测量方法，需要测8次。本例中，本发明提出的测量系统和方法效率是传统点频法、扫频法的4倍。同时，本实施例中用555定时器结合电阻阵列取代正弦波生成电路，具有成本更低的优点。
[0095]
由以上实施例可以看出，本发明提出的测量系统和方法，具有测量效率高、测量成本低的优势。
[0096]
可以理解，本发明是通过一些实施例进行描述的，本领域技术人员知悉的，在不脱离本发明的精神和范围的情况下，可以对这些特征和实施例进行各种改变或等效替换。另外，在本发明的教导下，可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此，本发明不受此处所公开的具体实施例的限制，所有落入本技术的权利要求范围内的实施例都属于本发明所保护的范围内。