[nC,dC] = tfdata(C, 'v');
sim('lab7sim');
grafiklar qurilishi (etalon i faktik o`tish jarayonlari) :
figure(1);
plot(out(:,1), out(:,2), out(:,1), out(:,3));
CHiqish va boshqarish bo`yicha integral kvadrat xatoliklar hisobidan ( trapetsiya usuli bilan
integrallash, trapz funktsiyasi):
uInf = 0;
% boshqarishning belgilangan qiymati
err2 = (out(:,2)-out(:,3)).^2;
% xatolik kvadrati
errOpt = trapz(out(:,1),err2)
% integral kvadrat xatolik
errU2 = (u(:,2) - uInf).^2;
% boshqarish bo`yicha xatolik
kvadrati
errU = trapz(u(:,1),errU2)
% integral kvadrat xatolik
22.
Modellashtirish vaqtini 100 s ga o`rnating (menyu Simulation-Simulation parameters-
Stop time). Skriptni ishga tushiring va nusxasini hisobotga ko`chiring.
23.
Sistema sifati egri chizig`ini qurish uchun skriptni o`zgartiramiz ( Pareto – optimal
rostlagichlari uchun boshqarish v chiqish bo`yicha integral kvadrat xatoliklar bog`liqligi):
kU = [0.2 0.25 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.2 1.5 1.7 2];
uInf = 0;
% boshqarishning belgilangan qiymati
errOut = [];
% massivlarni ozod qilmoq
errU = [];
for i=1:length(kU)
fprintf('\n%g', kU(i));
% kU ni chiqarish
sys = [ R -P
% standart sistema
0 -kU(i)
R -P];
[C,Poles] = h2opt ( sys );
% Optimal rostlagich sintezi
[nC,dC] = tfdata(C, 'v');
sim('lab7sim');
% modellashtirish
err2 = (out(:,2)-out(:,3)).^2;
u2 = (u(:,2) - uInf).^2;
errOut(i) = trapz(out(:,1),err2);
errU(i) = trapz(u(:,1),u2);
end;
24.
TSikldan keyin sifat egri chizig`ini qurish uchun buyruqlar kiriting
figure(1);
plot(errOut, errU, '-o');
xlabel('Error');
ylabel('Control error');
25.
Skriptni bajarishga ishga tushiring. Hosil qilingan sifat egri chizig`i nusxasini hisobotga
ko`chiring. Xulosalar chiqaring.
Do'stlaringiz bilan baham: |