Главная Нелинейные системы управления [0] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] [85] [86] [87] [88] [89] [90] [91] [92] [93] [94] [95] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] [107] [108] [109] [110] [111] [112] [113] [114] [115] [116] [117] [118] [119] [120] [121] [122] [123] [124] [125] [126] [127] [128] [129] [130] [131] [132] [133] [134] [ 135 ] [136] [137] [138] [139] [140] [141] [142] [143] [144] [145] [146] [147] [148] [149] [150] [151] [152] [153] [154] [155] [156] [157] [158] [159] [160] [161] [162] [163] [164] [165] [166] о наличии свойства экстремальности у У (х) и о необходимости соблюдения условия J (X*) < J (х), (11.4) где X - допустимая область изменения управляемого параметра. Отсутствие любой другой информации о свойствах функционала J {х} приводит к необходимости последовательно определять значения функционала качества внутри допустимой области изменения управляемого параметра. Если обозначить допустимый интервал изменения параметра X через А, а заданную точность в достижении экстремума - через 8 >0, то в результате сканирования определяют п значений функционала качества в точках Xi, х, .... х: Ji{x,). J,{x,), .... y„(n), (11.5) где п А /е + I. После перебора всех значений Ji (xt) выбирают максимальное или минимальное значение: J {x*)=min {JAXi)}. t-1,2,..., п. (11.6) Длительность процедуры поиска при сканировании в основном определяется задаваемой точностью е. Например, при допустимом интервале изменения А = 5 единиц и заданной точности 10 % необходимо сделать замеры в п = Л/е -1- 1 = 5/1 + 1 = 6 точках. (11.7) При заданной точности е = 1 % количество замеров существенно увеличивается: п = 5/0,1 -h 1 = 51. (11.8) Регулярность метода сканирования определяется заранее задаваемым порядком перебора значений. Чаще всего используют два способа при обходе точек: строчная развертка и спиральная развертка (соответственно рис. 11.3 и 11.4). Оба правила обхода обеспечивают просмотр всех допустимых точек без пропусков. Метод Гаусса-Зайделя. В методе Гаусса-Зайделя используют дополнительную информацию о виде функционала качества J (х), в частности предполагают, что J (х) является унимодальной функцией, т. е. функцией, имеющей один экстре- Ч 3 Z 2 3 Рис. 11.3 5 X/ 1 г J Рис. 11.4 мум. Условие унимодальности можно записать следующим образом (для поиска минимума): y(Xi)<y(x2) при J (Xi) > J (Хг) при Xi < Хг < Х„1„, (11.9) где Xmiij - положение минимума; Xj и Хг - произвольные положения относительно точки минимума. Условие унимодальности позволяет значительно сократить число просматриваемых точек по сравнению с полным перебором. В основу метода поиска положено исследование полной производной экстремизируемого функционала dJix)ldt [dJ{x)/dxt]xi, (11.10) (= 1 где X,- = 2- Он Ш {x)/dxi] г = 1, 2, .., п; /=1,2, /г; г=1 Пц- коэффициенты, характеризующие отклонение от экстремума. Отсюда dJ{x)/dt=!. ar,\dJ {x)ldxi\\dj{x)ldxi]. (11.И) В точке экстремума Xi - х, extr имеем [6J(x)/6xjpxtr = О-поэтому во всех точках, кроме Xextr. функция (П.П) должна удовлетворять условию монотонного приближения к экстремуму: d/(x)/d/>0 для максимума; (И 12) dJ {x)/dt <iO для минимума, f Пример 11.1. Пусть функция качества имеет вид, представленный па рис. 11.5. J [х,. х)=х1+х1-{-ахХг. (11.13) где а = const, о = 1,5 и поиск начинается из точки х, = 3, х - 3. Изменяем координату х, оставив х постоянной. Тогда функция качества ./(х,, 3)=Af-f9 + 4.5x,. (11.14) Минимум J (Xi, 3) находим, приравнивая нулю частную производную dJ (х,, 3)/dX], т.е. (Xi, 3)/axi = 2xi -f 4.5 = 0. (11.15) • Отсюда первый экстремум по х равен Xl = - 2,25 (11.16) и ему соответствует на рис. 11.5 точка с координатами Xi 2,25, Ха = 3. Теперь оставим координату х, на уровне - 2,25 и будем изменять !.55х,; I 2,55 = 0, X, =-1,27. j (11.19) / [(-2,2.5), XjJM5-fxi-3,4x,). (11.17) Находим минимум функции по х: dJ К-2,25), Xjj/dxa = 2xj -3,4 =0, Xj = 1,7. (11.18) Данному экстремуму соответствует точка И з на рис. 11.5 с координатами X, = - 2,25, Xj, = 1,7. Повторяем цикл вычислений лля координаты Xj, закрепляя найденное значение Xj = 1,7, т. е. J (X,, 1 ,7)=xf-b2.9-f2,55x,; dJ (X,. l,7)/ax, ==2A:,-f 2, В результате поисковых движений получаем ломаную линию, состоящую из взаимно перпендикулярных прямых, точки излома которой находятся в местах касания этих прямых с кривыми: J (Xl, Xj) s const. в методе Гаусса-Зайделя производится поочередное изменение координат Xl, Хг.....x„ и определяются частные экстремумы dJ (x)fdXi = О по каждой из координат, при этом все координаты, кроме выбранной, закрепляются. Взяв координату Xi, при Рис. п.,
[0] [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] [35] [36] [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] [49] [50] [51] [52] [53] [54] [55] [56] [57] [58] [59] [60] [61] [62] [63] [64] [65] [66] [67] [68] [69] [70] [71] [72] [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] [84] [85] [86] [87] [88] [89] [90] [91] [92] [93] [94] [95] [96] [97] [98] [99] [100] [101] [102] [103] [104] [105] [106] [107] [108] [109] [110] [111] [112] [113] [114] [115] [116] [117] [118] [119] [120] [121] [122] [123] [124] [125] [126] [127] [128] [129] [130] [131] [132] [133] [134] [ 135 ] [136] [137] [138] [139] [140] [141] [142] [143] [144] [145] [146] [147] [148] [149] [150] [151] [152] [153] [154] [155] [156] [157] [158] [159] [160] [161] [162] [163] [164] [165] [166] 0.0013 |