CN107830836B - 一种双冗余psd自准直仪吊索摆角的测量方法 - Google Patents

一种双冗余psd自准直仪吊索摆角的测量方法 Download PDF

Info

Publication number
CN107830836B
CN107830836B CN201711071977.2A CN201711071977A CN107830836B CN 107830836 B CN107830836 B CN 107830836B CN 201711071977 A CN201711071977 A CN 201711071977A CN 107830836 B CN107830836 B CN 107830836B
Authority
CN
China
Prior art keywords
psd
autocollimator
angle
sling
measuring
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201711071977.2A
Other languages
English (en)
Other versions
CN107830836A (zh
Inventor
谢文博
卜林卫
张健
许家忠
尤波
黄玲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin University of Science and Technology
Original Assignee
Harbin University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN201711071977.2A priority Critical patent/CN107830836B/zh
Publication of CN107830836A publication Critical patent/CN107830836A/zh
Application granted granted Critical
Publication of CN107830836B publication Critical patent/CN107830836B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C1/00Measuring angles

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

一种双冗余PSD自准直仪吊索摆角的测量方法,涉及PSD自准直仪的角度测量领域。本发明是为了解决现有的悬吊漂浮物随动控制系统中,使用单一PSD自准直仪测量吊索摆角时产生误差较的问题。本发明所述的一种双冗余PSD自准直仪吊索摆角的测量方法,首先建立了双冗余PSD自准值仪的测量结构,随后给出双冗余PSD自准直仪角度测量系统模型,在悬吊物铅锤并自旋状态时测试反射镜的安装情况,悬吊物铅锤并自旋运动时通过循环迭代算法求解数学模型中的各个参数,进而获得所需的测量角度。本发明能提高系统角度测量精度,较大程度上减小由于PSD镜面安装时水平度偏差导致的角度测量误差。

Description

一种双冗余PSD自准直仪吊索摆角的测量方法
技术领域
本发明属于PSD自准直仪的角度测量领域。
背景技术
近几十年来,航天卫星在人类的生产生活中占据了极重要的位置,并且对卫星的姿态和位置控制功能提出了更高的要求,但基于真实卫星的控制实验往往成本高昂,这就使得较低成本的地面低重力卫星模拟控制系统在航天工业中变的愈加重要。悬吊法作为地面低重力卫星模拟控制系统的一种主要方法,具有成本低、实验效果较贴近真实情况的优点,并且在悬吊漂浮物随动控制系统中,对于吊索的摆角测量精确程度成为影响系统控制精度的关键。
对于吊索摆角的测量方法,一般分为三种:机械式、电磁式和光学式。光学式测量方法的精度较高,所以近些年被较多的采用。PSD(位置检测)自准直仪是光学方式测量时非常重要的一种测量仪器,由于它具有较高的准确度和测量分辨力,因而得到了广泛的应用。现有的悬吊漂浮物随动控制系统吊索摆角的测量方法为使用单一PSD自准直仪进行测量,这种测量方法会因为悬吊索搭载的反射镜在真实环境下无法调整至完全水平和吊索自旋的问题,而导致测量结果产生较大的误差,进而会引发后续一系列的测量和控制误差。
当前对于PSD自准直仪测量角度的研究文献较少,文献《基于PSD的角度测量系统》后面的文献和专利也用相同方式处理阐述了一种利用PSD传感器和处理电路测量汇聚点位置进而得到入射光角度的方法。该方法不仅操作复杂,而且测量精度不高。文献《基于PSD的微小角度动态测量研究》提出了一种利用PSD传感器和高速旋转平面测量角度变化的方法,该方法适用于微小角度变化测试,不适用于长时间的动态角度测量,且装置安装复杂,无法保证运动过程中角度测量精度。发明专利《采用PSD的光电倾角测量装置》提出了一种采用PSD和半导体激光测量倾角的装置,该装置将分光折转镜作为固体摆,利用其中光的折射和反射最终达到测量的效果,该方法安装复杂,且测量精度不理想。
发明内容
本发明是为了解决现有的悬吊漂浮物随动控制系统中,使用单一PSD自准直仪测量吊索摆角时产生误差较大的问题,现提供一种双冗余PSD自准直仪吊索摆角的测量方法。
一种双冗余PSD自准直仪吊索摆角的测量方法,所述测量方法是基于PSD自准直仪角度测量系统实现的,所述系统包括:随动平台1、上位PSD自准直仪2、吊索3、下位反射镜4、悬吊物5、下位PSD自准直仪6和上位反射镜7;
上位反射镜7固定在随动平台1底面,上位反射镜7中心处固定有吊索3,下位反射镜4和悬吊物5由上至下依次悬吊在吊索3上,上位PSD自准直仪2和下位PSD自准直仪6均用于采集吊索摆角,上位PSD自准直仪2的出射光能够照射到下位反射镜4上,下位PSD自准直仪6的出射光能够照射到上位反射镜7上;
所述测量方法为:
使悬吊物5铅锤,分别记录上位PSD自准直仪2和下位PSD自准直仪6在各个时刻所采集的角度值,并判断角度值是否为零;
若角度值均为零,则下位反射镜4和上位反射镜7处于水平状态,当前上位PSD自准直仪2和下位自准直仪6所采集的角度值为吊索摆角;
若角度值不为零且呈现周期性变化,则建立系统模型构成的非线性方程,选取至少5个相邻时刻的角度值代入非线性方程中构成方程组,构造误差函数对方程组求解,将求解获得的参数值与tn时刻的上位PSD自准直仪2和下位PSD自准直仪6采集的角度值代入测量系统所测角度的数学模型,获得tn时刻的吊索摆角θx
所述求解获得的参数值为:下位反射镜4的倾斜角A1,上位反射镜7的倾斜角A2,下位PSD自准直仪2的初始相位
Figure GDA0002233073320000024
上位PSD自准直仪6的初始相位
Figure GDA0002233073320000025
吊索的转动角频率ω。
上位PSD自准直仪2在tn时刻采集的角度值y1n和下位PSD自准直仪6在tn时刻采集的角度值y2n分别为:
Figure GDA0002233073320000021
系统模型构成的非线性方程为:
Figure GDA0002233073320000022
其中,y1n为上位PSD自准直仪2在tn时刻采集的角度值,y2n为下位PSD自准直仪6在tn时刻采集的角度值。
测量系统所测角度的数学模型为:
Figure GDA0002233073320000023
其中,y1n为上位PSD自准直仪2在tn时刻采集的角度值,y2n为下位PSD自准直仪6在tn时刻采集的角度值。
构造误差函数对方程组求解的具体方法为:
选取5个相邻时刻的角度值构造误差函数:
Figure GDA0002233073320000031
令:
Figure GDA0002233073320000032
解出A1、A2
Figure GDA0002233073320000033
和ω。
本发明所述的一种双冗余PSD自准直仪吊索摆角的测量方法,使用双冗余PSD自准直仪对角度进行测量,利用双冗余PSD自准直仪和循环迭代算法减小反射镜倾斜和吊索自旋的误差。与现有技术相比,对于吊索摆角的测量准确度提高了70%以上,能够有效的减少由于人工安装反射镜不能完全水平问题导致的测量误差,将误差范围控制在±0.03°以内,有利于工程应用。
附图说明
图1为双冗余PSD自准直仪角度测量系统的结构示意图;
图2为两个反射物镜存在摆角时的示意图;
图3为误差参数曲线图,其中(a)表示上位反射镜的倾斜角,(b)表示下位反射镜的倾斜角,(c)表示下位PSD自准直仪的初始相位,(d)表示上位PSD自准直仪的初始相位,(e)表示吊索的转动角频率,其横坐标为时间;
图4为吊索摆角曲线图,其横坐标为时间;
图5为一种双冗余PSD自准直仪吊索摆角的测量方法的流程图。
具体实施方式
参照图1、2、3和5具体说明本实施方式,本实施方式所述的一种双冗余PSD自准直仪吊索摆角的测量方法。
所述测量方法是基于双冗余PSD自准直仪角度测量系统实现的,所述测量系统包括:随动平台1、上位PSD自准直仪2、吊索3、下位反射镜4、悬吊物5、下位PSD自准直仪6和上位反射镜7;
上位反射镜7固定在随动平台1底面,上位反射镜7中心处固定有吊索3,下位反射镜4和悬吊物5由上至下依次悬吊在吊索3上,上位PSD自准直仪2和下位PSD自准直仪6均用于采集吊索摆角,上位PSD自准直仪2的出射光能够照射到下位反射镜4上,下位PSD自准直仪6的出射光能够照射到上位反射镜7上。
所述方法为:
吊索3会因自身特性产生自旋运动,自旋角频率为ω。此时系统的测量角度会因自旋运动产生周期性的误差。根据双冗余PSD自准直仪角度测量系统的特性,建立测量系统所测角度的数学模型:
Figure GDA0002233073320000041
其中,y1为上位PSD自准直仪2的测量值,y2为下位PSD自准直仪6的测量值,A1为下位反射镜4的倾斜角,A2为上位反射镜7的倾斜角,
Figure GDA0002233073320000042
为上位PSD自准直仪2的初始相位,
Figure GDA0002233073320000043
为下位PSD自准直仪6的初始相位,θx为吊索摆角,ω为吊索的转动角频率。y1和y2为实验的真实测量值;A1、A2
Figure GDA0002233073320000044
和ω为待求解的未知误差量;θx为本申请要测量的吊索摆角。
对上式变形,获得系统模型构成的非线性方程:
Figure GDA0002233073320000045
当随动平台1处于静止状态时,由于地球引力的作用,悬吊物5铅锤并自旋,设定固定的采样周期T,分别记录上位PSD自准直仪2和下位PSD自准直仪6各个时刻所采集到的角度值。
上位PSD自准直仪2和下位PSD自准直仪6对应的角度值分别记为:
y10,y11,y12……y1n;y20,y21,y22……y2n;n为正整数;
记录角度值的时刻分别记为:t0,t1,t2……tn,其中,tn时刻与tn-1时刻之间的时间间隔为T。
若所采集到的角度值均为零,则下位反射镜4和上位反射镜7均处于水平状态,吊索摆角为零;
若所采集到的角度值不为零且呈现周期性变化,则下位反射镜4和上位反射镜7均处于非水平状态,即:两个反射镜安装有水平偏差。在确定了两个反射镜安装有水平偏差后,需要进行后续的误差消除处理,具体处理过程为:
通过计算机运用最小二乘法求解系统模型所构成的非线性方程组:
Figure GDA0002233073320000051
构造误差函数:
Figure GDA0002233073320000052
为使总误差Q取极小,则令:
Figure GDA0002233073320000053
Figure GDA0002233073320000054
Figure GDA0002233073320000055
Figure GDA0002233073320000056
Figure GDA0002233073320000057
式(1-2)-(1-6)为关于A1、A2
Figure GDA0002233073320000058
和ω的方程组,解出A1、A2
Figure GDA0002233073320000059
和ω,则有:
Figure GDA00022330733200000510
随后使用MATLAB软件基于最小二乘原理,利用fslove函数采用最小二乘法求解非线性方程组,循环迭代至第n次求解为:
Figure GDA0002233073320000061
此时可以求解出模型中的参数A1、A2
Figure GDA0002233073320000062
和ω的值,并使用上述参数值代入模型中得到更精确的摆角θx
Figure GDA0002233073320000063
具体实施例:
本实施例为一种双冗余PSD自准直仪吊索摆角的测量方法,能够针对悬吊漂浮物随动控制系统中的吊索摆角测量进行验证。
首先,令随动平台静止,吊索在悬吊物的作用下铅锤并自旋,此时校验上、下位反射镜是否完全水平,反射镜水平校验情况如图4所示。上、下位PSD自准直仪的测量结果都不为零,且呈现周期性变化。
随动平台静止,吊索处于铅锤并自旋的运动状态,建立测量系统模型如下:
Figure GDA0002233073320000064
利用MATLAB软件求解系统模型中的误差参数,误差参数示意图如图3所示,A1,A2
Figure GDA0002233073320000065
ω分别为:
A1=0.1134,A2=0.1834,
Figure GDA0002233073320000066
ω=0.0898。
根据求得的上述误差参数,求解x轴方向的吊索摆角θx,吊索摆角示意图如图4所示,此时吊索铅锤自旋,真实摆角为0度。明显可见,双冗余PSD自准值仪测量方法的精度较两个独立PSD有较大提高,测量误差明显减小,测量效果明显改善。

Claims (5)

1.一种双冗余PSD自准直仪吊索摆角的测量方法,所述PSD自准直仪为位置检测自准直仪,其特征在于,所述测量方法是基于PSD自准直仪角度测量系统实现的,所述系统包括:随动平台(1)、上位PSD自准直仪(2)、吊索(3)、下位反射镜(4)、悬吊物(5)、下位PSD自准直仪(6)和上位反射镜(7);
上位反射镜(7)固定在随动平台(1)底面,上位反射镜(7)中心处固定有吊索(3),下位反射镜(4)和悬吊物(5)由上至下依次悬吊在吊索(3)上,上位PSD自准直仪(2)和下位PSD自准直仪(6)均用于采集吊索摆角,上位PSD自准直仪(2)的出射光能够照射到下位反射镜(4)上,下位PSD自准直仪(6)的出射光能够照射到上位反射镜(7)上;
所述测量方法为:
使悬吊物(5)铅锤,分别记录上位PSD自准直仪(2)和下位PSD自准直仪(6)在各个时刻所采集的角度值,并判断角度值是否为零;
若角度值均为零,则下位反射镜(4)和上位反射镜(7)处于水平状态,当前上位PSD自准直仪(2)和下位自准直仪(6)所采集的角度值为吊索摆角;
若角度值不为零且呈现周期性变化,则建立系统模型构成的非线性方程,选取至少5个相邻时刻的角度值代入非线性方程中构成方程组,构造误差函数对方程组求解,将求解获得的参数值与tn时刻的上位PSD自准直仪(2)和下位PSD自准直仪(6)采集的角度值代入测量系统所测角度的数学模型,获得tn时刻的吊索摆角θx
所述求解获得的参数值为:下位反射镜(4)的倾斜角A1,上位反射镜(7)的倾斜角A2,上位PSD自准直仪(2)的初始相位
Figure FDA0002233073310000011
下位PSD自准直仪(6)的初始相位
Figure FDA0002233073310000012
吊索的转动角频率ω。
2.根据权利要求1所述的一种双冗余PSD自准直仪吊索摆角的测量方法,其特征在于,系统模型构成的非线性方程为:
Figure FDA0002233073310000013
其中,y1n为上位PSD自准直仪(2)在tn时刻采集的角度值,y2n为下位PSD自准直仪(6)在tn时刻采集的角度值。
3.根据权利要求1所述的一种双冗余PSD自准直仪吊索摆角的测量方法,其特征在于,测量系统所测角度的数学模型为:
Figure FDA0002233073310000021
其中,y1n为上位PSD自准直仪(2)在tn时刻采集的角度值,y2n为下位PSD自准直仪(6)在tn时刻采集的角度值。
4.根据权利要求1所述的一种双冗余PSD自准直仪吊索摆角的测量方法,其特征在于,上位PSD自准直仪(2)在tn时刻采集的角度值y1n和下位PSD自准直仪(6)在tn时刻采集的角度值y2n分别为:
Figure FDA0002233073310000022
5.根据权利要求1所述的一种双冗余PSD自准直仪吊索摆角的测量方法,其特征在于,构造误差函数对方程组求解的具体方法为:
选取5个相邻时刻的角度值构造误差函数:
Figure FDA0002233073310000023
令:
Figure FDA0002233073310000024
Figure FDA0002233073310000025
Figure FDA0002233073310000026
Figure FDA0002233073310000027
Figure FDA0002233073310000028
解出A1、A2
Figure FDA0002233073310000029
和ω。
CN201711071977.2A 2017-11-03 2017-11-03 一种双冗余psd自准直仪吊索摆角的测量方法 Expired - Fee Related CN107830836B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711071977.2A CN107830836B (zh) 2017-11-03 2017-11-03 一种双冗余psd自准直仪吊索摆角的测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711071977.2A CN107830836B (zh) 2017-11-03 2017-11-03 一种双冗余psd自准直仪吊索摆角的测量方法

Publications (2)

Publication Number Publication Date
CN107830836A CN107830836A (zh) 2018-03-23
CN107830836B true CN107830836B (zh) 2020-03-27

Family

ID=61654520

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711071977.2A Expired - Fee Related CN107830836B (zh) 2017-11-03 2017-11-03 一种双冗余psd自准直仪吊索摆角的测量方法

Country Status (1)

Country Link
CN (1) CN107830836B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2013130032A (ru) * 2013-07-02 2015-01-10 Эдуард Николаевич Кузин Горизонтальный маятниковый измеритель угла
CN105526906A (zh) * 2015-11-30 2016-04-27 上海卫星工程研究所 大角度动态高精度激光测角方法
CN106081173A (zh) * 2016-07-19 2016-11-09 哈尔滨工业大学 三维主动悬吊式空间飞行器微重力模拟装置
CN107179059A (zh) * 2016-03-11 2017-09-19 罗伯特·博世有限公司 用于确定角误差的方法和光发射设备

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2013130032A (ru) * 2013-07-02 2015-01-10 Эдуард Николаевич Кузин Горизонтальный маятниковый измеритель угла
CN105526906A (zh) * 2015-11-30 2016-04-27 上海卫星工程研究所 大角度动态高精度激光测角方法
CN107179059A (zh) * 2016-03-11 2017-09-19 罗伯特·博世有限公司 用于确定角误差的方法和光发射设备
CN106081173A (zh) * 2016-07-19 2016-11-09 哈尔滨工业大学 三维主动悬吊式空间飞行器微重力模拟装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A measuring instrument for angle based on invariant optical path and autocollimation principle;Yang Liangen;《International Conference on Electric Information》;20110527;全文 *
冗余惯组初始基准现场校准装置研究;崔桂利;《宇航计测技术》;20130430;第33卷(第2期);全文 *
悬吊式零重力模拟系统的水平位置随动技术研究;邹胜宇;《中国优秀硕士学位论文全文数据库》;20150215;全文 *

Also Published As

Publication number Publication date
CN107830836A (zh) 2018-03-23

Similar Documents

Publication Publication Date Title
CN104697747B (zh) 一种平台系统的光学瞄准棱镜安装精度偏差标定检测方法
WO2022174831A1 (zh) 一种连续寻北装置、方法、电子设备及存储介质
CN109459585A (zh) 一种加速度计零位偏置修正方法
RU2378616C1 (ru) Астронавигационная система
CN109974750B (zh) 一种基于模糊逻辑系统的环形激光器温度建模及补偿方法
CN111044077A (zh) 星敏感器测量坐标系与基准立方镜坐标系间的标定方法
CN109141868B (zh) 精密轴系误差运动的测量装置及测量方法
CN107830836B (zh) 一种双冗余psd自准直仪吊索摆角的测量方法
CN101082491A (zh) 一种摆式陀螺寻北仪测量方法
CN108318051B (zh) 一种闭环光纤陀螺阈值测试方法
Luo et al. Rotating shaft tilt angle measurement using an inclinometer
CN109443333B (zh) 一种陀螺阵列反馈加权融合方法
CN105444777B (zh) 一种摇摆条件下光纤陀螺仪误差测试方法
CN111102918B (zh) 一种立方镜坐标系的自动化测量系统
CN110954131B (zh) 一种光纤陀螺输入轴失准角标定工装
CN108489396B (zh) 一种二维转顶精度检测方法
CN108036756B (zh) 一种利用加速度计进行双轴旋转惯性测量装置相邻轴线垂直度检查的方法
CN105758340B (zh) 新型火炮身管内膛直线度检测设备
CN103293959A (zh) 空间激光干涉系统激光指向控制技术的模拟方法和装置
RU2534815C1 (ru) Способ измерения децентрировки оптической оси асферической поверхности и сферометр для осуществления способа
CN109188681B (zh) 经F-theta聚焦透镜后引起的畸变误差补偿方法
KR101742431B1 (ko) 개회로 및 폐회로 방식의 고 내구성 자체 검정 경위 태양광 추적방법 및 그를 위한 시스템
JP2017032453A (ja) 回折格子を用いた姿勢角計測方法及び姿勢角計測装置
CN105758341B (zh) 一种火炮身管内膛直线度检测设备
Toro Technology review of wind-tunnel angle measurement

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200327

Termination date: 20201103