Частотно-избирательный разветвитель-дециматор
Дипломная работа - Компьютеры, программирование
Другие дипломы по предмету Компьютеры, программирование
? идеального фильтра');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));
figure(7);(w*f_d/(2*pi),P);([0 60000000 min(P) max(P)]);('w');('ФЧХ искомого фильтра');
grid on
% 4. Прохождение сигнала через ЛДС
% 4.1. Через схему переноса спектра_c_f0=zeros(length(f0),K);
X_s_f0=zeros(length(f0),K);(8);c=1:length(f0)
[H_c_f0,w]=freqz(x_c_f0(c,:),1,K,'whole');
[H_s_f0,w]=freqz(x_s_f0(c,:),1,K,'whole');_c_f0(c,:)=H_c_f0';_s_f0(c,:)=H_s_f0';(2,4,c); plot(w*f_d/(2*pi),abs(X_c_f0(c,:)));([0, 60000000, 0 , abs(max(X_c_f0(c,:)))]);
xlabel('w');('Перенесенный спектр действ сост');
grid on(2,4,c+4); plot(w*f_d/(2*pi),abs(X_s_f0(c,:)));([0, 60000000, 0 , abs(max(X_s_f0(c,:)))]);
xlabel('w');('Перенесенный спектр мним сост');on
% 4.2 Формирование сигнала на выходах фильтра
y_c_filtered=zeros(length(f0),length(l)+N-1);_s_filtered=zeros(length(f0),length(l)+N-1);(9);c=1:length(f0) _c_filtered(c,:)=conv(h,x_c_f0(c,:));_s_filtered(c,:)=conv(h,x_s_f0(c,:)); (2,4,c); plot(y_c_filtered(c,:));
axis([0 500 -1.2 1.2]);('n');('Действ сигнал на вых ФНЧ');
grid on(2,4,c+4); plot(y_s_filtered(c,:));
axis([0 500 -1.2 1.2]);('n');('Мним сигнал на выходе ФНЧ');
grid on_c_filtered=zeros(length(f0),K);_s_filtered=zeros(length(f0),K);(10);c=1:length(f0) _c_filtered(c,:)=freqz(y_c_filtered(c,:),1,K,'whole');_s_filtered(c,:)=freqz(y_s_filtered(c,:),1,K,'whole');(2,4,c); plot(w*f_d/(2*pi),abs(Y_c_filtered(c,:)));([0, 60000000, 0 , abs(max(Y_c_filtered(c,:)))]);
xlabel('w(l)+N-1);_s_dmc=zeros(length(f0),length(l)+N-1);(11);c=1:length(f0) =0;b=1:(length(l)+N-1)i==3_c_dmc(c,b)=y_c_filtered(c,b);_s_dmc(c,b)=y_s_filtered(c,b);=0; =i+1;(2,4,c); plot(y_c_dmc(c,:));
axis([0 500 -1.2 1.2]);('n');('Действ сигнал на вых ДЦ');
grid on(2,4,c+4); plot(y_s_dmc(c,:));
axis([0 500 -1.2 1.2]);('n');('Мним сигнал на выходе ДЦ');
grid on_c_dmc=zeros(length(f0),K);_s_dmc=zeros(length(f0),K);(12);c=1:length(f0)_c_dmc(c,:)=freqz(y_c_dmc(c,:),1,K,'whole');_s_dmc(c,:)=freqz(y_s_dmc(c,:),1,K,'whole'); (2,4,c); plot(w*f_d/(2*pi),abs(Y_c_dmc(c,:)));
xlabel('w');('Спектр действ сигн на вых устройства');
grid on(2,4,c+4); plot(w*f_d/(2*pi),abs(Y_s_dmc(c,:)));
xlabel('w');('Спектр мним сигн на вых устройства');on