PhaseMatchingValue phase_start_norm;
50 bode(Fs_cont, 'k',opt); Původní systém
51 bode(Fs9, 'm'); Fit řád
52 bode(Fs11, 'r'); Fit 11. řád
65 impulse(Fs14, 5e7, 'b'); Fit 14. řád
64 impulse(Fs13, 5e7, 'c'); Fit 13. řád
66 grid on
Nejprve byl proveden převod diskrétního modelu, který přirozeně Model Linearizer
identifikuje, spojitý vypsán řád systému.30
31 Fitování různými řády přidání integrátoru
32 Fs14 zpk(tfest(sysgE, 14,14-RD,opttfest))/tf('s')
33 Fs13 zpk(tfest(sysgE, 13,13-RD,opttfest))/tf('s')
34 Fs11 zpk(tfest(sysgE, 11,11-RD,opttfest))/tf('s')
35 Fs9 zpk(tfest(sysgE, 9,9-RD,opttfest))/tf('s')
36
37 Zobrazení Bodeho diagramu pro všechny varianty
38 figure;
39 hold on;
40 Stanovení počáteční fáze odezvy
41 [~,phase,~] bode(E_cont,w);
42 phase_start phase(1);
43 phase_start_norm mod(phase(1) 180, 360) 180;
44
45 Posunutí fázového grafu
46 opt bodeoptions;
47 opt.PhaseMatching 'on';
48 opt. řád
54 bode(Fs14, 'b--'); Fit 14.
Tímto způsobem bylo dosaženo snížení řádu bez ovlivnění dynamiky systému, což
bylo ověřeno porovnáním amplitudových fázových charakteristik. Mimo model
43
.PhaseMatchingFreq 0;
49 opt. odstranění pólů nul, které jsou v
blízkosti navzájem vyruší, byla použita funkce minreal definovanou tolerancí. řád
53 bode(Fs13, 'c'); Fit 13.
Poté byla zahájena samotná redukce modelu. řád
55 legend('(originál)', 'Fs9', 'Fs11', 'Fs13', 'Fs14');
56 grid on;
57 xlim([wstart wstop])
58 hold off;
59 figure;
60 impulse(Fs_cont, 5e7, 'k'); Původní systém
61 hold on
62 impulse(Fs9, 5e7, 'm'); Fit řád
63 impulse(Fs11, 5e7, 'r'); Fit 11. Tento model následně převeden do
tvaru přenosové funkce typu ZPK (zero-pole-gain), tedy zápisu pomocí nul, pólů a
zesílení jeho chování porovnáno původním modelem pomocí Bodeho diagramu
