Strana 50 z 80
Poznámky redaktora
5 sqrt((T1 T2)^2 T2) 1));
44 --- Výpočet fázové bezpečnosti (Phase Margin) ---
45 phase_norm_bez_I' atan(T1*w)/pi*180 atan(T2*w)/pi*180 atan(e*w)/pi*180;
46
47 --- Vykreslení fázové bezpečnosti závislosti frekvenci řezu ---
48 figure
49 semilogx(w, PM)
50 hold on
51 semilogx(w, 45*ones(size(w)), "r") Referenční hodnota 45°
52 grid on
53 title('Fázová bezpečnost závislosti frekvenci řezu')
54
55 --- Určení optimální frekvence řezu (frekvence, kde největší) ---
56 [max_PM, max_PM_idx] max(PM);
57 wrez w(max_PM_idx); Optimální frekvence řezu
58 ampl_G_na_wrez ampl(max_PM_idx); Zesílení systému při této frekvenci
59
60 --- Výpočet zesílení regulátoru ---
61 ampl_G_na_wrez wrez .phase_start phase(1);
18 phase_start_norm mod(phase_start 180, 360) 180;
19 phase_diff phase_start_norm phase_start;
20 phase_norm phase phase_diff;
21
22 --- Vykreslení frekvenčních charakteristik ---
23 figure
24 tiledlayout(2,1)
25
26 nexttile
27 semilogx(w, 20*log10(ampl))
28 grid on
29 title('Amplitudová charakteristika systému')
30
31 nexttile
32 semilogx(w, phase_norm)
33 grid on
34 title('Fázová charakteristika systému')
35
36 --- Úprava fáze -90° pro systém bez integrační složky ---
37 phase_norm_bez_I phase_norm 90;
38
39 z;
40 T1; Obě nuly (T1, T2) jsou stejné
41
42 --- Výpočet realizačního pólu ---
43 (T1 T2)/2 0.
62 sqrt(e^2 wrez^2 ..
63 (sqrt(T1^2 wrez^2 sqrt(T2^2 wrez^2 1));
50
...
