Расследование и предупреждение техногенных катастроф. Научный детектив
Шрифт:
Если при изменении исходных данных расчета (например — диаметра круглой колонны) на 1% результат расчета (например — критическая нагрузка колонны) изменится в два раза, то такой расчет, разумеется, никакого практического смысла не имеет. Здание, построенное по такому нелепому расчету, разумеется, обязательно рухнет. Корректность решений для практики важна, очень важна. Поэтому корректность всегда тщательно проверяют. Но в 1987 году в СПбГУ было открыто, что существуют особые объекты, в математических моделях которых корректность изменяется при эквивалентных преобразованиях. Для таких особых объектов традиционные методы проверки корректности не достоверны,
Для того чтобы катастроф было меньше и наша жизнь стала безопаснее, надо уметь еще на стадии расчета и проектирования найти и выделить «особые» объекты. Об интереснейшей истории открытия особых объектов и разработки методов их распознавания и выделения мы далее расскажем, а пока приведем совсем простой числовой пример, который сразу прояснит суть дела. Никаких знаний, кроме школьной алгебры, для понимания примера не нужно.
Рассмотрим систему двух алгебраических уравнений:
(22 + 2)х = 2у (1)
(2+) = y (2)
с двумя переменными х и у и параметром .
Поскольку уравнения (1) и (2) однородны, то они, разумеется, имеют нулевое решение х = у = 0. Однако при некоторых значениях параметра система, состоящая из уравнений (1) и (2), имеет не нулевые решения. Значения параметра, при которых система однородных уравнений имеет не нулевые решения, называют собственными значениями (или собственными числами). Для системы (1) и (2) единственным собственным значением является = 1. Действительно, при подстановке в (1) и (2) значения = 1, система (1)-(2) переходит в систему:
4x = 2y (3)
2x = у (4)
и имеет, например, решения: х = 1; у = 2 или х = 2; у = 4 и многие другие. А вот при
Заметим сразу, что задача вычисления собственных значений (разумеется, для систем гораздо более сложных, чем простейшая система (1) и (2)) имеет очень важное значение в технике. От величин собственных значений зависит устойчивость того или иного технического объекта, здания, сооружения, зависят частоты его колебаний и т. п.
Поэтому задаче вычисления собственных значений, различным методам их расчета, посвящены целые книги (например книга: Х. Д. Икрамов. Несимметричная проблема собственных значений, издательство «Наука», 1991 г., 240 страниц или: Уилкинсон Д. Х. Алгебраическая проблема собственных значений, издательство «Наука», 1970 г., 564 страницы и многие другие). И все методы используют эквивалентные преобразования. А то, что может произойти при эквивалентных преобразованиях, мы покажем на простейшем примере системы (1)-(2).
Вместо громоздкого перебора всех возможных значений , собственное значение легко найти эквивалентным преобразованием — подстановкой. Подставив значение переменной у из уравнений (2) в уравнение (1), мы получим:
(22 + 2)х = 2(2 + ), (5)
Приведя подобные члены, получим:
(2 - 2) = 0. (6)
Из
Таким образом, эквивалентные преобразования позволили легко и просто найти (как и следовало ожидать) правильную величину собственных значений. Здесь все верно.
А теперь посмотрим, что получается при проверке корректности, при проверке зависимости собственных значений от малых изменений коэффициентов. После эквивалентных преобразований мы имеем дело с уравнением (6). В него входят два одинаковых коэффициента: двойка при и двойка как свободный член. Пусть свободный член изменился на 1% и стал равен 1,98. Тогда и собственное значение изменится на 1% и станет равным 1,01. То же самое произойдет, если на 1% изменится коэффициент при . Общий вывод: малым изменениям коэффициентов соответствуют малые изменения решения. Решение корректно.
А теперь (внимание!) исследуем корректность решения той же задачи до эквивалентных преобразований. Обратимся к исходным уравнениям (1) и (2) и посмотрим, что будет, если, например, коэффициент при 2 в уравнении (1) изменится на 1% и станет равным 1,98, а система (1)-(2) перейдет в систему:
(1,982 + 2)х = 2у (7)
(2+) = (8)
Отыскивая собственные значения для системы (7)-(8), мы убедимся, что их два: 1 = 0,99019; 2 =-100,99019 (с точностью до пятого знака после запятой). Можно проверить, что у системы (7)-(8) и при 1 =0,99019, и при 2 =—100,99019 действительно будут не нулевые решения х и у. Таким образом, для системы (1)-(2) решение задачи отыскания собственных значений — не корректно: уже при изменении одного из коэффициентов на 1% решение меняется коренным образом — вместо одного собственного значения появляется два, причем второе собственное значение резко отличается от первого (и даже имеет другой знак). Можно проверить, что если в системе (1)-(2) коэффициент при в уравнении (1) изменить не на 1%, а на 0,1% или даже на 0,01%, то все равно вместо одного собственного значения появятся два. Для системы (1)-(2) решение задачи о собственных значениях на самом деле не корректно — но мы не увидим этого, если будем исследовать корректность решения после эквивалентных преобразований системы, как это рекомендуют традиционные методы. Система (1)-(2) является особой системой — системой, у которой корректность решения изменяется после эквивалентных преобразований.
Для простой системы (1)-(2) все ясно и понятно: в уравнениях (1) и (2) коэффициенты при 2 после эквивалентного преобразования взаимно сокращаются и исчезают, хотя именно их малые изменения в исходной системе приводят к большим изменениям собственных значений. В более сложных системах все сложнее, распутать причины и следствия очень не просто, но главное заключается в другом: даже на примере очень простой системы (1)-(2) мы показали, что эквивалентные преобразования могут изменять многие важные свойства математических моделей. Могут изменять корректность решения, могут изменять запасы устойчивости и т. д. Впервые все это было опубликовано в 1987 году, в книге [1] (номер в квадратных скобках соответствует номеру в списке литературы в конце брошюры), а более подробно — в книгах [2], [3].