一种弹性飞翼布局飞机阵风减缓控制器参数优化方法
【专利摘要】本发明涉及一种弹性飞翼布局飞机阵风减缓控制器参数优化方法，属于飞机气动弹性设计领域。所述方法包括首先选择作为阵风减缓控制器输入的若干个机体运动量，以及确定阵风减缓控制器的通道个数；之后构建由滤波器阻尼系数构成的滤波器表达式，并由所述滤波器表达式与所述机体运动量输入信号的控制增益参数建立所述阵风减缓控制器的传递函数矩阵；最终通过多目标优化问题，设置若干个优化目标，优化待求参数，所述待求参数包括所述滤波器阻尼系数及所述控制增益参数，通过该方法保留对飞机阵风响应贡献较大的模态振动信息进入控制回路，在确保阵风减缓效果的同时，有效降低发生不利气动伺服弹性耦合的概率。
【专利说明】
一种弹性飞翼布局飞机阵风减缓控制器参数优化方法
技术领域
[0001] 本发明属于飞机气动弹性设计领域，属飞机气动弹性主动控制范畴，尤其涉及一 种弹性飞翼式布局飞机阵风减缓设计方法。
【背景技术】
[0002] 飞翼布局飞机翼载小，对阵风扰动非常敏感，特别当飞机尺寸和重量增加，结构固 有振动频率降低后，阵风扰动更易激励起低频弹性振动，致使结构所承受的阵风载荷进一 步增大，因此对于大型飞翼布局飞机，在进行阵风减缓设计时，必须考虑结构弹性的影响。
[0003] 飞翼布局飞机的升降副翼偏转不仅改变翼面升力分布，同时还会产生较大的俯仰 力矩变化，如何确定各控制通道的反馈增益，确保飞机姿态稳定的同时减缓阵风扰动产生 的结构附加载荷，对于飞翼布局实现阵风减缓至关重要。此外，当考虑结构弹性影响后，还 必须确保飞机与阵风减缓系统所组成的结构-控制耦合回路拥有合理的气动伺服弹性 (ASE)稳定裕度。
【发明内容】
[0004] 为了解决上述问题，本发明提供了一种弹性飞翼布局飞机阵风减缓控制器参数优 化方法，将阵风减缓设计转化为多目标优化问题，基于气动伺服弹性稳定裕度修正个体适 应度，建立改进的多目标非支配排序遗传-阵风减缓算法(NSGA_II_GLA算法)优化控制器反 馈增益参，确保飞机气动伺服弹性稳定裕度满足GJB67.7A-2008要求的同时，改善飞机的 阵风响应特性，实现阵风载荷减缓，提高控制器参数设计效率，其主要包括以下步骤：
[0005] S1、选择作为阵风减缓控制器输入的若干个机体运动量；
[0006] S2、确定阵风减缓控制器的通道个数；
[0007] S3、计算无控飞机的任一机体运动量对连续阵风输入的功率谱密度函数，并在所 述功率谱密度函数对应的功率谱密度曲线中自起始点开始选择至少一个峰值；
[0008] S4、对所述阵风减缓控制器的任一通道，构建在前一步所选峰值的峰值频率下由 滤波器阻尼系数构成的滤波器表达式；
[0009] S5、由所述滤波器表达式与所述机体运动量输入信号的控制增益参数建立所述阵 风控制减缓控制器的传递函数矩阵；
[0010] S6、根据所述传递函数矩阵设置若干个优化目标，优化待求参数，所述待求参数包 括所述滤波器阻尼系数及所述控制增益参数。
[0011 ]优选的是，所述机体运动量包括飞机的俯仰角速率《 y、重心过载气s以及翼尖过载 ft 0 wjip
[0012]在上述方案中优选的是，对进入任意第m个控制通道的任意第k个机体运动量rk， 选择nk个峰值频率.，其中j = 1，…，nk，设计nk个滤波器，所述滤波器表达式为：
[0014] 其中，s为拉普拉斯算子，a为滤波器阻尼系数。
[0015] 在上述方案中优选的是，在步骤S6中，所述优化目标包括使有控飞机翼根弯矩响 应的均方根值最小、有控飞机俯仰角速率响应的均方根值最小以及使最大舵偏角均方根值 最小情况下的优化参数组合。
[0016] 在上述方案中优选的是，在步骤S6中，所述优化过程包括：
[0017] T1、随机生成第一代染色体种群Po，染色体为优化变量X;
[0018] T2、计算Po中个体的优化目标的函数值，并根据个体的幅值和相位裕度设计罚函 数，修正个体适应度；
[0019] T3、将Po中个体按照其适应度进行非支配排序，计算同一非支配层级中个体拥挤 度；
[0020] T4、从Po中优先选择非支配层级编号小的个体生成新种群Pi，所述Pi种群内个体数 目为P〇内个体数目的一半；
[0021] T5、对种群丹中的个体进行交叉、变异操作，生成个体规模与Po相同的新种群Qo;
[0022] T6、计算种群Qo中个体的目标函数值，采用步骤T2的方法计算种群Qo中个体的适应 度；
[0023] T7、种群Po和Qo合并组成新种群&，对&中的个体进行非支配排序和并计算个体拥 挤度，并优先选择Q:中的个体的非支配层级编号小的个体生成新种群Q 2，所述Q2种群内个体 数目与P〇内的个体数目相同；
[0024] T8、重复步骤T3-T7，直至达到给定的最大进化代数，优化出待求参数。
[0025] 在上述方案中优选的是，在步骤T2中，根据个体X△对幅值和相位裕度的满足程度， 设计罚函数私胃和，公式包括：
[0028]其中，为若干个控制通道的幅值裕度的最小值，仏_为若干个控制通道的相位 裕度的最小值；
[0029] 之后，利用罚函数修正规范化处理后的优化目标函数，得到个体X△的适应度，修正 优化目标函数0〃i的公式为
[0030] 〇-；=〇-； + + Rfai\i = 1,2, ? ? ?,?
[0031] 所述(/ i为经规范化处理后的优化目标函数，所述规范化处理包括将所述优化目 标函数取倒数。
[0032] 在上述方案中优选的是，在步骤T5中，所述交叉包括二进制交叉，概率为0.9,所述 变异包括多项式变异。
[0033]在上述方案中优选的是，所述种群规模为80,迭代次数为100次。
[0034]本发明的有益效果在于：
[0035] a.在确定阵风减缓控制器结构时，根据无控飞机连续阵风响应功率谱密度曲线， 选择功率谱密度曲线的峰值频率作为滤波器设计参数，只保留对飞机阵风响应贡献较大的 模态振动信息进入控制回路，在确保阵风减缓效果的同时，有效降低发生不利气动伺服弹 性耦合的概率；
[0036] b .在进行控制器参数优化时，基气动伺服弹性稳定裕度对体的适应度进行修正， 能够确保优化结果实现阵风减缓的同时，满足气动伺服弹性稳定性设计要求，同时采用基 于随机联赛选择的精英保留策略，能够有效提升优化运行效率。
【附图说明】
[0037] 图1为本发明弹性飞翼布局飞机阵风减缓控制器参数优化方法的一优选实施例的 流程图。
[0038] 图2为图1所示实施例的通道控制指令解算示意图。
[0039] 图3为图1所示实施例的功率谱密度曲线示意图。
[0040] 图4为图1所示实施例的翼尖过载功率谱密度曲线示意图。
[0041] 图5为图1所示实施例的重心过载功率谱密度曲线示意图。
[0042] 图6为图1所示实施例的俯仰角过载功率谱密度曲线示意图。
[0043] 图7为图1所示实施例的迭代优化示意图。
[0044] 图8为图1所示实施例的优化结果示意图。
【具体实施方式】
[0045] 为使本发明实施的目的、技术方案和优点更加清楚，下面将结合本发明实施例中 的附图，对本发明实施例中的技术方案进行更加详细的描述。在附图中，自始至终相同或类 似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本发明 一部分实施例，而不是全部的实施例。下面通过参考附图描述的实施例是示例性的，旨在用 于解释本发明，而不能理解为对本发明的限制。基于本发明中的实施例，本领域普通技术人 员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。下 面结合附图对本发明的实施例进行详细说明。
[0046] 在本发明的描述中，需要理解的是，术语"中心"、"纵向"、"横向"、"前"、"后"、 "左"、"右"、"竖直"、"水平"、"顶"、"底" "内"、"外"等指示的方位或位置关系为基于附图所 示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装 置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明保护 范围的限制。
[0047] 下面通过实施例对本发明做进一步详细说明。
[0048] 如图1所示，本发明弹性飞翼布局飞机阵风减缓控制器参数优化方法主要包括以 下步骤：
[0049] S1、选择作为阵风减缓控制器输入的若干个机体运动量；
[0050] S2、确定阵风减缓控制器的通道个数；
[0051] S3、计算无控飞机的任一机体运动量对连续阵风输入的功率谱密度函数，并在所 述功率谱密度函数对应的功率谱密度曲线中自起始点开始选择至少一个峰值；
[0052] S4、对所述阵风减缓控制器的任一通道，构建在前一步所选峰值的峰值频率下由 滤波器阻尼系数构成的滤波器表达式；
[0053] S5、由所述滤波器表达式与所述机体运动量输入信号的控制增益参数建立所述阵 风控制减缓控制器的传递函数矩阵；
[0054] S6、根据所述传递函数矩阵设置若干个优化目标，优化待求参数，所述待求参数包 括所述滤波器阻尼系数及所述控制增益参数。
[0055] 下面对上述步骤进行详细说明，在步骤1中，飞机的机体运动量是指在飞机飞行或 停止状态中，对机体运动特性产生影响的参数值。比如，用于阵风减缓控制器输入的机体运 动量应包含但不限于飞机的俯仰角速率《 y、重心过载气s以及翼尖过载I# °
[0056] 需要说明的是，本实施例机体运动量的选取按实际需要选取，对于任意选取的若 干个机体运动量，通过本发明所述的步骤最终求得与这些机体运动量相关联的滤波器阻尼 系数及所述控制增益参数，以达到优化阵风减缓控制器参数的目的。
[0057]如图2所示，本实施例在选取上述三个机体运动量之后，用以进行通道控制指令解 算，因此，这里首先需要确定阵风减缓控制器的结构，具体的是，确定所述阵风减缓控制器 的通道个数，这是因为，对所述每一个阵风减缓控制器的通道，都需要通过所选机体运动量 来参与通道控制指令解算，假设通过对所述阵风减缓控制器的结构分析发现其包含m个通 道，那么针对k个机体运动量，当需要进行通道控制指令解算时，比如进行舵偏指令S m解算 之前，首先对这k个机体运动量进行多次滤波，该步骤称作滤波器设置。本实施例中，m为3， Si、s2和s3分别为外、中和内升降副翼偏角。
[0058] 需要说明的是，在进行滤波器设置时，首先要计算机体运动变量rk的功率谱密度 函数，从而绘制变量r k的功率谱密度曲线，如图3所示，从所述曲线中，选择前nk个峰值碎对 应的峰值频率?f，其中j = 1，…，nk，用于滤波器设计，可以理解的是，通过该步骤，将设计 有nk个滤波器，需要说明的是，所述峰值至少选择一个，也就是说对任一个机体运动变量r k 至少设计一个滤波器，即nk彡1。
[0059] 本实施例以前述三个机体运动量:俯仰角速率coy、重心过载以及翼尖过载 给出其功率谱密度曲线示意图，分别如图4、图5及图6所示，其中，图4为翼尖过载的功率谱 密度曲线，图5为重心过载的功率谱密度曲线，图6为俯仰角过载的功率谱密度曲线，图4中， 对翼尖过载设计一个滤波器，即其功率谱密度曲线中，选择第一响应峰值(纵坐标)及对应 的峰值频率(横坐标），图4中，对于翼尖过载，选取n 3 = 1，其功率谱密度曲线峰值频率为 < =9,4/W/'.v潮5中，对重心过载设计两个滤波器，即其功率谱密度曲线中，选择第一响应 峰值与第二响应峰值(纵坐标)及对应的峰值频率(横坐标），图6中，对俯仰角过载设计两个 滤波器，即其功率谱密度曲线中，选择第一响应峰值与第二响应峰值(纵坐标)及对应的峰 值频率(横坐标），图5及图6中，对于俯仰角速率《 y和重心过载，选取m=n2 = 2，其功率 谱密度曲线峰值频率分别为叫' =4: =9.4/m〃. s'。
[0060]进一步需要说明的是，计算无控飞机的阵风响应分析模型中，除了刚体俯仰和刚 体沉浮模态之外，应至少包含飞机的前10支对称弹性振动模态，本实施例中，选取飞机的刚 体俯仰、沉浮和前20支弹性模态，计算无控飞机的俯仰角速率《 y、重心过载\和翼尖过载 对冯卡门连续阵风输入的功率谱密度函数，该计算过程为现有技术，其计算公式可参阅 《飞机设计手册》北京:航空工业出版社，2001:1243-1248。
[0061 ] 本实施例对进入任意第m个控制通道的任意第k个机体运动量rk，选择nk个峰值频 率<??_，.其中i = l,…，nk，设计nk个滤波器，其公式如下。
[0063]上式中，s为拉普拉斯算子，^为滤波器阻尼系数。
[0064] 再次参考图2,将上述三个机体运动量代入以上求滤波器的公式中，确定出以下3 个控制通道5个峰值(可再次参考图4-图6中的共计5个峰值)的共计15个滤波器表达式：
[0074] 之后，根据上一步所设计的滤波器，建立阵风减缓控制器传递函数矩阵Gv的计算 表达式： ??? GnG'n---G\niKk
[0075] (7(^)= G2)G22'--G2nj/cll …啤與… G：Ltk2k G^l2-Gl,kml G；^2-G2Jm2 - Gl^2-Glkmh_
[0076] 如前所述，这里的m为通道数量，kmk为进入第m个控制通道中的第k个机体运动量输 入信号的控制增益。
[0077] 对本实施例给出的15个滤波器，代入上述公式，即 'G；1C?1Uu G^G^kn G；]ku
[0078] (?(s)= G'lk2i
[g;,g;^31 g^u32 0ik2i_
[0079] 之后，建立有控飞机连续阵风响应分析模型，将阵风减缓转换为多目标优化问题。
[0080] 计算有控飞机的阵风响应分析模型中，除了刚体俯仰和刚体沉浮模态之外，应至 少包含飞机的前10支对称弹性振动模态，优化目标如下。
[0082]其中：
[0087] %〕:为第j个控制通道的舵偏角均方根值、Gm/为第j个控制通道的幅值裕度、为 第j个控制通道的相位裕度，其中，j = l，2,…，m，即m个通道内的某一通道，其中优化目标 〇i，i = 1，2，…，n为有控飞机第i个响应量的均方根值;优化目标为m个舵面偏角均方根值 中的最大值;优化目标应包含但不限于翼根弯矩均方根值和最大舵偏角均方根值；代表 第m个控制通道，对第k个机体运动量rk输入信号所设置的第nk个滤波器的阻尼参数;为 /-I 选取的峰值个数。
[0088] 本实施例中仍以上述3个通道15个滤波器为例，即m为3,^]?/为5,其优化目标函数 1=1 如下： min =/,(X) min crc. = (X)
[0089] ^ min = maxiV(f (X),af (A'),rr(! (A )] > 6c/7i
[0090] 其中：
[0091] =mm(6；, (X),G,tl： (X),CrH< (X))
[0092] ^ =mm(\Pmi (X)|,|^ (X)|,|^ (X)|)
[0093] X=[Xo Xi X2 Xs]T
[0094] Xo= [kii ki2 ki3 k2i k22 k23 k3i k32 k33]T
[0095]
[0096]
[0097] x3=[aaciCiC0r
[0098] Xo(p,l)G[-3,3],p=l,2,---,9
[0099] Xi(qi,l)G[0.01,l],qi = l,2,---,5
[0100] X2(q2,l)G[0.01,l],q2 = l,2,---,5
[0101] X3(q3,l)G[0.01,l],q3 = l,2,---,5
[0102] 其中优化目标为有控飞机翼根弯矩响应的均方根值；
[0103] 其中优化目标〇q为有控飞机俯仰角速率响应的均方根值；
[0104] 其中优化目标os为3个舵偏角均方根值中的最大值。
[0105] 另外，对于闭环飞机连续阵风响应量均方根值、幅值和相位稳定裕度的计算，可参 阅《飞机设计手册》总编委会编.飞机设计手册第9册载荷、强度和刚度[M].北京:航空工业 出版社，2001。
[0106] 再之后，进行迭代优化，设定优化迭代代数、初始种群规模，改进多目标遗传优化 算法NSGA-II，从而求解上述步骤所建立的多目标优化问题。
[0107] 基气动伺服弹性稳定裕度建立罚函数，对个体适应度进行修正，采用规模为2的随 机联赛选择策略，有效加快了算法的收敛性。
[0108] 优化过程如图7所示，主要包括以下步骤：
[0109] Stepl:随机生成第一代染色体种群Po，染色体为优化变量X。本实施例中，设定优 化迭代代数gen = 100、初始种群规模为80。
[0110] Step2:计算P〇中个体目标函数值〇i，i = 1，2，…，n，并根据个体的幅值和相位裕度 设计罚函数，修正个体适应度；即在具体实施例中，计算上述Po中个体目标函数值^^~^ 和，并根据个体的幅值和相位裕度设计罚函数，修正个体适应度，具体的：
[0111] a.对个体Xa(下标A代表第A条染色体)所对应的每一个目标函数进行规范化处 理
[0115] b.根据个体Xa对幅值和相位裕度的满足程度，设计罚函数iC'和私'W
[0118] c.利用罚函数修正规范化处理后的目标函数，得到个体X△的适应度
[0119] erf = crj + if"1? + , i ~ 1,2, ? ?, n
[0120] S" = S' + R2l,g +RfttSL'
[0121] step3:将Po中个体按照其适应度进行非支配排序，计算同一非支配层级中个体拥 挤度；
[0122] Step4:采用规模为2的随机联赛选择方法，从Po中优先选择非支配层级编号小(第 1层最小）的个体，当个体的非支配层级相同时，再选择拥挤度较大的个体填充至种群Pi，直 至种群Pi中染色体数目等于Po中染色体数目的一半，本实施例中为40;
[0123] Step5 :对种群Pi中的个体进行交叉、变异操作，生成个体规模与Po相同的新种群 Q〇，需要说明的是，在本实施例中，模拟二进制交叉，概率为〇. 9，变异方式为多项式变异，概 率为1 /nchorm，nchorm = 24为单个染色体中的优化变量个数，Qo为80;
[0124] Step6:计算种群Qo中个体的目标函数值，采用Step2方法计算种群Qo中个体的适应 度；
[0125] Step7 :种群Po和Qo合并组成新种群&，对&中的个体进行非支配排序和并计算个 体拥挤度，并优先选择&中的个体的非支配层级编号小(第1层最小）的个体，对非支配层级 相同的个体，对同一非支配层级中的个体，优先选择拥挤度大的个体填充至种群Q2,直至种 群Q2中染色体数目与Po中染色体数目相同，本实施例中均为80;
[012 6 ] S t e p 8 :判断是否达给定的最大进化代数，达到则算法结束，否则令P 〇 = Q 2，转到 step3〇
[0127] 最终得到Pareto最优解集，参见附图8,选择Pareto前沿的下端面点所对应的优化 变量值作为阵风减缓控制器结构参数。
[0128] 本实施例中，阵风减缓控制器结构参数如图8所示，按给定的优化目标函数优化出 待求参数，结果如下所示。
[0129] X=[X。Xi X2 X3]t
[0130] X〇=[_0.5 -1.5 -0.7 0.9 -1.0 -0.1 -0.4 -0.1 0.6]t
[0131] Xi=[0.62 0.79 0.75 0.65 0.51] T
[0132] X2=[0.83 0.58 0.63 0.52 0.98] T
[0133] X3=[0.82 0.66 0.85 0.88 0.64] T
[0134] 本发明具有的优点和有益效果是：
[0135] a.在确定阵风减缓控制器结构时，根据无控飞机连续阵风响应功率谱密度曲线， 选择功率谱密度曲线的峰值频率作为滤波器设计参数，只保留对飞机阵风响应贡献较大的 模态振动信息进入控制回路，在确保阵风减缓效果的同时，有效降低发生不利气动伺服弹 性耦合的概率；
[0136] b.在进行控制器参数优化时，基于气动伺服弹性稳定裕度对体的适应度进行修 正，能够确保优化结果实现阵风减缓的同时，满足气动伺服弹性稳定性设计要求，同时采用 基于随机联赛选择的精英保留策略，能够有效提升优化运行效率。
[0137] 最后需要指出的是：以上实施例仅用以说明本发明的技术方案，而非对其限制。尽 管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解:其依然 可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替 换;而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精 神和范围。
【主权项】
1. 一种弹性飞翼布局飞机阵风减缓控制器参数优化方法，其特征在于，包括： 51、 选择作为阵风减缓控制器输入的若干个机体运动量； 52、 确定阵风减缓控制器的通道个数； 53、 计算无控飞机的任一机体运动量对连续阵风输入的功率谱密度函数，并在所述功 率谱密度函数对应的功率谱密度曲线中自起始点开始选择至少一个峰值； 54、 对所述阵风减缓控制器的任一通道，构建在前一步所选峰值的峰值频率下由滤波 器阻尼系数构成的滤波器表达式； 55、 由所述滤波器表达式与所述机体运动量输入信号的控制增益参数建立所述阵风控 制减缓控制器的传递函数矩阵； 56、 根据所述传递函数矩阵设置若干个优化目标，优化待求参数，所述待求参数包括所 述滤波器阻尼系数及所述控制增益参数。2. 如权利要求1所述的弹性飞翼布局飞机阵风减缓控制器参数优化方法，其特征在于， 所述机体运动量包括飞机的俯仰角速率《y、重心过载〃％以及翼尖过载气3. 如权利要求1所述的弹性飞翼布局飞机阵风减缓控制器参数优化方法，其特征在于， 对进入任意第m个控制通道的任意第k个机体运动量rk，选择nk个峰值频率，其中j = 1，…，nk，设计nk个滤波器，所述滤波器表达式为：其中，s为拉普拉斯算子，为滤波器阻尼系数。4. 如权利要求1所述的弹性飞翼布局飞机阵风减缓控制器参数优化方法，其特征在于， 在步骤S6中，所述优化目标包括使有控飞机翼根弯矩响应的均方根值最小、有控飞机俯仰 角速率响应的均方根值最小以及使最大舵偏角均方根值最小情况下的优化参数组合。5. 如权利要求1所述的弹性飞翼布局飞机阵风减缓控制器参数优化方法，其特征在于， 在步骤S6中，所述优化过程包括： T1、随机生成第一代染色体种群Po，染色体为优化变量X; T2、计算P〇中个体的优化目标的函数值，并根据个体的幅值和相位裕度设计罚函数，修 正个体适应度； T3、将P〇中个体按照其适应度进行非支配排序，计算同一非支配层级中个体拥挤度； T4、从P〇中优先选择非支配层级编号小的个体生成新种群Pi，所述Pi种群内个体数目为 P〇内个体数目的一半； T5、对种群丹中的个体进行交叉、变异操作，生成个体规模与Po相同的新种群Qo; T6、计算种群Qo中个体的目标函数值，采用步骤T2的方法计算种群Qo中个体的适应度； T7、种群Po和Qo合并组成新种群&，对&中的个体进行非支配排序和并计算个体拥挤度， 并优先选择Q:中的个体的非支配层级编号小的个体生成新种群Q2,所述出种群内个体数目 与Ρο内的个体数目相同； T8、重复步骤T3-T7，直至达到给定的最大进化代数，优化出待求参数。6. 如权利要求5所述的弹性飞翼布局飞机阵风减缓控制器参数优化方法，其特征在于， 在步骤Τ2中，根据个体Χλ对幅值和相位裕度的满足程度，设计罚函数$1胃和衣，公式包 括：其中，为若干个控制通道的幅值裕度的最小值，，.为若干个控制通道的相位裕度 的最小值； 之后，利用罚函数修正规范化处理后的优化目标函数，得到个体X△的适应度，修正优化 目标函数0〃i的公式为所述σ' i为经规范化处理后的优化目标函数，所述规范化处理包括将所述优化目标函数 取倒数。7. 如权利要求5所述的弹性飞翼布局飞机阵风减缓控制器参数优化方法，其特征在于， 在步骤T5中，所述交叉包括二进制交叉，概率为0.9，所述变异包括多项式变异。8. 如权利要求5所述的弹性飞翼布局飞机阵风减缓控制器参数优化方法，其特征在于， 所述种群规模为80，迭代次数为100次。
【文档编号】G05B13/04GK105929692SQ201610321700
【公开日】2016年9月7日
【申请日】2016年5月16日
【发明人】蒲利东, 张红波, 高怡宁, 程芳, 张庚庚
【申请人】中国航空工业集团公司西安飞机设计研究所