Быстрое преобразование Фурье — оптимизированный для генерации HDL-кода
Система БПФ HDL object™ обеспечивает две архитектуры, чтобы оптимизировать или пропускную способность или область. Используйте Основание потоковой передачи 2^2 архитектура для высокопроизводительных приложений. Эта архитектура поддерживает скалярные или векторные входные данные. Можно достигнуть пропускной способности выборки гига сэмплов в секунду (GSPS) с помощью векторного входа. Используйте пакетное Основание 2 архитектуры для минимизации ресурса размещения, особенно с большими размерами БПФ. Ваша система должна смочь допустить принятие пульсирующих данных. Эта архитектура поддерживает только скалярные входные данные. Объект принимает действительные или комплексные данные, обеспечивает благоприятные для оборудования управляющие сигналы и имеет дополнительные выходные сигналы управления кадром.
Вычислить быстрое преобразование Фурье:
Создайте dsp.HDLFFT объект и набор его свойства.
Вызовите объект с аргументами, как будто это была функция.
Чтобы узнать больше, как Системные объекты работают, смотрите то, Что Системные объекты?.
FFT_N = dsp.HDLFFTFFT_N, это выполняет быстрое преобразование Фурье.
FFT_N = dsp.HDLFFT(Name,Value)
fft128 = dsp.HDLFFT('FFTLength',128)[ возвращает быстрое преобразование Фурье (FFT) при использовании пакетного Основания 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)
[ возвращает БПФ когда 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)
Чтобы использовать объектную функцию, задайте Системный объект как первый входной параметр. Например, чтобы выпустить системные ресурсы Системного объекта под названием 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™. Можно сгенерировать HDL от этой функции.
Примечание: Этот объектный синтаксис запускается только в 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(:);
Отобразите данные о канале частоты на графике из БПФ. FFT выход находится в обратном битовом порядке. Переупорядочьте его перед графическим выводом.
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
Основание потоковой передачи 2^2 архитектура реализует архитектуру низкой задержки. Это сохраняет ресурсы по сравнению с Основанием потоковой передачи 2 реализации путем факторинга и группировки уравнения FFT. Архитектура имеет log4 (N) этапы. Каждый этап содержит двух бабочек обратной связи задержки одно пути (SDF) с контроллерами памяти. Когда вы используете векторный вход, каждый этап работает с меньшим количеством входных выборок, таким образом, некоторые этапы уменьшают до простой бабочки без SDF.
Первый этап SDF является обычной бабочкой. Второй этап умножает выходные параметры первой стадии –j. Чтобы избежать аппаратного множителя, блок подкачивает действительные и мнимые части входных параметров, и снова подкачивает мнимые части получившихся выходных параметров. Каждый этап округляет результат вертеть факторного умножения к входному размеру слова. Вертеть факторы имеют два целочисленных бита, и остальная часть битов используется для дробных битов. Вертеть факторы имеют ту же битную ширину как входные данные, WL. Вертеть факторы имеют два целочисленных бита и WL-2 дробных бита.
Если вы позволяете масштабироваться, алгоритм делит результат каждого этапа бабочки 2. Масштабирование на каждом этапе избегает переполнения, сохраняет размер слова тем же самым как вход и приводит к полному масштабному коэффициенту 1/N. Если масштабирование отключено, алгоритм избегает переполнения путем увеличения размера слова на 1 бит на каждом этапе. Схема показывает бабочек и внутренние размеры слова каждого этапа, не включая память.
Пакетное Основание 2 реализации архитектуры БПФ при помощи одного комплексного множителя бабочки. Алгоритм не может запуститься, пока он не сохранил целый входной кадр, и он не может принять следующую систему координат, пока расчеты не завершены. Порт выхода ready указывает, когда алгоритм готов к новым данным. Схема показывает пакетную архитектуру с конвейерными регистрами.
Входные данные процессов алгоритма только, когда порт входа valid равняется 1. Выходные данные допустимы только, когда порт выхода valid равняется 1.
Когда дополнительный порт входа reset равняется 1, алгоритм останавливает текущее вычисление и очищает все внутренние состояния. Алгоритм начинает новые вычисления, когда порт reset 0, и порт входа valid запускает новую систему координат.
Эта схема показывает значения порта valid ввода и вывода для непрерывных скалярных входных данных, передавая Основание потоком 2^2 архитектура, длина БПФ 1 024 и размер вектора 16.
Схема также показывает дополнительные значения порта start и end, которые указывают на контуры системы координат. Если вы включаете порт start, импульсы значения порта start для одного цикла с первым допустимым выходом системы координат. Если вы включаете порт end, импульсы значения порта start для одного цикла с последним допустимым выходом системы координат.
Если вы примените непрерывные входные кадры, выход также будет непрерывен после начальной задержки.
Порт входа valid может быть, состоящим из нескольких несмежных участков. Данные, сопровождаемые портом входа valid, обрабатываются, когда они прибывают, и получившиеся данные хранятся, пока система координат не заполнена. Затем алгоритм возвращает непрерывные выходные выборки в системе координат N (FFT length) циклы. Эта схема показывает вход, состоящий из нескольких несмежных участков, и непрерывный выход для длины БПФ 512 и размера вектора 16.
Когда вы используете пакетную архитектуру, вы не можете обеспечить следующую систему координат входных данных, пока пространство памяти не доступно. Порт ready указывает, когда алгоритм может принять новые входные данные.
Задержка меняется в зависимости от длины БПФ и размера вектора. Используйте getLatency функционируйте, чтобы найти задержку особой настройки. Задержка является количеством циклов между первым допустимым входом и первым допустимым выходом, принимая, что вход непрерывен.
При использовании пакетной архитектуры с непрерывным входом, если ваш проект ожидает ready к выходу 0 перед de-утверждением вход valid затем один дополнительный цикл данных прибывает во вход. Эта выборка данных является первой выборкой следующей системы координат. Алгоритм может сохранить одну выборку при обработке текущего кадра. Из-за этого демонстрационного усовершенствования, наблюдаемая задержка более поздних систем координат (от входа valid до выхода valid) является одним циклом короче, чем задержка, о которой сообщают. Задержка измеряется от первого цикла, когда вход valid 1 к первому циклу, когда выход valid равняется 1. Количество циклов между тем, когда порт ready 0 и порт выхода valid, равняется 1, всегда latency – FFTLength.
