Отправлено Artem 10 марта 2005 г. 05:32

1. Квантование 100 отсчетов на период
step = pi*2/100;

2. 10 периодов
x = 0:step:20*pi;

3. Входной сигнал.
inU = sin(x*1.2+0.2);

4. 2 опорных сигнала, для получения действительного и мнимого числа.
refRe = sin(x);
refIm = cos(x);

5. Находим разницу между частотой входного и порного сигнала.
for i=1:length(x)
errRe(i) = inU(i)* refRe(i);
errIm(i) = inU(i)* refIm(i);
err(i) = sqrt(errRe(i) + errIm(i));
end

на err(i)
имеем выходной сигнал, основаная гармоника которого соответсвует разницы между входным и выходным сигналом.

Интересует: Получить значение этой частоты, и в какую сторону она сдвинута.

(Практически уверен что это не реально, особенно без ФНЧ, но может чего общественность скажет)

• Не, ни хрена нереально, осчобенно, когда шумов нет:-) — -=ВН=-   (10.03.2005 13:13, 119 байт)
• Я так понимаю, что Вы решили заново изобрести ФАПЧ (PLL)? — AntZ   (10.03.2005 11:22, пустое)
  • О, ФАПЧ и никаких гвоздей. А без него снимать штаны и петь разлуку. Или застрелиться сразу. — -=ВН=-   (10.03.2005 13:16, пустое)
    • Ответ: Ой не надо про ВАПЧи, у меня от этих ЧФД и прочей пугани ... меряемой анализаторами спектра и измеряемыми в dBc уже тошнит :) — Microstrip   (16.03.2005 15:39, пустое)
  • Если честно, то да. — Artem   (10.03.2005 11:29, 277 байт)
    • Tutorials 18, 19 - Unlocking the Phase Locked Loop (PLL) - Parts 1,2 — AntZ   (10.03.2005 11:33, пустое, ссылка)
      • Спасибо. Короче облом. Без ФНЧ никак, так я и думал. А на ФНЧ у меня нет времени. — Artem   (10.03.2005 11:49, пустое)
        • Зачем слишком сложный фильтр городить в петле ФАПЧ - можно просто использовать ПИФ (два умножения одно сложение) — alostap   (10.03.2005 22:31, 145 байт)
        • Эффект ФНЧ может есть не только у стндартных equiripple FIR, частотная х-ка MA отчасти напоминает LP (иногда проходит как грубая апроксимация), думаю что IIR тоже можно попробовать, сильно дорогой не должен быть — AntZ   (10.03.2005 13:16, пустое)
          • Апроксимация это и есть ФНЧ. — Artem   (10.03.2005 14:54, пустое)
          • "Королева, мы в восхищении!" — -=ВН=-   (10.03.2005 13:42, пустое)
            • В крайности бросается! — Fat Robot   (10.03.2005 13:57, 11 байт)
        • Если у вас лично на него нет времени - поручите его процессору что используете и будет вам счасть. — советующий   (10.03.2005 12:11, пустое)
          • В том-то и дело что у проца. А ФНЧ придется делать очень длинный. Да и смысла нету тогда. Если у меня частота резко поменяется, простой пересчет лучше отреагирует. — Artem   (10.03.2005 12:27, пустое)
            • Вы лучше сформулируйте четко задачу - какой сигнал, как меняется и на чем планируется сделать — AntZ   (10.03.2005 13:13, пустое)
              • Все уже сделал. — Artem   (10.03.2005 13:23, 383 байт)