Исследованы процессы распространения и трансформации упругой и электромагнитной волновой энергии в пьезоэлектрических анизотропных средах и их матричная формулировка. Построены и проанализированы аналитический вид матриц коэффициентов и фундаментальных решений системы ОДУ, описывающие связанные электромагнитные и упругие волны в пьезоэлектрических средах разной сингонии для совместной системы уравнений Максвелла и уравнения движения упругой анизотропной среды. Получены инвариантные соотношения, определяющие законы сохранения для волновых процессов в среде из структуры матрицанта. Получено аналитическое решение задачи отражения и преломления, построены соответствующие графики коэффициентов отражения и преломления в зависимости от угла падения. Методом матрицанта аналитически получена зависимость угла сноса между групповой и фазовой скоростью, которая подтверждена экспериментом.*