Добрый день,
нужно найти функцию

, такую, что выражение

лежит в плоскости, заданной вектором нормали

.
При этом

можно подбирать как угодно, лишь бы компоненты были положительными (

менять нельзя). Такое же требование - на компоненты

. Также желательно, но не обязательно, удовлетворить ограничению на величину

. Матрицы

на самом деле диагональные.
Возможный вариант, причем вполне логичный - выкинуть производную, положив скорость изменения

малой, но добавить ограничение вроде

.
Подоплека задачи - построение ориентации спутника с помощью магнитных катушек, создающих дипольный момент. В результате механический момент имеет вид

. Он лежит в плоскости, перпендикулярной вектору геомагнитной индукции, так что система оказывается неуправляемой. Чтобы справиться с этим, используется скользящее управление - фактически, ищется программная траектория, которую управление обеспечить может, и которая приводит в нужное положение.
Такие дела.