Окно вычисления

Здравствуйте.

Прошу помочь с задачкой.
Есть выборка значений функции (гармонический сигнал). Выборка имеет, допустим, бесконечный объем (1000 отсчетов). Частота дискретизации 500 Гц.
Необходимо реализовать "окно" вычислений, где будет происходить какая-либо операция (например, нахождение среднеарифметического значения). Объем "окна" - 3 выборки. Сдвиг "окна" должен происходить на 1 отсчет каждый раз. Т.е. посчитал на первых 3х выборках - запомнил, сдвинулся на 2, 3 и 4ую выборки - посчитал - запомнил, и т.д.

В качестве исходного сигнала можно взять синусоиду.
Все вычисления необходимо провести на диапазоне угла от 0 до пи/2, от 0 до 100 Гц.

Миниатюры

 

0

Себя хватило только на это: (расширил окно до 5ти выборок)
Код Matlab M123456789101112131415time=[2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40]*10^-3;A=5;for fi=0i/360i/2for f=40:60for k=3length(time)-2)w=2*pi*f; Data=A*sin(w*time+fi);X=time((k-2)k+2));Y=Data((k-2)k+2));diff_eq=(Y(1)-8*Y(2)+4*Y(4)-Y(5))/12;;DIF(k-2)=diff_eq;endendend

Добавлено через 51 секунду
Поправьте, пожалуйста, код. У меня никак.

Добавлено через 15 секунд
time=[2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 38 40]*10^-3;
A=5;
for fi=0i/360i/2
for f=40:60
for k=3length(time)-2)
w=2*pi*f; Data=A*sin(w*time+fi);
X=time((k-2)k+2));
Y=Data((k-2)k+2));
diff_eq=(Y(1)-8*Y(2)+4*Y(4)-Y(5))/12;;
DIF(k-2)=diff_eq;
end
end
end


Добавлено через 13 минут
Ребята, посоветуйте как реализовать в циклах изменение параметров частоты f, угла fi. Хорошо бы еще добавить изменение амплитуды А.
Как организовать правильный вывод из всех циклов?

Добавлено через 11 часов 46 минут
Забыли про пред сообщение.
Немного поясню...
Изначально имеется 2 вектора: вектор времени Т и вектор какой-то функции Ф, которая представляет собой синусоиду A*sin(2*pi*f*t+ф), в эти моменты времени. Файл-программа должна представлять "окно" вычислений. Т.е. вычисляется diff_eq на основе 1, 2 и 3го отсчета и полученное значение будет соответствовать центральному для них отсчету, т.е. 2му. Записали в массив, например. Далее "окно" сдвигается. Теперь работаем с 2, 3 и 4м отсчетами. Результат записываем. И так далее в цикле до конца исходного вектора функции Ф.
Все действия необходимо повторить для углов ф в диапазоне от 0 до 100 градусов, также для частоты от 0 до 100 Гц.
В итоге получается 3 цикла: цикл движения окна, цикл изменения угла и цикл изменения частоты.
Вопрос: как организовать вывод всех результатов? Для каждой точки, для каждого угла, для каждой частоты.

0