Вычленение огибающей частоты на фоне несущей, Расчет в Матлаб Опции

[Stefan1]

Доброго времени суток!
Подскажите: как можно вычлененить 1-ю гармонику несущего сигнала на фоне модулирующего?
Т.е. при разложении сигнала в ряд Фурье, модулированного низкой частотой и выделении из него первой гармоники - у меня получается первая гармоника этого модулирующего сигнала, а мне нужна 1-я гармоника несущего сигнала.

Считаю в Матлабе так:
u=fft(U1); - раскладываю в ряд Фурье сигнал U1.
u([1 3:length(u)-1])=0; - обнуляю все элементы кроме второй гармоники.
idi0f=ifft(idi1f); - вывожу его первую гармонику.

Сообщение отредактировал Stefan1 - Nov 23 2012, 18:58

[V_G]

модулированного низкой частотой и выделении из него первой гармоники - у меня получается первая гармоника этого модулирующего сигнала, а мне нужна 1-я гармоника несущего сигнала.

Насколько я понял, это тема не про демодуляцию.
Уточнитесь со способом модуляции. Для правильно модулированного сигнала в спектре ВООБЩЕ не должно быть НЧ составляющих. Для АМ нужно искать надо не первую гармонику, а максимальную спектральную составляющую.
При ЧМ мне непонятем смысл 1-й гармоники несущего сигнала. Он там весь несущий, и его частота меняется по закону модуляции.

[Stefan1]

Уточнитесь со способом модуляции...

Сигнал выглядит так:

У этого сигнала провал посередине. Период сигнала = 0.77 с. Мне нужно вывести на график 1-ю гармонику с тем же периодом 0.77 с.
При этом частот модуляции я не знаю, знаю только частоту несущей.

Если сигнал комплексный - то sqrt(I*I + Q*Q) для каждого отсчета + ФНЧ с соответствующей частотой среза

Что-то не совсем понятно. Можно чуть подробнее.
Что-то вроде этого надо записать?
z = ddemod(Idi, Fc, Fd, Fs, 'qask/arb', M, opt1, opt2, opt3, opt4);
Но тут я не знаю откуда взять параметры: Fd, Fs, M, opt1, opt2...
И как мне затем отсюда выразить 1-ю гармонику?

Сообщение отредактировал Stefan1 - Nov 23 2012, 19:02

[DRUID3]

Сигнал выглядит так:
...
Знаю я только частоту несущей.
...
Что-то не совсем понятно.

Да уж... Действительно непонятно. И непонятен сам закон модуляции - а как тогда можно что-то реконструировать? Что это у Вас на картинке? Огибающая? Или сам сигнал? И что тогда значит "несущая"? Первая гармоника?

Можно чуть подробнее.

... не имеет смысла. Я думал у Вас нечто вроде радиоимпульса. А во вских

ddemod()

я не разбираюсь.

[Stefan1]

Да уж... Действительно непонятно. И непонятен сам закон модуляции - а как тогда можно что-то реконструировать? Что это у Вас на картинке? Огибающая? Или сам сигнал? И что тогда значит "несущая"? Первая гармоника?

Это сам сигнал. Несущая - имею ввиду частоту с которой следуют эти импульсы с провалом посередине (т.е. на графике период равен 0.77 с). Мне надо отфильтровать все высшие гармоники и низкочастотную составляющую модулирующего сигнала, оставить только 1-ю гармонику несущего сигнала (а не 1-ю гармонику модулирующего сигнала).
Если обнулить все высшие гармоники, кроме второй:
idi([1 3:length(idi)-1])=0;
то у меня будет 1-я гармоника НЧ модулирующего сигнала:

А как раз ее мне и надо убрать.

Сообщение отредактировал Stefan1 - Nov 23 2012, 08:29

[DRUID3]

Это сам сигнал. Несущая - имею ввиду частоту с которой следуют эти импульсы с провалом посередине (т.е. на графике период равен 0.77 с).

...это тактовая - частота первой гармоники.

Мне надо отфильтровать все высшие гармоники и низкочастотную составляющую модулирующего сигнала, оставить только 1-ю гармонику несущего сигнала (а не 1-ю гармонику модулирующего сигнала).

...ну так за чем дело стало?

Если обнулить все высшие гармоники, кроме второй:
idi([1 3:length(idi)-1])=0;
то у меня будет 1-я гармоника НЧ модулирующего сигнала:
А как раз ее мне и надо убрать.

...ну так выделите 1-ую. или в вашей нумерации это уже постоянная составляющая? Так значит Вы уже выделили ту, что нужно. А колеблется она потому как это никакая не АМ, а обычный спектр сигнала и это самая обычная его первая гармоника - в разные периоды доля энергии сигнала на нее выпадает разная оттуда и колебания.

[Stefan1]

...ну так выделите 1-ую. или в вашей нумерации это уже постоянная составляющая? Так значит Вы уже выделили ту, что нужно. А колеблется она потому как это никакая не АМ, а обычный спектр сигнала и это самая обычная его первая гармоника - в разные периоды доля энергии сигнала на нее выпадает разная оттуда и колебания.

Значит мне не она нужна. Мне нужна 1-я гармоника (имею ввиду обычную синусоиду, без высших гармоник) у которой должен быть период тот же, что и у несущей - 0.77 с.

Сообщение отредактировал Stefan1 - Nov 23 2012, 09:04

[DRUID3]

Значит мне не она нужна. Мне нужна 1-я гармоника (имею ввиду обычную синусоиду, без высших гармоник) у которой должен быть период тот же, что и у несущей - 0.77 с.

Я понял. Тем графиком(последним) Вы отрисовали гармонику(период показали) а я думал изменение этой гармоники от фрейма к фрейму. Ясно... Но... я так и не понял Ваших затруднений. Я не знаю синтаксис матлаба. Но в чем затруднение то? Сделайте FFT. Пройдитесь пороговой функцией(этож полноценный скриптовый язык) и выделите гармонику с номером 1 превышающую определенный порог.

[Stefan1]

Пройдитесь пороговой функцией(этож полноценный скриптовый язык) и выделите гармонику с номером 1 превышающую определенный порог.

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

Сообщение отредактировал Stefan1 - Nov 23 2012, 16:17