И не найдете

А нафиг Вам спектр синуса?
"Прога" не дла этого
Надо "дискретный стохастический ряд, имеющий зашумлённые гармонические компоненты"
С "низкочастотным трендом и скачками в потоке данных"

Шо есть трЕнд, кстати?

Мне ? Ни на фига



"Плавное" изменение среднего в частности по ср. с общим или локально-генеральная тенденция.
На самом деле - еще та ересь.

Да все та же сумма низкочастотных гармоник.

И кто-то выдает это за планы "Партии и Правительства".


Формируем сигнал, как и задумано авторами:

for i:= 1 to 1000 do..
x := step*i;
rnd := 0.5*(2*Random-1);
if i mod 133 = 0 then shift := shift + 0.3;
if i mod 33 = 0 then pulse := 3*(2*Random-1);
y := Sin(x) + rnd + pulse + shift + 0.3*Sin(3*x) + 0.2*Cos(5*x) + 0.3*Cos(20*x);

Лажа полная:




Я даже боюсь представить, что покажет spectrand, если скормлю ему стохастический процесс с заданными корреляционными свойствами , сформированный из белого шума методом формирующего фильтра

P.S.
Это "чистый" синус-композиция
y := Sin(x) + 0.3*Sin(3*x) + 0.2*Cos(5*x) + 0.3*Cos(20*x);

Да ну... Прикольно так все. Разноцветно...

Надо поддержать авторов, зря что-ли они на многих ворумах тычат пальчиком в свой уникальный алгоритм

TSerg,
"Не поленился, сгенерил чистый синус в один период и поимел вот такой спектр"

Да, Вы получили честный спектр - проверьте другим спектроанализатором.
Если не понимаете, что искать спектр по одному периоду - смешно.
Вы бы ещё пол-периода прогнали, или четверть.
Чтобы получить выраженный пик, нужно несколько периодов.
Ваши дальнейшие примеры страдают тем же. Подучитесь, что ли.
Неужели Вы думаете, что программа не тестировалась?

Вы господа совершаете ошибку пытаясь впарить свое поделие на этом форуме. Тут знаете ли люди разбираются в цифровой обработке и ваша программа вызывает в лучшем случае улыбку. Идите премьер-министру это показывайте он любит всякие нанотехнологии. И да еще если уж вы и говорите что изобрели чудо-юдо алгоритм, то дайте ссылку где про него почитать.

ЗЫ спектр одного периода синусоиды есть одна палка а не как у вас там нарисовано.

Сообщение отредактировал bahurin - Oct 30 2012, 16:15

bahurin,
я ничего не впариваю, а просто предлагаю - потому что есть чем поделиться.
Так что смените тон, пожалуйста.
Ссылка на статью, где описан алгоритм, дана в Описании.

"спектр одного периода синусоиды есть одна палка а не как у вас там нарисовано"

Это спектр бесконечной синусоиды есть одна палка. Я смотрю, Вы тоже крутейший
специалист. По понтам.

2 палки. В смысле, 2 дельта-функции.

ps
А апликуха ваша, увы, абсолютно unusable...

ТС на заметку:

Генерация сигнала:
a := 0.1;
b := 0.5;
c := 0.1;
d := 0.05;
alpha := 2*cPi*t; // t =0..1.0
Value := a + b*Sin(alpha) + c*Sin(3*alpha) + d*Sin(17*alpha);
P.S.
Число точек: 1030

Для анализа спектра имеем один период основной частоты и пару гармоник ( 3 и 17 ):



После честного и всем понятного FFT на 1024 точки ( без какой-либо оконной обработки ):




P.S.
Что-то примерно в этом роде должна показывать Ваша прога, а иначе - см. посты выше.

Для ТС

В статье приведен алгоритм. Обоснования законности (корректности, работоспособности и другие синонимы) таких преобразований там нет, где можно уведеть? Каковы ограничения на применимость алгоритма?

FFT подразумевает периодичность сигнала. Может быть растекание спектра, как показал TSerg, если укладывается не целое число периодов на обрабатываемой выборке.

ТС оперирует прикладной темой связанной с движением спутников, планет, галактик.. etc.
Наблюдаемые параметры движения относятся к непрерывным процессам., поэтому ТС-ру приводить контрпримеры связанные с ограниченной длительностью наблюдаемого сигнала по меньшей мере некорректно.
Если бы он изучал спектры одиночных импульсных сигналов, то в этом случае поправки на спектральное описание были бы уместны, да и то, чисто теоретически.

Открываем матлаб вставляем и запускаем скрипт:



смотрим график:



Очевидно что на графике только одна палка S(1) все остальные равны нулю (надеюсь не надо объяснять что такое вторая палка).

После этого извиняемся перестаем дерзить и идем допиливать свою программу до рабочего состояния или перестаем называть
свою программу анализатором спектра поскольку к спектральному анализу она в своем текущем виде не имеет никакого отношения.