Предложения по функциональности программы для SDR

[Relayer]

слона то ты и не приметил ссылки ведут вот сюда:
http://dttsp.sourceforge.net
а там написано вот чего:

DttSP is an open source project started by Dr.Frank Brickle and Dr. Robert McGwier of the DTTS Microwave Society to provide code to be used in various DSP projects with an emphasis on Software Defined and Cognitive Radio.
DttSP is the digital signal processing powering the SDR-1000: Flex Radio and we are contributors active contributors to the open source project associated with it.

[NOP]

Нда-с... по идеи остаеться только собрать в кучу этот SDR/DSP процессер, в нормальный си-шный билд, для винды (походу в линуксе с этой прогой, давно полный порядок). И научиться им управлять. А то вон оболочка есть даже на java
http://microsat.homelinux.org/dttsp/softrock.jar
для GTK+ гдет видел, лихо, однако... :~\

[Serg]

To Serg
Попробуйте у товарисча исходники его работ выпросить, хотелось бы взглянуть как у него сделана фильтрация.

На сколько я помню, там был использован модуль ffts.pas готовый и найденный в инете, а фильтра он сам делал полосовые, складывая ФНЧ и ФВЧ. Насчет поделиться чем-то - не гарантирую, но спрошу.

[RA9YTJ]

Вопрс спецам:
есть комплекс. число. типа

typedef struct
{
double real;
double imag;

}complex;

и число double Alpha;
это угол поворота фазы, может -360 до +360,
какой алгоритм поворота?
Ни как не могу сделать правильно, помогите....

[Relayer]

какой алгоритм поворота?
Ни как не могу сделать правильно, помогите....

вики говорит что надо комплексно домножить на cos(φ)+i*sin(φ)
http://en.wikipedia.org/wiki/Rotation_(mathematics)

[DRUID 3]

Вопрс спецам:
есть комплекс. число. типа

typedef struct
{
double real;
double imag;
}complex;

и число double Alpha;
это угол поворота фазы, может -360 до +360,
какой алгоритм поворота?
Ни как не могу сделать правильно, помогите....

))) домножить на нормированное (sqrt(I*I+Q*Q) == 1) комплексное число, причем искомый угол определяется fi = atan2(Q, I). Другими словами масштабировать I и Q коэффициенты, что думаю, очевидно...

Вобщем вот ссылка http://mat.1september.ru/2001/10/no10_1.htm, почитать в низу страницы и немного подумать...

Но еще вопрос какие это числа (временные отсчеты или комплексный спектр) . "Подкрутить" I Q каналы для корректировки, вполне реально, но избавиться от преобразования Гильберта или DFT для получения SSB абсолютно нереально.

[Sergey_gh]

Угол +- 180

Верно написал Relayer.

complex z, result;

result.real = z.real * cos(φ) - z.imag * sin(φ);
result.imag = z.real * sin(φ) + z.imag * cos(φ);

DRUID 3 не искал легких путей и показал что теорема Пифагора и тут верна .

[RA9YTJ]

Так, большая проблема:
Вот упрощенный алгоритм:
1. копируем вх. буф. отчетов АЦП в массив БПФ, в real, imag обнуляем.
2 Прямое БПФ.
3. По формулам сверху поворачиваем фазу.(например на 90гр, что эквивалентно FFTnew(i).real=-FFTold(i).imag FFTnew(i).imag=FFTol d(i).real
4. Обратное БПФ.
5. Копируем из FFT.real в выходной буф. АЦП.
так вот во выходном буф. 0!!! если во вых. буф. копировать корень суммы квадратов FFT.real FFTold(i).imag(Пифаг ор) то получим поворот 180гр!!!!!
Я построил графики на всех этапах.
Что получилось, после поворота на 180гр, мы на вых в real получаем поворот противофазный(так и должно быть)!, при 90 и 270- 0!, при 45 0.5 амплитуды.
И вообще если строить графики real и imag, то он друг относительно друга повернуты на 180гр.
Уменя вообще создалось впечатление, что так поворачивать фазу нельзя. НО! Если перед прямым БПФ во FFT массиве не обнулять imag, а скопировать туда real но повернутый на 90гр.(пример делался с синусойдой), то весь этот алгоритм ДЕЙСТВУЕТ ПРАВИЛЬНО!, даже если усложнить вход до например перемножения разных частот, а в imag. тоже самое только sin поменять на cos, то все как надо.
Кто-нибуть сможет это объяснить?

[Sergey_gh]

To RA9YTJ

Все правильно. Вы получили то, что и должны были получить. Попробую объяснить.

1. копируем вх. буф. отчетов АЦП в массив БПФ, в real, imag обнуляем.

На вход БПФ Вы подали действительный сигнал. Для упрощения допустим что это cos(wt). Представим входной действительный сигнал cos(wt) в виде:

cos(wt) = 0.5(exp(iwt) + exp(-iwt))

т.е. в виде суммы двух частот w и -w. Графически это можно представить на еденичной окружности как два вращающихся в разные стороны еденичных вектора. Векторная сумма которых - это вектор лежащий на действительной оси. Т.е. наш исходный cos(wt). Далее БПФ суммы = сумме БПФ. Отсюда видно, что спектр действительного сигнала комплексный. Имеет отрицательные и положительные частоты (для нашего случая w и -w). При этом действительная часть и модуль симметрична относительно нуля частоты, а фаза и мнимая часть антисеммитрична.

3. По формулам сверху поворачиваем фазу.(например на 90гр, что эквивалентно FFTnew(i).real=-FFTold(i).imag FFTnew(i).imag=FFTol d(i).real

Осмелюсь предположить, что по формулам сверху вы крутанули в одну и ту же сторону и отрицательные и положительные частоты на 90гр. Т.е. умножили в частотной области всё на j (H(w)=exp(jφ) где φ = pi/2 => H(w)=j). Это эквивалентно свертке исходного сигнала с функцией:

h(n) = j исли n=0
h(n) = 0 исли n!=0

Проще говоря, Вы развернули сигнал вдоль мнимой оси и убили Re часть. Нужно было поворачивать отрицательные и положительные частоты в разные стороны. Для преобразователя Гильберта:

H(w) = -j если w>=0
H(w) = j если w<0

... если во вых. буф. копировать корень суммы квадратов FFT.real FFTold(i).imag(Пифаг ор) то получим поворот 180гр!!!!!

Должно быть просто |cos(wt)|. Никакого смысла в этой операции не вижу.

Я построил графики на всех этапах.
Что получилось, после поворота на 180гр, мы на вых в real получаем поворот противофазный(так и должно быть)!, при 90 и 270- 0!, при 45 0.5 амплитуды.

На 180гр не важно в какую сторону поворот, поэтому и получилось верно. При 45гр должно было быть 1/sqrt(2).

Если перед прямым БПФ во FFT массиве не обнулять imag, а скопировать туда real но повернутый на 90гр.(пример делался с синусойдой), то весь этот алгоритм ДЕЙСТВУЕТ ПРАВИЛЬНО!

Таким образом Вы подаете на вход аналитический сигнал не имеющий отрицательной частоты. Поэтому тут все работает.

P.S. Впечатлён Вашей внимательностью и аккуратностью. Посоветую классику "Л. Рабинер. Б. Гоулд. Теория и применение цифровой обработки сигналов". На русском языке, в интернете легко находится и скачивается.

[Relayer]

к примеру можно найти вот здесь:
http://library.tstu.edu.ua/resources...eering/Signals