CN108108523A - 一种飞机荷兰滚等效拟配初值选取方法 - Google Patents
一种飞机荷兰滚等效拟配初值选取方法 Download PDFInfo
- Publication number
- CN108108523A CN108108523A CN201711230947.1A CN201711230947A CN108108523A CN 108108523 A CN108108523 A CN 108108523A CN 201711230947 A CN201711230947 A CN 201711230947A CN 108108523 A CN108108523 A CN 108108523A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msup
- mfrac
- initial value
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种飞机荷兰滚等效拟配初值选取方法，属于飞机操稳系统设计技术领域。包括以下步骤：步骤一、获得侧滑角对脚蹬力高阶系统在0.1～10rad/s范围内的幅值和相位数据；步骤二、基于典型的二阶系统获取与步骤一中高阶传递函数等效的二阶系统；步骤三、依据幅值曲线的极大值与ω＝0.1处的幅值确定荷兰滚阻尼比初值；步骤四、计算荷兰滚频率初值ωd；步骤五、利用ω＝0.1处的幅值和频率初值计算增益K初值；步骤六、确定时间延迟值τ；本发明依据高阶系统的幅值和相位数据，获得等效拟配的初值，该初值与精确解接近，在后续求解中能保证满足等效拟配要求。该方法避免了反复调试初值的工作，提高工作效率，简化计算难度。
Description
技术领域
本发明属于飞机操稳系统设计技术领域，具体涉及一种飞机荷兰滚等效拟配初值选取方法。
背景技术
当前多数飞机设计时采用了控制增稳系统。控制增稳后的飞机模态远远多于未增稳飞机，直接应用成熟的飞行品质规范比较困难。目前比较典型的方法是等效拟配，即针对某一高阶系统，找到等效的低阶系统，两个系统在一定的频域范围内，相应输出量的差值满足一定的指标条件。等效的低阶系统可利用品质规范，进行模态分析、飞行品质评定等工作。
等效拟配原理是飞机高阶系统与低阶等效系统的幅值、相位差的指标函数取极小值，且拟配点的幅值差、相位差处于限制范围内。指标函数如下式：
式中：J为失配度；k为拟配点个数，常取20；ωi为拟配点，在0.1～10rad/s的对数坐标取平均点；K为加权系数，常取0.0175；G为对数幅值，Φ为相位；HOS为高阶系统，LOES为低阶系统。
在荷兰滚等效拟配时，低阶等效系统的表达式如下：
式中，需拟配的未知参数包括增益Kd，荷兰滚阻尼比ξd,荷兰滚频率ωd，时间延迟τ。失配度函数表达式复杂，使用最小二乘、牛顿迭代等算法时，计算初值会影响等效拟配的准确度，初值选取不当会得到不合理的结果，无法反应高阶系统的特性。
初值选取在等效拟配中是一项关键工作，目前初值选取并无统一的方法，主要依赖工程人员的经验。对于缺少经验的工程人员，如果初值选取不合适时，需要对初值多次调整试验，保证等效拟配的结果符合要求。
发明内容
本发明的目的：为了解决上述问题，本发明提出了一种飞机荷兰滚等效拟配初值选取方法，基于高阶系统的幅频和相频特性曲线，估计出计算需要的初值，用于荷兰滚等效拟配，能够快速获得与准确结果接近的一组初值，提高对高阶系统等效拟配的效率，减少初值调试的工作量。
本发明的技术方案：一种飞机荷兰滚等效拟配初值选取方法，包括以下步骤：
步骤一、获得侧滑角对脚蹬力高阶系统在0.1～10rad/s范围内的幅值和相位数据；
其高阶系统传递函数公式如下：
其中，a0-am是传递函数分子表达式的系数；b0-bn是传递函数分母表达式的系数；s是拉普拉斯变换算子；
步骤二、基于典型的二阶系统获取与步骤一中高阶传递函数等效的二阶系统；
式中，为固定增益，只平移幅值；第二部分为典型的二阶系统，第三部分e-τs为时延环节；
步骤三、依据幅值曲线的极大值与ω＝0.1处的幅值确定荷兰滚阻尼比初值；
若存在极大值，利用公式(3)及公式(4)计算极大值与ω＝0.1处的幅值Aω＝0.1之比Amax/Aω＝0.1，并计算阻尼比初值ξd；
若不存在极大值，幅值曲线单调递减，则为阻尼比初值；
步骤四、计算荷兰滚频率初值ωd；
当存在极大值时，根据极大值Amax求解其对应的频率ωAmax并计算荷兰滚频率初值ωd：
若不存在极大值时，寻找相角90°对应的频率ω90，并计算荷兰滚频率初值ωd：
步骤五、利用ω＝0.1处的幅值和频率初值计算增益K初值；
步骤六、确定时间延迟值τ＝0.1；
本发明技术方案的有益效果：本发明依据高阶系统的幅值和相位数据，获得等效拟配的初值，该初值与精确解接近，在后续求解中能保证满足等效拟配要求。该方法避免了反复调试初值的工作，提高工作效率，简化计算难度。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚，对本发明实施例中的技术方案进行更加详细的描述。所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明一种飞机荷兰滚等效拟配初值选取方法，快速获得与准确结果接近的一组初值，提高对高阶系统等效拟配的效率，减少初值调试的工作量。
本发明的工作原理如下：
如果高阶系统和低阶系统满足等效拟配的要求，高低阶系统的幅频、相频曲线应是基本一致的，可认为高阶系统具有低阶系统的幅频和相频特性。
荷兰滚低阶系统可表示为如下式。第一部分为固定增益，只平移幅值；第二部分为典型的二阶系统，第三部分e-τs为时延环节，只改变相角，时间延迟较小时影响可忽略。
典型的二阶系统表达式如下：
上述二阶系统据有如下特性：
当ωd＝0时，幅值Aωd＝0＝1；
当时，在幅值有极大值当时，幅值单调递减无峰值；
当ω＝ωn时，相角为90°。
依据等效拟配的要求，一般情况下，可认为高阶系统应有如下特性：
当低阶系统时，对应的高阶系统在附近存在与接近的峰值；
近似认为ω＝0.1与ω＝0的幅值相等，即
不考虑时延影响时，高阶系统ω＝ωn对应的相角为90°。
基于高阶系统上述特性，首先求解侧滑角对脚蹬力高阶系统的幅值相位曲线；然后依据幅值曲线的极大值与ω＝0.1处的幅值确定荷兰滚阻尼比初值；之后依据幅值曲线极大值、相位90°对应的频率确定荷兰滚频率初值；再利用ω＝0.1处的幅值和频率初值计算增益K初值；最后依据飞行品质规范要求给出时间延迟τ初值。
因此本发明一种飞机荷兰滚等效拟配初值选取方法，包括以下步骤：
步骤一、获得侧滑角对脚蹬力高阶系统在0.1～10rad/s范围内的幅值和相位数据；
其高阶系统传递函数公式如下：
其中，a0-am是传递函数分子表达式的系数；b0-bn是传递函数分母表达式的系数；s是拉普拉斯变换算子；
步骤二、基于典型的二阶系统获取与步骤一中高阶传递函数等效的二阶系统；
式中，为固定增益，只平移幅值；第二部分为典型的二阶系统，第三部分e-τs为时延环节；
步骤三、依据幅值曲线的极大值与ω＝0.1处的幅值确定荷兰滚阻尼比初值；
若存在极大值，利用公式(3)及公式(4)计算极大值与ω＝0.1处的幅值Aω＝0.1之比Amax/Aω＝0.1，并计算阻尼比初值ξd；
若不存在极大值，幅值曲线单调递减，则为阻尼比初值；
步骤四、计算荷兰滚频率初值ωd；
当存在极大值时，根据极大值Amax求解其对应的频率ωAmax并计算荷兰滚频率初值ωd：
若不存在极大值时，寻找相角90°对应的频率ω90，并计算荷兰滚频率初值ωd：
步骤五、利用ω＝0.1处的幅值和频率初值计算增益K初值；
步骤六、确定时间延迟值τ；
在俯仰轴和滚转轴飞行品质规范中，τ＝0.1s是等级1的上限，荷兰滚拟配时，时间延迟初值τ＝0.1。
至此，实现荷兰滚等效拟配的初值选取。
针对荷兰滚等效拟配初值选取的问题，本方法原理直观，计算简单，仅依据高阶系统的幅值和相位数据，获得等效拟配的初值，该初值与精确解接近，在后续求解中能保证满足等效拟配要求。该方法避免了反复调试初值的工作，提高工作效率，简化计算难度。
最后需要指出的是：以上实施例仅用以说明本发明的技术方案，而非对其限制。尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (1)
1.一种飞机荷兰滚等效拟配初值选取方法，其特征在于：包括以下步骤：
步骤一、获得侧滑角对脚蹬力高阶系统在0.1～10rad/s范围内的幅值和相位数据；
其高阶系统传递函数公式如下：
<mrow>
<mfrac>
<mi>&beta;</mi>
<msub>
<mi>F</mi>
<mi>r</mi>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<msup>
<mi>s</mi>
<mi>m</mi>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msup>
<mi>s</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>+</mo>
<mo>...</mo>
<mo>+</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>a</mi>
<mi>m</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<msup>
<mi>s</mi>
<mi>n</mi>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<msup>
<mi>s</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>+</mo>
<mo>...</mo>
<mo>+</mo>
<msub>
<mi>b</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>b</mi>
<mi>n</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，a0-am是传递函数分子表达式的系数；b0-bn是传递函数分母表达式的系数；s是拉普拉斯变换算子；
步骤二、基于典型的二阶系统获取与步骤一中高阶传递函数等效的二阶系统；
<mrow>
<mfrac>
<mi>&beta;</mi>
<msub>
<mi>F</mi>
<mi>r</mi>
</msub>
</mfrac>
<mo>=</mo>
<mfrac>
<mi>K</mi>
<msubsup>
<mi>&omega;</mi>
<mi>d</mi>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>&CenterDot;</mo>
<mfrac>
<msubsup>
<mi>&omega;</mi>
<mi>d</mi>
<mn>2</mn>
</msubsup>
<mrow>
<msup>
<mi>s</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&xi;</mi>
<mi>d</mi>
</msub>
<msub>
<mi>&omega;</mi>
<mi>d</mi>
</msub>
<mi>s</mi>
<mo>+</mo>
<msubsup>
<mi>&omega;</mi>
<mi>d</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mfrac>
<mo>&CenterDot;</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>&tau;</mi>
<mi>s</mi>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，为固定增益，只平移幅值；第二部分为典型的二阶系统，第三部分e-τs为时延环节；
步骤三、依据幅值曲线的极大值与ω＝0.1处的幅值确定荷兰滚阻尼比初值；
若存在极大值，利用公式(3)及公式(4)计算极大值与ω＝0.1处的幅值Aω＝0.1之比Amax/Aω＝0.1，并计算阻尼比初值ξd；
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>/</mo>
<msub>
<mi>A</mi>
<mrow>
<mi>&omega;</mi>
<mo>=</mo>
<mn>0.1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>2</mn>
<msub>
<mi>&xi;</mi>
<mi>d</mi>
</msub>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msubsup>
<mi>&xi;</mi>
<mi>d</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>&omega;</mi>
<mo>=</mo>
<mn>0.1</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mrow>
<mo>|</mo>
<mfrac>
<mrow>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>+</mo>
<mo>...</mo>
<mo>+</mo>
<msub>
<mi>a</mi>
<mrow>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>a</mi>
<mi>m</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mi>n</mi>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>+</mo>
<mo>...</mo>
<mo>+</mo>
<msub>
<mi>b</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>b</mi>
<mi>n</mi>
</msub>
</mrow>
</mfrac>
<mo>|</mo>
</mrow>
<mrow>
<mi>&omega;</mi>
<mo>=</mo>
<mn>0.1</mn>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
若不存在极大值，幅值曲线单调递减，则为阻尼比初值；
步骤四、计算荷兰滚频率初值ωd；
当存在极大值时，根据极大值Amax求解其对应的频率ωAmax并计算荷兰滚频率初值ωd：
<mrow>
<msub>
<mi>&omega;</mi>
<mi>d</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>&omega;</mi>
<mrow>
<mi>A</mi>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<msup>
<mi>&xi;</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
若不存在极大值时，寻找相角90°对应的频率ω90，并计算荷兰滚频率初值ωd：
<mrow>
<msub>
<mi>&omega;</mi>
<mi>d</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>&omega;</mi>
<mn>90</mn>
</msub>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mn>2</mn>
<msup>
<mi>&xi;</mi>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤五、利用ω＝0.1处的幅值和频率初值计算增益K初值；
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>&omega;</mi>
<mo>=</mo>
<mn>0.1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mi>K</mi>
<msubsup>
<mi>&omega;</mi>
<mi>d</mi>
<mn>2</mn>
</msubsup>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤六、确定时间延迟值τ＝0.1。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711230947.1A CN108108523B (zh) | 2017-11-29 | 2017-11-29 | 一种飞机荷兰滚等效拟配初值选取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711230947.1A CN108108523B (zh) | 2017-11-29 | 2017-11-29 | 一种飞机荷兰滚等效拟配初值选取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108108523A true CN108108523A (zh) | 2018-06-01 |
| CN108108523B CN108108523B (zh) | 2021-07-09 |
Family
ID=62208891
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711230947.1A Active CN108108523B (zh) | 2017-11-29 | 2017-11-29 | 一种飞机荷兰滚等效拟配初值选取方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108108523B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109711034A (zh) * | 2018-12-24 | 2019-05-03 | 西北工业大学 | 一种飞机起落架健康监测方法 |
| CN110598295A (zh) * | 2019-09-03 | 2019-12-20 | 中国航空工业集团公司西安飞机设计研究所 | 飞机纵向短周期拟配方法及飞机 |
| CN111428311A (zh) * | 2020-03-24 | 2020-07-17 | 中国人民解放军空军工程大学 | 基于平均差分进化算法的等效系统拟配方法 |
| CN113128745A (zh) * | 2020-01-16 | 2021-07-16 | 中车青岛四方机车车辆股份有限公司 | 一种控制参数优化方法及装置 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5669582A (en) * | 1995-05-12 | 1997-09-23 | The Boeing Company | Method and apparatus for reducing unwanted sideways motion in the aft cabin and roll-yaw upsets of an airplane due to atmospheric turbulence and wind gusts |
| US20060237594A1 (en) * | 2002-05-01 | 2006-10-26 | Kioumars Najmabadi | Aircraft multi-axis modal suppression system |
| CN106507927B (zh) * | 2005-05-20 | 2010-03-10 | 成都飞机工业（集团）有限责任公司 | 用时域气动等效确定飞机运动模态的方法 |
| US20140114505A1 (en) * | 2012-02-27 | 2014-04-24 | Textron Innovations Inc. | Yaw Damping System And Method For Aircraft |
| CN104573177A (zh) * | 2014-11-28 | 2015-04-29 | 中国飞行试验研究院 | 运输机试飞员飞行品质评价训练的模拟方法 |
| CN105868893A (zh) * | 2016-03-24 | 2016-08-17 | 四川汉科计算机信息技术有限公司 | 固定翼无人机横航向动稳定性飞行品质评估系统 |
-
2017
- 2017-11-29 CN CN201711230947.1A patent/CN108108523B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5669582A (en) * | 1995-05-12 | 1997-09-23 | The Boeing Company | Method and apparatus for reducing unwanted sideways motion in the aft cabin and roll-yaw upsets of an airplane due to atmospheric turbulence and wind gusts |
| US20060237594A1 (en) * | 2002-05-01 | 2006-10-26 | Kioumars Najmabadi | Aircraft multi-axis modal suppression system |
| CN106507927B (zh) * | 2005-05-20 | 2010-03-10 | 成都飞机工业（集团）有限责任公司 | 用时域气动等效确定飞机运动模态的方法 |
| US20140114505A1 (en) * | 2012-02-27 | 2014-04-24 | Textron Innovations Inc. | Yaw Damping System And Method For Aircraft |
| CN104573177A (zh) * | 2014-11-28 | 2015-04-29 | 中国飞行试验研究院 | 运输机试飞员飞行品质评价训练的模拟方法 |
| CN105868893A (zh) * | 2016-03-24 | 2016-08-17 | 四川汉科计算机信息技术有限公司 | 固定翼无人机横航向动稳定性飞行品质评估系统 |
Non-Patent Citations (4)
| Title |
|---|
| MOHAMMAD HOSSEIN MOHAJERANI ET AL.: "Identification of Low Order Equivalent Transfer Function Model of Trex-700E Helicopter from Flight Test Data", 《AIAA MODELING AND SIMULATION TECHNOLOGIES CONFERENCE》 * |
| 支真莉: "虚拟试飞中的横航向品质研究", 《航空技术推广》 * |
| 江飞鸿等: "一种新的等效拟配初值计算方法", 《飞行力学》 * |
| 谭文倩等: "飞机横航向等效拟配模型的对比研究", 《飞行力学》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109711034A (zh) * | 2018-12-24 | 2019-05-03 | 西北工业大学 | 一种飞机起落架健康监测方法 |
| CN109711034B (zh) * | 2018-12-24 | 2022-04-15 | 西北工业大学 | 一种飞机起落架健康监测方法 |
| CN110598295A (zh) * | 2019-09-03 | 2019-12-20 | 中国航空工业集团公司西安飞机设计研究所 | 飞机纵向短周期拟配方法及飞机 |
| CN113128745A (zh) * | 2020-01-16 | 2021-07-16 | 中车青岛四方机车车辆股份有限公司 | 一种控制参数优化方法及装置 |
| CN113128745B (zh) * | 2020-01-16 | 2022-12-13 | 中车青岛四方机车车辆股份有限公司 | 一种控制参数优化方法及装置 |
| CN111428311A (zh) * | 2020-03-24 | 2020-07-17 | 中国人民解放军空军工程大学 | 基于平均差分进化算法的等效系统拟配方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108108523B (zh) | 2021-07-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108108523A (zh) | 一种飞机荷兰滚等效拟配初值选取方法 | |
| US10334032B2 (en) | Energy-saving deployment method of virtual CDN | |
| CN104199295B (zh) | 基于神经网络的机电伺服系统摩擦补偿和变结构控制方法 | |
| CN104376581B (zh) | 一种采用自适应重采样的高斯混合无迹粒子滤波算法 | |
| CN106802660A (zh) | 一种复合强抗扰姿态控制方法 | |
| CN110417011A (zh) | 一种基于互信息与迭代随机森林的在线动态安全评估方法 | |
| CN105512502B (zh) | 一种基于残差归一化的权函数最小二乘状态估计方法 | |
| CN104316049A (zh) | 高精度低信噪比椭圆化星点光斑细分定位方法 | |
| CN110032706A (zh) | 一种低阶时滞系统的两阶段参数估计方法及系统 | |
| CN107290742A (zh) | 一种非线性目标跟踪系统中平方根容积卡尔曼滤波方法 | |
| CN110687781A (zh) | 一种二阶非线性系统的精确不确定性和控制增益估计方法 | |
| CN109540459A (zh) | 一种气动特性数值计算结果修正方法 | |
| CN115688288B (zh) | 飞行器气动参数辨识方法、装置、计算机设备及存储介质 | |
| CN114326405B (zh) | 一种基于误差训练的神经网络反步控制方法 | |
| CN104965506B (zh) | 一种基于可调参数分布式飞行控制系统实时故障诊断方法 | |
| CN105892297A (zh) | 一种自适应分数阶动态滑模控制算法 | |
| CN111669150A (zh) | 一种基于采样点调整的无迹卡尔曼滤波方法 | |
| CN108120452A (zh) | Mems陀螺仪动态数据的滤波方法 | |
| CN114063458A (zh) | 一种非三角结构系统不依赖初始条件的预设性能控制方法 | |
| CN114611420A (zh) | 非定常气动力计算精度评估及修正方法 | |
| CN106599381A (zh) | 通过调整风轮桨距角及转速提高风力机效率的方法 | |
| CN103473470B (zh) | 一种地面效应风洞试验数据处理方法 | |
| CN109904863B (zh) | 一种附加阻尼控制器、水轮机调速器及水轮机控制系统 | |
| CN110598295B (zh) | 飞机纵向短周期拟配方法及飞机 | |
| Xu et al. | Synchronization time and energy consumption for multiweighted complex networks |
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 |