В кратце, желательно, чтобы в массиве дискретизации было целое количество периодов входного сигнала, иначе необходимо накладывать окно, с целью минимизации краевых эффектов. Далее сделать преобразование Гильберта, найти ортогональный сигнал, после чего найти огибающую: A[iTд]=sqrt(Uг[iТд]^2+Uф[[iТд]^2), где Uг[iТд] ортогональный сигнал, а Uф[[iТд]- ваш исходный.