Товарищи! я новичок в этом деле, очень прошу помочь.
Мой вопрос состоит в следующем. Необходимо построить цифровой детектор. Пусть есть сигнал s(t)=А*cos(wt+). В моем понимании квадратурной обработки я должен умножить его сначала на сигнал с выхода генераторов cos(wt+) и sin(wt+). Далее пропустить результаты произведения через ФНЧ.
Какие следующие операции необходимо для детектирования сигнала фазовой манипуляции?
Встречал другой способ детектирования, где происходит умножение принятого сигнала на exp(-I*w*t), какой формы должен быть принятый сигнал и что происходит при умножении на exp(-I*w*t)?

edmund Гм, прикольно. Книжки есть для новичков и нетолько. Тяжев, Окунев, Сергеенко, Марпл младший, Прокис и в таком разрезе. Есть тут на сайте, шукайте. :)

в книгах я смотрел, там пишут, что после фильтрации сигналов с выхода опорных генераторов для фазовой демодуляции необходимо вычислить фазовый угол. что мне это дает?
а по поводу умножения на exp(-I*w*t) необходимо сначала представить сигнал в комплексной форме, для чего я провожу преобразование гилберта. после перемножения опять ФНЧ, а далее вычисление фазового угла?

edmund
в книгах я смотрел, там пишут, что после фильтрации сигналов с выхода опорных генераторов для фазовой демодуляции необходимо вычислить фазовый угол. что мне это дает?
по фазовому углу принимается решение о переданном символе

то есть к примеру при двоичной фазовой манипуляции угол 0 и 180. именно такие значения фазового угла я и получу?

edmundто есть к примеру при двоичной фазовой манипуляции угол 0 и 180. именно такие значения фазового угла я и получу?в идеале да

я строю детектор в матлаб. беру сигнал s(t)=А*cos(wt+phi) умножаю его на стробирующий сигнал. после чего полученный результат перемножаю с s1(t)=s(t)*cos(wt + phi) и s2(t)=s(t)*(-sin(wt + phi)). фильтрую и получаю Real и Im составляющие. Вычисляю фазовый угол angle(s1(t)-1i*s2(t)) по переменной t?

ошибся***
Вычисляю фазовый угол angle(s1(t)-1i*s2(t)) и строю график по переменной t? или для восставновления импульсов мне необходимо получить амплитуду?

edmund А сигналик который вы там мучаете положите а? Без метерьяла оно можно говорить сколько угодно и о чём угодно.

А сигналик который вы там мучаете положите а
да сигналика то вобщем то особенного никакого нет.

Делаю так
A0 =1;
w0 = 50e6
st = 1e-8;
t = 0 : st : 5e-5 - st;
Strobe = 0.5 * square(1e6 * t, 70) + 0.5;
S0 = A0 * cos(w0 * t + phi0).*Strobe ;

потом перемножаю S0 на синус и косинус

edmund А, сигналика нет значит. Ну и ладно, значит не судьба. Не люблю теоризысков, можно обговорится, а на реал сигналах встрять наглухо. Так что пардоне, я пас.

edmund
А зачем это Вам? Если что-то типа курсовика - это один разговор, если для декодирования реальных сигналов- совсем другой. Начать хотя бы с того, что на реальных, если идти вашим методом, придётся ещё искать истинную частоту сигнала. + компенсировать многолучёвость.

то есть к примеру при двоичной фазовой манипуляции угол 0 и 180
В реале без дополнительной обработки 0 и 180 таким способом Вы не получите.

А возможно ли компенсировать многолучёвсоть, если сигнал уже принят и записан в файл, и при его приёме многолучёвость не компенсировалась?

Smerch66 Запись положите, и там видно будет. А так, вариант мерятся линейками и говорить ни об чём гарантирован. :)

Smerch66
А возможно ли компенсировать многолучёвсоть, если сигнал уже принят и записан в файл, и при его приёме многолучёвость не компенсировалась?
В общем от сигнала зависит. Если известны парамметры- один разговор, если ничего неизвестно- другой. В большинстве случаев, если сигнал правильно записан, вопрос решаемый.

Подскажите пожалуйста, а если у меня на микросхеме 2-I (положительный и отрицательный) и 2-Q выхода, то как можно это все загнать в компьютер?
Поможет ли, если на один канал стереовхода SB дать 2-I а на другой 2-Q (не повлияют ли они друг на друга)?
Сигнал беру с выхода м/сх TDA8263.
Если кто знает, то подскажите пожалуйста.

