Помощь - Поиск - Пользователи - Календарь
Полная версия этой страницы: Распределение фаз.
Форум разработчиков электроники ELECTRONIX.ru > Цифровая обработка сигналов - ЦОС (DSP) > Алгоритмы ЦОС (DSP)
ivan219
Есть 255 частот которые кратны целому числу т.е. 1 частота 1 полный периуд 2 два периуда 3 три и.тд. нужно сложить из фазы так чтобы значения в выходном сигнале были как можно ближе друг другу т.е. разница между максимальным и минимальным значением была как можно меньше.

Если такой алгоритм распределения фаз???

Можно конечно воспользоватся Рандомом но может всётаки есть более универсальный smile.gif

Код
for I1 := 1 to 255 do
  begin
   R := Random * 2;
   for I := 0 to 511 do
    SignalOut[I]:= SignalOut[I] +  600 * Sin(2 * Pi * I * I1 / 512 + Pi * R);
  end;
petrov
ifft(ones(255,1)) все поднесущие сложились в фазе и на выходе получаем дельта-импульс
alex_os
[quote name='ivan219' date='Aug 11 2009, 16:09' post='633530']
..
Если такой алгоритм распределения фаз???
quote]

OFDM PARP reduction?
petrov
Наверное неправильно понял. Если имеется ввиду минимальный пик-фактор на выходе то закон распределения фаз квадратичный, на выходе получится ЛЧМ импульс.
ivan219
А по подробнее можно.
ЛЧМ это помойму плавно нарастающая частота а мне нужно все 255 частот звучали одновременно.
petrov
Цитата(ivan219 @ Aug 12 2009, 15:09) *
А по подробнее можно.
ЛЧМ это помойму плавно нарастающая частота а мне нужно все 255 частот звучали одновременно.


Подробнее смотрите например в учебниках по радиолокации, амплитудный спектр ЛЧМ константа - содержит все частоты, а фазовый меняется по квадратичному закону. Вообще непонятно что вам нужно.
ivan219
ЛЧМ это хорошо но вот как подобрать коэфиценты А и В (Phase[k] = A * k + B * k * k) ???
petrov
Цитата(ivan219 @ Aug 13 2009, 13:57) *
ЛЧМ это хорошо но вот как подобрать коэфиценты А и В (Phase[k] = A * k + B * k * k) ???


angle(exp(i*pi*[0:256-1].^2/256))
ivan219
Цитата
angle(exp(i*pi*[0:256-1].^2/256))


angle это что???

[0:256-1] это я так понял от нуля и до 255???

^2/256 2 делёное на 256???
petrov
Цитата(ivan219 @ Aug 13 2009, 15:24) *
angle это что???

[0:256-1] это я так понял от нуля и до 255???

^2/256 2 делёное на 256???


Это матлабовская нотация, копируете в коммандную строку и получается массив фаз. angle - эта функция возвращает аргумент комплексного числа.
[0:256-1] - это вектор строка 0...255. [0:256-1].^2 - это поэлементное возведение в квадрат вектора строки 0, 1^2, 2^2,...255^2. / - это деление возведённого в квадрат вектора на 256, приоритет операций в матлабе такой же как и в стандартной математической нотации.
ivan219
С матлабом не в ладах можно как нибудь на язык высокого уровня перевести smile.gif Delphi, C++
petrov
Цитата(ivan219 @ Aug 13 2009, 15:41) *
С матлабом не в ладах можно как нибудь на язык высокого уровня перевести smile.gif Delphi, C++


Матлаб де факто язык инженерных вычислений, он более высокого уровня чем Delphi, C++, на нём например можно свободно обращаться с матрицами как с обычными числами, чем быстрее вы начнёте использовать матлаб тем быстрее получите наглядный результат того о чём пишут в книжках по ЦОС.
ivan219
Блин ну я то программу не на матлабе пишу sad.gif мне надо чисто только фазы.

А вы может выложить данные которые получились у вас???
petrov
Цитата(ivan219 @ Aug 13 2009, 20:32) *
Блин ну я то программу не на матлабе пишу sad.gif мне надо чисто только фазы.


Так я тоже в конечном итоге на VHDL пишу, в матлабе, симулинке, маткаде проверка идей, приходит понимание сути.


Цитата(ivan219 @ Aug 13 2009, 20:32) *
А вы может выложить данные которые получились у вас???


Какие данные то?
ivan219
Цитата
Какие данные то?


Значение фаз сигнала.
petrov
Цитата(ivan219 @ Aug 13 2009, 22:51) *
Значение фаз сигнала.


сумма гармоник с низким пик-фактором ifft(exp(i*pi*[0:256-1].^2/256))

фазы поднесущих меняются pi*0^2/256, pi*1^2/256, pi*2^2/256, ... pi*255^2/256
ivan219
Цитата
фазы поднесущих меняются pi*0^2/256, pi*1^2/256, pi*2^2/256, ... pi*255^2/256


То что нужно biggrin.gif

Спасибо.
Для просмотра полной версии этой страницы, пожалуйста, пройдите по ссылке.
Invision Power Board © 2001-2025 Invision Power Services, Inc.