вот так будет работать?
Что такое 2*pi*f*t?
Получатся две боковые по обе стороны от частоты sin/cos

http://www.dsplib.ru/content/hilbertquad/hilbertquad.html#r0

Увеличить

rst59 Ну собсно sea уже ответил. :) Придётся с полосой определится заранее, или только бла-бла с микрофона, или если в планах музЫку или чего пошире, то сразу лучше прикидывать.

спасибо за ответы
и ещё вопрос
где будет располагаться спектр сверху или снизу? как перекинуть на другую сторону? cos и sin местами поменять?

sea
cos(pi/2*n)
т.е. эта функция принимает только значения -1, 0 и 1?

rst59 как перекинуть на другую сторону? cos и sin местами поменять? Да.
т.е. эта функция принимает только значения -1, 0 и 1? Блин наваждение какое-то с частными случаями. :-) В общем случае нет, в каких-то частных да там будут константы, парит чекить те что вы напсали или чего другое, и в каких случаях именно, это имхо несложно прикинуть.

То есть мы переносим сигнал на промежуточную частоту, отфильтовываем одну боковую и вроде как переносим обратно в ноли кажный канал (0 и 90 град).

Есть еще похожий вариант. Например при сигнале 3 КГц перемножаем на 1.5 ГКц так, чтобы сигнал был вокруг 0 Гц. И так далее, прямо как по книге ППП, только в цифровом виде.
http://www.microtelecom.it/ssbdex/ssbdex-e.htm

Еще можно попробовать сделать Фазовый фильтр(All pass filter) но мне это показалось сложнее для понимания. Много гуглил по поводу Hilbert filter design нашел немного информации. По русски еще меньше. И как понял КИХ фильтр даст нелинейную АЧХ, а БИХ я чего-то не смог синтезировать.

Вот тут делают, как я понял, на базе цепочки из нескольких БИХ фильтров:

http://www.jurnal.org/articles/2008/izmer10.html

Вот еще несколько файлов полистать
10.1.1.26.1054.pdf
10.1.1.41.5876.pdf
10.1.1.139.1109.pdf
Еще
Че то про SDR
Эхх, устал грузить. Остальное тут http://rghost.ru/27484101

Вот еще ссылки
http://www.dsprelated.com/showmessage/88412/1.php
http://www.katjaas.nl/hilbert/hilbert.html
http://www.student.oulu.fi/~oniemita/dsp/hilbert/

Приятного прочтения! ;)

Придётся с полосой определится заранее, или только бла-бла с микрофона, или если в планах музЫку или чего пошире, то сразу лучше прикидывать.

А если с микроконтроллера генерить, например CW, BPSK и прочее, на базе несущих, то их сразу же можно генерировать в 0 и 90 и ничего фильтровать не надо.

Mesh, sea
спасибо! попробую реализовать

rst59 Есть еще похожий вариант. Например при сигнале 3 КГц перемножаем на 1.5 ГКц так, чтобы сигнал был вокруг 0 Гц. И так далее, прямо как по книге ППП, только в цифровом виде. Вот если взять дискрет для такого примеру 6000 килогерц это и будет один из частных случаев. Имхо. :) 6000 килогерцев это я конечно маху дал, пардоне.

Таки да ) Получается умножение можно заменить за изменение знака, потому, что умножать надо будет один канал на 0, 1, 0, -1, второй, соответственно на 1, 0, -1, 0. Но тут все равно, что использовать ключевой смеситель, т.е. получается перемножение на меандр. В данном случае получается, что после ЦАП с дискретизацией в 6 КГц нужно поствить соотв. ФНЧ в 6 КГц, иначе сигнал много раз повторится по частоте (гармоники). Физика - штука злая. Если же при 48 КГц умножать на меандр 1.5 КГц, то гармоники вылезут уже в цифровом виде. Может я конечно и ошибаюсь, надо глядеть на условия реализации, но мне кажется, что лучше сделать 48 КГц дискретизацию, умножать на нечто похожее на синус с косинусом и тогда хотя бы в цифре гармоник не будет до 24 КГц, дальше на тот же ФНЧ КГц, и он подавит то, что будет выше 24 КГц намного лучше, чем то, что останется выше 6 КГц при дискретизации в 6 КГц. Но надо смотреть вычислительную мощь, применяемые ЦАП/АЦП. Все ИМХО, конечно.

sea Гм, как бы оно б мне не очень надо, разберётесь конечно, но вроде как то после ЦАП с дискретизацией в 6 КГц нужно поствить соотв. ФНЧ в 6 КГц, что-то тут не то? :) Что может отфильтровать ФНЧ со срезом 6 килогерц при дискрете 6 килогерц? Как бы особо дискутировать не очень хоца, тема для меня не очень волнующая, но вроде косяк какой в прикидах?

Ну может быть... На два поделить... Просто после ЦАП надо ФНЧ. А если ЦАП на 6, а ФНЧ на 3 зачем вообще цифровая обработка? :)

спасибо! попробую реализовать
Интересно на чем

ps. Я в БПФ не специалист, но думаю, что его тоже можно использовать для поворота фазы на 90 градусов. Там же вроде как раз и получаются амплитуды и фазы. Но это я так... Слышал звон.

sea А, ЦАП... ну да, чёт не просёк сходу, думал о чём-то о своём. :) Да юзают и БПФ конешно. Тут от платформы зависит изделия ну и цели какие нужно достигнуть.

