0 引 言

风能作为一种重要的自然资源被广泛应用于风能发电、 风力助航和气象测量等领域。因此风速和风向的准确测量一直是国内外学者研究的热点[1,2]。目前已经开发了多种类型的测风仪器, 其中常用的仪器是机械风速仪, 但由于其具有旋转部件, 存在机械磨损、 寿命短、 维修费用高、 精度低等缺点, 影响其应用[3]。超声波测风技术中使用最广泛的方法是时间差分理论[4],它利用顺风和逆风之间的传播时间差测量风速。这种方法结构简单, 但没有考虑噪声对风速传播的影响。

阵列信号处理方法可有效抑制复杂环境中存在的噪声[5,6]。目前很少有人将信号处理方法与风速测量相结合。Han等[7]提出了一种利用超声波传感器阵列的风速和风向测量方法, 采用传感器位于上层发射超声波信号, 其他4个传感器位于下层接收超声波信号的阵列结构, 并采用时差法作为测风算法。然而其测量精度完全取决于超声波传输时间差的测量, 特别是在强电磁干扰或低信噪比的情况下, 超声波传输的时间测量将变得困难甚至失效。Li等[8]提出了一种基于超声波传感器阵列的MUSIC(Multiple Signal Classification)风速风向的测量方法。笔者提出了一种基于弧形的超声波传感器阵列, 并且对弧形阵列结合MUSIC算法估计风速风向的方法展开了相关问题的分析和讨论。其采用的MUSIC算法是一种高分辨估计算法[9,10], 但需对阵列接收数据的协方差进行奇异值分解, 增加了算法的复杂性。

波束形成算法是阵列信号处理理论中的一种重要方法[11], 一方面具有较强的噪声抑制能力,另一方面可得到高精度的目标分辨力[12]。值得注意的是, MVDR(Minimum Variance Distortioniess Response)功率估计是作为一个凸优化问题的解决方案而获得的, 这使MVDR估计更加稳健[13]。

笔者将超声波传感器阵列结构与MVDR波束形成算法相结合, 应用于风速测量系统, 能在低信噪比条件下能准确测量风速和风向, 仿真结果验证了笔者方法的有效性。

1 信号模型

对于由M个相同的传感器组成的阵元任意分布在同一平面上的传感器阵列, 假设所有阵元都是各向同性的, 且无信道不一致或互耦现象, 采样后的传感器 M 维向量输出模型为

2)加强ADSS光缆设计、施工和验收等环节的工程管理。在光缆敷设时不仅要选择适宜悬挂点，还应对其进行力学校验，进行必要加固后方可加挂。ADSS光缆与输电线路同塔(杆)架设，周围存在高压电场，受到电腐蚀侵害在所难免，建议组织技术攻关，研究优化耐电腐外护套和防腐方案，在源头上预防断缆事故的发生。

X(t)=AS(t)+N(t)

(1)

其中t=1,…,T为采样时间, 阵列流型矩阵

A=[e-j2πfτ1,e-j2πfτ2,…,e-j2πfτM]T

(2)

待估计参数隐藏在时间延迟τ中, 只要确定阵元间的时间延迟表达式τ, 就可得到空间阵列的阵列流型。

  

图1 弧形超声波传感器阵列结构示意图Fig.1 Diagram of arc ultrasonic sensor array structure

笔者采用高伟[14]提出的一种弧形超声波传感器阵列, 具体结构如图1所示。传感器0为发送超声波传感器。传感器1～4均匀排列在以传感器0为圆心, R为半径的圆弧上, 组成接收超声波传感器阵列。传感器0～4共同构成超声波阵列测风系统。

图1中, V为待测风速, θ为待测风向角(风的来向为以竖直正方向顺时针偏θ方向)。V1～V4分别为风速V在传感器0与传感器1～4连线方向上的分量, α为两个相邻阵元与发射阵元连线夹角, 这里α取30°。

超声波传感器阵列流型推导过程如下。发射信号传播到传感器i的时间为

1/1-1d淹水处理的产量较对照显著增加，可能的原因一是分蘖期水稻耐淹能力最强，经过短历时(1 d)全淹处理后水稻植株体内抗氧化酶含量增加，使水稻表现出一定的抵抗淹水逆境的能力[14]，二是淹水胁迫去除后，稻株出现超补偿效应[15]。具体表现为水稻叶片的光合功能和活力提高，促进光合作用的进行及光合产物的生产和积累，1/1-1d处理的净光合速率(Pn值)较对照增加可以说明这一点。

ti=R/(c+Vi)

(3)

则图1所示的弧形阵列的阵列流型矢量为

式中：P1、P2为断面Ⅰ—Ⅰ、Ⅱ—Ⅱ处的静压，Pa；ρ1、ρ2为断面Ⅰ—Ⅰ、Ⅱ—Ⅱ处的空气密度，kg/m3；v1、v2为断面Ⅰ—Ⅰ、Ⅱ—Ⅱ处的风速，m/s；Z1、Z2为断面Ⅰ—Ⅰ、Ⅱ—Ⅱ处的标高，m；g为重力加速度，取9.8 m/s2；hr为风流从断面Ⅰ—Ⅰ到断面Ⅱ—Ⅱ过程中的阻力损失，Pa。

