1. czas = 0:0.1:15;
2. k_t = [0.7, 1.6];
3. T_t = [1.3, 0.9];
4. kolor = ['-g', '-r'];
5. tlumienie1 = 3;
6. tlumienie2 = 0.3;
7.  
8. for i=1:2
9.     k = k_t(i);
10.     T = T_t(i);
11.    
12.     licznik = [0, k];
13.     mianownik = [T, 1];
14.    
15.     [y_s, x_s, czas_s] = step(licznik, mianownik, czas);
16.     [y_i, x_i, czas_i] = impulse(licznik, mianownik, czas);
17.    
18.     figure(1)
19.     subplot(1, 2, 1)
20.     plot(czas_s, y_s, kolor(i))
21.     hold on
22.    
23.     subplot(1, 2, 2)
24.     plot(czas_i, y_i, kolor(i))
25.     hold on
26. end
27.  
28. for i=1:2
29.     k = k_t(i);
30.     T1 = T_t(i);
31.     T2 = T1 + 1;
32.    
33.     licznik = [0, 0, k];
34.     mianownik = [T1*T2, T1+T2, 1];
35.    
36.     [y_s, x_s, czas_s] = step(licznik, mianownik, czas);
37.     [y_i, x_i, czas_i] = impulse(licznik, mianownik, czas);
38.    
39.     figure(2)
40.     subplot(1, 2, 1)
41.     plot(czas_s, y_s, kolor(i))
42.     hold on
43.    
44.     subplot(1, 2, 2)
45.     plot(czas_i, y_i, kolor(i))
46.     hold on
47. end
48.  
49. for i=1:2
50.     k = k_t(i);
51.     T = T_t(i);
52.    
53.     licznik = [0, 0, k];
54.     mianownik = [T^2, 2*tlumienie1*T, 1];
55.    
56.     [y_s, x_s, czas_s] = step(licznik, mianownik, czas);
57.     [y_i, x_i, czas_i] = impulse(licznik, mianownik, czas);
58.    
59.     figure(3)
60.     subplot(1, 2, 1)
61.     plot(czas_s, y_s, kolor(i))
62.     hold on
63.    
64.     subplot(1, 2, 2)
65.     plot(czas_i, y_i, kolor(i))
66.     hold on
67. end
68.  
69. for i=1:2
70.     k = k_t(i);
71.     T = T_t(i);
72.    
73.     licznik = [0, 0, k];
74.     mianownik = [T^2, 2*tlumienie2*T, 1];
75.    
76.     [y_s, x_s, czas_s] = step(licznik, mianownik, czas);
77.     [y_i, x_i, czas_i] = impulse(licznik, mianownik, czas);
78.    
79.     figure(4)
80.     subplot(1, 2, 1)
81.     plot(czas_s, y_s, kolor(i))
82.     hold on
83.    
84.     subplot(1, 2, 2)
85.     plot(czas_i, y_i, kolor(i))
86.     hold on
87. end