Курсовая работа: Анализ системы автоматического регулирования угловой скорости вращения турбины

то коэффициенты , , , … называют коэффициентами ошибок. Их можно определить по известным формулам

Затем величину ошибки можно рассчитать по формуле

В задании требуется определить ошибку при x_ВХ (t) = 1; t; t² .

Программа расчёта величины ошибки:

disp (' *** Определение коэффициентов ошибок ***');

syms p Wd x dx t % p, Wd, x, dx - символьные переменные

a(2)=252.5*1.0e+8;

disp (' Передаточная функция по ошибке ');

Wd=(a(7)*p^6+a(6)*p^5+a(5)*p^4+a(4)*p^3+a(3)*p^2+a(2)*p+a(1))/(a(7)*p^6+a(6)*p^5+a(5)*p^4+

+a(4)*p^3+a(3)*p^2+(a(2)+0.075*1.0e+8)*p+a(1));

pretty(Wd) % выводвудобочитаемомвиде

disp (' Коэффициенты ошибок ');

S0=subs(diff(Wd,p,0),p,0)

S1=subs(diff(Wd,p,1),p,0)

S2=subs(diff(Wd,p,2),p,0)

S3=subs(diff(Wd,p,3),p,0)

S4=subs(diff(Wd,p,4),p,0)

S5=subs(diff(Wd,p,5),p,0)

S6=subs(diff(Wd,p,6),p,0)

disp (' ');

disp ('Определение ошибки при различных функциях входного сигнала');

x=1

dx=eval(S0*x+S1*diff(x,t)+S2*diff(x,t,2)+S3*diff(x,t,3)+S4*diff(x,t,4)+S5*diff(x,t,5)+S6*diff(x,t,6))

%pretty(dx)

x=t

dx=eval(S0*x+S1*diff(x,t)+S2*diff(x,t,2)+S3*diff(x,t,3)+S4*diff(x,t,4)+S5*diff(x,t,5)+S6*diff(x,t,6));

pretty(dx)