G_D_licznik=1;
G_D_mianownik=[4,1];
G_D=tf(G_D_licznik, G_D_mianownik); %filtr dolnoprzepustowy

G_G_licznik=[0.4,0];
G_G_mianowni=[0.4,1];
G_G=tf(G_G_licznik, G_G_mianownik); %filtr górnoprzepustowy

G_P_licznik=[1,0];
G_P_mianownik=[1,2,1];
G_P=tf(G_P_licznik, G_P_mianownik); %filtr pasmoprzepustowy

G_2_licznik=1;
G_2_mianownik=[1,0,1,1];
G_2=tf(G_2_licznik, G_2_mianownik); %filtr drugiego rzędu

G_liczniki=[G_D_licznik, G_G_licznik, G_P_licznik, G_2_licznik];
G_mianowniki=[G_D_mianownik, G_G_mianownik, G_P_mianownik, G_2_mianownik];
G=[G_D, G_G, G_P, G_2];
G_titles={'Filtr dolnoprzepustowy'
    'Filtr górnoprzepustowy'
    'Filtr pasmoprzepustowy'
    'Filtr drugiego rzędu'};

A_values=[2, 3, 1, 2, 4];
omega_values=[0.001, 1, 30, 0.01, 5];
fi_values=[1, 5, 2, 1, 0];

%charakterystyki filtrów
for i=1:4
    bode(G(i));
    title(G_titles(i));
end;

%dla każdego sygnału
for i=1:5
    figure(i);
    A=A_values(i);
    omega=omega_values(i);
    fi=fi_values(i);
    switch i
        case 1
            A1=A;
            omega1=omega;
            fi1=fi;
        case 2
            A2=A;
            omega2=omega;
            fi2=fi;
        case 3
            A3=A;
            omega3=omega;
            fi3=fi;
        case 4
            A4=A;
            omega4=omega;
            fi4=fi;
        case 5
            A5=A;
            omega5=omega;
            fi5=fi;
    end;
    hold on;
    %dla każdego filtra
    for j=i:4
        G_licznik=G_liczniki(j);
        G_mianownik=G_mianowniki(j);
        sim('model3');
        plot(y);
    end;
    hold off;
end;
sim('model3sum');
hold on;
%dla każdego filtra
for j=i:4
    G_licznik=G_liczniki(j);
    G_mianownik=G_mianowniki(j);
    sim('model3');
    plot(y1);
end;
hold off;