Генерировать реальные или сложные синусоидальные сигналы
Осциллятор с числовым управлением или NCO объект генерирует реальные или сложные синусоидальные сигналы. Амплитуда генерируемого сигнала всегда равна 1.
Для генерации действительных или сложных синусоидальных сигналов:
Создать dsp.NCO и задайте его свойства.
Вызовите объект с аргументами, как если бы это была функция.
Дополнительные сведения о работе системных объектов см. в разделе Что такое системные объекты?.
nco = dsp.NCONCO object™ системы, nco, который генерирует многоканальный реальный или комплексный синусоидальный сигнал с независимой частотой и фазой в каждом выходном канале.
nco = dsp.NCO( возвращает Name,Value)NCO Системный объект, nco, каждое указанное свойство имеет заданное значение.
Y = nco()PhaseIncrementSource и PhaseOffsetSource оба свойства имеют значение 'Property'.
Y = nco(phInc,OFFSET)Y, с приращением фазы, phIncи фазовое смещение, OFFSET, когда PhaseIncrementSource и PhaseOffsetSource свойства оба 'Input port'. phInc и OFFSET оба должны быть векторами строк одинаковой длины, где длина определяет количество каналов в выходном сигнале.
Чтобы использовать функцию объекта, укажите объект System в качестве первого входного аргумента. Например, для освобождения системных ресурсов объекта System с именем obj, используйте следующий синтаксис:
release(obj)
Примечание.Этот пример выполняется только в R2016b или более поздних версиях. При использовании более ранней версии замените каждый вызов функции эквивалентным step синтаксис. Например, myObject () становится step (myObject).
Проектирование источника NCO в соответствии с заданными спецификациями.
df = 0.05; % Frequency resolution = 0.05 Hz minSFDR = 96; % Spurious free dynamic range >= 96 dB Ts = 1/8000; % Sample period = 1/8000 sec dphi = pi/2; % Desired phase offset = pi/2;
Вычислите количество битов-накопителей, необходимых для заданного частотного разрешения.
Nacc = ceil(log2(1/(df*Ts)));
Достигнуто фактическое разрешение частоты.
actdf = 1/(Ts*2^Nacc);
Вычислить требуемое количество квантованных битов накопителя из требования SFDR
Nqacc = ceil((minSFDR-12)/6);
Расчет смещения фазы
phOffset = 2^Nacc*dphi/(2*pi);
Спроектируйте источник NCO.
nco = dsp.NCO('PhaseOffset', phOffset,... 'NumDitherBits', 4, ... 'NumQuantizerAccumulatorBits', Nqacc,... 'SamplesPerFrame', 1/Ts, ... 'CustomAccumulatorDataType', numerictype([],Nacc)); san = dsp.SpectrumAnalyzer('SampleRate', 1/Ts, ... 'PlotAsTwoSidedSpectrum', false);
Просмотрите выходные данные источника NCO на анализаторе спектра. Измените выходную частоту в середине моделирования с 510 Гц на 1520 Гц.
tic; while toc < 10 if toc < 5 F0 = 510; else F0 = 1520; end % Calculate the phase increment phIncr = int32(round(F0*Ts*2^Nacc)); y = nco(phIncr); san(y) end
Информация признака объекта NCO определяется следующими полями:
NumPointsLUT - Количество точек данных в таблице подстановки.
SineLUTSize - Размер таблицы поиска sine в четверти волны в байтах.
TheoreticalSFDR - Теоретический паразитный свободный динамический диапазон (SFDR) в dBc.
FrequencyResolution - Частотное разрешение NCO.
Получение вышеуказанных признаков для конкретного NCO объект, вызовите info функция на объекте.
nco = dsp.NCO
nco =
dsp.NCO with properties:
PhaseIncrementSource: 'Input port'
PhaseOffsetSource: 'Property'
PhaseOffset: 0
Dither: true
NumDitherBits: 4
PhaseQuantization: true
NumQuantizerAccumulatorBits: 12
PhaseQuantizationErrorOutputPort: false
Waveform: 'Sine'
SamplesPerFrame: 1
OutputDataType: 'Custom'
Show all properties
info(nco)
ans = struct with fields:
NumPointsLUT: 1025
SineLUTSize: 2050
TheoreticalSFDR: 84
FrequencyResolution: 1.5259e-05
Поля и соответствующие им значения изменяются в зависимости от настроек объекта. Например, если PhaseQuantization свойство имеет значение false, TheoreticalSFDR поле не отображается.
nco.PhaseQuantization = false; info(nco)
ans = struct with fields:
NumPointsLUT: 16385
SineLUTSize: 32770
FrequencyResolution: 1.5259e-05
Этот объект реализует алгоритм, входы и выходы, описанные на странице ссылок блока NCO. Свойства объекта соответствуют свойствам блока, за исключением того, что нет свойства объекта, соответствующего параметру Sample time block. Время выборки объектов составляет одну секунду.