CN113218379B - 一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法 - Google Patents

一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法 Download PDF

Info

Publication number
CN113218379B
CN113218379B CN202110501762.XA CN202110501762A CN113218379B CN 113218379 B CN113218379 B CN 113218379B CN 202110501762 A CN202110501762 A CN 202110501762A CN 113218379 B CN113218379 B CN 113218379B
Authority
CN
China
Prior art keywords
lithium battery
electromagnetic interference
angle
coordinate system
magnetic compass
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
Application number
CN202110501762.XA
Other languages
English (en)
Other versions
CN113218379A (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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN202110501762.XA priority Critical patent/CN113218379B/zh
Publication of CN113218379A publication Critical patent/CN113218379A/zh
Application granted granted Critical
Publication of CN113218379B publication Critical patent/CN113218379B/zh
Active 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
    • G01C17/00Compasses; Devices for ascertaining true or magnetic north for navigation or surveying purposes
    • G01C17/02Magnetic compasses
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/04Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means
    • G01C21/06Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means involving measuring of drift angle; involving correction for drift
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/04Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means
    • G01C21/08Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by terrestrial means involving use of the magnetic field of the earth
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B11/00Automatic controllers
    • G05B11/01Automatic controllers electric
    • G05B11/36Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
    • G05B11/42Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P.I., P.I.D.
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Abstract

本发明涉及一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法,属于磁传感技术领域。本发明针对无人机常用的航模锂电池,首先对描述锂电池和飞行轨迹的导航坐标系、地磁坐标系、机体坐标系、大地坐标系,以及方位角与欧拉角进行定义,并推导方位角与欧拉角之间的转换关系;根据锂电池的供能机理,结合锂电池与磁罗盘之间的空间位置关系,建立锂电池的电磁干扰机理模型;将多片锂电池电磁干扰机理模型的三轴分量投影到磁罗盘所在位置的地磁矢量上;对模型进行线性化,优化得到锂电池在机舱内的安置倾角,并采用反馈控制实现锂电池安置倾角对最优倾角进行实时跟踪,从而实现锂电池在磁罗盘处电磁干扰的自抵消。

Description

一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法
技术领域
本发明涉及一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法,适用于小型固定翼锂电池无人机磁罗盘处的电磁干扰抑制,属于磁传感技术领域。
背景技术
锂电池动力无人机具有振动低、噪声低、能源清洁环保、系统简单可靠的优势,可执行复杂多样的任务,因而得到了广泛的应用,具有广阔的应用前景。无人机依靠磁罗盘进行航向导航,磁罗盘测量无人机所在位置的地磁场强度大小,来判断无人机的航向。但是,随着任务需求的增多,锂电池无人机携带了大量的电子设备,狭小的机舱内电子设备密度激增,各种设备及电源产生的电磁干扰在时域上持续不断,在频域上互相重叠,使得锂电池无人机的磁罗盘面临着恶劣的电磁干扰环境。这种电磁干扰将导致磁罗盘对地磁场的测量误差增大,导致飞行控制器对航向误判,对飞行安全造成极大危害,使得无人机发生速度漂移,甚至飞行失稳、坠毁。
为降低电磁干扰对磁罗盘的干扰,需对复杂的电磁干扰进行分析。根据电磁干扰的产生来源,将无人机的电磁干扰分为与轨迹无关的电磁干扰和与轨迹相关的电磁干扰。开关电路和电机的高频电磁干扰为与轨迹无关的电磁干扰,其对磁罗盘的电磁干扰可通过低通滤波器进行滤波消除;通电导线和锂电池产生的电磁干扰是低频电磁干扰,与无人机的飞行需求功率相关,因而与轨迹相关,且无法通过滤波器消除。通电导线的电磁干扰可以通过导线双绞的方式消除。因此,锂电池的电磁干扰是对磁罗盘的主要影响因素,抑制锂电池产生的电磁干扰,对磁罗盘的正常工作具有重要意义,直接关系到无人机的航向测量和飞行安全。
在磁传感技术领域的实际应用中,根据电磁干扰产生、传播和影响的环节不同,抑制锂电池无人机中通电导线和锂电池的电磁干扰的方法主要分为2类:
第一类是通过电磁屏蔽的方法,通过设计合适的电磁屏蔽装置,将锂电池包裹包覆,隔断电磁干扰的传播途径。电磁屏蔽装置中,通常采用由磁屏蔽材料,如坡莫合金、镀锡铜丝,编织成屏蔽网或制成壳体结构,并通过一层或多层的组合搭配,既抑制电磁干扰对外辐射,又避免外部干扰对所包覆设备产生影响。但是,锂电池在使用过程中发热,电磁屏蔽装置的包覆导致装置内部温度升高,不利于锂电池工作,也不利于锂电池使用寿命。
第二类是通过电路设计的方法,对电池内部开关电路图进行特殊设计,并对电路板进行抗干扰印刷制作,在电磁干扰源头降低电磁干扰。电路设计和电路板制作不仅需要专业的电路设计知识,对广泛的无人机使用人员不够简单、友好,而且不能完全消除电磁干扰,对磁罗盘的影响依然存在。
综上所述,上述方法从电磁干扰的源头和传播途径入手,对影响磁罗盘的电磁干扰进行干预,但均缺乏对无人机内电磁干扰的模型描述,无法完全消除电磁干扰。同时,电磁屏蔽设备制作成本较高,电路设计需要专业的人员与知识。
近年来,通过建模描述系统行为得到了广泛的关注。建模方法可分为3类,第一类方法通过理论分析建立机理模型,第二类方法通过实验数据建立数据模型,第三类方法结合前两种方法建立混合模型。通过建立锂电池在磁罗盘处的电磁干扰机理模型,依据模型描述,在电磁干扰的影响对象处,精准降低电磁干扰,使对磁罗盘的影响最小,是降低电磁干扰影响的新思路。
发明内容
本发明的目的是为了解决锂电池对磁罗盘的电磁干扰,提供一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法,该方法通过对锂电池进行与飞行轨迹耦合的机理建模,从模型的数学机理角度,降低锂电池对无人机磁罗盘的电磁干扰。具有下述优点:成本低,操作简单,易实现,效果显著等优点。此外,本发明即使在缺乏专业知识背景的情况下也能进行操作,能够在磁传感技术工程应用中广泛应用于降低无人机电磁干扰的影响。
本发明的目的是通过下述技术方案实现的。
本发明公开的一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法,针对无人机常用的航模锂电池,首先对描述锂电池和飞行轨迹的导航坐标系、地磁坐标系、机体坐标系、大地坐标系,以及方位角与欧拉角进行定义,并推导方位角与欧拉角之间的转换关系;根据锂电池的供能机理,结合锂电池与磁罗盘之间的空间位置关系,建立锂电池的电磁干扰机理模型;将多片锂电池电磁干扰机理模型的三轴分量投影到磁罗盘所在位置的地磁矢量上;对模型进行线性化,优化得到锂电池在机舱内的安置倾角,并采用反馈控制实现锂电池安置倾角对最优倾角进行实时跟踪,从而实现锂电池在磁罗盘处电磁干扰的自抵消。
一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法,包括以下步骤:
步骤1:对描述锂电池和飞行轨迹的导航坐标系、地磁坐标系、机体坐标系、大地坐标系,以及方位角与欧拉角进行定义,并推导方位角与欧拉角之间的转换关系。
步骤1具体实现方法包括如下步骤:
步骤1.1:定义导航坐标系、地磁坐标系、机体坐标系和大地坐标系,如图2所示。导航坐标系采用北东地导航系,记作Σe,其xe轴指向地理正北,ye轴指向地理正东,ze轴指向地面。地磁坐标系,记作Σm,通过地磁场强度Be、磁偏角D和磁倾角I描述地磁矢量,磁偏角是地磁正北与地理正北的夹角,磁倾角是地磁矢量与地磁正北的夹角;地磁坐标系xm轴指向地磁正西,ym轴指向地磁正北,zm轴指向地面。机体坐标系与无人机机体固连,记作Σb,其xb轴指向左翼,yb轴指向机头,zb轴指向机腹。大地坐标系,记作Σg,其yg轴为航向的水平分量,zg轴指向地面,xg轴与yb和zg满足右手螺旋定则。无人机航向,记作θ,为yg轴与xe轴之间的夹角。无人机磁航向,记作
Figure GDA0003415163340000032
为yg轴与ym轴之间的夹角。那么,可以得到
Figure GDA0003415163340000033
步骤1.2:推导方位角与欧拉角的角度转换关系。地磁矢量与机体坐标系三轴的夹角为方位角,分别记作X,Y,Z,如图3所示。其中,X为机体坐标系xb轴与地磁矢量的夹角;Y为机体坐标系yb轴与地磁矢量的夹角;Z为机体坐标系zb轴与地磁矢量的夹角。无人机欧拉角为俯仰角λ,偏航角Ω和滚转角ψ。由于无人机飞行时仅能测量欧拉角,无法直接得知方位角,因此,需要进行方位角与欧拉角之间的转换,得到
Figure GDA0003415163340000031
其中,为了简化公式的表达形式,c代表cos,s代表sin。
步骤2:根据锂电池的供能机理,结合锂电池与磁罗盘之间的空间位置关系,建立锂电池的电磁干扰机理模型。
步骤2具体实现方法包括如下步骤:
步骤2.1:建立锂电池与磁罗盘之间的空间位置关系。
测量锂电池的尺寸,包括锂电池的宽度l,每一片锂电池的厚度d。如图4所示,测量锂电池和磁罗盘之间的相对距离,包括第i片锂电池中心与磁罗盘的距离Di,第i片锂电池中心与磁罗盘距离的水平分量l1,第i片锂电池中心与磁罗盘距离的垂直分量l2i,Λi为第i片锂电池中心到磁罗盘的方位角。那么,存在几何关系
Figure GDA0003415163340000041
其中,l21为第1片锂电池与磁罗盘距离的垂直分量。
步骤2.2:建立每一片锂电池的电磁干扰机理模型。
根据毕奥萨伐尔定理,锂电池内部长等于锂电池宽度l的电流Ib产生的电磁干扰BI
Figure GDA0003415163340000042
其中,μ0=4π×10-7V·s/(A·m)为真空磁导率,π为圆周率。
在磁罗盘处建立与机体坐标系固连的电磁干扰坐标系,记作Σc,xc轴指向左翼,与xb轴平行;yc轴指向机尾,与yb轴平行;zc轴指向机背,与zb轴平行,如图4所示。电磁干扰坐标系用于描述电磁干扰在机体坐标系xb轴、yb轴和zb轴的分量。
BI在磁罗盘处的三轴分量为BIx,BIy和BIz,表示为
Figure GDA0003415163340000043
其中,Θ为锂电池的安置倾角。
步骤2.3:建立多片锂电池的电磁干扰机理模型。
将多片锂电池在磁罗盘处产生的电磁干扰进行求和,得到多片锂电池在磁罗盘处的电磁干扰分量Bx,By和Bz
Figure GDA0003415163340000044
其中,N为锂电池总片数;C1,C2和C3为公式化简后的常数项;(-1)i表示第i片锂电池在磁罗盘处产生电磁干扰的方向。
步骤3:根据多片锂电池的电磁干扰机理模型,将多片锂电池在磁罗盘处的电磁干扰分量Bx,By和Bz投影到磁罗盘所在位置的地磁矢量上,此投影量BI为磁罗盘对地磁场的测量影响量,投影后,BI与航向和欧拉角耦合。
Figure GDA0003415163340000051
步骤4:由于无人机飞行时欧拉角变化较小,因此,将欧拉角视为小角度,将步骤3得到的电磁干扰投影量BI线性化。
步骤4具体实现方法包括如下步骤:
步骤4.1:将式线性化,得到
Figure GDA0003415163340000052
步骤4.2:将式化简,得到
Figure GDA0003415163340000053
步骤5:根据辅助角公式,优化得到最优锂电池安置倾角。
Figure GDA0003415163340000054
其中,
Figure GDA0003415163340000055
当无人机航向固定,欧拉角均为0时,
Figure GDA0003415163340000056
那么,当锂电池安置倾角满足
Figure GDA0003415163340000057
此时,BI=0,即,锂电池在磁罗盘处产生的电磁干扰对磁罗盘的影响实现自抵消,与电流大小无关。此时,锂电池的安置倾角为最优倾角,最优倾角记作Θ*。
步骤6:根据反馈控制实时调整锂电池的安置倾角,使无人机在飞行的过程中,其锂电池倾角对式获得的最优倾角实时跟踪,实现锂电池对磁罗盘的电磁干扰最小。
PID控制是工程上最常用、最成熟的反馈控制方法。为了实现锂电池对最优倾角的跟踪,作为优选,通过PID反馈控制实时跟踪锂电池最优倾角。
有益效果:
1、本发明公开的一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法,建立了锂电池在磁罗盘处电磁干扰的机理模型,在机理上实现锂电池在磁罗盘处产生的电磁干扰模型的定量描述,揭示了锂电池在磁罗盘处产生的电磁干扰与电流之间的关系。
2、本发明公开的一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法,在机理上实现锂电池在磁罗盘处产生的电磁干扰与飞行轨迹的耦合描述,揭示了锂电池在磁罗盘处产生的电磁干扰与飞行轨迹之间的关系,并通过模型线性化简化了模型的表述,易于工程应用,具有广泛的工程应用前景。
3、本发明公开的一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法,可仅通过调整锂电池在无人机机舱内的安置倾角,实现锂电池在磁罗盘处产生的电磁干扰影响的自抵消,无需其他复杂的电磁屏蔽装置或电子电路设计,有模型理论支持,易于工程实现,具有广泛的应用前景。
附图说明
图1为本发明公开的一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法的流程图;
图2为坐标系定义中的导航坐标系、地磁坐标系、大地坐标系;
图3为坐标系定义中的机体坐标系和方位角;
图4为锂电池在磁罗盘处生成的电磁干扰的几何关系示意图;
图5为具体实施方式中无人机的飞行轨迹;
图6为具体实施方式中采用和未采用本发明中方法时锂电池在磁罗盘处产生的电磁干扰的对比。
具体实施方式
为了更好地说明本发明的目的和优点,下面通过锂电池无人机点到点飞行的实施例,并结合附图与表格对本发明做出详细解释。
实施例:
本实施例的无人机采用6S锂电池,电池容量22Ah,放电倍率25C,宽9cm,单片锂电池厚度0.01m;锂电池距离磁罗盘l1为0.1m,l21为0.15m。无人机从预定的原点位置,坐标(0,0,100)m,飞往坐标(3500,1500,110)m处。如图1所示,本实施例公开的一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法,包括如下步骤:
步骤1:对描述锂电池和飞行轨迹的导航坐标系、地磁坐标系、机体坐标系、大地坐标系,以及方位角与欧拉角进行定义,并推导方位角与欧拉角之间的转换关系。
步骤1的实现方法如下:
步骤1.1:根据选取的飞行任务区域,选取当地的地理北向、地理东向和地面分别作为导航坐标系的xe轴,ye轴和ze轴。
步骤1.2:测量选取的飞行任务区域的地磁场信息,包括磁偏角D和磁倾角I,并以此标定地磁坐标系,指向地磁正西为xm轴,指向地磁正北位ym轴,指向地面为zm轴。飞行任务区域内,磁偏角D=π/36,磁倾角I=π/6。
步骤1.3:标定无人机,选取无人机质心位置为机体坐标系原点,指向左机翼为xb轴,指向机头为yb轴,指向机腹为zb轴。
步骤1.4:依据式建立方位角与欧拉角之间的转换关系。
步骤2:根据锂电池的供能机理,结合锂电池与磁罗盘之间的空间位置关系,建立无人机上的锂电池电磁干扰机理模型。
步骤2的实现方法如下:
步骤2.1:建立每一片锂电池在磁罗盘处产生的电磁干扰机理模型
Figure GDA0003415163340000071
考虑到
Figure GDA0003415163340000072
将BI的三轴分量化简得到
Figure GDA0003415163340000073
步骤2.2:建立6片锂电池的电磁干扰机理模型
Figure GDA0003415163340000074
步骤3:根据锂电池的电磁干扰机理模型,将其投影到磁罗盘所在位置的地磁矢量上,此投影量BI为磁罗盘对地磁场的测量影响量,投影后,BI与航向和欧拉角耦合。
Figure GDA0003415163340000081
步骤4:由于无人机飞行时欧拉角变化较小,因此,可将欧拉角视为小角度,将步骤3得到的电磁干扰投影量线性化,如式所示。并将式代入式,得到
BI=[cXC1sinΘ-(cYC2+cZC3)cosΘ]·Ib (15)
步骤5:根据辅助角公式,优化得到最优锂电池安置倾角。
Figure GDA0003415163340000082
步骤6:采用PID控制方法实现锂电池倾角Θ对最优倾角Θ*的实时跟踪。
步骤7:设计无人机的飞行轨迹,使无人机按照预定轨迹飞行,测量无人机在采用本发明方法和未采用本发明方法时磁罗盘测得的地磁线的电磁干扰量。
表1锂电池在磁罗盘处产生的电磁干扰BI对比
锂电池安置倾角 跟踪最优倾角
平均电磁干扰 3947nT 0nT
从本实施例的计算结果可以看出,本发明提出的一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法能够有效降低锂电池在磁罗盘处产生的电磁干扰,在本实施例中能够完全消除锂电池对磁罗盘的电磁干扰影响。在获取锂电池和磁罗盘相对位置关系、设计的飞行任务轨迹后,仅通过优化锂电池在机舱内的安置倾角,即可极大地降低锂电池对磁罗盘的干扰。
根据本实例公开的一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法,能够得到锂电池在磁罗盘处产生的电磁干扰的定量描述;能够得到锂电池在磁罗盘处产生的电磁干扰与飞行轨迹的耦合关系,具有较高的应用价值。此外,根据本实例公开的一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法,可以仅通过优化锂电池在机舱的安置角度实现锂电池在磁罗盘处产生的电磁干扰的自抵消,应用规则简单,符合无人机应用需求,对无人机应用人员友好,易于工程实现,具有广泛的应用前景。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例,用于解释本发明,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (1)

1.一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法,其特征在于:包括以下步骤:
步骤1:对描述锂电池和飞行轨迹的导航坐标系、地磁坐标系、机体坐标系、大地坐标系,以及方位角与欧拉角进行定义,并推导方位角与欧拉角之间的转换关系;
步骤1具体实现方法包括如下步骤:
步骤1.1:定义导航坐标系、地磁坐标系、机体坐标系和大地坐标系;导航坐标系采用北东地导航系,记作Σe,其xe轴指向地理正北,ye轴指向地理正东,ze轴指向地面;地磁坐标系,记作Σm,通过地磁场强度Be、磁偏角D和磁倾角I描述地磁矢量,磁偏角是地磁正北与地理正北的夹角,磁倾角是地磁矢量与地磁正北的夹角;地磁坐标系xm轴指向地磁正西,ym轴指向地磁正北,zm轴指向地面;机体坐标系与无人机机体固连,记作Σb,其xb轴指向左翼,yb轴指向机头,zb轴指向机腹;大地坐标系,记作Σg,其yg轴为航向的水平分量,zg轴指向地面,xg轴与yb和zg满足右手螺旋定则;无人机航向,记作θ,为yg轴与xe轴之间的夹角;无人机磁航向,记作
Figure FDA0003415163330000011
为yg轴与ym轴之间的夹角;那么,可以得到
Figure FDA0003415163330000012
步骤1.2:推导方位角与欧拉角的角度转换关系;地磁矢量与机体坐标系三轴的夹角为方位角,分别记作X,Y,Z;其中,X为机体坐标系xb轴与地磁矢量的夹角;Y为机体坐标系yb轴与地磁矢量的夹角;Z为机体坐标系zb轴与地磁矢量的夹角;无人机欧拉角为俯仰角λ,偏航角Ω和滚转角ψ;由于无人机飞行时仅能测量欧拉角,无法直接得知方位角,因此,需要进行方位角与欧拉角之间的转换,得到
Figure FDA0003415163330000013
其中,为了简化公式的表达形式,c代表cos,s代表sin;
步骤2:根据锂电池的供能机理,结合锂电池与磁罗盘之间的空间位置关系,建立锂电池的电磁干扰机理模型;
步骤2具体实现方法包括如下步骤:
步骤2.1:建立锂电池与磁罗盘之间的空间位置关系;
测量锂电池的尺寸,包括锂电池的宽度l,每一片锂电池的厚度d;测量锂电池和磁罗盘之间的相对距离,包括第i片锂电池中心与磁罗盘的距离Di,第i片锂电池中心与磁罗盘距离的水平分量l1,第i片锂电池中心与磁罗盘距离的垂直分量l2i,Λi为第i片锂电池中心到磁罗盘的方位角;那么,存在几何关系
Figure FDA0003415163330000021
其中,l21为第1片锂电池与磁罗盘距离的垂直分量;
步骤2.2:建立每一片锂电池的电磁干扰机理模型;
根据毕奥萨伐尔定理,锂电池内部长等于锂电池宽度l的电流Ib产生的电磁干扰BI
Figure FDA0003415163330000022
其中,μ0=4π×10-7V·s/(A·m)为真空磁导率,π为圆周率;
在磁罗盘处建立与机体坐标系固连的电磁干扰坐标系,记作Σc,xc轴指向左翼,与xb轴平行;yc轴指向机尾,与yb轴平行;zc轴指向机背,与zb轴平行;电磁干扰坐标系用于描述电磁干扰在机体坐标系xb轴、yb轴和zb轴的分量;
BI在磁罗盘处的三轴分量为BIx,BIy和BIz,表示为
Figure FDA0003415163330000023
其中,Θ为锂电池的安置倾角;
步骤2.3:建立多片锂电池的电磁干扰机理模型;
将多片锂电池在磁罗盘处产生的电磁干扰进行求和,得到多片锂电池在磁罗盘处的电磁干扰分量Bx,By和Bz
Figure FDA0003415163330000024
其中,N为锂电池总片数;C1,C2和C3为公式化简后的常数项;(-1)i表示第i片锂电池在磁罗盘处产生电磁干扰的方向;
步骤3:根据多片锂电池的电磁干扰机理模型,将多片锂电池在磁罗盘处的电磁干扰分量Bx,By和Bz投影到磁罗盘所在位置的地磁矢量上,此投影量BI为磁罗盘对地磁场的测量影响量,投影后,BI与航向和欧拉角耦合;
Figure FDA0003415163330000031
步骤4:由于无人机飞行时欧拉角变化较小,因此,将欧拉角视为小角度,将步骤3得到的电磁干扰投影量BI线性化;
步骤4具体实现方法包括如下步骤:
步骤4.1:将式线性化,得到
Figure FDA0003415163330000032
步骤4.2:将式化简,得到
Figure FDA0003415163330000033
步骤5:根据辅助角公式,优化得到最优锂电池安置倾角;
Figure FDA0003415163330000034
其中,
Figure FDA0003415163330000035
当无人机航向固定,欧拉角均为0时,
Figure FDA0003415163330000036
那么,当锂电池安置倾角满足
Figure FDA0003415163330000037
此时,BI=0,即,锂电池在磁罗盘处产生的电磁干扰对磁罗盘的影响实现自抵消,与电流大小无关;此时,锂电池的安置倾角为最优倾角,最优倾角记作Θ*;
步骤6:根据反馈控制实时调整锂电池的安置倾角,使无人机在飞行的过程中,其锂电池倾角对式获得的最优倾角实时跟踪,实现锂电池对磁罗盘的电磁干扰最小。
CN202110501762.XA 2021-05-08 2021-05-08 一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法 Active CN113218379B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110501762.XA CN113218379B (zh) 2021-05-08 2021-05-08 一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110501762.XA CN113218379B (zh) 2021-05-08 2021-05-08 一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法

Publications (2)

Publication Number Publication Date
CN113218379A CN113218379A (zh) 2021-08-06
CN113218379B true CN113218379B (zh) 2022-03-22

Family

ID=77093848

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110501762.XA Active CN113218379B (zh) 2021-05-08 2021-05-08 一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法

Country Status (1)

Country Link
CN (1) CN113218379B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113984032A (zh) * 2021-11-04 2022-01-28 北京理工大学 一种与飞行轨迹耦合的燃料电池无人机电磁干扰抑制方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4919142B2 (ja) * 2005-10-19 2012-04-18 アイチ・マイクロ・インテリジェント株式会社 磁気コンパス
CN106683378A (zh) * 2016-12-30 2017-05-17 厦门安胜网络科技有限公司 一种无人机干扰方法及系统
CN110672092B (zh) * 2019-09-24 2021-05-25 北京理工大学 一种降低固定翼无人机平台磁干扰的航迹生成方法
CN112525156A (zh) * 2020-11-29 2021-03-19 陈启渊 一种抗干扰性能好的导轮式固定测斜仪

Also Published As

Publication number Publication date
CN113218379A (zh) 2021-08-06

Similar Documents

Publication Publication Date Title
CN102722184B (zh) 单脉冲天线角度模拟跟踪方法
CN108563243A (zh) 一种基于改进rrt算法的无人机航迹规划方法
CN110672092B (zh) 一种降低固定翼无人机平台磁干扰的航迹生成方法
CN105891821B (zh) 一种机载下视测量目标的自动跟踪方法
CN113489560B (zh) 一种无人机非平稳空地mimo信道的几何随机建模方法
CN107861522A (zh) 无人靶机控制系统
CN113218379B (zh) 一种与飞行轨迹耦合的锂电池无人机电磁干扰抑制方法
CN108469838A (zh) 基于特高压线区电场等势面的无人机自主导航系统及方法
CN109708639A (zh) 飞行器平飞跟踪直线和圆弧路径的侧向制导指令生成方法
CN108254630A (zh) 一种短波天线方向图和增益的测量系统及方法
CN108254716A (zh) 一种基于粒子群算法的观测平台轨迹优化方法
CN108170154A (zh) 一种无人机多传感器正向摄影倾斜飞控调试方法
CN105119683A (zh) 一种基于实时嵌入式控制系统的无人机通信干扰对抗方法
CN204557220U (zh) 一种小型智能化无人机飞控导航系统
CN105676170B (zh) 一种用于多目标无源定位与跟踪的虚假点去除方法
CN115857538A (zh) 三维空间下满足落角约束的多飞行器协同制导方法
CN109917373B (zh) 运动补偿搜索的动平台雷达的动态规划检测前跟踪方法
CN115113651A (zh) 一种基于椭圆拟合的无人机长僚协同覆盖寻优方法
WO2021217341A1 (zh) 避障方法、可移动平台、控制设备和存储介质
CN109581080A (zh) 用于评估短波天线性能的空中测试设备
CN112241589B (zh) 一种基于系统可观测度的多水下无人艇协同运动编队设计方法
CN105510907B (zh) 一种基于强散射点目标探测的弱散射点目标跟踪逼近方法
CN113984032A (zh) 一种与飞行轨迹耦合的燃料电池无人机电磁干扰抑制方法
CN116206490A (zh) 一种用于跨域协同导航的导航信息源群优化交互方法
CN111539151B (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