[1351] | 1 | function [ztrata1 ztrata2 ztrata3]=srovnej(system,apriori)
|
---|
| 2 | kresli=0;
|
---|
| 3 | p=0;
|
---|
| 4 | seed=randn(6,system.horizont_rizeni);
|
---|
| 5 | [ztrata1 omega1 rizeni1]=jednoduche_rizeni(system,apriori,seed);
|
---|
| 6 | [ztrata2 omega2 rizeni2]=dualni_rizeni(system,apriori,seed);
|
---|
| 7 | [ztrata3 omega3 rizeni3]=nic(system,apriori,seed);
|
---|
| 8 | if (kresli)
|
---|
| 9 | figure
|
---|
| 10 | plot(omega1(1,:),'r-'); hold on
|
---|
| 11 | plot(omega1(2,:),'r.-.');
|
---|
| 12 | plot(omega2(1,:),'b-');
|
---|
| 13 | plot(omega2(2,:),'b.-.');
|
---|
| 14 | legend('SKUTECNE_simple','ODHAD_simle','SKUTECNE_dual','ODHAD_dual');
|
---|
| 15 | TITLE([ztrata1;ztrata2;ztrata3]);
|
---|
| 16 | figure
|
---|
| 17 | if (p)
|
---|
| 18 | plot(rizeni1(1,:),'r-'); hold on
|
---|
| 19 | plot(rizeni1(2,:),'r.-.');
|
---|
| 20 | plot(rizeni2(1,:),'b-');
|
---|
| 21 | plot(rizeni2(2,:),'b.-.');
|
---|
| 22 | legend('simple_1','simple_1','dual_1','dual_1');
|
---|
| 23 | else
|
---|
| 24 | plot(rizeni1(1,:).^2+rizeni1(2,:).^2-(rizeni2(1,:).^2+rizeni2(2,:).^2),'b-');
|
---|
| 25 | end
|
---|
| 26 | end
|
---|
| 27 | end |
---|