Базисный вейвлет

[Legato]

Базисный, или материнский вейвлет ψ образует по средствам растяжений и сдвигов семейство ψ = w((x - b))/a)
где a и b — параметры, определяющие соответственно масштаб и смещение функции ψ.
Процедура преобразования стартует с масштаба а=1 и продолжается при увеличивающихся значениях а, т.e. анализ начинается с высоких частот и проводится в сторону низких частот.
Для полного включения в обработку всех точек входного сигнала требуется задание начальных (и конечных) условий преобразования (определенных значений входного сигнала при t<0 и t>tmax на полуширину окна вейвлета).
Затем вейвлет масштаба а=1 сдвигается вправо на значение b и процедура повторяется. Получаем значение, соответствующее t=b в строке а=1 на частотно-временном плане. Процедура повторяется до тех пор, пока вейвлет не достигнет конца сигнала. Таким образом получаем строку точек на масштабно-временном плане для масштаба а=1.

Необходимо реализовать функцию Y = f(t)* w ((x - b))/a)

Взял в качестве базового вейвлета w и функции f(t) - синусойды, то есть Y= sin(t)* sin((x - b))/a). Попытался построить, ерунда какая то вышла...что то отдаленно похожее на правду,но не то.
Алгоритм был такой
for(b=0,1; b<10; b=b+0.1)
for(y=10; y<30; y=y+1)
for(x=0; x<20; x=x+1)
f = sin(x-y)/b*sin(x);
Ну только в лаб вью конечно реализован....если непонятно то,что здесь написано - вот номер аськи - 233749253, могу выслать программу.
В любом случае спасибо,что хоть прочитали)

[Eugen Graf]

Выложи прогу в форум, посмотрим что там за ерунда вышла, ну или хотя бы скриншоты блок диаграммы и результатов.

[Konstantin Sumenko]

Y = f(t)* w ((x - b))/a)

t и х как на сетке связаны? И, да, нужна программа.

[Legato]

Файл прикрепил...
Еще нужно контролировать число периодов,а я пока не совсем понимаю как...в файле "Работает" вроде задал функцию для периода,но тоже не уверен...чтобы он на выходе точки выдавал, так сказать.

[Legato]

Crowbar прошу прощения - это из теории выхвачено было,напутал,когда печатал...сейчас поясню..

Y = f(t)* w ((x - b))/a) Это :

Y = f(t)* w ((t - b))/a)

b - сдвиг, а - масштаб

[Legato]

На фото,что я прикрепил f(t) = sin(250t) + sin(300t)....в моем случае должно получиться примерно что то похожее, на вот эти вершины,если конечно правильно все написать...а у меня в программе синусойда кривая какая то или непонятно что...
На еще одном фото как раз формула для Непрерывное вейвлет-преобразование которую и нужно реализовать...