Witam! Mam problem z filtrem polifazowym decymacyjnym.. Model ten oparłem na modelu komutatorowym, który zamieszczony jest w załączniku.. Chodzi mi o to, w jaki sposób przedstawić poniższy zestaw procedur, by odnosił się do decymacji o czynnik M, a nie tylko o czynnik 4.. Z góry dzękuję za pomoc!Pozdrawiam!
(h- wsólczynniki wcześniej zadeklarowanego filtru)
4- współczynnik decymacji
% Polyphase decomposition
P = reshape(h,4,length(h)/4);
% Polyphase downsampling and filtering
x0 = x(1:4:length(x));
x1 = [0,x(4:4:length(x))];
x2 = [0,x(3:4:length(x))];
x3 = [0,x(2:4:length(x))];
y0 = filter(P(1,: ),1,x0);
y1 = filter(P(2,: ),1,x1);
y2 = filter(P(3,: ),1,x2);
y3 = filter(P(4,: ),1,x3);
sygdec =y0+y1+y2+y3; % Decimated signal
[Sygdec,F2] = freqz(sygdec,1,512);
subplot(4,1,3), plot(F2/pi,abs(Sygdec), 'r'), legend('Widmo sygnału zdecymowanego'), ylabel('|Y_d_e_c(F)|'),
axis([0,1,0,60])
(h- wsólczynniki wcześniej zadeklarowanego filtru)
4- współczynnik decymacji
% Polyphase decomposition
P = reshape(h,4,length(h)/4);
% Polyphase downsampling and filtering
x0 = x(1:4:length(x));
x1 = [0,x(4:4:length(x))];
x2 = [0,x(3:4:length(x))];
x3 = [0,x(2:4:length(x))];
y0 = filter(P(1,: ),1,x0);
y1 = filter(P(2,: ),1,x1);
y2 = filter(P(3,: ),1,x2);
y3 = filter(P(4,: ),1,x3);
sygdec =y0+y1+y2+y3; % Decimated signal
[Sygdec,F2] = freqz(sygdec,1,512);
subplot(4,1,3), plot(F2/pi,abs(Sygdec), 'r'), legend('Widmo sygnału zdecymowanego'), ylabel('|Y_d_e_c(F)|'),
axis([0,1,0,60])