Главная страница Метод функций Ляпунова [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] ~ L Г {у,..... Задача В. При заданных векторах найти систему функций yi{t), у{t), тп, минимизирующих 8. Прежде всего следует заметить, что выбор значений функций Уа(0 в данный момент t не влияет на выбор значений этих функций в любой другой момент времени. Поэтому задача сводится к минимизации величины я т Теперь уже в качестве гильбертова пространства Н следует рассматривать пространство л-мерных векторов с, а в качестве подпространства будет фигурировать подпространство, порожденное векторами .....Сщ. Согласно (8.8) получаем систему для определения (f): 2 yk it) ic„ Cj) = (г (О, Cj), / = 1, 2, ..., /л. ft=i Здесь выражения (с, Cj), {r{t), cj) означают скалярные произведения соответствующих векторов в обычном смысле. Согласно (8.9) получим 8=Д(01Р=-Г Г(„ r{t))dt. Задача имеет простую геометрическую интерпретацию. Обозначим через плоскость, порожденную векторами с,.....с. В л-мерном пространстве Н задана кривая z = r{t)\ требуется найти кривую г= сУкЦ), которая имеет наименьшее среднеквадратичное отклонение Квадрат ошибки приближения в данном случае будет равен от заданной кривой. Решение задачи состоит в том, что каждую точку искомой кривой мы получаем как проекцию на плоскость соответствующей точки первой кривой z = r(t). Отсюда следует, что решение задачи не изменится, если в качестве порождающей системы векторов Cj, ..., с„ возьмем эквивалентную ей ортонормированную систему. В этом случае получим yj(t) = ir(t), cj), (8.13) = 5 (О dt - J]ir (t). с,Г dt. (8.14) о ;=1 Задача С. Найти оптимальную систему векторов Cft и функций {t) минимизирующих 8. Рассмотрим матрицу В, элементы которой имеют вид biu\ri{t)ru it)dt, где г,-(О - проекции вектора r{t). Так как матрица В симметрична, то ее собственные числа вещественны. Более того, собственные числа матрицы В неотрицательны. В самом деле, матрица В является матрицей знакоположительной квадратичной формы /(X) = 2 ft,.,x,x, = Ц 2 г, (О X, у. i,k=i о ft=l Отсюда согласно § 10 первой главы выводим, что наименьшее собственное число Х„, как минимум формы 1{х) на единичной сфере, должно быть неотрицательным. Расположим собственные числа матрицы В в порядке убывания Xi 2 ... Х„ 0. Каждому собственному числу Х соответствует по крайней мере один собственный вектор с- Известно ([5], стр. 20), что систему собственных векторов можно ортонормировать. Таким образом получим ортонормированную систему собственных векторов матрицы В. Теорема 8.1. Ортонормированная система собственных векторов Ci, матрицы В является оптимальной системой управляющих векторов, причем, если система управляющих функций у {t) выбрана согласно (8.13), то 5= 2 К. (8.15) В самом деле, если система векторов Cj, с„ орто-нормирована, то согласно формуле (8.14) имеем b=\r\t)dt- 2 [(с, nmt. о а=1 О Необходимо систему векторов подобрать таким образом, чтобы сумма J 2 \{с„ r{t))4t ft = I о была наибольшей. Рассмотрим первое слагаемое этой суммы: т h = I{cu) = \{Ci, r{t))4t = я 7 я = S (5 ri{t)rk{t)dtci,cu)= S -i/Cift- i, ft = l О i, = l Здесь - проекции вектора с,-. Задача сводится к тому, чтобы найти максимум квадратичной формы /(и) при условии i = l Согласно § 10 первой главы искомый максимум равен наибольшему собственному числу X, и достигается для собственного вектора Cj, отвечающего этому числу. Чтобы сделать максимальным следующее слагаемое, т k = \{c„ r{t))4t = I{c.J о при дополнительных условиях (Cj, c<i) = 0, с=1, нужно взять собственный вектор Cj, отвечающий следующему собственному числу. Согласно экстремальной теории квадратичных форм ([5]), максимум будет равен Х. Если \i = \, то в качестве с следует выбрать среди бесконечного множества собственных векторов кратного корня X, вектор. [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] 0.0142 |