CN110990772A - 星载天线指向精度分析方法 - Google Patents
星载天线指向精度分析方法 Download PDFInfo
- Publication number
- CN110990772A CN110990772A CN201910986674.6A CN201910986674A CN110990772A CN 110990772 A CN110990772 A CN 110990772A CN 201910986674 A CN201910986674 A CN 201910986674A CN 110990772 A CN110990772 A CN 110990772A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- pointing
- antenna
- theoretical
- normal vector
- 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.)
- Granted
Links
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 (7)
1.一种星载天线指向精度分析方法,其特征在于,包括:
设定星体坐标系O0X0Y0Z0和Z0轴上两个坐标转换参考点a、b,其中,O0X0水平向右,O0Z0为星体安装面法向量,O0Y0垂直于XZ平面并符合右手定则;
根据所述星体坐标系O0X0Y0Z0和Z0轴上两个坐标转换参考点a、b,并通过齐次坐标变换求得天线顶端的理论坐标系O3X3Y3Z3以及含误差的实际坐标系O3’X3’Y3’Z3’,得到理论参考点c、d及实际参考点c’、d’的全局坐标;
利用matlab数学分析软件计算所述机构指向精度的表达式。
2.如权利要求1所述的星载天线指向精度分析方法,其特征在于,设定Z0轴上两个坐标转换参考点a、b,包括:
在全局坐标系Z0轴上设定两个坐标转换参考点a=[0 0 0.5]T和b=[0 0 1]T。
3.如权利要求1所述的星载天线指向精度分析方法,其特征在于,根据所述星体坐标系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轴的加工误差;
天线实际指向法向量为:
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 true CN110990772A (zh) | 2020-04-10 |
CN110990772B 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) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111611742A (zh) * | 2020-06-11 | 2020-09-01 | 上海卫星工程研究所 | 平面变形绝对-相对欧拉角计算方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107799875A (zh) * | 2017-09-18 | 2018-03-13 | 西安空间无线电技术研究所 | 一种星载天线三轴指向机构精度补偿方法 |
CN108225239A (zh) * | 2018-01-25 | 2018-06-29 | 哈尔滨工业大学 | 一种星载天线的伸展臂末端位姿静态误差域的计算方法 |
CN109781060A (zh) * | 2019-01-08 | 2019-05-21 | 上海卫星工程研究所 | 星载点波束天线对地指向精度评估方法 |
CN109781059A (zh) * | 2019-01-08 | 2019-05-21 | 上海卫星工程研究所 | 星载点波束天线对地指向精度评估系统 |
-
2019
- 2019-10-14 CN CN201910986674.6A patent/CN110990772B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107799875A (zh) * | 2017-09-18 | 2018-03-13 | 西安空间无线电技术研究所 | 一种星载天线三轴指向机构精度补偿方法 |
CN108225239A (zh) * | 2018-01-25 | 2018-06-29 | 哈尔滨工业大学 | 一种星载天线的伸展臂末端位姿静态误差域的计算方法 |
CN109781060A (zh) * | 2019-01-08 | 2019-05-21 | 上海卫星工程研究所 | 星载点波束天线对地指向精度评估方法 |
CN109781059A (zh) * | 2019-01-08 | 2019-05-21 | 上海卫星工程研究所 | 星载点波束天线对地指向精度评估系统 |
Non-Patent Citations (4)
Title |
---|
张锋等: "星载天线指向机构误差分析与建模", 《电子机械工程》 * |
张锋等: "星载天线指向机构误差分析与建模", 《电子机械工程》, no. 01, 15 February 2010 (2010-02-15) * |
潘博等: "星载天线指向精度建模与分析", 《航天器工程》 * |
潘博等: "星载天线指向精度建模与分析", 《航天器工程》, no. 05, 15 September 2011 (2011-09-15) * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111611742A (zh) * | 2020-06-11 | 2020-09-01 | 上海卫星工程研究所 | 平面变形绝对-相对欧拉角计算方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN110990772B (zh) | 2023-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107738254B (zh) | 一种机械臂坐标系的转换标定方法与系统 | |
US10241496B2 (en) | Multi-axis motor synchronization control system and method thereof | |
CN110802585B (zh) | 机械臂末端传感器补偿方法及接触力/力矩测量方法 | |
CN103792885B (zh) | 一种数控弯管加工仿真方法及装置 | |
CN109163675B (zh) | 一种基于激光跟踪仪检测角摆轴位置精度的方法 | |
US20190285396A1 (en) | Method for extracting gear tooth profile edge based on engagement-pixel image edge tracking method | |
CN113580148B (zh) | 基于等效运动链的并联机器人运动学标定方法 | |
CN113211445B (zh) | 一种机器人参数标定方法、装置、设备及存储介质 | |
CN108062071B (zh) | 参数曲线轨迹伺服轮廓误差的实时测定方法 | |
CN107992647A (zh) | 一种受几何结构影响的局部并联尺寸链误差获取方法 | |
CN111489399A (zh) | 一种视觉跟踪组件安装参数标定装置和方法 | |
CN110990772A (zh) | 星载天线指向精度分析方法 | |
CN108182330A (zh) | 一种基于b样条计算柔性矩形薄板刚柔耦合动力学响应的方法 | |
CN113009883B (zh) | 基于3-rpsr并联机构的自由弯曲成形设备运动位姿解析方法 | |
CN108108555B (zh) | 一种装配尺寸链计算方法、装置及设备 | |
CN110569582B (zh) | 一种用于天线的二维指向机构的指向精度计算方法及装置 | |
CN110888123A (zh) | 一种基于旋转矩阵的雷达坐标转换方法 | |
Srinivas et al. | A novel method to synthesize a single topology for dynamically loaded members | |
CN111044289B (zh) | 基于闭环动态测量的大型高速回转装备对准误差测量方法 | |
CN106934108B (zh) | 一种并联构型主轴头的动态误差传递特性评价方法 | |
Duong et al. | Robot Control Using Alternative Trajectories Based on Inverse Errors in the Workspace | |
Wang et al. | A generating method for digital gear tooth surfaces | |
CN111504276B (zh) | 一种基于视觉投影尺度因子集的联合目标函数多推进器姿态角获取方法 | |
Fan et al. | Improvement in hole-pose error for aerospace drilling applications based on Hermite surface reconstruction and manifold error similarity | |
CN112659121A (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 |