三千零一刀！短短一炷香的工夫，天葬师一共割出了三千零一刀，尸体的皮肉、内脏、骨骼，每一寸都是细碎而整齐的。这种强度的工作，让年迈的天葬师有些吃不消，他面色黑红，头顶热气氤氲，就连呼吸都打了颤。他知道自己老了，年轻的时候，他最多曾一连给四具尸身施行割礼，一口气割出了整整一万一千五百二十六刀！即使那样，他也没有觉得像如今这般劳累。

τ1=0

(4)

(5)

仪器设备在使用过程中应注意保养维护，做好仪器设备使用的档案记录，确保设备的完好安全和有效使用，避免损坏，造成浪费。对已达到或超过使用年限的仪器设备，按有关规定申请报废或降级使用[2]。

(6)

由Lagrange乘子法可以求得上述带约束条件的最小化问题的解为

(7)

其中i=1,2,3,4； c为超声波传播速度； R为弧形半径。

A(θ,V)=[e-j2πfτ1,e-j2πfτ2,e-j2πfτ3,e-j2πfτ4]

(8)

设超声波发送信号为s(t), 则第i个阵元接收到的信号为

由图1中的几何关系, 运用矢量分解的方法可得：以传感器1为参考阵元, 则各阵元与基准阵元的时延τi为

xi(t)=s(t)e-j2πfτi+ni(t)

(9)

其中f为发射的超声波频率； τi为发射信号到达第i个阵元相对于基准阵元的时延； ni(t)为第i个阵元上的噪声。

2 基于MVDR波束形成的风速风向测量方法

权力寻租理论，曾经是西方经济学家提出的一经济学上的概念。它是指掌握一定权力的人运用其手中的权力为资本，参与商品交换或市场竞争，来为自身谋取经济利益的现象。权力寻租是一种非生产性的活动，其自身并不会产生经济利益。在权力寻租的过程中，都是先有租金的诱惑，即丰厚的利润或者是较大的经济利益，然后再是掌权者将权力出租。根据权力寻租理论，掌权者为了获得经济利益往往会找机会将手中的权力出租出去，而那些拥有金钱而无权力的人则需要用金钱来租取一定的权力。权力寻租所带来的利益，是导致权力腐败的原动力。在市场经济在商业活动中，当出现法律不健全，权力制约机制薄弱的情况，就会不可避免的出现商业贿赂。

笔者采用经典的由Capon[15]提出的最小方差无畸变响应即MVDR波束形成算法(以下均称MVDR波束形成算法)估计风速V和风向θ。MVDR波束形成算法是一种主瓣约束的自适应方法, 其在所需方向有信号通过的情况下, 尽可能地抑制其他方向信号, 其实现手段就是调节阵列的权矢量。设MVDR波束形成的权为ω(θ,V) , 有

min E[]=min ωHRω, s.t. ωHa(θ,V)=1

(10)

 

(11)

其中R为接收信号的相关矩阵。则阵列的输出功率谱为

 

(12)

通过二维谱峰搜索即可实现对风速V和风向θ的估计。

3 仿 真

风向的均方根误差表示为 风速的均方根误差表示为

实验条件： 超声波传感器阵列的形式如图1所示, 超声波频率f=40 kHz； 声速c=340 m/s； 弧形半径R=0.1 m； 阵元噪声为加性高斯白噪声； 风速扫描范围为0～60 m/s, 步长为0.1 m/s(风速分辨率)； 风向角扫描范围为0°～359°, 步长为1°(风向分辨率)； 快拍数为512。信噪比范围为-10～10 dB, 间隔为 2 dB。对每个信噪比, 做100次Monte-Carlo实验。分别用于估计下列两组的参数： 1) V=2.7 m/s, θ=178°; 2) V=43.5 m/s, θ= 178°。

图2和图3中实线分别是基于MVDR算法的风速和风向均方根误差图, 图2和图3中虚线分别是基于MUSIC算法的风速和风向均方根误差图。可见两种算法对于风速和风向的均方根误差基本相同, 均是随信噪比的增加而逐渐减小, 当RSNR≥-6 dB(或-8 dB)时, 风速(或风向)的均方根误差较小, 估计精度很高, 可实现对风速(或风向)的估计。

有无活动性出血或明显渗血，填塞组48h后逐渐拔除高膨胀海绵，根据出血情况分为无明显出血、少量渗血、较多出血。不填塞组明胶海绵可吸收无需取出。

  

图2 MVDR与MUSIC风速均方根误差对比图 图3 MVDR与MUSIC风向均方根误差对比图 Fig.2 Comparison of the RMSE of wind velocity between MVDR and MUSIC Fig.3 Comparison of the RMSE of wind angle between MVDR and MUSIC

4 结 语

