一种全角度无限旋转单电位器及高精度风速检测方法

1.本发明涉及气象事业中的风速测量技术领域，提供了一种全角度无限旋转单电位器及高精度风速检测方法。


背景技术：

2.空气处于不停的运动之中，空气运动产生的气流称为风，风的运动既有速度又有方向，它在天气预报中有着重要的作用，常作为良好的预报指标而被广泛应用。目前，广泛应用于气象传感器数字化网络化系统中的风速采集器，对测量出的风速数据的精度和分辨率有较高要求，而且由于应用于气象服务中，风速采集器一般安装于户外，要求传输距离远，测量精度高。按照气象台站要求，瞬时风速即0.25秒的平均风速。
3.测风仪器在我国所有的气象部门中占有不可或缺的地位，然而针对绝大多数的气象单位来说，对于测风仪器的选择仍然是传统的机械旋转式测风传感器，这些传感器就是我们平时生活中常见的风杯风速传感器和单翼风速传感器，以此来进行对于风速和风向的测量。目前风速测量通常采用三杯式光电型风速传感器，其主要工作原理是：当风杯转动时，会带动安置于光电耦合器的发光二极管和光电三极管之间的多齿截光盘转动，发光二极管发射的光束被齿盘上的齿度切割，就会使下面的光敏三极管有时接收到上面的发光二极管发射的光线而导通，有时接收不到上面发光二极管照射来的光线而截止，这样，光电三极管即产生脉冲输出，每个齿遮住光束时间表现为低电平，光束通过两齿盘间隙之间的时间表现为高电平。转动轴转动一圈，多次切割光束而输出一高一低的脉冲信号，因此风速传感器输出的均是频率与风速成正比的脉冲信号，输出信号：1～1221hz方波信号，所以，测出信号频率就可得到风速。
4.风杯式风速传感器，是一种十分常见的风速传感器，最早由英国鲁宾孙发明。感应部分是由三个或四个圆锥形或半球形的空杯组成。空心杯壳固定在互成120
°
的三叉星形支架上或互成90
°
的十字形支架上，杯的凹面顺着一个方向排列，整个横臂架则固定在一根垂直的旋转轴上。风杯1的压力在垂直于风杯轴方向上的分力近似为零。风杯2与3同风向成60度角相交，对风杯2而言，其凹面迎着风，承受的风压最大；风杯3其凸面迎风，风的绕流作用使其所受风压比风杯2小，由于风杯2与风杯3在垂直于风杯轴方向上的压力差，而使风杯开始顺时针方向旋转，风速越大，起始的压力差越大，产生的加速度越大，风杯转动越快。
5.当风杯转动时，带动同轴的多齿截光盘或磁棒转动，通过电路得到与风杯转速成正比的脉冲信号，该脉冲信号由计数器计数，经换算后就能得出实际风速值。
6.因此，现有风向风速指标分别为：
7.风速：测量范围：0～60m/s；测量精度：
±
0.3m/s；分辨率：0.1m/s
8.风向：测量范围：0～360
°
；测量精度：
±3°
；分辨率：3
°
。


技术实现要素：

