Главное меню

Децимация

Уменьшение частоты дискретизации сигнала в m раз (m- коэффициент децимации) осуществляется дециматором, формирующим сигнал путем взятия только каждого m-го отсчета из последовательности . Реализуем этот алгоритм с помощью Matlab. Полученная последовательность представляет собой сигнал на выходе устройства.

Рисунок 16 - Сигналы на выходах дециматоров

В ходе работы при помощи схемы переноса спектра, фильтров низкой частоты и дециматора был спроектирован частотно-избирательный разветвитель-дециматор, который имеет несколько выходов, сигнал на которых появляется в случае, если спектр входного сигнала содержит поддиапазон, сопоставленный этому выходу. Также осуществляется фильтрация сигнала и понижение частоты дискретизации от 120 МГц до 30 МГц.

С помощью среды Matlab был рассчитан ФНЧ, а также была смоделировна работа ФНЧ и дециматора на примере единичного импульса, так как это простой сигнал, имеющий бесконечный спектр.

Для расчета ФНЧ использовалось окно Кайзера, так как оно является параметрическим в отличие от других окон и наиболее оптимальным.

Приложение

Расчет ФНЧ, модель ФНЧ и дециматора, испытанная на пробном сигнала в среде

Matlab

clear all;;=23;=-0.5*(N-1):1:0.5*(N-1);

K=2048;_d= 120000000;

% 1. Формирование пробного сигнала

l=1:1:1000;=35000000;=15000000;= -5000000;= -25000000;=exp(i*2*pi*f1*l/f_d)+exp(i*2*pi*f2*l/f_d)+exp(i*2*pi*f3*l/f_d)+exp(i*2*pi*f4*l/f_d);_c= real(x);_s= imag(x);

[Hx_c,w]=freqz(x_c,1,K,'whole');

[Hx_s,w]=freqz(x_s,1,K,'whole');(1);(2,2,1); plot (x_c);('n');('Действ составл вх сигнала');on(2,2,2); plot(w*f_d/(2*pi),abs(Hx_c));([0 60000000 0 abs(max(Hx_c))]);

xlabel('w');('Спектр действит составл');

grid on(2,2,3); plot(x_s);('n');('Мним составл вх сигнала');on(2,2,4); plot(w*f_d/(2*pi),abs(Hx_s));([0 60000000 0 abs(max(Hx_s))]);

xlabel('w');('Спектр мним составл');on

% 2. Перенос спектра=[-30000000 -10000000 10000000 30000000];

x0=zeros(length(f0),length(l));_c_f0=zeros(length(f0),length(l));_s_f0=zeros(length(f0),length(l));c=1:length(f0)

x0(c,:)=x.*exp(i*2*pi*f0(c)*l/f_d); _c_f0(c,:)= real(x0(c,:));_s_f0(c,:)= imag(x0(c,:));

end

% 3. Расчёт характеристик ФНЧ

% 3.1.Построение ИХ идеального фильтра

h_u=sin(5*pi*n/24)./(pi*n);_u(N/2+0.5)=5/24;(2);(n, h_u);

xlabel('n');('ИХ идеального фильтра');on

% 3.2. Построение ИХ и АЧХ окна Кайзера

w_k=kaiser(N, 0);(3);(n, w_k);('n');('ИХ окна Кайзера');on

[W_k,w]=freqz(w_k,1,K,'whole');=abs(max(W_k));_kdB=20*log10(abs(W_k)/M);(4);(w*f_d/(2*pi), W_kdB );([0 60000000 -90 0]);('w');('АЧХ окна Кайзера, дБ');

grid on

% 3.3. Построение ИХ и АЧХ искомого фильтра

h=h_u.*w_k';(5);(n, h);('n');

ylabel('ИХ искомого фильтра');on

[H,w]=freqz(h,1,K,'whole');=abs(max(H));=20*log10(abs(H)/M);(6);(w*f_d/(2*pi), HdB );([0 60000000 -90 0]);('w');('АЧХ искомого фильтра, дБ');on= unwrap(angle(H));

Перейти на страницу: 1 2

Другое по теме:

Расчет управляемого выпрямителя и СИФУ
В настоящее системы преобразования переменного синусоидального напряжения и тока в постоянные практически полностью представлены полупроводниковыми выпрямителями. Очень часто также необходимо регулировать величину полученного постоянного напряже ...

Copyright © www.techproof.ru