Извлеките неотвернутую фазу комплексного входа
The dsp.PhaseExtractor Системная object™ извлекает неотвернутую фазу действительного или комплексного входа.
Чтобы извлечь неотвернутую фазу входного сигнала:
Создайте dsp.PhaseExtractor Объекту и установите его свойства.
Вызывайте объект с аргументами, как будто это функция.
Дополнительные сведения о работе системных объектов см. в разделе «Что такое системные объекты?».
Чтобы использовать функцию объекта, задайте системный объект в качестве первого входного параметра. Например, чтобы освободить системные ресурсы системного объекта с именем 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:
Объект действует на эту систему координат и формирует этот выход:
где:
Здесь я бегу от 1 до N. The angle функция возвращает угол фазы в радианах.
Если входной сигнал состоит из нескольких систем координат:
Если вы задаете TreatFramesIndependently на trueобъект обрабатывает каждую систему координат независимо. Поэтому в каждой системе координат объект вычисляет фазу с помощью предыдущей формулы, где:
равен 0.
равен 1.
Если вы задаете TreatFramesIndependently на falseобъект игнорирует контуры между системами координат. Поэтому в каждой системе координат step метод вычисляет фазу по предыдущей формуле, где:
- последняя неотвернутая фаза из предыдущей системы координат.
- последняя выборка из предыдущей системы координат.