Заключение

Оптимальным коэффициентом передачи линейного звена является такое значение, которому соответствует наименьшее значение параметра M. С настоящими исходными данными минимума для функции М нет, так как зависимость монотонно возрастающая.

СПИСОК ИСПОЛЬЗОВАННЫХ ИСТОЧНИКОВ

1. Бесекерский В.А., Попов Е.П. Теория систем автоматического управления. - СПб: Профессия, 2003.

2. Емельянов В. Ю. Методы моделирования стохастических систем управления: учеб. пособие для ВУЗов. – СПб: БГТУ, 2004.

3. Андриевский Б.Р., Емельянов В.Ю., Коротков Б.Ф. Теория управления: Лабораторный практикум в среде MATLAB/SIMULINK. – СПб: БГТУ, 2001.

Приложение а

function kurs;

t1=0.5; t2=3; g=1; a=5; d=10;

mx=1; sigx=1;

k0=2; k1=1; e=1;

i=1; k=0.5

while k<=20

mx=1; sigx=1;

k0=2; k1=1; e=1; mx1=g/(1+k0*k);

while e>0.0001

mx1=g/(1+k0*k); sigx1=sqrt(k*k*k1*k1*2*a*d*(t1*t2*(a*t1*t2+t1+t2))/2/t1/t2/(a*(1+k1*k)*((a*t1*t2+t1+t2)*(a*t1+a*t2+1+k1*k)-t1*t2*(a+a*k1*k))));

k01=2/mx*erf(mx/sigx/sqrt(2))/2;

kk1=1/sigx*sqrt((1-4*erf(mx/sigx/sqrt(2))/2*erf(mx/sigx/sqrt(2))/2));

kk2=2/(sigx*sqrt(2*3.14))*exp(-(mx*mx/2/sigx/sigx));

k11=(kk1+kk2)/2;

e=abs(mx-mx1)+abs(sigx-sigx1)+abs(k0-k01)+abs(k1-k11);

mx=mx1;

sigx=sigx1;

k0=k01;

k1=k11;

end;

mas1(i)=k0;

akx1(i)=k;

mas2(i)=k1;

m = fzero(@(m) p(m,mx,sigx),1.0);

mx2(i)=mx;

sigx2(i)=sigx;

m1(i)=m;

i=i+1;

k=k+0.1

end;

sigx=sigx1

plot(akx1,m1,'.'); grid;