Разместил статью про квадратурной обработке сигнала с использованием программы Матлаб. Может быть, она для кого-нибудь окажется полезной.

Я разместил статью о квадратурной обработке сигналов в процессе их генерации в цифровых модуляторах CPFSK и PSK.

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

Vanadi Большое спасибо вам за статью. Но к сожалению математику уже подзабыл, и перевести эти формулы в целочисленные алгоритмы для меня уже затруднительно. Хотя суть дела описана действительно хорошо и понятно.

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

Vanadi К стати, это ваша задумка отображать вещественный сигнал одной стрелкой, а комплексный двойной? Нахожу это чрезвычайно мощной и наглядной вещью, иногда реально в голове пухнет, где что, очень удачный ход, если ваш, то вдвойне удачный. Ну и статья само собой отличная.

Programmist
xman
Mesh
Я благодарен Вам за ваше поощряющее одобрение и признание, что в эпохе картиночных коммуникаций стоит приподнимать информативность блок-схем (между прочим) увеличивая число разновидностей стрелок соответствующих разным классам сигналов. Такое различие сигналов: вещественный / комплексный я позаимствовал с учебника: Трахтман А.М. "Введение в обобщенную спектральную теорию сигналов", Москва 1972, стр. 30, рис. 1.17 и 1.18, несколько лет тому назад.

Vanadi У меня вопрос по последней статье. В формуле 3, в правой части появился параметер альфа? Нет ли ошибки опечатки и там должно быть d[n/L], а не как написано d[a/L]? Собственно это не так уж важно, понятно всё из контектста, но всё ж. И с функцией fra неясно в вычислениях pwr, формула 9, вроде round есть в МатЛабе с аналогичным резалтом? По описанию в статье, fra это округление до ближайшего целого, но по резалтам это извлечение дробной части c "перескоком через целое" если можно так выразится, тож скорее всего опечатка в описании функционала fra. Вобще теория вами преподнесена шикарно, имхо.

Mesh
Я Вам искренне благодарен за вашу тщательную внимательность. Вы дважды правы. Вместо [альфа/L] в формуле (3) должно быть [n/L] и после (10) вместо "округление" должно быть "дробная часть (fraction)". Всех читателей прошу прощения.

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

Разместил очередную статью на тему квадратурной обработки радиосигналов а именно про цифровую фазоразностную модуляцию при нецелом отношении скорости дискретизации Fs к скорости модуляции (бодовой скорости) FB на примере моды SDPSK. Эта статья есть продолжением всех трёх статей предлагаемого цикла.

Vanadi

Допустим на входе 3.(3) отсчёта на символ, увеличиваем полифазным фильтром в два раза до 6.(6) и берём ближай к точке отстутствия МСИ отсчёт, чем ваш интерполятор лучше?

petr0v Уточните, пожалуйста, ваш вопрос. Если на входе (интерполятора?) имеется 3,33... Sa/symbol а на выходе (интерполятора) хотите иметь 6,66 Sa/symbol то ваша проблема интерполяции с целым шагом L = 2 вне темы моей статьи.
Если Вас интересует сравнение моего двускоростного упрощенного интерполятора с нецелым L = 20/3 с, упомянутым в вашем вопросе, трехскоростным каскадным интерполятором с шагом L = L1 x L2 = L2 x L1; L1 = 2, L2 = 10/3, то без каких либо оговорок интерполяция кратности L=20/3 (при сохранении её качества) по схеме с Рис. 1 и формуле (3) вычислительно самая дешёвая.

Vanadi

Уточните, пожалуйста, ваш вопрос.

Практическая задача, в приёмнике АЦП работает с 3.33... отсчёта на символ, в идеале мне нужно интерполировать один отсчёт на символ в точке отсутствия МСИ синхронно с символьной частотой передатчика. Допустим устраивает ближайший отсчёт на частоте 6.66... отсчётов на символ. Берём схему с рисунка 2, при полифазной структуре фильтр разделяется на два параллельных с вдвое меньшим количеством коэффициентов в каждом и работает на исходной частоте 3.33..., с выхода фильтра мультиплексируем ближайший отсчёт. Как при вашем подходе к интерполяции мне при той же МСИ получить меньшие вычислительные затраты?