Применить формирование импульсов путем интерполяции сигнала с помощью КИХ-фильтра повышенной косинусности
comm.RaisedCosineTransmitFilter Система object™ применяет формирование импульсов путем интерполяции входного сигнала с использованием фильтра с увеличенной косинусной конечной импульсной характеристикой (КИХ). Фильтр FIR имеет (FilterSpanInSymbols × OutputSamplesPerSymbol + 1) коэффициенты отводов.
Для применения формирования импульса путем интерполяции входного сигнала с использованием КИХ-фильтра с увеличенным косинусом:
Создать comm.RaisedCosineTransmitFilter и задайте его свойства.
Вызовите объект с аргументами, как если бы это была функция.
Дополнительные сведения о работе системных объектов см. в разделе Что такое системные объекты?.
txfilter = comm.RaisedCosineTransmitFilter
txfilter = comm.RaisedCosineTransmitFilter(Name,Value)comm.RaisedCosineTransmitFilter('FilterSpanInSymbols',15) конфигурирует восходящий объект системы косинусного фильтра передачи с диапазоном фильтра, равным 15 символам.
Чтобы использовать функцию объекта, укажите объект System в качестве первого входного аргумента. Например, для освобождения системных ресурсов объекта System с именем obj, используйте следующий синтаксис:
release(obj)
Интерполяция сигнала с использованием объекта фильтра передачи с возведением квадратного корня в косинус (SRRC) и отображение спектра отфильтрованного сигнала.
Создание случайных биполярных символов со скоростью 1e6 символов в секунду.
data = 2*randi([0 1],1e6,1) - 1;
Создайте объект фильтра передачи SRRC. По умолчанию фильтру присваивается форма квадратного корня, а число выборок на символ - 8.
txfilter = comm.RaisedCosineTransmitFilter
txfilter =
comm.RaisedCosineTransmitFilter with properties:
Shape: 'Square root'
RolloffFactor: 0.2000
FilterSpanInSymbols: 10
OutputSamplesPerSymbol: 8
Gain: 1
Фильтрация данных с помощью фильтра SRRC.
filteredData = txfilter(data);
Создайте объект анализатора спектра с частотой выборки 8e6. Эта частота дискретизации соответствует частоте дискретизации отфильтрованного сигнала.
spectrumAnalyzer = dsp.SpectrumAnalyzer('SampleRate',8e6);Просмотрите спектр отфильтрованного сигнала с помощью объекта анализатора спектра.
spectrumAnalyzer(filteredData)
Создание интерполированных сигналов из фильтра с квадратным корнем и увеличенным косинусом (SRRC) с различными диапазонами фильтра. Проверьте амплитудную характеристику различных конструкций фильтров.
Создание объектов фильтра SRRC, задающих различные области фильтра. Используйте coeffs объектная функция для получения коэффициентов фильтра.
txfilt2 = comm.RaisedCosineTransmitFilter('FilterSpanInSymbols',2); txfilt4 = comm.RaisedCosineTransmitFilter('FilterSpanInSymbols',4); txfilt6 = comm.RaisedCosineTransmitFilter('FilterSpanInSymbols',6); txfilt8 = comm.RaisedCosineTransmitFilter('FilterSpanInSymbols',8); txfilt16 = comm.RaisedCosineTransmitFilter('FilterSpanInSymbols',16); taps2 = coeffs(txfilt2).Numerator; taps4 = coeffs(txfilt4).Numerator; taps6 = coeffs(txfilt6).Numerator; taps8 = coeffs(txfilt8).Numerator; taps16 = coeffs(txfilt16).Numerator;
Запустите инструмент визуализации фильтра для отображения импульсной характеристики. Укажите частоту дискретизации 1 кГц. Отображение двухсторонней центрированной реакции.
h = fvtool(taps2,1,taps4,1,taps8,1,taps16,1); h.Fs = 1e3; h.FrequencyRange = '[-Fs/2, Fs/2)'; legend('Span 2 symbols','Span 4 symbols','Span 8 symbols','Span 10 symbols') title('Magnitude Response (dB) for Various Filter Spans')
Создайте объект фильтра передачи с квадратным корнем и увеличенным косинусом (SRRC). Используйте FVTool для построения графика отклика фильтра. Результаты показывают, что линейный коэффициент усиления фильтра больше единицы. В частности, коэффициент усиления полосы пропускания составляет более 0 дБ.
txfilter = comm.RaisedCosineTransmitFilter; fvtool(txfilter)
Используйте coeffs объектная функция для получения коэффициентов фильтра и регулировки коэффициента усиления фильтра в соответствии с единичной энергией.
b = coeffs(txfilter);
Поскольку фильтр с единичным коэффициентом усиления полосы пропускания должен иметь коэффициенты фильтра, которые суммируются в 1, установите линейный коэффициент усиления фильтра в обратную сторону от суммы коэффициентов отвода фильтра, b.Numerator.
txfilter.Gain = 1/sum(b.Numerator);
Убедитесь, что результирующие коэффициенты фильтра равны 1.
bNorm = coeffs(txfilter); sum(bNorm.Numerator)
ans = 1.0000
Постройте график частотной характеристики фильтра. Результаты теперь показывают, что коэффициент усиления полосы пропускания равен 0 дБ, что является единичным коэффициентом усиления.
fvtool(txfilter)
