Mimo model
43
. řád
53 bode(Fs13, 'c'); Fit 13.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
66 grid on
Nejprve byl proveden převod diskrétního modelu, který přirozeně Model Linearizer
identifikuje, spojitý vypsán řád systému. řád
64 impulse(Fs13, 5e7, 'c'); Fit 13. 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.
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.PhaseMatchingFreq 0;
49 opt.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.
Poté byla zahájena samotná redukce modelu. řád
65 impulse(Fs14, 5e7, 'b'); Fit 14. řád
54 bode(Fs14, 'b--'); Fit 14.PhaseMatching 'on';
48 opt. řá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. odstranění pólů nul, které jsou v
blízkosti navzájem vyruší, byla použita funkce minreal definovanou tolerancí
