Подскажите как в Matlab пропустить через рассчитанный фильтр такой сигнал: 101011001101001?

close all; clear all; clc;
f_d = 15e7; w_d = 2*pi*f_d; %Частота дискретизации
F_p = 4000000; w_p = 2*pi*F_p; %Частота среза
R_pzad = 0.4; R_szad = 43; %Заданные допуски на ЛАЧХ
F_s = 5400000; w_s = 2*pi*F_s; %Граничная частота
                                     
w_0 = (w_p+w_s)/2; %Частота разрыва (среза) идеального АЧХ
 
delta1 = (10^(0.05*R_pzad - 1))/(10^(0.05*R_pzad + 1));
delta2 = 10^(-0.05*R_szad);
delta = min(delta1,delta2);
R_s = -20*log10(delta);
if R_s <= 21
    alfa = 0;
elseif (R_s > 21)&(R_s <= 50)
    alfa = 0.5842*(R_s-21)^0.4+0.07886*(R_s-21);
elseif R_s > 50
    alfa = 0.1102*(R_s-8.7);
end;
if R_s <= 21
    D = 0.9222;
else
    D = (R_s-7.95)/14.36;
end;
N = ceil(w_d*D/(w_s-w_p)+1);                %Выбор порядка фильтра R=N-1
if mod(N,2) == 0                            %Если N четное прибавляем 1
    N = N+1;
end;
h_imp((N-1)/2+1) = w_0/(f_d*pi);            %Нулевой отсчет БИХ
for n = -(N-1)/2:-1                         %Вычисление фрагмента БИХ
    h_imp((N-1)/2+1+n) = w_0/(f_d*pi)*sin(w_0*n/f_d)*f_d/(w_0*n);
end;
for n = 1:(N-1)/2
    h_imp(n+(N-1)/2+1) = w_0/(f_d*pi)*sin(w_0*n/f_d)*f_d/(w_0*n);
end;
%Построение фрагмента бесконечной импульсной хар-ки
figure('name','Фрагмент бесконечной импульсной х-ки'); 
stem(1:N, h_imp);
grid on;
%Построение импульсной хар-ки реального фильтра
w_n = kaiser(N, alfa);                      
figure('name','Отсчеты функции окна');  
stem(1:N, w_n);
grid on;
ylabel('w(n)');
xlim([0 N+1]);
%Расчет имп. х-ки реального фильтра
for n = 0:N-1                             
    h_n(n+1) = h_imp(n+1)*w_n(n+1);
