GNU Octave: как вычислить верхнюю и нижнюю огибающую набора сигналов?

Question

У меня есть набор сигналов (ряды Фурье, изменяющие параметры синуса и косинуса), и я хотел бы найти верхнюю и нижнюю огибающую этого набора.

Работоспособный пример ряда Фурье:

%interval length and number of points
L = 2*pi;
N = 60;
%data points
x = (linspace(0,L,N+1))';
s1 = 1*sin(x) + 2*cos(x) + 3*sin(2*x) + 4*cos(2*x);
s2 = 4*sin(x) + 1*cos(x) + 2*sin(2*x) + 3*cos(2*x);
s3 = 3*sin(x) + 4*cos(x) + 1*sin(2*x) + 2*cos(2*x);
#
clf;
plot(x, s1, 'k', x, s2, 'b',  x, s3, 'r' )
#

Как я могу вычислить нужные мне конверты?

Answer 1

Благодаря @Tasos Papastylianou я нашел решение.Достаточно расположить значения сигналов в матрице и взять минимальное и максимальное значения для каждой строки:

#To find envelope: arrange the function values in  matrix
#And take max and min by row
foo = [s1, s2, s3];
env_hi = max(foo,[], 2);
env_low = min(foo, [], 2);
#
clf;
hold on;
plot(x, s1, 'g', x, s2, 'b',  x, s3, 'r' )
plot(x, env_low, 'k^', x, env_hi, 'kv')
legend('s1', 's2', 's3', 'env-low', 'env-hi', 'location', 'north');
#