单经纬仪测风.docVIP

单经纬仪测风实验报告

姓名：刘浩

学号：111170052

2013/5/8

一、实验目的

掌握GYR1电子式光学经纬仪的操作使用方法，学会冲灌、捕捉气球方法和准确读数的技术，编程处理观测数据。

二、实验原理

按固定升速数值求取氢气的净举力值，冲罐并释放固定升速气球，用经纬仪观测气球，保证气球一直在观测中心，导出数据，用单经纬仪测风计算程序计算出高空中不同高度的风向风速。

实验步骤

测风经纬仪安装、调试和使用

经纬仪应设置在开阔区域。

拉开并牢固的架设三脚架，使其顶面尽量调整水平，使顶面与观测员的第二颗衣扣同高。

从箱内取出经纬仪，用三脚架顶面中心螺钉把经纬仪固定。

转动望远镜，使水准仪与三个水平调节螺钉中的任意两个平行。调节这两个螺钉，使气泡居中央，然后将水准仪转过90度再调节另一个螺钉使气泡居中，这时望远镜转动到任意位置上，水准仪气泡偏离中心不得大于1格，否则重复上述步骤以至符合要求为止。

指北针定方位：打开磁针制动螺钉，松开磁针；转动镜筒，使指北针指向正北，并转动方位调整螺旋，使指北针的尖端对转磁针盒两条刻线，拨动刻度盘，使游尺“0”与方位盘上“0”对齐。使用指南针定方位，步骤同上，只是磁针的指向为正南。定好方位后，严禁转动调零手轮。

瞄准远方的任意目标物，调焦距使目标物最清晰。

操作变倍手轮，熟练掌握主、辅望远镜的转换方法。

灌球和观测

按天气条件确定气球型号、颜色。

10号球采用100米/分升速。

用当时的气温、气压查取标准密度升速值，用气球和附加物重量查取净举力。冲灌气球。

施放气球，跟踪观测。每隔30秒读取仰角和方位角，在读数时应使气球正好聚中在十字线中间，并停止转动仰角和方位角手轮。读取10组。

四、观测资料处理

1、列出两次观测的仰角和方位角，并用观测数据计算出不同高度的风向、风速。

编号

第一次

第二次

仰角

方位角

仰角

方位角

1

11.87

240.09

23.42

230.20

2

11.15

243.60

25.71

235.08

3

10.73

245.16

26.33

238.88

4

12.56

245.47

26.03

242.29

5

13.78

244.24

26.43

246.72

6

14.69

242.75

27.78

252.53

7

16.98

243.10

29.33

259.32

8

18.50

245.22

30.97

265.95

9

20.59

249.03

32.56

272.93

10

22.66

254.57

33.54

279.91

鉴于第一次观测为尝试，现使用第二次观测数据分析处理。

取南北方向为x轴（北为正），东西方向为y轴（东为正），气球水平面投影点的坐标分量为：

xn=Lncosαn.

yn=Lnsinαn.

其中

Ln=Hnctgβn。

αn为方位角，βn为仰角。

则相应层的风速:

Vn=[(xn-xn-1)2+(yn-yn-1)2]1/2/t.

其中t为相邻两次观测的时间间隔。

高度（m）

风速(m/s)

方位角

25

3.85

129.8

75

3.11

118.9

125

3.23

113.0

175

3.62

108.1

225

3.33

94.8

275

2.86

71.1

325

2.95

51.1

375

2.88

37.6

425

3.06

25.1

475

3.39

22.7

2、利用实测风速采用最小二乘法拟合出风速幂指数。

在近地层(250米以下)，风速随高度的变化可用幂指数表示，其表达式为：

V2=V1(Z2/Z1)p……①

其中，V1、V2分别为Z1、Z2高度处的风速（m/s），p为风廓线指数。

最小二乘法公式：

对于平面上的直线：y=a+bx如有一组数据（xi,yi）

利用最小二乘法计算a,b时公式如下：

a=y0-bx0

b=∑(xi-x0)(yi-y0)/∑(xi-x0)2

将①式取对数并编程（见附录）计算可得到：

风廓线指数P=0.16878

近地层风速表达式V=1.76272*Z^0.16878

画出拟合曲线与观测数据差异对比图：

3、误差分析。

实验误差来源可能有以下几个：

（1）实验观测存在误差，观测时气球不断运动，很难保证每次观测都在视野正中央。

（2）假设条件的理论误差，在计算气球升速时所用气压等条件可能与实际条件有出入，导致气球的升速并不是100m/s。

实验总结

通过此次单经纬仪测风实验，了解了经纬仪测风的原理，初步掌握了电子式光学经纬仪的操作使用方法，学会了编程处理数据的方法。对近地层风速分布有了更深的认识。

附Fortran程序：

programmain

implicitnone

integer,parameter::n=10,f=4055,m=