使用数字示波器DS6104测量交流信号的幅值和相位
01简介
使用普通的万用表测量交流信号的时候,通常会遇到 万用表的频率响应 的问题。使用可以联网的示波器可以获得它采集到的数据,进而可以计算出所测量的交流信号的有效值和相位。
这里通过实验来确定使用示波器测量交流信号的参数方法中的以下问题:
- 测量的精度有多高? 测量的频率范围有多宽?
- 获得高精度的测量参数,对于示波器参数如何设置?
▲ 实验所使用的示波器:ds6104
实验所使用的DS6104是一款可以通过联网读取各通道采集数据的示波器,为自动化测量提供了很好的支持。对它的 DS6104局域网络编成接口 可以从网络中下载相关的文档。
DS6104对于每个通道采样数据点个数N=1400N = 1400N=1400,采样的有效位数为8bit。下面显示了DS6104采集一个有效值为E=2.39VE = 2.39VE=2.39V,频率为f=1000Hzf = 1000Hzf=1000Hz的正弦波波形。该信号是由 DS345 信号发生器产生。该信号发生器同时还可以产生同相的同步方波信号,可以便于示波器用于同步采集信号的波形。
▲ 通过DS6104读取的信号发生器所产生的1000Hz信号及其同步方波信号
02测量方法
假设被测量信号的频率为fff,有效值为EEE,那么信号的表达式为:f(t)=2⋅E⋅sin(2πf⋅t+θ)f\left( t \right) = \sqrt 2 \cdot E \cdot \sin \left( {2\pi f \cdot t + \theta } \right)f(t)=2⋅E⋅sin(2πf⋅t+θ)示波器采样时间间隔为TsT_sTs,采样点个数NNN。所得到的采样数据为:f[n]=2⋅E⋅sin(2πf⋅nTs+θ),n=0,1,...,N−1f\left[ n \right] = \sqrt 2 \cdot E \cdot \sin \left( {2\pi f \cdot nT_s + \theta } \right),\,\,\,\,n = 0,1,...,N - 1f[n]=2⋅E⋅sin(2πf⋅nTs+θ),n=0,1,...,N−1
从数据中计算信号的E,θE,\thetaE,θ,采用以下方法。首先计算序列与cos[nTs],sin[nTs]\cos \left[ {nT_s } \right],\sin \left[ {nT_s } \right]cos[nTs],sin[nTs]的投影(内积):
Ys=1N∑n=0N−1f[n]⋅sin(2πf⋅nTs)Y_s = {1 \over N}\sum\limits_{n = 0}^{N - 1} {f\left[ n \right] \cdot \sin \left( {2\pi f \cdot nT_s } \right)}Ys=N1n=0∑N−1f[n]⋅sin(2πf⋅nTs)Yc=1N∑n=0N−1f[n]⋅cos(2πf⋅nTs)Y_c = {1 \over N}\sum\limits_{n = 0}^{N - 1} {f\left[ n \right] \cdot \cos \left( {2\pi f \cdot nT_s } \right)}Yc=N1n=0∑N−1f[n]⋅cos(2πf⋅nTs)
将上述两个公众是的f[n]f\left[ n \right]f[n]代入,并进行化简。对于YsY_sYs可以化简如下:
Ys=1N∑n=0N−12⋅E⋅sin(2πf⋅nTs+θ)⋅sin(2πf⋅nTs)Y_s = {1 \over N}\sum\limits_{n = 0}^{N - 1} {\sqrt 2 \cdot E \cdot \sin \left( {2\pi f \cdot nT_s + \theta } \right) \cdot \sin \left( {2\pi f \cdot nT_s } \right)}Ys=N1n=0∑N−12⋅E⋅sin(2πf⋅nTs+θ)⋅sin(2πf⋅nTs)=2⋅EN∑n=0N−112[cos(θ)−cos(4πf⋅nTs+θ)]= {{\sqrt 2 \cdot E} \over N}\sum\limits_{n = 0}^{N - 1} {{1 \over 2}\left[ {\cos \left( \theta \right) - \cos \left( {4\pi f \cdot nT_s + \theta } \right)} \right]}=N2⋅En=0∑N−121[cos(θ)−cos(4πf⋅nTs+θ)]=E2cos(θ)+E2⋅N∑n=0N−1cos(4πf⋅nTs+θ)= {E \over {\sqrt 2 }}\cos \left( \theta \right) + {E \over {\sqrt 2 \cdot N}}\sum\limits_{n = 0}^{N - 1} {\cos \left( {4\pi f \cdot nT_s + \theta } \right)}=2Ecos(θ)+2⋅NEn=0∑N−1cos(4πf⋅nTs+θ)
根据 等差角度cos,sin序列之和公式
:
∑n=0N−1cos(a+nd)=R⋅cos[a+(N−1)d2]\sum\limits_{n = 0}^{N - 1} {\cos \left( {a + nd} \right)} = R \cdot \cos \left[ {a + \left( {N - 1} \right){d \over 2}} \right]n=0∑N−1cos(a+nd)=R⋅cos[a+(N−1)2d]∑n=0N−1sin(a+nd)=R⋅sin[α+(N−1)d2]\sum\limits_{n = 0}^{N - 1} {\sin \left( {a + nd} \right)} = R \cdot \sin \left[ {\alpha + \left( {N - 1} \right){d \over 2}} \right]n=0∑N−1sin(a+nd)=R⋅sin[α+(N−1)2d]
其中:R=sin(N⋅d2)sin(d2)R = {{\sin \left( {{{N \cdot d} \over 2}} \right)} \over {\sin \left( {{d \over 2}} \right)}}R=sin(2d)sin(2N⋅d)
当sin(d2)=0\sin \left( {{d \over 2}} \right) = 0sin(2d)=0时,R=NR = NR=N。
对于前序列采样公式,d=4πf⋅Tsd = 4\pi f \cdot T_sd=4πf⋅Ts,N=1400。下面假设f=1,对ts取不同值绘制出R对应的曲线:
▲ ts的不同取值下的R的取值,其中N=1400
可以看出,当ts取(0.5)整数倍数的时候,对应的R值很大。当ts偏离(0.5)整数倍数时,R迅速减少。下面给出Ts取0~0.之间的值对应的R的曲线:
▲ ts取值0~0.1之间的数值时对应的R取值
▲ ts取值0.45~0.55之间时,对应的R的取值
从上面分析可以看出,当ts取值范围在0.11f0.1{1 \over f}0.1f1到0.451f0.45{1 \over f}0.45f1之间时,对应正弦波每个周期采样点在2.5~10之间,R的取值都会远远小于N。因此,YsY_sYs就可以近似为:
Ys≈E2cos(θ)Y_s \approx {E \over {\sqrt 2 }}\cos \left( \theta \right)Ys≈2Ecos(θ)
同样的分析可以知道:Yc≈E2sin(θ)Y_c \approx {E \over {\sqrt 2 }}\sin \left( \theta \right)Yc≈2Esin(θ)
然后在计算信号的幅值和相位:
E=2⋅Ys2+Yc2E = \sqrt 2 \cdot\sqrt {Y_s^2 + Y_c^2 }E=2⋅Ys2+Yc2
θ=arctan(YcYs)\theta = \arctan \left( {{{Y_c } \over {Y_s }}} \right)θ=arctan(YsYc)
def data2ap(tt, y, freq):sint = [sin(t * 2 * pi * freq) for t in tt]cost = [cos(t * 2 * pi * freq) for t in tt]ys = sum([a * b for a,b in zip(y, sint)]) / len(tt)yc = sum([a * b for a,b in zip(y, cost)]) / len(tt)amp = sqrt(ys**2 + yc**2)*sqrt(2)phase = atan2(ys,yc)return amp,phase
通过对一个1000Hz测试正弦波的测量,所获得的测量结果为:E=2.4382V, θ=1.3108\theta = 1.3108θ=1.3108。
示波器的设置:时基时间TTb=500μsTb = 500\mu sTb=500μs,幅值尺度Vb=2VV_b = 2VVb=2V。
使用FLUKE45数字万用表的交流档测量该交流信号的有效值为:Ef=2.394VE_f = 2.394VEf=2.394V。可以看出与示波器所测得得到的结果之间有1.85%的误差。
03测量性能
1.测量数值精度
测量100次信号的幅值和相位,统计相应的平均值和方差:
- 幅值:Emean=2.4383V,Var(E)=3.51×10−7E_{mean} = 2.4383V,\,\,\,\,Var\left( E \right) = 3.51 \times 10^{ - 7}Emean=2.4383V,Var(E)=3.51×10−7
- 相位:θmean=1.3112,Var(θ)=6.2134×10−8\theta _{mean} = 1.3112,\,\,Var\left( \theta \right) = 6.2134 \times 10^{ - 8}θmean=1.3112,Var(θ)=6.2134×10−8
▲ 测量100个数据的分布
2.示波器时基长度对测量结果的影响
(1)直接采集信号进行测量
设置示波器的时基参数,可以改变数据采集时间间隔,同时也影响了每帧数据(对于DS6104,数据长度为1400)对应信号的周期个数。
▲ 不同时基下采集到的信号波形
下面显示了对于时基信号从30us变化到5000us过程中,200个信号的幅值和相位的测量数据。
▲ 不同的时下测量幅度和相位的变化
在上面分析由于不同时基的变化,对应采样时间的变化对于测量结果的影响,在ts小于0.002秒之前似乎还能够通过理论分析得到这个结果。但是对于ts大于0.002秒之后的这种变化,感觉就无法从理论上来分析了。
在博文 对示波器测量正弦波幅值和相位仿真实验 中通过理论仿真来验证不同时基ts下的测量误差,可以看到在前面ts大于0.002秒之后的部分的确不在理论可以分析的范围之内。所以这部分误差需要在实际测量中加以避免。
a=[2.13,2.84,2.27,2.40,2.60,2.32,2.45,2.52,2.35,2.47,2.48,2.36,2.48,2.45,2.37,2.48,2.44,2.39,2.48,2.42,2.41,2.48,2.41,2.43,2.47,2.40,2.41,2.45,2.40,2.44,2.43,2.38,2.44,2.44,2.41,2.42,2.41,2.40,2.45,2.42,2.42,2.43,2.40,2.42,2.45,2.42,2.40,2.43,2.40,2.43,2.45,2.38,2.42,2.43,2.41,2.44,2.40,2.39,2.43,2.42,2.42,2.45,2.39,2.40,2.44,2.42,2.43,2.41,2.40,2.42,2.43,2.42,2.40,2.41,2.40,2.43,2.43,2.43,2.41,2.40,2.24,2.40,2.24,2.37,2.24,2.36,2.20,2.37,2.18,2.34,2.42,2.32,2.42,2.32,2.40,2.31,2.40,2.28,2.41,2.27,2.38,2.25,2.38,2.22,2.38,2.21,2.36,2.20,2.34,2.17,2.34,2.41,2.32,2.41,2.30,2.41,2.30,2.39,2.28,2.39,2.25,2.39,2.24,2.37,2.23,2.36,2.20,2.36,2.18,2.34,2.17,2.32,2.41,2.32,2.40,2.30,2.40,2.27,2.40,2.27,2.38,2.25,2.37,2.22,2.38,2.21,2.35,2.19,2.34,2.16,2.34,2.41,2.32,2.41,2.29,2.41,2.29,2.39,2.27,2.39,2.25,2.39,2.23,2.37,2.22,2.36,2.19,2.35,2.17,2.33,2.15,2.32,2.41,2.31,2.40,2.29,2.40,2.27,2.40,2.26,2.38,2.24,2.37,2.21,2.37,2.20,2.35,2.18,2.34,2.15,2.33,2.13,2.31,2.41,2.29,2.41,2.29,2.38,2.26,2.40]
p=[1.06,1.32,1.15,1.20,1.28,1.17,1.25,1.23,1.20,1.26,1.29,1.28,1.26,1.23,1.31,1.28,1.22,1.31,1.28,1.36,1.31,1.28,1.36,1.31,1.28,1.36,1.36,1.32,1.25,1.22,1.46,1.40,1.36,1.31,1.25,1.50,1.45,1.40,1.36,1.31,1.26,1.49,1.45,1.40,1.36,1.30,1.55,1.50,1.45,1.40,1.36,1.59,1.38,1.33,1.28,1.23,1.46,1.41,1.37,1.31,1.25,1.21,1.45,1.39,1.34,1.29,1.24,1.48,1.42,1.37,1.32,1.27,1.51,1.46,1.40,1.36,1.30,1.25,1.49,1.37,1.76,1.40,1.81,1.47,1.85,1.52,1.90,1.56,1.95,1.61,1.28,1.64,1.32,1.71,1.37,1.75,1.42,1.81,1.47,1.85,1.51,1.90,1.56,1.95,1.59,2.00,1.66,2.03,1.71,2.09,1.75,1.41,1.80,1.46,1.83,1.51,1.89,1.54,1.94,1.61,1.99,1.65,2.04,1.70,2.07,1.75,2.14,1.79,2.18,1.84,2.22,1.48,1.14,1.53,1.18,1.56,1.22,1.61,1.27,1.65,1.31,1.69,1.35,1.73,1.39,1.77,1.43,1.84,1.47,1.86,1.52,1.17,1.55,1.22,1.60,1.26,1.64,1.30,1.68,1.35,1.73,1.38,1.77,1.47,1.80,1.46,1.85,1.51,1.89,1.55,1.95,1.59,1.25,1.63,1.29,1.67,1.34,1.71,1.37,1.75,1.42,1.79,1.45,1.86,1.49,1.88,1.54,1.92,1.58,1.96,1.62,2.01,1.66,1.33,1.70,1.36,1.85,1.41,1.79,1.45]
tbase=[0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.01]
(2)通过加窗方法进行测量
对采集数据通过加窗处理,可以避免由于对数据的截断所引起的频率特性的波动,进而提高对信号测量的精度。
常用到的窗口包括有三角波窗口,升余弦(Hanning)窗口等。
▲ 三角窗口与升余弦(Hanning)窗口
下面是统一正弦波在三角窗口和Hanning窗口加权下的波形。可以看出,经过加窗处理之后,数据在在时域内不再会因为截取而不连续了。
▲ 加窗后的数据波形(左:三角窗;右:Haning窗口)
由于加窗对数据幅度有了改变,所以会直接影响测量的结果。下面显示不同加窗方法处理后的幅度和相位。幅度会因为增加三角窗口和Hanning窗口而减少1倍左右,相位测量不受窗口方法的影响。
窗口 | 幅度 | 相位 |
---|---|---|
矩形窗 | 2.4400 | 1.4182 |
三角窗 | 1.2209 | 1.4177 |
汉宁窗 | 1.2226 | 1.4175 |
下面对数据经过加窗处理之后,看示波器的时基大小对测量结果的影响。
A.使用三角窗测量信号的幅度和相位
测量过程中信号经过三角波加窗。
▲ 使用三角窗测量信号的幅值和相位
测量的幅值和相位的结果如下:
▲ 使用三角窗测量信号的幅度和相位
a=[1.71,2.31,2.42,2.36,2.42,2.43,2.41,2.44,2.43,2.42,2.44,2.43,2.43,2.44,2.43,2.43,2.44,2.43,2.43,2.43,2.43,2.44,2.43,2.43,2.44,2.44,2.42,2.43,2.43,2.43,2.42,2.42,2.43,2.43,2.43,2.41,2.42,2.43,2.42,2.43,2.43,2.42,2.42,2.43,2.43,2.43,2.42,2.42,2.43,2.43,2.45,2.41,2.42,2.43,2.43,2.43,2.41,2.42,2.42,2.43,2.43,2.43,2.42,2.42,2.43,2.43,2.43,2.41,2.42,2.42,2.43,2.43,2.41,2.42,2.42,2.43,2.43,2.43,2.42,2.41,2.33,2.40,2.33,2.40,2.32,2.39,2.31,2.39,2.30,2.38,2.41,2.37,2.41,2.37,2.41,2.36,2.41,2.35,2.41,2.34,2.40,2.34,2.40,2.32,2.39,2.31,2.39,2.30,2.38,2.29,2.37,2.41,2.37,2.41,2.36,2.41,2.35,2.41,2.34,2.40,2.33,2.40,2.33,2.40,2.32,2.39,2.31,2.38,2.30,2.38,2.29,2.37,2.41,2.36,2.41,2.36,2.41,2.35,2.41,2.34,2.40,2.33,2.40,2.32,2.39,2.31,2.39,2.30,2.38,2.29,2.37,2.41,2.37,2.41,2.36,2.41,2.35,2.41,2.34,2.41,2.33,2.40,2.32,2.39,2.31,2.39,2.30,2.38,2.29,2.37,2.28,2.37,2.41,2.36,2.41,2.35,2.41,2.34,2.40,2.34,2.40,2.33,2.39,2.32,2.39,2.31,2.38,2.29,2.38,2.28,2.37,2.27,2.36,2.41,2.35,2.41,2.35,2.40,2.34,2.42]
p=[0.66,1.15,1.20,1.17,1.23,1.22,1.24,1.20,1.22,1.25,1.27,1.30,1.25,1.22,1.32,1.27,1.22,1.32,1.27,1.37,1.32,1.27,1.37,1.32,1.27,1.37,1.36,1.31,1.26,1.21,1.45,1.41,1.36,1.31,1.26,1.50,1.45,1.40,1.35,1.31,1.26,1.49,1.45,1.40,1.35,1.31,1.55,1.50,1.45,1.40,1.36,1.59,1.38,1.33,1.28,1.23,1.46,1.41,1.37,1.31,1.26,1.21,1.45,1.39,1.34,1.29,1.24,1.47,1.42,1.37,1.32,1.28,1.51,1.45,1.40,1.36,1.30,1.25,1.49,1.37,1.76,1.40,1.81,1.47,1.85,1.52,1.90,1.56,1.95,1.61,1.28,1.64,1.32,1.71,1.37,1.76,1.42,1.80,1.47,1.85,1.51,1.89,1.56,1.95,1.59,1.99,1.66,2.04,1.70,2.09,1.75,1.42,1.80,1.46,1.82,1.51,1.89,1.54,1.94,1.61,1.99,1.65,2.04,1.70,2.08,1.75,2.13,1.79,2.18,1.84,2.23,1.48,1.14,1.53,1.18,1.56,1.22,1.62,1.27,1.65,1.31,1.69,1.35,1.73,1.39,1.77,1.43,1.84,1.47,1.86,1.52,1.18,1.56,1.22,1.60,1.26,1.64,1.30,1.69,1.35,1.73,1.38,1.77,1.47,1.80,1.46,1.85,1.51,1.89,1.55,1.95,1.59,1.25,1.63,1.29,1.67,1.34,1.71,1.37,1.75,1.42,1.80,1.45,1.86,1.50,1.88,1.54,1.92,1.58,1.96,1.62,2.01,1.66,1.32,1.70,1.36,1.85,1.41,1.79,1.45]
tbase=[0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.01]
B.使用Hanning窗口加权:
▲ 使用Hnning窗测量信号的幅值和相位
▲ 使用Hanning窗测量信号的幅值和相位
a=[1.67,2.19,2.49,2.45,2.43,2.45,2.44,2.44,2.44,2.44,2.44,2.44,2.44,2.44,2.44,2.44,2.44,2.44,2.44,2.44,2.44,2.44,2.44,2.44,2.44,2.44,2.43,2.43,2.43,2.44,2.42,2.43,2.43,2.43,2.44,2.42,2.43,2.43,2.43,2.43,2.44,2.42,2.43,2.43,2.43,2.44,2.42,2.43,2.43,2.43,2.45,2.42,2.43,2.43,2.43,2.43,2.42,2.42,2.43,2.43,2.43,2.44,2.42,2.43,2.43,2.43,2.43,2.42,2.43,2.43,2.43,2.43,2.42,2.42,2.43,2.43,2.43,2.43,2.42,2.41,2.34,2.40,2.35,2.40,2.34,2.40,2.33,2.39,2.33,2.39,2.42,2.38,2.42,2.38,2.41,2.37,2.41,2.36,2.42,2.36,2.41,2.36,2.40,2.35,2.40,2.34,2.40,2.33,2.39,2.32,2.38,2.42,2.38,2.42,2.38,2.41,2.37,2.41,2.36,2.41,2.35,2.40,2.35,2.40,2.34,2.40,2.33,2.39,2.32,2.39,2.32,2.38,2.42,2.38,2.41,2.37,2.41,2.37,2.41,2.36,2.41,2.35,2.40,2.34,2.40,2.34,2.40,2.33,2.39,2.32,2.39,2.42,2.38,2.42,2.37,2.41,2.37,2.41,2.36,2.41,2.35,2.40,2.35,2.40,2.34,2.40,2.33,2.39,2.32,2.39,2.31,2.38,2.42,2.38,2.41,2.37,2.41,2.36,2.41,2.36,2.40,2.35,2.40,2.34,2.40,2.33,2.39,2.32,2.39,2.31,2.38,2.30,2.38,2.41,2.37,2.41,2.37,2.41,2.36,2.42]
p=[0.54,1.09,1.22,1.21,1.23,1.22,1.25,1.20,1.23,1.25,1.27,1.30,1.25,1.22,1.32,1.27,1.22,1.32,1.27,1.37,1.32,1.27,1.37,1.32,1.27,1.37,1.36,1.31,1.26,1.21,1.45,1.41,1.36,1.31,1.26,1.50,1.45,1.40,1.35,1.31,1.26,1.49,1.45,1.40,1.35,1.31,1.55,1.50,1.45,1.40,1.36,1.59,1.38,1.33,1.28,1.23,1.46,1.41,1.37,1.31,1.26,1.21,1.45,1.39,1.34,1.29,1.24,1.47,1.42,1.37,1.32,1.28,1.51,1.45,1.40,1.36,1.30,1.25,1.49,1.37,1.76,1.40,1.81,1.47,1.85,1.52,1.90,1.56,1.95,1.61,1.28,1.64,1.32,1.71,1.37,1.75,1.42,1.80,1.47,1.85,1.51,1.89,1.56,1.95,1.59,1.99,1.66,2.04,1.70,2.09,1.75,1.42,1.80,1.46,1.82,1.51,1.89,1.54,1.94,1.61,1.99,1.65,2.04,1.70,2.08,1.75,2.13,1.79,2.18,1.84,2.23,1.48,1.14,1.53,1.18,1.56,1.22,1.62,1.27,1.65,1.31,1.69,1.35,1.73,1.39,1.77,1.43,1.84,1.47,1.86,1.52,1.18,1.56,1.22,1.60,1.26,1.64,1.30,1.69,1.35,1.73,1.38,1.77,1.47,1.80,1.46,1.85,1.51,1.89,1.55,1.95,1.59,1.25,1.63,1.29,1.67,1.34,1.71,1.37,1.75,1.42,1.79,1.45,1.86,1.50,1.88,1.54,1.92,1.58,1.96,1.62,2.01,1.66,1.32,1.70,1.36,1.85,1.40,1.79,1.45]
tbase=[0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.00,0.01]
C.对比矩形窗、三角窗、Hanning窗的效果
通过对比可以看出,经过加窗处理之后,时基长度小于2ms之内,测量的结果收到时基长度的影响就大大减弱了。
▲ 对比矩形窗、三角窗、Hanning窗对于幅度测量的影响
特别是在时基设置为0.25~0.7毫秒左右(对应1000Hz的信号,示波器显示大约3.5至10个周期)使用Hanning窗口所测量的幅度谱基本上不受时基长短的 影响。
▲ 对比矩形窗、三角窗、Hanning窗对于幅度测量的影响(局部)
#------------------------------------------------------------
def winddata(len, mode = 0):if mode == 0:return [1] * lenelif mode == 1: # Triangle Windowwin = [1] * lenfor i in range(len):if i < len / 2:win[i] = i / ((len + 1) / 2)else: win[i] = (len-i) / ((len + 1) / 2)return winelif mode == 2:win = [1] * lenfor i in range(len):theta = (i - (len - 1) / 2) * pi * 2 / lenwin[i] = (1 + cos(theta)) / 2return winelse: return [1] * len
#------------------------------------------------------------
def data2ap(tt, y, freq):sint = [sin(t * 2 * pi * freq) for t in tt]cost = [cos(t * 2 * pi * freq) for t in tt]ys = sum([a * b for a,b in zip(y, sint)]) / len(tt)yc = sum([a * b for a,b in zip(y, cost)]) / len(tt)amp = sqrt(ys**2 + yc**2)*sqrt(2)phase = atan2(ys,yc)return amp,phase
#------------------------------------------------------------
def ds6104data(freq, win=0):x,y = ds6104readcal(1)if win == 1:wind = winddata(len(x), 1)y = [d * w * 2 for d,w in zip(y, wind)]elif win == 2:wind = winddata(len(x), 2)y = [d*w*2 for d,w in zip(y, wind)]a,p = data2ap(x,y,freq)return a, p, x, y
04示波器采样模式对测量精度的影响
示波器采样模式包括有四种:正常采样模式、平均采样模式、峰值采样模式以及高分辨率采样模式。下面通过对比实验,查看这些采样模式对于测量精度的影响。、
下面实验中,采样信号为1000Hz,有效值为2.3V左右的正弦波信号。示波器的时基设为0.7ms,处理数据使用Hanning窗加窗。
1.正常采样模式
测量结果:
幅值:2.44283V;相位:1.29056
▲ 正常采样模式下的信号的波形
2.平均采样模式
测量结果:
幅值:2.44128;相位:1.2282
▲ 平均采样模式下信号的波形
3.峰值采样模式
测量结果:
幅值:2.44415V;相位:1.2903
▲ 峰值采样模式下的信号波形
4.高分辨率采样模式
测量结果:
幅值:2.4243V;相位:1.4068
▲ 高分辨率采样模式下的信号波形
对比以上四种示波器采样模式,对于信号的幅度没有太多的影响,对于信号的相位,高分辨率采样则出现了与其它三类较大的差异。根据示波器实际波形显示,采用高分辨率采样模式所得到的结果与实际相位是最接近的。
采样模式 | 幅值 | 相位 | |
---|---|---|---|
1 | 正常采样 | 2.44283 | 1.29056 |
2 | 平均采样 | 2.44128 | 1.2282 |
3 | 峰值采样 | 2.44415 | 1.2903 |
4 | 高分辨率采样 | 2.4243 | 1.4068 |
05不同频率对测量结果的影响
对于信号频率从10Hz变化到10MHz,使用相同的方法进行测量它的幅值和相位,结果如下:
频率 | 幅值 | 相位 |
---|---|---|
10.0 | 2.3879 | 1.3380 |
100.0 | 2.4073 | 1.2545 |
1000.0 | 2.4385 | 1.3111 |
10000.0 | 2.5401 | 1.2231 |
100000.0 | 2.5516 | 1.2740 |
1000000.0 | 2.5952 | 1.0833 |
10000000.0 | 2.3971 | -0.6503 |
#------------------------------------------------------------
def ds6104settimebase(freq, period):timebase = period / freq / 14ds6104timebase(timebase)
06结论
通过对DS6104联网示波器采集到的数据进行加窗计算,可以获得被测量信号的幅值和相位。这种方法所得到的数值精度在1000Hz时精度很高。同样,这种方法所能适应的频率范围可以是示波器的有效显示频率范围。
为了提高测量精度,示波器在设置参数的需要尽可能满足:
- 设置示波器通道增益,尽可能使得信号波形达到屏幕满幅值,这样可以最大利用示波器采样8bit的精度;
- 设置示波器时基,使得能够采样到7个周期左右的数据;
- 示波器的采样模式使用“高分辨率”采样模式;
在处理数据时,建议使用增加Hanning窗口,来减少对数据截断所带来的幅值波动。
使用数字示波器DS6104测量交流信号的幅值和相位相关推荐
- 不用傅里叶变换,提取某一频率的幅值和相位
不用傅里叶变换,提取某一频率的幅值和相位 摘要: 本文从实际工程问题入手,探寻解决办法,为引入信号正交分解,和广义傅里叶级数做铺垫. 转子做周期性旋转时,不平衡质量所产生的周期性惯性离心力会引起转子产 ...
- matlab复数的相位,复数的幅值和相位
[i]); x[i]=x[i]*180/PI-90; printf("第%d次谐波的相位为为%f \\n", i , x[i]) ; } printf("第%d次谐波的幅 ...
- matlab 求电流幅值,输出信号的幅值与相位.ppt
输出信号的幅值与相位 1.什么叫频率分析?2.频率分析的问题引入?3.频率分析的仿真.4.频率分析的计算方法.5.频率分析的作用与意义6.幅频特性与相频特性 王选择 1 频率分析的定义 1 频率分析的 ...
- fft之后求模值和相位_50Hz交流信号经ADC在一个周期采样有限个点后,怎么用FFT变换求得有效值、幅值和相位等?...
FFT是离散傅立叶变换的快速算法,可以将一个信号变换 到频域.有些信号在时域上是很难看出什么特征的,但是如 果变换到频域之后,就很容易看出特征了.这就是很多信号 分析采用FFT变换的原因.另外,FFT ...
- 幅值测量c语言程序,幅值和电平测量程序块
幅值和电平测量程序块可以通过LabVIEW提供的测量Ⅵ的测量得到周期性波形的幅值和电平大小.幅值可 以通过高低电平位置相减得到.这3个测量参数同样反映出周期性波形的基本特性,在周期性波形的测量 中常常 ...
- matlab幅值解调,DSB解调后信号的幅值与信号的幅值为啥差别那么大
做了一个DSB调制和解调程序,在信号相干过滤波器的时候,我自己写了一个先将它们频域相乘再逆FFT得到滤波后的信号,但是它的幅度比我预期的小一半左右,然后我又使用filter函数,直接对时域信号进行滤波 ...
- 脉冲信号是数字信号吗?为什么时间离散的信号在幅值上可以是连续的?
先来看看信号的分类: 连续信号=连续时间信号, 分为以下两类: 时间连续, 幅值连续. 定义为模拟信号. 时间连续, 幅值离散. 比如脉冲信号. 离散信号=离散时间信号, 分为以下两类: 时间离散 ...
- 幅值、相位与强度之间的关系
"相位"(phase)是描述周期信号中某点在特定时刻位于循环中的相对位置的标度,其是描述信号波形变化的度量,通常以弧度(角度)作为单位,也称作相角. 在波动光学中,单色相干光波场由 ...
- 傅里叶变换的理解----计算幅值和相位
先推荐表格文章 如果看了此文你还不懂傅里叶变换,那就过来掐死我吧[完整版] 傅里叶原理表明:任何连续测量的时序或信号,都可以表示为不同频率的正弦波或余弦信号的无限叠加. 1.FT的理论就会告诉你可以 ...
最新文章
- java spark 读取json_apache-spark - 与aws-java-sdk链接时,在读取json文件时发生Spark崩溃 - 堆栈内存溢出...
- Java日常错误及需要注意细节,持续更新......
- 分布式存储绝不简单 —— UCan下午茶-武汉站纪实
- Python 语言程序设计(3-1)字符串处理函数和相关功能
- 写于Silverlight整装待发之际【瞿杰】
- Ubuntu 安装deb格式的安装包(收集)
- mongodb,Mysql,redis基础教程
- Java开发不得不会!java私塾初级模拟银源代码
- Androidstudio svn 使用 版本控制 详细步骤(一)
- matlab绘制符号函数的ezplot函数
- mac上 网易mumu模拟器在打开代理后无法使用键盘 解决办法
- 计算机制作贺卡教案,WORD电子贺卡制作的优秀教案_电子贺卡
- Kafka性能监控与优化
- ANS.1的基础总结
- 陈睿竟也是个“大龄二次元宅”?还不知道的集美看过来
- msgbox php,HTML_MsgBox函数语言参考,描述 在对话框中显示消息 - phpStudy
- VC6.0实现网络编程弹出一个消息框
- android hdmi解码器,如何使用HDMI高清编码器进行直播
- 学生学籍管理系统_登陆界面设计
- 新用户启动Anaconda时,一直卡在loading applications