Уменьшение частоты дискретизации сигнала 
в 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);
=freqz(x_c,1,K,'whole');
=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));();
xlabel('w');('Спектр действит составл');
grid on(2,2,3); plot(x_s);('n');('Мним составл вх сигнала');on(2,2,4); plot(w*f_d/(2*pi),abs(Hx_s));();
xlabel('w');('Спектр мним составл');on
% 2. Перенос спектра=;
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
=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 );();('w');('АЧХ окна Кайзера, дБ');
grid on
% 3.3. Построение ИХ и АЧХ искомого фильтра
h=h_u.*w_k';(5);(n, h);('n');
ylabel('ИХ искомого фильтра');on
=freqz(h,1,K,'whole');=abs(max(H));=20*log10(abs(H)/M);(6);(w*f_d/(2*pi), HdB );();('w');('АЧХ искомого фильтра, дБ');on= unwrap(angle(H)); Перейти на страницу: 1 2
Другое по теме:
Ситсема Секам моделирование видеосигнал телевидение секам Приведение содержания подготовки современных специалистов в соответствии с требованиями времени и достигнутым уровнем развития техники особенно остро стоит перед техническими учебными заведениями из-за ...