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