Быстрое преобразование Фурье - оптимизировано для генерации кода HDL
HDL FFT System object™ предоставляет две архитектуры для оптимизации пропускной способности или области. Используйте потоковую архитектуру Radix 2 ^ 2 для высокопроизводительных приложений. Эта архитектура поддерживает скалярные или векторные входные данные. С помощью векторного ввода можно достичь пропускной способности giga-sample-per-second (GSPS). Используйте пакетную архитектуру Radix 2 для минимальной реализации ресурсов, особенно с большими размерами FFT. Система должна быть в состоянии выдерживать пакетные данные и более высокую задержку. Эта архитектура поддерживает только скалярные входные данные. Объект принимает реальные или сложные данные, обеспечивает аппаратные управляющие сигналы и имеет дополнительные выходные управляющие сигналы кадра.
Для вычисления быстрого преобразования Фурье:
Создать dsp.HDLFFT и задайте его свойства.
Вызовите объект с аргументами, как если бы это была функция.
Дополнительные сведения о работе системных объектов см. в разделе Что такое системные объекты?.
FFT_N = dsp.HDLFFTFFT_N, который выполняет быстрое преобразование Фурье.
FFT_N = dsp.HDLFFT(Name,Value)
fft128 = dsp.HDLFFT('FFTLength',128)[ возвращает быстрое преобразование Фурье (FFT) при использовании пакетной архитектуры Radix 2. Y,validOut,ready] = FFT_N(X,validIn)ready сигнал указывает, когда объект может принимать входные выборки.
Чтобы использовать этот синтаксис, задайте для свойства Architecture значение 'Burst Radix 2'. Например:
FFT_N = dsp.HDLFFT(___,'Architecture','Burst Radix 2'); ... [y,validOut,ready] = FFT_N(x,validIn)
[ также возвращает сигналы управления кадрами Y,startOut,endOut,validOut] = FFT_N(X,validIn)startOut и endOut. startOut является true на первом образце кадра выходных данных. endOut является true для последней выборки кадра выходных данных.
Чтобы использовать этот синтаксис, задайте для свойств StartOutputPort и EndOutputPort значение true. Например:
FFT_N = dsp.HDLFFT(___,'StartOutputPort',true,'EndOutputPort',true); ... [y,startOut,endOut,validOut] = FFT_N(x,validIn)
[ возвращает FFT, когда Y,validOut] = FFT_N(X,validIn,resetIn)validIn является true и resetIn является false. Когда resetIn является trueобъект останавливает текущий расчет и очищает все внутренние состояния.
Чтобы использовать этот синтаксис, задайте для свойства ResetInputPort значение true. Например:
FFT_N = dsp.HDLFFT(___,'ResetInputPort',true); ... [y,validOut] = FFT_N(x,validIn,resetIn)
Чтобы использовать функцию объекта, укажите объект System в качестве первого входного аргумента. Например, для освобождения системных ресурсов объекта System с именем obj, используйте следующий синтаксис:
release(obj)
Создайте спецификации и входной сигнал.
N = 128; Fs = 40; t = (0:N-1)'/Fs; x = sin(2*pi*15*t) + 0.75*cos(2*pi*10*t); y = x + .25*randn(size(x)); y_fixed = sfi(y,32,24);
Запишите функцию, которая создает и вызывает системный object™. Из этой функции можно создать ЛПВП.
Примечание.Синтаксис этого объекта выполняется только в R2016b или более поздних версиях. При использовании более ранней версии замените каждый вызов объекта эквивалентным step синтаксис. Например, заменить myObject(x) с step(myObject,x).
function [yOut,validOut] = HDLFFT128(yIn,validIn) %HDLFFT128 % Processes one sample of FFT data using the dsp.HDLFFT System object(TM) % yIn is a fixed-point scalar or column vector. % validIn is a logical scalar value. % You can generate HDL code from this function. persistent fft128; if isempty(fft128) fft128 = dsp.HDLFFT('FFTLength',128); end [yOut,validOut] = fft128(yIn,validIn); end
Вычислите БПФ, вызвав функцию для каждой выборки данных.
Yf = zeros(1,3*N); validOut = false(1,3*N); for loop = 1:1:3*N if (mod(loop, N) == 0) i = N; else i = mod(loop, N); end [Yf(loop),validOut(loop)] = HDLFFT128(complex(y_fixed(i)),(loop <= N)); end
Отменить недопустимые выборки данных. Затем постройте график результатов частотного канала из БПФ.
Yf = Yf(validOut == 1); Yr = bitrevorder(Yf); plot(Fs/2*linspace(0,1,N/2), 2*abs(Yr(1:N/2)/N)) title('Single-Sided Amplitude Spectrum of Noisy Signal y(t)') xlabel('Frequency (Hz)') ylabel('Output of FFT (f)')
Создание спецификаций и входного сигнала. В этом примере используется 128-точечный БПФ и одновременно вычисляется преобразование по 16 выборкам.
N = 128; V = 16; Fs = 40; t = (0:N-1)'/Fs; x = sin(2*pi*15*t) + 0.75*cos(2*pi*10*t); y = x + .25*randn(size(x)); y_fixed = sfi(y,32,24); y_vect = reshape(y_fixed,V,N/V);
Запишите функцию, которая создает и вызывает системный object™. Функция не должна знать размер вектора. Объект сохраняет размер входного сигнала при первом вызове.
Примечание.Синтаксис этого объекта выполняется только в R2016b или более поздних версиях. При использовании более ранней версии замените каждый вызов объекта эквивалентным step синтаксис. Например, заменить myObject(x) с step(myObject,x).
function [yOut,validOut] = HDLFFT128V16(yIn,validIn) %HDLFFT128V16 % Processes 16-sample vectors of FFT data % yIn is a fixed-point column vector. % validIn is a logical scalar value. % You can generate HDL code from this function. persistent fft128v16; if isempty(fft128v16) fft128v16 = dsp.HDLFFT('FFTLength',128); end [yOut,validOut] = fft128v16(yIn,validIn); end
Вычислите БПФ, передав объекту 16-элементные векторы. Используйте getLatency функция для определения того, когда будет готова первая выборка выходных данных. Затем добавьте длину кадра, чтобы определить, сколько раз вызывать объект. Поскольку переменная объекта находится внутри функции, используйте второй объект для вызова getLatency. Используйте счетчик цикла для разворота validIn кому false после N входных выборок.
tempfft = dsp.HDLFFT; loopCount = getLatency(tempfft,N,V)+N/V; Yf = zeros(V,loopCount); validOut = false(V,loopCount); for loop = 1:1:loopCount if ( mod(loop,N/V) == 0 ) i = N/V; else i = mod(loop,N/V); end [Yf(:,loop),validOut(loop)] = HDLFFT128V16(complex(y_vect(:,i)),(loop<=N/V)); end
Отменить недопустимые выходные выборки.
C = Yf(:,validOut==1); Yf_flat = C(:);
Постройте график данных частотного канала из БПФ. Выход БПФ выполнен в битовом порядке. Измените порядок перед построением графика.
Yr = bitrevorder(Yf_flat); plot(Fs/2*linspace(0,1,N/2),2*abs(Yr(1:N/2)/N)) title('Single-Sided Amplitude Spectrum of Noisy Signal y(t)') xlabel('Frequency (Hz)') ylabel('Output of FFT (f)')
Задержка объекта зависит от длины БПФ и размера вектора. Используйте getLatency для поиска задержки конкретной конфигурации. Задержка - это количество циклов между первым действительным входом и первым действительным выходом, предполагая, что вход является непрерывным.
Создание нового dsp.HDLFFT и запросить задержку.
hdlfft = dsp.HDLFFT('FFTLength',512);
L512 = getLatency(hdlfft)L512 = 599
Запрос информации о гипотетической задержке для аналогичного объекта с другой длиной БПФ. Свойства исходного объекта не изменяются.
L256 = getLatency(hdlfft,256)
L256 = 329
N = hdlfft.FFTLength
N = 512
Запрос информации о гипотетической задержке аналогичного объекта, который принимает ввод вектора из восьми выборок.
L256v8 = getLatency(hdlfft,256,8)
L256v8 = 93
Разрешить масштабирование на каждом этапе БПФ. Задержка не изменяется.
hdlfft.Normalize = true; L512n = getLatency(hdlfft)
L512n = 599
Запросите тот же порядок вывода, что и порядок ввода. Задержка увеличивается, поскольку объект должен собрать выходные данные перед переупорядочиванием.
hdlfft.BitReversedOutput = false; L512r = getLatency(hdlfft)
L512r = 1078
Потоковая архитектура Radix 2 ^ 2 реализует архитектуру с низкой задержкой. Это экономит ресурсы по сравнению с потоковой реализацией Radix 2 за счет факторинга и группирования уравнения FFT. Архитектура имеет этапы log4 (N). Каждая ступень содержит две кнопки обратной связи с однопутевой задержкой (SDF) с контроллерами памяти. При использовании векторного ввода каждый этап работает на меньшем количестве входных выборок, поэтому некоторые этапы сводятся к простой бабочке, без SDF.
Первая ступень SDF - обычная бабочка. Второй каскад умножает выходы первого каскада на -j. Чтобы избежать аппаратного множителя, блок заменяет действительную и мнимую части входов и снова заменяет мнимые части результирующих выходов. Каждый этап округляет результат умножения множителя до длины входного слова. Сдвоенные коэффициенты имеют два целых бита, а остальные биты используются для дробных битов. Сдвоенные коэффициенты имеют такую же битовую ширину, как и входные данные WL. Множители twiddle имеют два целых бита и WL-2 дробные биты.
При включении масштабирования алгоритм делит результат каждого этапа бабочки на 2. Масштабирование на каждом этапе позволяет избежать переполнения, сохранить длину слова на уровне входных данных и получить общий масштабный коэффициент 1/N. Если масштабирование отключено, алгоритм избегает переполнения, увеличивая длину слова на 1 бит на каждом этапе. На диаграмме показаны бабочки и внутренние длины слов каждого этапа, не включая память.
Пакетная архитектура Radix 2 реализует FFT с использованием единого комплексного умножителя бабочек. Алгоритм не может запуститься, пока не сохранит весь входной кадр, и не сможет принять следующий кадр до завершения вычислений. Выходной порт готовности указывает, когда алгоритм готов к новым данным. На схеме показана архитектура пакета с конвейерными регистрами.
Алгоритм обрабатывает входные данные только в том случае, если входной допустимый порт равен 1. Выходные данные действительны, только если выходной допустимый порт равен 1.
Если дополнительный входной порт сброса равен 1, алгоритм останавливает текущее вычисление и сбрасывает все внутренние состояния. Алгоритм начинает новые вычисления, когда порт сброса равен 0, а входной допустимый порт запускает новый кадр.
Эта диаграмма показывает входные и выходные допустимые значения портов для смежных скалярных входных данных, потоковой архитектуры Radix 2 ^ 2, длины FFT 1024 и размера вектора 16.
На схеме также показаны необязательные значения начального и конечного портов, указывающие границы кадра. Если начальный порт включен, значение начального порта будет импульсным в течение одного цикла с первым допустимым выходом кадра. Если конечный порт включен, значение начального порта будет импульсным в течение одного цикла с последним допустимым выходом кадра.
При применении непрерывных входных кадров выходные данные также будут непрерывными после начальной задержки.
Входной допустимый порт может быть несмежным. Данные, сопровождаемые входным допустимым портом, обрабатываются по мере их поступления, и полученные данные сохраняются до тех пор, пока не будет заполнен кадр. Затем алгоритм возвращает непрерывные выходные выборки в кадре из N циклов (длины БПФ). Эта диаграмма показывает несмежный вход и смежный выход для длины БПФ 512 и размера вектора 16.
При использовании архитектуры пакета невозможно предоставить следующий кадр входных данных до тех пор, пока не будет доступно пространство памяти. Порт готовности указывает, когда алгоритм может принимать новые входные данные.
Задержка зависит от длины БПФ и размера вектора. Используйте getLatency для поиска задержки конкретной конфигурации. Задержка - это количество циклов между первым действительным входом и первым действительным выходом, предполагая, что вход является непрерывным.
При использовании пакетной архитектуры с непрерывным входом, если проект ожидает готовности к выходу 0 перед отменой действительности ввода на вход поступает один дополнительный цикл данных. Этот образец данных является первым образцом следующего кадра. Алгоритм может сохранить один образец при обработке текущего кадра. Благодаря этому опережению одной выборки наблюдаемая задержка более поздних кадров (от действительного входного до действительного выходного) на один цикл короче, чем сообщаемая задержка. Задержка измеряется от первого цикла, когда действительность ввода равна 1, до первого цикла, когда действительность вывода равна 1. Число циклов между, когда порт готовности равен 0, а выходной допустимый порт равен 1, всегда является задержкой - FFTLength.