end;
%Построение импульсной хар-ки реального фильтра
figure('name','Импульсная хар-ка реального фильтра'); 
stem(1:N, h_n);
ylabel('hи(n)');
grid on;
%Построение АЧХ фильтра
figure('name', 'АЧХ фильтра');             
[H, W] = freqz(h_n,[1]);
plot(W*f_d/(2*pi),abs(H),'g');
grid on;
hold on;
%Вычисление БИХ с помощью Matlab-функци
[R_st, W_n_st, alfa_st, ftype] = kaiserord([F_p F_s], [1 0], [delta1 delta2], f_d, 'cell');  
b = fir1(R_st, W_n_st, ftype, kaiser(R_st+1,alfa_st), 'noscale');                         
[H1, W1] = freqz(b, 1); 
%Построение АЧХ фильтра сравнение с ранее полученным результатом
plot(W1*f_d/(2*pi),abs(H1), '--k');
legend('Bez ispolz stand funkcii','S ispolz stand funkcii');
xlabel('f, MHz');
ylabel('Amplitude');
%Построение ЛАЧХ фильтра в dB
figure('name', 'ЛАЧХ фильтра (dB)');       
plot(W*f_d/(2*pi),20*log10(abs(H)), 'b');
hold on;
%построение границы полос пропускания и задержки
plot([F_s, F_s], [-140, -R_s], 'r--');
plot([F_s, max(W*f_d/(2*pi))], [-R_s, -R_s], 'r--');
plot([F_p, F_p], [20, -R_pzad], 'r--');
plot([0, F_p], [-R_pzad, -R_pzad], 'r--');
hold off;
grid on;
xlabel('f, MHz');
ylabel('Amplitude(dB)');
%построение ФЧХ фильтра в градусах
FCH = angle(H);
FCH = unwrap(FCH);   
figure('name', 'ФЧХ фильтра (градус)');    
plot(W*f_d/(2*pi), FCH*180/pi);
grid on;
title('FCH');
xlabel('f, MHz');
ylabel('Fase(gradus)');
%Построение диаграммы нулей
figure('name', 'Диаграмма нулей фильтра');  
zplane(h_n, [1]);
grid on;
xlabel('Real');
ylabel('Imag');
ylim([-2 2]);
%Построение переходной и импульсной х-ки
figure('name', 'Переходная и импульсная характеристики'); 
pereh = filter(h_n, [1], ones(1, 300));
subplot(2,1,1);
stem((0:N-1)./f_d, h_n);
grid on;
title('Impulsnaja');
xlabel('Time(s)');
ylabel('Amplitude');
subplot(2,1,2);
stem((0:length(pereh)-1)./f_d, pereh);
grid on;
title('Perehodnaja');
xlabel('Time(s)');
ylabel('Amplitude');
%Построение реакции фильтра на различные входные сигналы
f1 = 2000000;
f2 = 3800000;
f3 = 6000000;
t_sign = 0:1/f_d:0.000009;
L = length(t_sign);
signal_1 = 10*sin(2*pi*f1*t_sign);
signal_2 = 10*sin(2*pi*f2*t_sign);
signal_3 = 10*sin(2*pi*f3*t_sign);
%--------------------------------------------------------------------------
smp = [1 0 1 0 1 1 0 0 1 1 0 1 0 0 1];
ones_vect = ones(round(length(t_sign)/length(smp)), 1);
sig = ones_vect*smp;
signal_4 = sig(:)';
if length(signal_4)<L
    signal_4(1,end:L) = sig(end);
end
f_nz = 10e6;
noize = 0.5*cos(2*pi*f_nz*t_sign);
signal_5 = signal_4+noize;
%--------------------------------------------------------------------------
vyx1 = filtfilt(h_n, [1], signal_1);
vyx2 = filtfilt(h_n, [1], signal_2);
vyx3 = filtfilt(h_n, [1], signal_3);
%--------------------------------------------------------------------------
vyx4 = filtfilt(h_n, [1], signal_4);
vyx5 = filtfilt(h_n, [1], signal_5);
%--------------------------------------------------------------------------
%% Графики
figure('name', 'Реакция фильтра на различные входные сигналы');
subplot(2,1,1);
plot(t_sign, signal_1, 'r-');
grid on;
title('2 MHz');
subplot(2,1,2);
plot(t_sign, vyx1, 'k-');
grid on;
figure('name', 'Реакция фильтра на различные входные сигналы');
subplot(2,1,1);
plot(t_sign, signal_2, 'r-');
grid on;
title('3,8 MHz');
subplot(2,1,2);
plot(t_sign, vyx2, 'k-');
grid on;
figure('name', 'Реакция фильтра на различные входные сигналы');
subplot(2,1,1);
plot(t_sign, signal_3, 'r-');
grid on;
title('6 MHz');
subplot(2,1,2);
plot(t_sign, vyx3, 'k-');
grid on;
%--------------------------------------------------------------------------
figure('name', 'Реакция фильтра на различные входные сигналы');
subplot(2,1,1);
plot(t_sign, signal_4, 'g-');
grid on;
title('Битовая последовательность');
subplot(2,1,2);
plot(t_sign, vyx4, 'm-');
grid on;
 
figure('name', 'Реакция фильтра на различные входные сигналы');
subplot(2,1,1);
plot(t_sign, signal_5, 'g-');
grid on;
title(['Битовая последовательность+шум (' num2str(f_nz/1e6) ' Мгц)']);
subplot(2,1,2);
plot(t_sign, vyx5, 'm-');
grid on;
%--------------------------------------------------------------------------
[sos,g] = tf2sos(b,[1]);
[r,p,k] = sos2zp(sos);