вот заюзал тут фильтр TIirLPF из исходников Eralgen и Mesh
но как-то он не по нужной частоте срезает, а выше процентов на 40-50
кто вкурсе в чём дело?

вот заюзал тут фильтр TIirLPF из исходников Eralgen и Mesh
А можно чуть подробнее? Просто вариантов LPF в виде IIR вагон и маленькая тележка.

действительно, забыл упомянуть
вот отсюда брал Конвертер из IQ в wav исходники на Delphi и прога, добавил пропуск одного отсчёта, и игнорирование не стерео файла, ещё добавил реверс, или инверсию, что одно и тоже, конечного файла.

собственно хочу сделать конвертер wav в IQ :)

rst59 Я так думаю, там просто вы попарились с частотой дискрета судя по тому что 40-50 эт как раз в два раза. Вся софтина и все классы это eralgen, я тока там чуток совего дописал. Дискрет меняется при перегонах туда суда.

собственно хочу сделать конвертер wav в IQ :)
Ага, а потом нарежете на болвалку и сделаете безкомпьютерный SDR маяк.

вот файл до и после фильтра
http://db.tt/JdtjEmNI
http://dl.dropbox.com/u/9597146/bfr.jpg
http://dl.dropbox.com/u/9597146/aft.jpg

фильтр задаю так
TIirLPF.Create(12000, 48000, 0.8, 12, 0.9)

нарежете на болвалку и сделаете безкомпьютерный SDR маяк
можно и так :)
но желательно риалтайм обработка

rst59 Да я не могу проверить, у меня дельфи давно уже нет. Вот думкаю может скачать да поставить? А другой какой фильтр, там вроде их куча классов была, да и порядок с ходу 12 может не надо, так попробовать-то можно и на меньшем потом увеличить. Блин я уж не помню чего там и как.

странно... с другими фильтрами (TIIRFilter, TIIRFilter1) тоже не получилось отсечь частоты выше 12кГц

rst59 Сгенерите шум во всей полосе в том ж Куле, тока с амплитудой там не переборщите, прогоните через фильтр и будет видно в реале чего там где, может оно всё и нормал и что б получить то что вы хотите нужно 120 порядку фильтры или ещё круче. Срез вобще вроде по уровню 3 дб задаётся, там ни как не будет что сразу после 12 полный мрак типо всё исчезло.

пробовал и больший порядок - бестолку

rst59 Положите листинг кода где инизилизируете и используете фильтр, так не бывает, где-то в консерватории косяк.

Я вот на дельфи делаю белый шум рандомом. Потом фильтр. Потом все это в файл dat. В SV старом открываю и смотрю что получается. Дискретизацию сейчас считаю и ставлю на 48 КГц. Фильтры БИХ. Порядок 6-8. Тогда получается нормально. Если КИХ, то там и 120 выглядят как-то... Фильтры считаю в Filter Solutions там есть сразу же вывод текста в с++, чтобы понять алгоритм хватает.

По поводу подавления заметил такую интересную весчь, что БИХ ФНЧ чем выше частота среза, тем меньше подавление, а ФВЧ практически ничего не меняют (странно). Полосовые работают.

sea
можете выложить выши делфяные исходники?

Mesh
function Norm(r: real): SmallInt;
begin
if r > 65535 then r := 65535;
if r < -65536 then r := -65536;
Result := Round(r);
end;

LPF1 := TIirLPF.Create(12000, 48000, 0.8, 12, 0.9);
SetLength(waveout, Length(wavein));
for i := Low(wavein) to High(wavein) do
waveout[i] := Norm(LPF1.Exec(wavein[i]));

далее пишем в файл заголок wav и массив waveout

вот пример обработки с белым шумом
http://dl.dropbox.com/u/9597146/noise.rar

function Norm(r: real): SmallInt;
Нельзя так делать, Smallint: -32768..32767
Сначала надо вычислять максимальное значение массива, потом приводить к Smallint.

з.ы.
Совсем прааильно будет, если исходный массив RealArray привести к SingleArray, чтобы максимальные значения не вылетали за -+1.0, потом умножить на 32767, поучится нормальная 16-ти битная Data.
Все это делать надо после фильтра, потому, как он нормально работает с плавающей точкой.
Во избежание непонятных результатов тип Real, без особой на то необходимости использовать не надо, Single вполне достаточно, если мало, везде ставим Double.

rst59 Ну да Programmist дело говорит. Так-то косяков особых в консерватории не видать, но с типами промашка конкретная. Запросто может быть причиной чего угодно. Уж хотя б фаншен в рамки вогнать для приличия что ли, где так

function Norm(r: real): SmallInt;
const
max = 32760;
min = -32760;
begin
if r > max then r := max;
if r < min then r := min;
Result := Round(r);
end;

з.ы. Ессно wavein должен быть smallint array если читаете из 16 битных файлов, иначе тож бардак будет.

Mesh
function Norm(r: real): SmallInt
Не совсем в курсе, что тут хотят, но, имхо эта функция все испортит. Она только ограничит диапазон до нормальных значений. Для правильного результата, нужно приводить весь массив к SmallInt.

з.ы.
Ессно wavein должен быть smallint array если читаете из 16 битных файлов, иначе тож бардак будет.
Ну да, цикл нормализации нужен. Не знаю пока стандартов на 32-х битную дату, но похоже, что все стараются сейчас не выходить за единицу.