9.本发明的目的在于解决现有技术是通过脉冲计数来测量风速时存在的误差较大
的问题。
10.为了实现上述目的本发明提供了如下方案：
11.本发明提供了一种全角度无限旋转单电位器，包括底座、顶盖，其特征在于，在底座上设置有内碳膜圈和外碳膜圈，内碳膜圈设置有内碳膜圈开口，外碳膜圈设置有外碳膜圈开口，内碳膜圈开口与外碳膜圈开口错开，电极座中心设置有转轴，转轴上设置有电极片，电极片包括与内碳膜圈接触的内碳膜圈触点和与外碳膜圈接触的外碳膜圈触点。
12.上述技术方案中，外碳膜圈开口处的一端为外碳膜圈电极引出端，外碳膜圈电极引出端通过电极引出，另一端为外碳膜圈跳线端；
13.内碳膜圈开口处的一端为内碳膜圈电极引出端，内碳膜圈电极引出端通过电极引出，另一端为内碳膜圈跳线端；
14.外碳膜圈跳线端与内碳膜圈跳线端通过跳线电连接。
15.上述技术方案中，内碳膜圈跳线端与外碳膜圈跳线端的错位角度为θ，θ即电器角度相位差。
16.上述技术方案中，电位器旋转角度α为：
[0017][0018]
其中vc为电极片的电压，α为电位器旋转角度，i为恒流源电流，r1内碳膜的电阻、r2为外碳膜的电阻，γ1为内碳膜圈开口应的机械角度、γ2为外碳膜圈开口应的机械角度，θ为内碳膜圈跳线端与外碳膜圈跳线端的错位角度，γ1＜θ≤π且γ2＜θ≤π。
[0019]
本发明还提供了一种高精度风速检测方法，包括以下步骤：
[0020]
步骤1、cpu中央处理器启动adc每毫秒转换一次数据，得到时长m1毫秒，连续n次电位器中心抽头的毫秒电压瞬时数据，并存储为瞬时电压v1、v2、
……
、v4、vn，判断瞬时电压vn是否经过拐点并进行数据处理得到电位器旋转角度αn：
[0021]
步骤2：间隔时间m1毫秒，重复第1步n 1次得到集合{α1，α2，
……
，αn，α
n 1
}，其中n大于等于4；
[0022]
步骤3：采用步骤2得到的集合{α1，α2，
……
，αn，α
n 1
}计算角速度存储角速度集合{ω1，ω2，......ωn}；
[0023]
第7步：构建插值表；
[0024]
通过风洞可以测试出风速与角速度的对应数据，构成一个风速角速度查值表，
[0025]
得到风速值数据集{v
b1
、v
b2
、λλ、v
b14
、v
bn
}与对应的角速度值数据集{ω
b1
、ω
b2
、λλ、ω
b14
、ω
bn
}；
[0026]
第8步：插值算法计算风速。在插值表给出的数据中，将相邻两点间的风速看成线性增加，将前面计算存储的角速度值集合{ω
b1
、ω
b2
、λλ、ω
b14
、ω
bn
}分别按照下面计算处理：
[0027]
判断ωn值：ωn等于{ω
b1
、ω
b2
、λλ、ω
b14
、ω
bn
}中的某个值ω
bn
，则vn＝v
bn
，否则，ωn将落在{ω
b1
、ω
b2
、λλ、ω
b14
、ω
bn
}中的两个数之间，即ω
bn
＜ωn＜ω
bn-1
，则按如下公式
计算：vn为瞬时风速，得到瞬时风速集合{v1，v2，
……
，vn}，平均风速
[0028]
上述技术方案中，判断瞬时电压vn是否经过拐点包括以下步骤：
[0029]
a、电极片的内碳膜圈触点顺时针经过内碳膜圈跳线端：
[0030]
前n-1个瞬时电压值连续增加或相邻值相等，且瞬时电压值小于或等于vs，第n个值接近且小于或等于v
l
，处理办法：第n个值用vs替换；
[0031]
或者前n-2个瞬时电压值连续增加或相邻值相等，瞬时电压值接近且小于或等于vs，第n-2，n-1个瞬时电压值接近且小于或等于v
l
，且第n个值小于等于第n-1个值，处理办法：第n-1和第n个值用vs替换；
[0032]
或者第n-4个和第n-3个瞬时电压值依次增大或相等，瞬时电压值接近且小于或等于vs，第n-2、n-1、n个值依次减小或相等，瞬时电压值接近且小于或等于v
l
，处理办法：第n-4个和第n-3个值均用v
l
替换；
[0033]
或者第n-4个值接近或等于vs，后第n-3，n-2，n-1，n个值依次减小或相邻值相等，接近且小于或等于v
l
，处理办法：第n-4个值用v
l
替换。
[0034]
前面四种情形可判断内碳膜圈触点顺时针经过内碳膜圈跳线端，其中vs＝ir1，
[0035]
b、电极片的内碳膜圈触点逆时针经过内碳膜圈跳线端：
[0036]
第n-4，n-3，n-2，n-1个瞬时电压值连续增加或相邻值相等，接近且小于或等于v
l
，第n个值接近且小于或等于vs，处理办法：第n个值用v
l
替换；
[0037]
或者第n-4，n-3，n-2个瞬时电压值连续增加或相邻值相等，瞬时电压值接近且小于或等于v
l
，第n-1、n个值接近且小于或等于vs，第n个值小于等于第n-1个值，处理办法：第n-1和第n个值用v
l
替换；
[0038]
或者第n-4，n-3个瞬时电压值依次增大或相等，瞬时电压值接近且小于或等于v
l
，第n-2，n-1、n个值依次减小或相邻值相等，瞬时电压值接近且小于或等于vs，处理办法：第n-4，n-3个值均用vs替换；
[0039]
或者第n-4个值接近且小于或等于v
l
，后n-3、n-2、n-1、n个值依次减小或相邻值相等，瞬时电压值接近且小于或等于vs，处理办法：第n-4个值用vs替换。
[0040]
前面四种情形可判断内碳膜圈触点逆时针经过内碳膜圈跳线端。
[0041]
c.内碳膜圈触点顺时针经过内碳膜圈电极引出端：
[0042]
第n-4，n-3，n-2，n-1个瞬时电压值连续减小或相邻值相等，且瞬时电压值接近且大于或等于vm，第n个值接近或等于0，处理办法：第n个值用vm替换，
[0043]
或者前第n-4，n-3，n-2个瞬时电压值连续减小或相邻值相等，且瞬时电压值接近且大于或等于vm，第n-1、n个瞬时电压值接近或等于0，且第n个值大于或等于第n-1个值，处理办法：第n-1和第n个值用vm替换；
[0044]
或者第n-4，n-3个瞬时电压值依次减小或相等，接近且大于或等于vm，第n-2、n-1、
n个瞬时电压值依次增加或相邻值相等，并接近或等于0，处理办法：第n-4，n-3个瞬时电压值均用0替换；
[0045]
或者第n-4个瞬时电压值接近且大于或等于vm，后第n-3，n-2，n-1，n个值依次增加或相邻值相等，并接近或等于0，处理办法：第n-4个值用0替换。
[0046]
前面四种情形可判断内碳膜圈触点顺时针经过内碳膜圈电极引出端。
[0047]
d.内碳膜圈触点逆时针经过内碳膜圈电极引出端。
[0048]
第n-4，n-3，n-2，n-1个瞬时电压值连续减小或相邻值相等，且接近或等于0，第n个值接近且大于或等于vm，处理办法：第n个值用0替换；
[0049]
或者第n-4，n-3，n-2个瞬时电压值连续减小或相邻值相等，且接近或等于0，第n-1、n个值接近且大于或等于vm，且第n个值大于或等于第n-1个值，处理办法：第n-1和第n个值用0替换；
[0050]
或者第n-4，n-3个值依次减小或相等并接近0，第n-2，n-1，n个值依次增加或相邻值相等，接近且大于或等于vm，处理办法：第n-4、n-3个值均用vm替换；
[0051]
或者第n-4个值接近或等于0，第n-3，n-2，n-1，n个值依次增加或相邻值相等，接近且大于或等于vm，处理办法：第n-4个值用vm替换。
[0052]
前面四种情形可判断内碳膜圈触点逆时针经过内碳膜圈电极引出端。
[0053]
上述技术方案中，进行数据处理得到电位器旋转角度αn包括以下步骤：
[0054]
经过拐点判断处理后，n个值的相邻电压值是否有超过m2毫伏的，若有去掉该值，插值代替，即用该值前后值的均值代替该值，若最后一个值vn超过m2伏，则用2v
n-1-v
n-2
代替vn，去掉最大和最小两个值，求剩下的值的算术平均值，得到的均值即为m1毫秒电压瞬时值v
ss
，判断如0≤v
ss
≤vs，否则存储αn。
[0055]
本发明还提供了一种全角度无限旋转单电位器的参数自动检测方法，其特征在于，内碳膜圈开口应的机械角度γ1，外碳膜圈开口应的机械角度γ2，内碳膜圈跳线端与外碳膜圈跳线端的错位角度θ的测量方法，包括以下步骤：通过plc程控步进电机(步进角度0.01
°
)带动满角度电位器向一个方向旋转(避免回差)，其旋转角度可控可测，每旋转一次停留3秒钟，待电路稳定，输出一组角度电压数据；
[0056]
第1步：初测，控步进电机带动电位器每次转动1
°
，测出v
s1
、v
l1
、v
m1
电位器对应的角度，α
s1
、α
l1
、α
m1
；
[0057]
具体方法：步进电机带动全角度无限旋转单电位器向一个方向旋转，每次旋转1
°
，每旋转一次停留3秒钟，待电路稳定再检测，电位器从0
°
开始旋转359
°
，记录一组360个电压角度数据，电压角度数据一一对应；
[0058]
寻找判断v
s1
、v
l1
、v
m1
对应的α
s1
、α
l1
、α
m1
，方法：
[0059]
步进电机旋转，检测到0电压时，开始计算步进电机的步进次数，步进359次，记录0
°
到359
°
对应的电压值，构成电压角度表，表中找到最大电压值即为v
l1
，该值对应的角度即为α
l1
，比α
l1
小1
°
的即为α
s1
，对应的电压值即为v
s1
，α
m1
＝359
°
，对应的电压即为v
m1
，则有：
[0060][0061]
其中γ
11
对应为内碳膜圈开口应的机械角度γ1，θ1对应为内碳膜圈跳线端与外碳
膜圈跳线端的错位角度θ，γ
21
对应为外碳膜圈开口应的机械角度；
[0062]
第2步：细测，控步进电机带动电位器每次转动0.01
°
，测出v
s2
、v
l2
、v
m2
电位器对应的角度，α
s2
、α
l2
、α
m2
；
[0063]
具体方法：步进电机带动满角度电位器向一个方向快速旋转，每次旋转0.01
°
，检测到电压为0时，停留3秒钟，确保电路稳定，再次检测确认电压为0，启动步进电机开始计数，快速转动到α
s1-2
°
，开始3秒钟步进一次，每次步进0.01
°
，直到α
s1
1
°
，记录该3
°
的电压角度数据300个，构成电压角度表，表中找到最大电压值即为v
l2
，该值对应的角度即为α
l2
，比α
l2
小0.01
°
的即为α
s2
，对应的电压值即为v
s2
，继续快速旋转到359.99
°
，停留检测输出电压值，该值即v
m2
，α
m2
＝359.99
°
，此步得到的电压和角度精度高，代入前面公式可计算出得到：
[0064]
γ1＝γ
12
＝2π-α
s2
1，
[0065][0066][0067]
因为本发明采用上述技术方案，因此具备有以下有益效果：
[0068]
1、本发明通过无限旋转的电位器配合风速测量方法，能够精确的测量到风速。
[0069]
2、本发明采用上电极片和下电极片，分别与内碳膜圈和外碳膜圈滑动相接，可无限圈旋转，且不存在测量死角。
[0070]
3、本发明提供的方法能够进行数据质量控制，由于电路受到雷击、闪电和强磁场干扰，这些干扰虽然短暂，但也势必造成产品质量下降，我们采用数据质量控制后的高精度角度校测方法来处理。解决雷击、闪电和强磁场干扰。
附图说明
[0071]
图1为一种全角度无限旋转单电位器示意图；
[0072]
图2为一种全角度无限旋转单电位器顶盖剖掉部分后示意图；
[0073]
图3为一种全角度无限旋转单电位器俯视方向透视图；
[0074]
图4为一种全角度无限旋转单电位器等效原理图；
[0075]
图5为数据坐标图。
[0076]
附图标记说明：
[0077]
1-底座、2-顶盖，3-转轴，4-电极片，5-外碳膜圈，5-1外碳膜圈开口，6-内碳膜圈6-1-内碳膜圈开口，a-外碳膜圈跳线端，b-外碳膜圈电极引出端，c-内碳膜圈电极引出端，d-内碳膜圈跳线端。
[0078]
其中vc为电极片的电压，α为电位器旋转角度，i为恒流源电流，r1内两碳膜的电阻、r2为外碳膜的电阻，γ1为内碳膜圈开口应的机械角度、γ2为外碳膜圈开口应的机械角度，θ为内碳膜圈跳线端d与外碳膜圈跳线端a的错位角度。
具体实施方式
[0079]
本发明涉及气象事业中的风速测量技术领域，通过引入全角度电位器测量方法可
以很好的改善风速风向的测量精度，实现了高精度风速风向的检测。
[0080]
全角度无限旋转单电位器(简称全角度电位器)如图4，可用于风速、风向、风力发电等0-360度转动的旋转角度检测和转速检测场合。
[0081]
说明：两导电膜半径内小外大，其电阻也不一样，四个端分别为a、b和d、c，ad端通过导线相连，恒流源接到b端，经过第一个导电膜电流流到a端，再经过第二个导电膜流到c端，c端接地。
[0082]
所述全角度电位器碳膜片内设有两个独立的内碳膜圈和外碳膜圈，内碳膜圈和外碳膜圈半径不一样电阻也不一样，碳膜片用导线将其首尾相连，余下两端作为电位器两个连接端；所述设计一个滑动接触片分内外接触头，分别与内碳膜圈和外碳膜圈滑动相接，可无限圈旋转，滑片导线引出为电位器的中心抽头端；所述电位器上下碳膜圈电器角度错位，相位差最大180
°
，确保滑动接触片上下接触头旋转到任何位置至少有一个接触头接触到碳膜圈。两导电膜四个端分别为a、b和c、d，ad端通过导线相连，给电位器加一个恒流源，恒流源从b端流入，经第一个导电膜电流流到a端，再经过第二个导电膜流到c端，c端接地。
[0083][0084][0085]
其中，v
c-电位器中心抽头电压，α-电位器旋转角度，i-恒流源电流，r1、r
2-分别是两碳膜的电阻(通常r1小于r2)，γ1、γ
2-分别是两碳膜圈缺口对应的机械角度，一般不相等，θ-相位差，电位器滑动触头从内圈a端开始顺时针旋转到外圈a端的角度，γ1＜θ≤π且γ2＜θ≤π。本算法结合图4，更方便理解，解出得到α，则有：
[0086][0087]
针对已选定全角度电位器和设计调试好的电路，i、γ1、γ1、r2、r2、θ均为已知数，通过有关方法可以测得。令vs＝ir1，，则有
[0088][0089]
瞬时角速度公式：
[0090][0091]
ω一瞬时角速度，通过软件处理可以计算出瞬时角速度。基本方法：通过前面的计算得到α，每5毫秒取样一个值，测出电压值，计算出α，单位时间(单位时间取250毫秒，则有50个α值的)内的α构成一个数组，取某时刻相邻前后两时刻点的角度差即为δα，δt＝0.01秒，则可计算出ω。
[0092]
瞬时角速度公式：
[0093][0094]
λ-瞬时角加速度，通过软件处理可以计算出瞬时加角速度。基本方法：通过前面的计算得到ω，每50毫秒秒计算出一个ω，单位时间(单位时间取250毫秒，则有50个ω值的)内的ω构成一个数组，取某时刻相邻前后两时刻点的角度差即为δω，δt＝0.01秒，则可计算出λ。
[0095]
由(2)可知，当α＝2π-γ1，
[0096]vs
＝ir1ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0097]
γ1＝2π-α
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0098][0099]
当α从大变小，无限接近2π-γ1，即α＝2π-γ1 1/∞≈2π-γ1，∞-无穷大，则有无穷大，则有
[0100]
物理意义：电位器内圈从0
°
旋转到2π-γ1时，再转动无穷小，实际电路中考虑为步进电机转动最小步进(比如本测试电路考虑步进电机步进为0.01
°
，1/∞＝0.01/360≈0.000028，单位为弧度)。
[0101]
当α从小变大，无限接近2π，即α＝2π-1/∞≈2π，则有
[0102][0103]
由(10)和(11)解出由(10)和(11)解出
[0104]
风速检测方法：现代气象业务要求提供的瞬时风速为0.25秒平均风速。
[0105]
(风速检测：风速传感器工作原理风速传感器的感应元件为三杯式风杯组件和全角度无限旋转单电位器组成。风吹动风杯旋转，通过轴带动全角度无限旋转单电位器旋转，风速与单位时间转动的角度成单调增加。测量单位时间内转动的角度，可计算风杯的角速度，风速与传感器风杯的角速度成线性单调增加，但不是线性增加，角速度需处理后才能得到高精度风速。
[0106]
数据质量控制：由于电路受到雷击、闪电和强磁场干扰，这些干扰虽然短暂，但也势必造成产品质量下降，我们采用数据质量控制的办法来处理。雷击、闪电和强磁场干扰是瞬间发射，影响的数据是极个别的，根据各机械旋转，其角速度和角度变化不能超过一定值，针对大多数实际情况，假定每秒钟旋转小于一周(特殊业务超过一周，可调整标准)，则有每毫秒角度变化不超过一定值，设计ad参考电压3.3伏，对应电压约10毫伏，取10毫伏作为标准，即相邻两数据差不超过10毫伏，若超过即可判断有干扰。具体方法如下：
[0107]
角度、角速度、角加速度及风速数据处理方法如下：
[0108]
第1步：电压数据质量控制与风杯角度计算：cpu中央处理器启动adc每毫秒转换一次数据，得到时长5毫秒，连续5次电位器中心抽头的毫秒电压瞬时数据，并存储为v1、v2、v3、
v4、v5。
[0109]
(1)判断是否经过拐点及数据处理：
[0110]
a.滑动触头顺时针经过内圈d点。前4个值连续增加(或相邻值相等)接近且小于或等于vs，第5个值接近且小于或等于v
l
，处理办法：第5个值用vs替换；
[0111]
或者前3个值连续增加(或相邻值相等)接近且小于或等于vs，第4、5个值接近且小于或等于v
l
，且第5个值小于等于第4个值，处理办法：第4和第5个值用vs替换；
[0112]
或者前2个值依次增大(或相等)且小于或等于vs，第3、4、5个值依次减小(或相等)接近且小于或等于v
l
，处理办法：第1、2个值均用v
l
替换；
[0113]
或者第1个值接近或等于vs，后4个值依次减小(或相邻值相等)并接近且小于或等于v
l
，处理办法：第1个值用v
l
替换。前面四种情形可判断滑动触头顺时针经过内圈a点(内圈a拐点)。
[0114]
b.滑动触头逆时针经过内圈d点。前4个值连续增加(或相邻值相等)并接近且小于或等于v
l
，第5个值接近或等于vs，处理办法：第5个值用v
l
替换：
[0115]
或者前3个值连续增加(或相邻值相等)并接近且小于或等于v
l
，第4、5个值接近且小于或等于vs，第5个值小于等于第4个值，处理办法：第4和第5个值用v
l
替换；
[0116]
或者前2个值依次增大(或相等)并接近且小于或等于v
l
，第3、4、5个值依次减小(或相邻值相等)并接近且小于或等于vs，处理办法：第1、2个值均用vs替换；
[0117]
或者第1个值接近且小于或等于v
l
，后4个值依次减小(或相邻值相等)并接近且小于或等于vs，处理办法：第1个值用vs替换。
[0118]
前面四种情形可判断滑动触头逆时针经过内圈c点(内圈c拐点)。
[0119]
c.滑动触头顺时针经过内圈c点。前4个值连续减小(或相邻值相等)并接近大于或等于vm，第5个值接近或等于0，处理办法：第5个值用vm替换；
[0120]
或者前3个值连续减小(或相邻值相等)并接近大于或等于vm，第4、5个值接近或等于0，且第5个值大于或等于第4个值，处理办法：第4和第5个值用vm替换；
[0121]
或者前2个值依次减小(或相等)并接近vm，第3、4、5个值依次增加(或相邻值相等)并接近或等于0，处理办法：第1、2个值均用0替换；
[0122]
或者第1个值接近或等于vm，后4个值依次增加(或相邻值相等)并接近或等于0，处理办法：第1个值用0替换。
[0123]
前面四种情形可判断滑动触头顺时针经过内圈c点(内圈c拐点)。
[0124]
d.滑动触头逆时针经过内圈c点。前4个值连续减小(或相邻值相等)并接近或等于0，第5个值接近且大于或等于vm，处理办法：第5个值用0替换；
[0125]
或者前3个值连续减小(或相邻值相等)并接近或等于0，第4、5个值接近且大于或等于vm，且第5个值大于或等于第4个值，处理办法：第4和第5个值用0替换；
[0126]
或者前2个值依次减小(或相等)并接近0，第3、4、5个值依次增加(或相邻值相等)接近且大于或等于vm，处理办法：第1、2个值均用vm替换；
[0127]
或者第1个值接近或等于0，后4个值依次增加(或相邻值相等)接近且大于或等于vm，处理办法：第1个值用vm替换。
[0128]
前面四种情形可判断滑动触头逆时针经过内圈a点(内圈c拐点)。
[0129]
(2)进一步判断是否受到干扰：经过拐点判断处理后，5个值的相邻电压值是否有
超过35毫伏的，若有去掉该值，插值代替，即用该值前后值的均值代替该值，(比如v3超过，用代替v3)，若最后一个值v5超过35毫伏，则用2v
4-v3代替v5，去掉最大和最小两个值，求3个值的算术平均值，得到的均值即为5毫秒电压瞬时值v
ss5
。判断如0≤v
ss5
≤vs，否则存储α1；
[0130]
第2步：间隔时间5毫秒，重复第1步得到α2，存储α1、α2；
[0131]
第3步：计算角速度存储ω1；
[0132]
第4步：间隔5毫秒再重复第1步，第2步和第3步，循环后得到并存储角度α1、α2、α3，角速度ω1、ω2；
[0133]
第5步：计算角加速度，存储λ1；
[0134]
第6步：循环步骤1至步骤5，得到时长250毫秒，每个数据间隔5毫秒的50个连续的角度瞬时值α1、α2、λλ、α
50
，连续49个角速度瞬时值ω1、ω2、、λλ、ω
49
，连续48个角加速度瞬时值λ1、λ2、、λλ、λ
48
，全部存储，α9，ω8，每5毫秒循环更新一次，经过n 2次更新后得到的数据为：角度瞬时值为α
n 3
、α
n 4
、λλ、α
n 52
，角速度瞬时值为ω
n 2
、ω
n 3
、λλ、ω
n 51
，角加速度为λ
n 1
、λ
n 2
、λλ、λ
n 50
。
[0135]
第7步：构建插值表。
[0136]
风速与传感器风杯的角速度成线性单调增加，但不是线性增加。采用插值法处理：传感器设计加工好后，风速与风杯的角速度保持不变，通过风洞可以测试出风速与角速度的对应数据，构成一个风速角速度查值表，通过分析和实际测量，选取风速为0.1m/s，0.3m/s，0.5m/s，1m/s，1.5m/s，2，5m/s，10m/s，15m/s，20m/s，30m/s，40m/s，50m/s，60m/s，70m/s，80m/s，分别计为v
b1
、v
b2
、λλ、v
b14
、v
b15
，通过风洞提供上述风速分别测出对应风杯的角速度值，计为ω
b1
、ω
b2
、λλ、ω
b14
、ω
b15
，
[0137]
第8步：插值算法计算风速。在插值表给出的数据中，将相邻两点间的风速看成线性增加，将前面计算存储的50个角速度值分别按照下面计算处理：
[0138]
判断ω值：ω等于ω
b1
、ω
b2
、λλ、ω
b14
、ω
b15
中的某个值ω
bn
，则v＝v
bn
，否则，ω将落在ω
b1
、ω
b2
、λλ、ω
b14
、ω
b15
中的两个数之间，即ω
bn
＜ω＜ω
bn-1
，则按如下公式计算：v为瞬时风速(5毫秒平均风速)。
[0139]
将ω1、ω2、λλ、ω
14
、ω
50
依次全部代入处理计算得到相应的250毫秒时长的50个风速值，即为：v1、v2、λλ、v
14
、v
50
，平均即得到250毫秒平均风速v
ss250
：
[0140][0141]
第9步：信号输出，一种是spi数字输出，一种将电压转化为4-20毫秒信号输出，一种是0-3.3伏电压输出。
[0142]
电位器参数自动检测方法：已选定的全角度电位器和设计调试好的电路，i、γ1、γ2、r1、r2、θ均已确，但其参数还需测量，i、r1、r2很容易测量，γ1、γ2、θ均是角度不易测量，针对这3个参数编写一个检测程序，通过plc程控步进电机(步进角度0.01
°
)带动满角度电
位器向一个方向旋转(避免回差)，其旋转角度可控可测，每旋转一次停留3秒钟，待电路稳定，输出一组角度电压数据。
[0143]
第1步：初测。步进1
°
，测出v
s1
、v
l1
、v
m1
电位器对应的角度，α
s1
、α
l1
、α
m1
。
[0144]
具体方法：步进电机带动满角度电位器向一个方向旋转，每次旋转1
°
，每旋转一次停留3秒钟，待电路稳定再检测，电位器从0
°
开始旋转359
°
，记录一组360个电压角度数据，电压角度数据一一对应。寻找判断v
s1
、v
l1
、v
m1
对应的α
s1
、α
l1
、α
m1
，方法：步进电机旋转，检测到0电压时，开始计算步进电机的步进次数，步进359次，记录0
°
到359
°
对应的电压值，构成电压角度表，表中找到最大电压值即为v
l1
，该值对应的角度即为α
l1
，比α
l1
小1
°
的即为α
s1
，对应的电压值即为v
s1
，α
m1
＝359
°
，对应的电压即为v
m1
。则有：
[0145][0146]
第2步：细测。步进0.01
°
，测出v
s2
、v
l2
、v
m2
电位器对应的角度，α
s2
、α
l2
、α
m2
。
[0147]
具体方法：步进电机带动满角度电位器向一个方向快速旋转，每次旋转0.01
°
，检测到电压为0时，停留3秒钟，确保电路稳定，再次检测确认电压为0，启动步进电机开始计数，快速转动到α
s1-2，开始3秒钟步进一次(每次步进0.01
°
)，直到α
s1
1，记录该3
°
的电压角度数据300个，构成电压角度表，表中找到最大电压值即为v
l2
，该值对应的角度即为α
l2
，比α
l2
小0.01
°
的即为α
s2
，对应的电压值即为v
s2
，继续快速旋转到359.99
°
，停留检测输出电压值，该值即v
m2
，α
m2
＝359.99
°
。此步得到的电压和角度精度高。代入前面公式可计算出准确的(误差允许范围内)γ1、γ2、θ。
[0148][0149]
说明：文中角度单位为度(
°
)，角速度单位为弧度/秒(rand/s)，在进行角速度计算时自动换算为弧度，转换关系为1弧度＝180/π度，1π＝π/180弧度。具体应用中要注意各参数的量纲。