CN110990772B - 星载天线指向精度分析方法 - Google Patents
星载天线指向精度分析方法 Download PDFInfo
- Publication number
- CN110990772B CN110990772B CN201910986674.6A CN201910986674A CN110990772B CN 110990772 B CN110990772 B CN 110990772B CN 201910986674 A CN201910986674 A CN 201910986674A CN 110990772 B CN110990772 B CN 110990772B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- pointing
- axis
- antenna
- theoretical
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000013598 vector Substances 0.000 claims abstract description 49
- 239000011159 matrix material Substances 0.000 claims abstract description 36
- 230000009466 transformation Effects 0.000 claims abstract description 35
- 230000014509 gene expression Effects 0.000 claims abstract description 18
- 230000007246 mechanism Effects 0.000 claims description 23
- 238000006243 chemical reaction Methods 0.000 claims description 22
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 238000003754 machining Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 7
- 238000013461 design Methods 0.000 abstract description 7
- 238000004422 calculation algorithm Methods 0.000 abstract description 5
- 238000004458 analytical method Methods 0.000 abstract description 3
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明提供了一种星载天线指向精度分析方法,本发明基于齐次坐标变换矩阵对天线指向法向量进行表示,增加了算法的准确性;通过matlab数值分析软件对指向精度代数表达式进行化简,保证了表达式的简洁性,提高了指向精度分析计算精度和效率,为天线指向精度设计提供了可靠的理论依据。
Description
技术领域
本发明涉及一种星载天线指向精度分析方法。
背景技术
指向机构是实现天线实时跟踪、定位等功能的重要空间机构,其指向精度是评判星载天线指向机构能力的重要指标。由于星载天线运动日益复杂,对指向机构的设计要求不断提高,因此需要综合考虑安装面加工误差、驱动轴回转误差、热变形误差等因素,同时需要考虑到计算效率的问题。
发明内容
本发明的目的在于提供一种星载天线指向精度分析方法。
为解决上述问题,本发明提供一种星载天线指向精度分析方法,包括:
设定星体坐标系O0X0Y0Z0和Z0轴上两个坐标转换参考点a、b,其中,O0X0水平向右,O0Z0为星体安装面法向量,O0Y0垂直于XZ平面并符合右手定则;
根据所述星体坐标系O0X0Y0Z0和Z0轴上两个坐标转换参考点a、b,并通过齐次坐标变换求得天线顶端的理论坐标系O3X3Y3Z3以及含误差的实际坐标系O3’X3’Y3’Z3’,得到理论参考点c、d及实际参考点c’、d’的全局坐标;
利用matlab数学分析软件计算所述机构指向精度的表达式。
进一步的,在上述方法中,设定Z0轴上两个坐标转换参考点a、b,包括:
在全局坐标系Z0轴上设定两个坐标转换参考点a=[0 0 0.5]T和b=[0 0 1]T。
进一步的,在上述方法中,根据所述星体坐标系O0X0Y0Z0和Z0轴上两个坐标转换参考点a、b,并通过齐次坐标变换求得天线顶端的理论坐标系O3X3Y3Z3,包括:
确定参考坐标系O1X1Y1Z1相对于星体坐标系O0X0Y0Z0的齐次坐标变换矩阵为:
式中,X1和Z1分别为参考坐标系O1X1Y1Z1在星体坐标系X0和Z0方向的位移;
确定参考坐标系O2X2Y2Z2相对于参考坐标系O1X1Y1Z1的齐次坐标变换矩阵为:
式中,θ为参考坐标系O2X2Y2Z2相对于参考坐标系O1X1Y1Z1,以Y1为转轴的偏转角度;
确定理论坐标系O3X3Y3Z3相对于参考坐标系O2X2Y2Z2的齐次坐标变换矩阵为:
式中,Z3为理论坐标系(参考坐标系)O3X3Y3Z3在参考坐标系O2X2Y2Z2的Z2方向的位移。
天线理论指向法向量为:
用绕X、Y、Z三个坐标轴旋转δx、δy和δz来等价绕任意轴k转动微量角dθ,绕X轴旋转的微分变换矩阵为:
绕Y轴旋转的微分变换矩阵为:
绕Z轴旋转的微分变换矩阵为:
微分平移矩阵为:
驱动组件的回转误差,输出轴摆动误差,其误差微分变换矩阵为:
式中,δx、δy和δz为驱动组件相对于参考坐标系O1X1Y1Z1的X轴、Y轴和Z轴的偏转误差;
星体安装面加工误差,其误差微分变换矩阵为:
式中,dx、dy和dz为星体安装面相对于全局坐标系的X轴、Y轴和Z轴的加工误差;
天线实际指向法向量为:
进一步的,在上述方法中,所述机构指向精度为:
进一步的,在上述方法中,利用matlab数学分析软件计算所述机构指向精度的表达式,包括:
与现有技术相比,本发明基于齐次坐标变换矩阵对天线指向法向量进行表示,增加了算法的准确性;通过matlab数值分析软件对指向精度代数表达式进行化简,保证了表达式的简洁性,提高了指向精度分析计算精度和效率,为天线指向精度设计提供了可靠的理论依据。
附图说明
图1是本发明一实施例的星载指向机构示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
本发明是通过利用齐次坐标变换所得的实际指向法向量和理论指向法向量点积求得指向精度的方法,如图1所示,首先设定全局坐标系上的坐标转换参考点,以及各个部件的局部坐标系,再而通过齐次坐标变换矩阵将参考点坐标逐步转换至天线坐标系,通过转换后的参考点坐标作差得到实际和理论指向法向量,最后利用两者的点积求得指向精度;为了使结论简洁、清晰,利用matlab对指向精度表达式进行化简,提取主要的误差产生因素,使该方法的计算效率大幅提升。
本发明提供一种星载天线指向精度分析方法,包括:
步骤S1,设定星体坐标系O0X0Y0Z0和Z0轴上两个坐标转换参考点a、b,其中,O0X0水平向右,O0Z0为星体安装面法向量,O0Y0垂直于XZ平面并符合右手定则;
步骤S2,根据所述星体坐标系O0X0Y0Z0和Z0轴上两个坐标转换参考点a、b,并通过齐次坐标变换求得天线顶端的理论坐标系O3X3Y3Z3以及含误差的实际坐标系O3’X3’Y3’Z3’,得到理论参考点c、d及实际参考点c’、d’的全局坐标;
步骤S5,利用matlab数学分析软件计算所述机构指向精度的表达式。
在此,由于星载天线运动日益复杂,对指向机构的设计要求不断提高,因此需要综合考虑安装面加工误差、驱动轴回转误差、热变形误差等因素,得出指向精度精确分析的一种方法,为指向机构的指向精度设计提供理论依据。同时考虑到计算效率,需对算法进行简化,本发明提供一种准确、高效、灵活的指向精度分析方法。
本发明采用了齐次坐标变换矩阵进行坐标变换,简化了繁琐的逐步推导过程,便于利用matlab进行代数表达式的计算。本发明在理论指向法向量推导过程的基础上,考虑实际制造工艺上产生的误差、工作时的热变形等因素,利用齐次坐标变换矩阵的运算法则,增加误差矩阵,得到实际指向法向量。
本发明利用matlab对代数式进行化简,可得到天线指向精度与各轴回转误差、偏摆误差的代数关系;该方法将天线理论法向量与实际法向量比较,有效地提高了天线指向精度计算精度和效率,为天线指向精度设计提供了可靠的理论依据。
本发明的星载天线指向精度分析方法一实施例中,设定Z0轴上两个坐标转换参考点a、b,包括:
在全局坐标系Z0轴上设定两个坐标转换参考点a=[0 0 0.5]T和b=[0 0 1]T,以便于坐标系转换之后的指向法向量的计算和表示。
本发明的星载天线指向精度分析方法一实施例中,根据所述星体坐标系O0X0Y0Z0和Z0轴上两个坐标转换参考点a、b,并通过齐次坐标变换求得天线顶端的理论坐标系O3X3Y3Z3,包括:
确定参考坐标系O1X1Y1Z1相对于全局坐标系O0X0Y0Z0的齐次坐标变换矩阵为:
式中,X1和Z1分别为参考坐标系O1X1Y1Z1在星体坐标系(全局坐标系)X0和Z0方向的位移;
确定参考坐标系O2X2Y2Z2相对于参考坐标系O1X1Y1Z1的齐次坐标变换矩阵为:
式中,θ为参考坐标系O2X2Y2Z2相对于参考坐标系O1X1Y1Z1,以Y1为转轴的偏转角度;
确定理论坐标系(参考坐标系)O3X3Y3Z3相对于参考坐标系O2X2Y2Z2的齐次坐标变换矩阵为:
式中,Z3为理论坐标系(参考坐标系)O3X3Y3Z3在参考坐标系O2X2Y2Z2的Z2方向的位移。
天线理论指向法向量为:
用绕X、Y、Z三个坐标轴旋转δx、δy和δz来等价绕任意轴k转动微量角dθ,绕X轴旋转的微分变换矩阵为:
绕Y轴旋转的微分变换矩阵为:
绕Z轴旋转的微分变换矩阵为:
微分平移矩阵为:
驱动组件的回转误差,输出轴摆动误差,其误差微分变换矩阵为:
式中,δx、δy和δz为驱动组件相对于参考坐标系O1X1Y1Z1的X轴、Y轴和Z轴的偏转误差;
星体安装面加工误差,其误差微分变换矩阵为:
式中,dx、dy和dz为星体安装面相对于全局坐标系的X轴、Y轴和Z轴的加工误差;
天线实际指向法向量为:
本发明的星载天线指向精度分析方法一实施例中,所述机构指向精度为:
本发明的星载天线指向精度分析方法一实施例中,利用matlab数学分析软件计算所述机构指向精度的表达式,包括:
在此,应用时,只需了解指向机构零部件各向的误差即可求得指向精度。
综上所述,本发明可以在了解各轴偏转误差的前提下,快速计算出天线的指向精度,具有准确、高效、灵活的特点。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
专业人员还可以进一步意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,为了清楚地说明硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
显然,本领域的技术人员可以对发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包括这些改动和变型在内。
Claims (1)
1.一种星载天线指向精度分析方法,其特征在于,包括:
设定星体坐标系O0X0Y0Z0和Z0轴上两个坐标转换参考点a、b,其中,O0X0水平向右,O0Z0为星体安装面法向量,O0Y0垂直于XZ平面并符合右手定则;
根据所述星体坐标系O0X0Y0Z0和Z0轴上两个坐标转换参考点a、b,并通过齐次坐标变换求得天线顶端的理论坐标系O3X3Y3Z3以及含误差的实际坐标系O3’X3’Y3’Z3’,得到理论参考点c、d及实际参考点c’、d’的全局坐标;
利用matlab数学分析软件计算所述机构指向精度的表达式;
设定Z0轴上两个坐标转换参考点a、b,包括:
在全局坐标系Z0轴上设定两个坐标转换参考点a=[0 0 0.5]T和b=[0 0 1]T;
根据所述星体坐标系O0X0Y0Z0和Z0轴上两个坐标转换参考点a、b,并通过齐次坐标变换求得天线顶端的理论坐标系O3X3Y3Z3,包括:
确定参考坐标系O1X1Y1Z1相对于星体坐标系O0X0Y0Z0的齐次坐标变换矩阵为:
式中,X1和Z1分别为参考坐标系O1X1Y1Z1在星体坐标系X0和Z0方向的位移;
确定参考坐标系O2X2Y2Z2相对于参考坐标系O1X1Y1Z1的齐次坐标变换矩阵为:
式中,θ为参考坐标系O2X2Y2Z2相对于参考坐标系O1X1Y1Z1,以Y1为转轴的偏转角度;
确定理论坐标系O3X3Y3Z3相对于参考坐标系O2X2Y2Z2的齐次坐标变换矩阵为:
式中,Z3为理论坐标系(参考坐标系)O3X3Y3Z3在参考坐标系O2X2Y2Z2的Z2方向的位移;
天线理论指向法向量为:
用绕X、Y、Z三个坐标轴旋转δx、δy和δz来等价绕任意轴k转动微量角dθ,绕X轴旋转的微分变换矩阵为:
绕Y轴旋转的微分变换矩阵为:
绕Z轴旋转的微分变换矩阵为:
微分平移矩阵为:
驱动组件的回转误差,输出轴摆动误差,其误差微分变换矩阵为:
式中,δx、δy和δz为驱动组件相对于参考坐标系O1X1Y1Z1的X轴、Y轴和Z轴的偏转误差;
星体安装面加工误差,其误差微分变换矩阵为:
式中,dx、dy和dz为星体安装面相对于全局坐标系的X轴、Y轴和Z轴的加工误差;
天线实际指向法向量为:
所述机构指向精度为:
利用matlab数学分析软件计算所述机构指向精度的表达式,包括:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910986674.6A CN110990772B (zh) | 2019-10-14 | 2019-10-14 | 星载天线指向精度分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910986674.6A CN110990772B (zh) | 2019-10-14 | 2019-10-14 | 星载天线指向精度分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110990772A CN110990772A (zh) | 2020-04-10 |
CN110990772B true CN110990772B (zh) | 2023-07-07 |
Family
ID=70082093
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910986674.6A Active CN110990772B (zh) | 2019-10-14 | 2019-10-14 | 星载天线指向精度分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110990772B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111611742B (zh) * | 2020-06-11 | 2022-03-18 | 上海卫星工程研究所 | 平面变形绝对-相对欧拉角计算方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109781060A (zh) * | 2019-01-08 | 2019-05-21 | 上海卫星工程研究所 | 星载点波束天线对地指向精度评估方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107799875B (zh) * | 2017-09-18 | 2020-07-14 | 西安空间无线电技术研究所 | 一种星载天线三轴指向机构精度补偿方法 |
CN108225239B (zh) * | 2018-01-25 | 2019-10-22 | 哈尔滨工业大学 | 一种星载天线的伸展臂末端位姿静态误差域的计算方法 |
CN109781059B (zh) * | 2019-01-08 | 2021-05-11 | 上海卫星工程研究所 | 星载点波束天线对地指向精度评估系统 |
-
2019
- 2019-10-14 CN CN201910986674.6A patent/CN110990772B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109781060A (zh) * | 2019-01-08 | 2019-05-21 | 上海卫星工程研究所 | 星载点波束天线对地指向精度评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110990772A (zh) | 2020-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11371821B2 (en) | Method for extracting gear tooth profile edge based on engagement-pixel image edge tracking method | |
CN108692688B (zh) | 一种机器人测量-加工系统扫描仪坐标系自动标定方法 | |
US10241496B2 (en) | Multi-axis motor synchronization control system and method thereof | |
US20110218780A1 (en) | Method and Apparatus for Estimating Error in Multi-Axis Controlled Machine | |
CN106289641B (zh) | 挠性航天器质心位置和转动惯量参数联合辨识方法 | |
CN109163675B (zh) | 一种基于激光跟踪仪检测角摆轴位置精度的方法 | |
CN105015800A (zh) | 航天器舱段在地面上的自动化装配系统 | |
CN109781059A (zh) | 星载点波束天线对地指向精度评估系统 | |
CN110990772B (zh) | 星载天线指向精度分析方法 | |
CN113211445A (zh) | 一种机器人参数标定方法、装置、设备及存储介质 | |
CN108108555B (zh) | 一种装配尺寸链计算方法、装置及设备 | |
Avizzano et al. | An optimal geometric model for clavels delta robot | |
Deng et al. | An improved posture evaluation method for cylindrical intersecting holes on large aerospace components based on monocular vision | |
CN113009883B (zh) | 基于3-rpsr并联机构的自由弯曲成形设备运动位姿解析方法 | |
CN110569582B (zh) | 一种用于天线的二维指向机构的指向精度计算方法及装置 | |
Srinivas et al. | A novel method to synthesize a single topology for dynamically loaded members | |
Chu et al. | A novel calibration method for ball joint position in a posture adjustment system | |
CN110702145A (zh) | 一种二维导航系统的地图误差补偿方法及系统 | |
CN111504276B (zh) | 一种基于视觉投影尺度因子集的联合目标函数多推进器姿态角获取方法 | |
CN112066943B (zh) | 一种矢量调节转动机构转动角度的计算方法及装置 | |
CN110427686B (zh) | 一种考虑约束条件的大刚体位移参量计算方法 | |
CN106707962B (zh) | 时变温度场中快插机构热变形误差补偿法 | |
CN111612847A (zh) | 机器人抓取操作的点云数据匹配方法及系统 | |
CN114647893B (zh) | 适应蒙皮装夹位姿及形变偏差的自动钻铆在线编程方法 | |
CN119845306B (zh) | 一种多源误差融合的二维转台在轨标定方法 |
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 |