Извлечение развернутой фазы сложного входа
dsp.PhaseExtractor Система object™ извлекает развернутую фазу вещественного или сложного ввода.
Для извлечения развернутой фазы входного сигнала:
Создать dsp.PhaseExtractor и задайте его свойства.
Вызовите объект с аргументами, как если бы это была функция.
Дополнительные сведения о работе системных объектов см. в разделе Что такое системные объекты?.
Чтобы использовать функцию объекта, укажите объект System в качестве первого входного аргумента. Например, для освобождения системных ресурсов объекта System с именем obj, используйте следующий синтаксис:
release(obj)
Примечание.Этот пример выполняется только в R2016b или более поздних версиях. При использовании более ранней версии замените каждый вызов функции эквивалентным step синтаксис. Например, myObject (x) становится шагом (myObject, x).
Создать dsp.SineWave object™ системы. Укажите, что объект создает экспоненциальный вывод со сложной экспонентой.
sine = dsp.SineWave('Frequency',10,... 'ComplexOutput',true,'SamplesPerFrame',128);
Создать dsp.PhaseExtractor Системный объект?. Укажите, что объект игнорирует границы кадра при возврате развернутой фазы.
phase = dsp.PhaseExtractor('TreatFramesIndependently',false);
Извлеките развернутую фазу синусоидальной волны. Постройте график зависимости фазы от времени с помощью timescope Системный объект.
timeplot = timescope('PlotType','Line','SampleRate',1000,... 'TimeSpanSource','Property','TimeSpan',1.5,'YLimits',[0 80],... 'ShowGrid',true,... 'YLabel','Unwrapped Phase (rad)'); for ii = 1:10 sineOutput = sine(); phaseOutput = phase(sineOutput); timeplot(phaseOutput) end
Примечание.При использовании R2016a или более ранней версии замените каждый вызов объекта эквивалентным step синтаксис. Например, obj(x) становится step(obj,x).
Создать dsp.TransferFunctionEstimator object™ системы.
tfe = dsp.TransferFunctionEstimator('FrequencyRange','centered');
Создать dsp.PhaseExtractor object™ системы. Укажите, что объект должен обрабатывать каждый кадр данных независимо.
phase = dsp.PhaseExtractor('TreatFramesIndependently',true);
Создать dsp.IIRFilter object™ системы. Вычислите передаточную функцию БИХ-фильтра третьего порядка. Используйте butter для формирования коэффициентов для фильтра.
[b,a] = butter(3,.3); iir = dsp.IIRFilter('Numerator',b,'Denominator',a);
Извлеките фазовый отклик передаточной функции. Печать с использованием dsp.ArrayPlot object™ системы.
sampleRate = 1e3; phaseplot = dsp.ArrayPlot('PlotType','Line','XOffset',-sampleRate/2,... 'YLimits',[-15 0],... 'YLabel','Phase Response (rad)',... 'XLabel','Frequency (Hz)',... 'Title','System Phase response'); for ii = 1:100 % Generate input input = 0.05*randn(1000,1); % Pass through IIR filter filterOutput = iir(input); % Estimate transfer function transferFunction = tfe(input,filterOutput); % Plot transfer function phase phaseOutput = phase(transferFunction); phaseplot(phaseOutput); end
Рассмотрим входной кадр длиной N:
Объект действует на этот кадр и выдает следующие выходные данные:
где:
xi − 1 * xi)
Здесь я пробегаю от 1 до N. angle функция возвращает фазовый угол в радианах.
Если входной сигнал состоит из нескольких кадров:
Если установить TreatFramesIndependently кому trueобъект обрабатывает каждый кадр независимо. Поэтому в каждом кадре объект вычисляет фазу, используя предыдущую формулу, где:
равно 0.
равно 1.
Если установить TreatFramesIndependently кому falseобъект игнорирует границы между кадрами. Поэтому в каждом кадре step способ вычисляет фазу по предыдущей формуле, где:
- последняя развернутая фаза из предыдущего кадра.
- последняя выборка из предыдущего кадра.