В работе предложен конструктивный алгоритм вычисления матриц исключения L и дублирования D для операции векторизации произведения P⊗P при P=P^T. Матрица L, получаемая в соответствии с алгоритмом, позволяет формировать из упомянутого произведения вектор, содержащий только уникальные элементы. Матрица D, в свою очередь, позволяет выполнять обратное преобразование. Предложена программная реализация процедуры вычисления матриц L и D. На основе отмеченных результатов предложена новая операция vecu(.), определенная над произведением P⊗P при P=P^T и описаны ее свойства. Показаны отличия и преимущества разработанной операции от известных операций vec(.) и vech(.) (vecd(.)) в случае их применения для векторизации произведения P⊗P при P=P^T. На примере параметризации алгебраического уравнения Риккати продемонстрирована эффективность операции vecu(.) для снижения перепараметризации задачи идентификации неизвестных параметров.