Примените формирование импульса путем интерполяции сигнала с помощью фильтра конечной импульсной характеристики приподнятого косинуса
The comm.RaisedCosineTransmitFilter Система object™ применяет формирование импульса путем интерполяции входного сигнала с помощью фильтр с приподнятой косинусоидальной конечной импульсной характеристикой (КИХ). Фильтр конечной импульсной характеристики имеет (FilterSpanInSymbols × OutputSamplesPerSymbol + 1) коэффициенты касания.
Для применения формирования импульсов путем интерполяции входного сигнала с помощью фильтр конечной импульсной характеристики приподнятого косинуса:
Создайте comm.RaisedCosineTransmitFilter Объекту и установите его свойства.
Вызывайте объект с аргументами, как будто это функция.
Дополнительные сведения о работе системных объектов см. в разделе «Что такое системные объекты?».
txfilter = comm.RaisedCosineTransmitFilter
txfilter = comm.RaisedCosineTransmitFilter(Name,Value)comm.RaisedCosineTransmitFilter('FilterSpanInSymbols',15) конфигурирует фильтр передачи приподнятого косинуса с Системным объектом диапазона фильтра, установленным на 15 символов.
Чтобы использовать функцию объекта, задайте системный объект в качестве первого входного параметра. Например, чтобы освободить системные ресурсы системного объекта с именем 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)
