Извлеките развернутую фазу комплексного входа
Система dsp.PhaseExtractor object™ извлекает развернутую фазу действительного или комплексного входа.
Извлекать развернутую фазу входа сигнала:
Создайте объект dsp.PhaseExtractor и установите его свойства.
Вызовите объект с аргументами, как будто это была функция.
Чтобы узнать больше, как Системные объекты работают, смотрите то, Что Системные объекты? MATLAB.
phase = dsp.PhaseExtractorphase = dsp.PhaseExtractor(Name,Value)Для версий ранее, чем R2016b, используйте функцию step, чтобы запустить алгоритм Системного объекта. Аргументы к step являются объектом, который вы создали, сопровождаемый аргументами, показанными в этом разделе.
Например, y = step(obj,x) и y = obj(x) выполняют эквивалентные операции.
p = phase(input)Чтобы использовать объектную функцию, задайте Системный объект как первый входной параметр. Например, чтобы выпустить системные ресурсы Системного объекта под названием obj, используйте этот синтаксис:
release(obj)
Примечание: Этот пример запускается только в R2016b или позже. Если вы используете более ранний релиз, заменяете каждый вызов функции с эквивалентным синтаксисом step. Например, myObject (x) становится шагом (myObject, x).
Создать Системный объект dsp.SineWave?. Укажите, что объект генерирует экспоненциальный вывод с комплексной экспонентой.
sine = dsp.SineWave('Frequency',10,... 'ComplexOutput',true,'SamplesPerFrame',128);
Создать Системный объект dsp.PhaseExtractor?. Укажите, что объект игнорирует контуры кадра при возврате развернутой фазы.
phase = dsp.PhaseExtractor('TreatFramesIndependently',false);
Извлеките развернутую фазу синусоиды. Построить фазу по сравнению со временем с помощью Системного объекта dsp.TimeScope?.
timeplot = dsp.TimeScope('PlotType','Line','SampleRate',1000,... 'TimeSpan',1.5,'YLimits',[0 80],... 'ShowGrid',true,... 'YLabel','Unwrapped Phase (rad)'); for ii = 1:10 sineOutput = sine(); phaseOutput = phase(sineOutput); timeplot(phaseOutput) end
Примечание: Этот пример запускается только в R2016b или позже. Если вы используете более ранний релиз, заменяете каждый вызов функции с эквивалентным синтаксисом step. Например, myObject (x) становится шагом (myObject, x).
Создать Системный объект dsp.TransferFunctionEstimator?.
tfe = dsp.TransferFunctionEstimator('FrequencyRange','centered');
Создать Системный объект dsp.PhaseExtractor?. Укажите, что объект должен обработать каждый кадр информационно-независимо.
phase = dsp.PhaseExtractor('TreatFramesIndependently',true);
Создать Системный объект dsp.IIRFilter?. Вычислите передаточную функцию БИХ-фильтра третьего порядка. Используйте функцию butter, чтобы сгенерировать коэффициенты для фильтра.
[b,a] = butter(3,.3); iir = dsp.IIRFilter('Numerator',b,'Denominator',a);
Извлеките фазовый отклик передаточной функции. График с помощью Системного объекта dsp.ArrayPlot?.
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:
Метод step действует на этот кадр и производит этот вывод:
где:
Здесь, i запускаюсь от 1 до N. Функция angle возвращает угол фазы в радианах.
Если входной сигнал состоит из нескольких кадров:
Если вы устанавливаете TreatFramesIndependently на true, метод step обрабатывает каждый кадр независимо. Поэтому в каждом кадре, метод step вычисляет фазу с помощью предыдущей формулы где:
0.
1.
Если вы устанавливаете TreatFramesIndependently на false, метод step игнорирует контуры между кадрами. Поэтому в каждом кадре, метод step вычисляет фазу с помощью предыдущей формулы где:
последняя развернутая фаза от предыдущего кадра.
последняя выборка от предыдущего кадра.