Скалярное произведение – простой, но чрезвычайно полезный математический инструмент. Он кодирует отношение между величинами и направлениями двух векторов в единственное скалярное значение. Его можно использовать для вычисления проекции, отражения, расчета тени и постановки освещения. Из этого руководства вы узнаете:
- Геометрический смысл скалярного произведения.
- Как спроецировать один вектор на другой.
- Как измерить размер объекта вдоль произвольной оси.
- Как отразить вектор относительно плоскости.
- Как создать эффект отскока мяча от наклонной поверхности.
Представьте два вектора a и b . Вектор характеризуется только направлением и величиной (длиной), поэтому не имеет значения, в каком месте плоскости он расположен. Допустим, оба вектора начинаются в одной точке:
Два вектора a и b , выходящие из одной точки плоскости
Проекцию одного вектора на другой можно представить как отбрасывание тени от первого вектора в направлении, перпендикулярном второму:
Проекция одного вектора на другой
Как вы знаете, операция вычисления произведения векторов записывается так:
Далее в статье мы будем использоваться запись a * b .
Если между векторами острый угол, то длина проекции будет положительной величиной, если больше – то отрицательной.
Реверсирование порядка векторов не изменяет результат скалярного произведения
Если b – единичный вектор, то величина проекции a на b – это просто произведение a * b .
На схеме изображен прямоугольный треугольник. Угол между векторами a и b равен θ .
Два вектора образуют прямоугольный треугольник
Для начала требуется рассчитать величину проекции вектора a на вектор b – это нижний катет в нарисованном нами треугольнике. Длину катета стороны можно найти, умножив длину гипотенузы треугольника на косинус прилежащего угла.
Проекция вектора a на вектор b – это катет прямоугольного треугольника
Итак, длина проекции равна произведению модуля вектора a на косинус угла θ . Скалярное произведение можно выразить следующим образом:
Формула вычисления скалярного произведения через косинус
Эта формула лишний раз подтверждает, что порядок умножения не важен – в результат входят беззнаковые длины обоих векторов. Если оба вектора – единичные, правая часть формулы упрощается до cos(θ) . А если угол равен 90° (векторы перпендикулярны), то их произведение равно 0 .
Если угол острый (меньше 90°), результат будет положительным, так как косинус такого угла больше 0. Аналогично для тупого угла получится отрицательный результат. Таким образом, знак скалярного произведения дает нам некоторое представление о направлениях векторов.
Функция косинуса монотонно убывает на промежутке от 0 до 180° (от 1 до -1 ). Следовательно, чем ближе направления двух векторов, тем больше их скалярное произведение и наоборот.
- Направления совпадают, угол θ равен 0°, произведение равно |a| * |b| .
- Направления противоположны, угол θ равен 180°, произведение – -1* |a| * |b| .
Если наши векторы расположены в 3D-пространстве и имеют по три координаты каждый, не совсем понятно, где тот угол, косинус которого нужно вычислить. К счастью, существует другой способ расчета скалярного произведения – без всякой тригонометрии! Для начала нужно разложить каждый вектор на компоненты:
Разложение векторов на компоненты
Вычисление скалярного произведения векторов через сумму произведений соответствующих компонент
Намного проще и без всяких косинусов! В Unity есть встроенный метод Vector3.Dot для вычисления скалярного произведения двух векторов:
Его реализация выглядит следующим образом:
Нам известно, как найти длину вектора по его координатам:
Но ее можно выразить и через скалярное произведение вектора на себя:
Вернемся к формуле a * b=|a| * |b| * cos θ . При известных длинах векторов мы можем вычислить угол между ними с помощью функции арккосинуса:
Вычисление величины угла между векторами
Если оба вектора являются единичными, мы можем упростить формулы:
Упрощенные формулы для единичных векторов
Проекция вектора
Теперь, когда нам известно геометрическое значение скалярного произведения векторов (произведение длины со знаком первого проецируемого вектора и длины второго вектора), мы можем перейти к практическому применению этого знания. Например, спроецируем один вектор на другой.
Пусть вектор с=projectba – это проекция вектора a на вектор b .
Вектор c – проекция вектора a на вектор b
Возьмем единичный вектор в направлении вектора b . Он будет равен b / |b| . Если мы возьмем величину проекции a на b со знаком и умножим на этот единичный вектор, то получим вектор c . Cкалярное произведение a * b – это результат умножения длины b на длину проекции a на b . Отсюда получаем, что длину c со знаком можно найти, разделив скалярное произведение a * b на длину b :
Вычисление длины проекции одного вектора на другой
Умножив полученное значение на единичный вектор b / |b| , получаем формулу для нахождения проекции вектора:
Вычисление проекции одного вектора на другой
Вспомним теперь, что квадрат длины вектора равен его скалярному произведению на самого себя, и перепишем формулу:
Вычисление проекции одного вектора на другой
Если b – единичный вектор, то можно упростить еще больше:
Вычисление проекции вектора на единичный вектор
В Unity для вычисления проекции одного вектора на другой есть специальная функция Vector3.Project :
Вот так выглядит ее реализация:
Следует остерегаться возможного вырожденного случая, когда вектор, на который происходит проекция, – нулевой или имеет малую величину. При этом произойдет численный «взрыв» из-за деления на 0 или близкое к нему значение. Один из способов решить проблему – заранее вычислять величину вектора и при необходимости использовать резервный вариант (единичный вектор).
Упраженение #1. Линейка
В качестве упражнения для закрепления знаний о векторной проекции сделаем линейку, которая измеряет размер объекта вдоль произвольной оси.
Линейка характеризуется базовой позицией (точка) и осью (единичный вектор):
Как спроецировать какую-либо точку ( Point ) на линейку? Прежде всего, найдем относительный вектор от базовой позиции линейки ( Base ) до этой точки. Затем спроецируем его на ось линейки ( Axis ). Проекция точки ( Projection ) – это базовое положение линейки, смещенное на проекцию относительного вектора.
Проекция точки на линейку
Промежуточное значение relativeDot в коде отражает, насколько далеко проекция точки находится от базового положения линейки – в направлении оси линейки, если она положительная, или в противоположном направлении, если отрицательная.
Чтобы найти размер объекта вдоль оси линейки, нужно провести такие измерения для каждой вершины меша (Mesh) и найти минимальное и максимальное значение. Ответ будет равен разнице между ними.
Отражение вектора
Еще одно практическое применение скалярного произведения – отражение вектора относительно плоскости. Рассмотрим вектор v и плоскость с нормальным вектором (перпендикуляром) n .
Отражение вектора v от плоскости
Мы можем разложить отражаемый вектор на параллельную и перпендикулярную к плоскости составляющие:
Разложение отражаемого вектора на составляющие
Сам вектор является суммой параллельной и перпендикулярной составляющих:
Разложение отражаемого вектора на составляющие
Перпендикулярный компонент, как вы уже догадались – это проекция вектор на нормаль плоскости. Параллельный компонент можно получить вычитанием этой проекции из самого вектора:
Получение параллельного плоскости компонента вектора
Теперь, чтобы получить искомый отраженный вектор, нужно перевернуть перпендикулярную составляющую и прибавить ее к параллельной:
Получение отраженного от плоскости вектора
Параллельный компонент можно заменить разностью самого вектора и перпендикулярной составляющей, тогда получим следующий вид формулы:
То есть разность самого вектора и его удвоенной проекции на нормаль плоскости.
В Unity, конечно же, есть встроенная функция для расчета отраженного вектора – Vector3.Reflect :
Так выглядит реализация согласно первой выведенной нами формуле:
Упражнение #2. Отскок мяча от наклонной плоскости
Теперь, когда мы знаем все об отражении вектора от плоскости, можно переходить на следующий уровень сложности. Попробуем сымитировать движение мяча, отскакивающего от склона.
Для моделирования траектории движения шара под действием силы тяжести мы будем использовать метод Эйлера.
Чтобы определить, когда мяч ударяется о склон, нужно как-то определить момент, когда мяч проникает в плоскость.
Сфера может быть определена центром ( C ) и радиусом ( R ). Плоскость определяется нормальным вектором ( n ) и точкой на плоскости ( P ). Вектор от P до С обозначим u .
Сфера и плоскость
Если сфера НЕ проникает в плоскость, перпендикулярный плоскости компонент вектора u , должен иметь то же направление, что и вектор n , а также длину не менее R .
Перпендикулярный плоскости компонент вектора
Другими словами, сфера не проникает в плоскость, если скалярное произведение векторов u и n больше R . В противном случае величина проникновения составляет R – u * n , и положение сферы нужно исправить.
Чтобы это сделать, можно просто переместить сферу в направлении нормали плоскости n на величину проникновения. Это лишь приближенное решение, которое не является физически правильным, но для упражнения оно хорошо подходит.
Добавим логику для коррекции позиции:
Нам также нужно отразить скорость, с которой движется сфера, относительно наклона, чтобы она правильно отскакивала.
Эта анимация демонстрирует идеальное отражение и выглядит неестественной. Мы ожидаем, что с каждым отскоком скорость мяча будет уменьшаться.
Это поведение обычно моделируется значением реституции (восстановления) между двумя сталкивающимися объектами. При 100% реституции мяч идеально отскакивает от плоскости. При 50% – величина перпендикулярной к плоскости составляющей скорости мяча будет уменьшена вдвое.
Величина реституции – это отношение величин перпендикулярного к плоскости компонента скорости мяча до и после отскока.
Вот пересмотренный с учетом коэффициента восстановления скорости вариант функции отражения:
Вот так выглядит обновленная функция SphereVsPlane :
Логика корректировки позиции заменяется логикой полноценного отскока:
Теперь мы можем устанавливать разные коэффициенты реституции для разных шариков:
Заключение
Мы ответили на все вопросы, заданные в начале этого руководства.
- Скалярное произведение двух векторов – это произведение проекции первого вектора на второй (с учетом знака) и модуля второго вектора.
- Существует две формулы вычисления скалярного произведения: через косинус угла и через компоненты векторов.
- Скалярное произведение имеет множество полезных практических применений. Например, оно позволяет рассчитать проекцию вектора на другой вектор.
- С помощью скалярного произведения можно найти отраженный от плоскости вектор. На этой основе строятся различные физические модели, например, имитация отскока шарика от плоскости.
Расскажите в комментариях – разрабатываете ли вы сейчас какую-нибудь игру и если да, то о чем она. Каких уроков по Unity вам не хватает?