Курсовая работа: Синтез системы автоматического управления непрерывным объектом
Таким образом, рассчитать компенсатор можно следующим образом:
Рассчитаем дискретный компенсатор с помощью Matlab’а.
|
W1=tf([0.9],[20 1],'td', 1) % задаем передаточную функцию W2=tf([1],[500 100 1],'td', 15) % задаем передаточную функцию Wf=tf([0.7],[10 1]) % задаем передаточную функцию Wob=W1*W2 % общая передаточная двух последовательных частей системы T=1 % время квантования Wdiskr=c2d(Wob,T,'zoh') % передаточная в дискретной области W1d=c2d(W1,T,'zoh') % передаточная в дискретной области W2d=c2d(W2,T,'zoh') % передаточная в дискретной области Wfd=c2d(Wf,T,'zoh') % передаточная в дискретной области [NumerDenom]=tfdata(Wdiskr, 'v') % находим числитель и знаменатель m=length (Numer) Denom1=Denom(2:m) Numer1=Numer(2:m) q0=1/sum(Numer1) for i=1:(m-1) q(i)=q0*Denom1(i) |
|
p(i)=q0*Numer1(i) end Q=[q0 q] % матрица коэффициентов числителя P=[1 -p] % матрица коэффициентов знаменателя Wr=tf(Q, P, T) % передаточная функция регулятора Wkomp=(Wfd)/(Wr*W1d) % передаточная функция компенсатора [NkDk]=tfdata(Wkomp, 'v') % коэффициенты числителя и знаменателя [NfDf]=tfdata(Wfd, 'v') % коэффициенты числителя и знаменателя [N1 D1]=tfdata(W1d, 'v') % коэффициенты числителя и знаменателя [N2 D2]=tfdata(W2d, 'v') % коэффициенты числителя и знаменателя |