笔者提出了一种基于四元弧形超声波传感器及MVDR波束形成算法的风速风向测量方法。该方法利用超声波传感器阵列的冗余空间采样及阵列信号处理思想, 在较低信噪比的条件下达到了较高的测量精度, 并在相同信噪比环境下与基于MUSIC算法的风速风向测量方法进行统计性能的比较, 结果显示随着信噪比的增加, 两种算法可达到较相似的估计性能, 但计算量相对较少, 因此在实时性方面则更加优越。

参考文献：

[1]VAN B J S, RICHMOND A D. Radar Interferometry Technique: Three-Dimensional Wind Measurement Theory [J]. Radio Science, 2016, 26(5): 1209-1218.

[2]褚志刚, 段云炀, 沈林邦, 等. 奇异值分解波束形成声源识别方法 [J]. 电子测量与仪器学报, 2014(11): 1177-1184.

CHU Zhigang, DUAN Yunyang, SHEN Linbang, et al. Singular Value Decomposition Beamforming Method for Sound Source Identification [J]. Journal of Electronic Measurement and Instrumentation, 2014(11): 1177-1184.

[3]HYSON P. Cup Anemometer Response to Fluctuating Wind Speeds [J]. Journal of Applied Meteorology, 2010, 11(5): 843-848.

[4] 李晶, 莫德举. 基于改进时差法的超声流量计 [J]. 中国仪器仪表, 2004(12): 10-12.

LI Jing, MO Deju. Design of Ultrasonic Flow Meter Based on Improved Transit-time Method [J]. China Instrumentation, 2004(12): 10-12.

[5]王玉祥, 孙晓东, 单泽彪, 等. 海杂波背景下基于四元数电磁矢量阵列DOA估计 [J]. 吉林大学学报: 信息科学版, 2016, 34(1): 23-28.

WANG Yuxiang, SUN Xiaodong, SHAN Zebiao, et al. DOA Estimation Using Quaternion Based on Electromagnetic Vector Sensor Array on Real Clutter Data [J]. Journal of Jilin University: Information Science Edition, 2016, 34(1): 23-28.

[6]黄志强, 王树勋, 王波. 基于四阶累积量的近场源三维参数联合估计 [J]. 吉林大学学报: 信息科学版, 2006, 24(3): 225-229.

HUANG Zhiqiang, WANG Shuxun, WANG Bo. Jointly Estimating 3-D Parameters of Near-Field Sources Based on Fourth-Order Cumulant [J]. Journal of Jilin University: Information Science Edition, 2006, 24(3): 225-229.

[7]HAN D, KIM S, PARK S. Two-Dimensional Ultrasonic Anemometer Using the Directivity Angle of an Ultrasonic Sensor [J]. Microelectronics Journal, 2008, 39(10): 1195-1199.

[8]LI X, SUN H, GAO W, et al. Wind Speed and Direction Measurement Based on Arc Ultrasonic Sensor Array Signal Processing Algorithm [J]. ISA Transactions, 2016, 65: 437-444.

[9]石宇, 王波, 王树勋. 基于四阶累积量Root-Music的近场信源二维参数估计 [J]. 吉林大学学报: 信息科学版, 2007, 25(5): 465-470.

SHI Yu, WANG Bo, WANG Shuxun. Algorithm for 2-D Near-Field Sources Parameter Estimate Based on Fourth-Order Cumulant Root-Music [J]. Journal of Jilin University: Information Science Edition, 2007, 25(5): 465-470.

[10]姚林宏, 高鹰, 石宇, 等. 基于均匀线阵的改进MUSIC算法 [J]. 吉林大学学报: 信息科学版, 2010, 28(6): 558-562.

YAO Linhong, GAO Ying, SHI Yu, et al. Research on Improved Algorithm of MUSIC Based on ULA [J]. Journal of Jilin University: Information Science Edition, 2010, 28(6): 558-562.

[11]于涧, 胡亮, 王剑辉. 3G智能天线算法优化 [J]. 吉林大学学报: 信息科学版, 2008, 26(1): 16-20.

YU Jian, HU Liang, WANG Jianhui. Optimizing of Smart Antenna Algorithm in 3G [J]. Journal of Jilin University: Information Science Edition, 2008, 26(1): 16-20.

[12]何子述, 黄振兴, 向敬成. 一种有效的MVDR波束形成器 [J]. 信号处理, 2000, 16(4): 302-305.

HE Zishu, HUANG Zhenxing, XIANG Jingcheng. An Effective MVDR Beamformer [J]. Signal Processing, 2000, 16(4): 302-305.

[13]AUBRY A, CAROTENUTO V, MAIO A D. A New Optimality Property of the Capon Estimator [J]. IEEE Signal Processing Letters, 2017(99): 1.

[14]高伟. 基于超声波传感器阵列及多重信号分类算法的风速风向测量方法 [D]. 长春： 吉林大学通信工程学院, 2017.

GAO Wei. Wind Speed and Direction Measurement Method Based on Ultrasonic Sensor Array and Multiple Signal Classification Algorithm [D]. Changchun: College of Communication Engineering, Jilin University, 2017.

[15]CAPON J. High-Resolution Frequency-Wavenumber Spectrum Analysis [J]. Proceedings of the IEEE, 2005, 57(8): 1408-1418.