Курсовые / kurstay

% Курсовая работа по ТАУ
%Pn-мощность
%Un-номинальное напряжение
%In-номинальный ток
%R-сопротиление,L-индуктивность
%J-момент инерции
%n-номинальная угловая скорость
----------------------------------------------------------------------------------
L=360*10^(-3);
R=74;
Pn=20;
Un=110;
In=0.5;
J=3.4*0.01^2;
Kp=1;
M=0.2;
n=1000;
w=2*pi*n/60;
----------------------------------------------------------------------------------
%нахождение передаточной функции
K1=(Un-In*R)/w
K2=M/In
T1=L*J/(K1*K2)
T2=R*J/(K1*K2)
K=1/K1
----------------------------------------------------------------------------------
%нахождение Tp
v1=-T2/(2*T1)+sqrt((T2/T1)^2-4/T1)/2
v2=-T2/(2*T1)-sqrt((T2/T1)^2-4/T1)/2
a=1/(v1*v2*K1*T1)
c=1/((v2-v1)*v2*K1*T1)
b=-a-c
%Tp=0.3
----------------------------------------------------------------------------------
%уравнения состояния непрерывного объекта
A=[0 Kp 0;0 0 (K2/J);0 (-K1/L) (-R/L)]
B=[0;0;(1/L)]
C=[1 0 0]
Su=[B A*B A*A*B]
Sh=[C;C*A;C*A*A]
----------------------------------------------------------------------------------
%уравнения состояния дискретной модели объекта
Ad=[1 (0.1*Tp*Kp) (0.005*(Tp^2)*K2*Kp/J);
0 (1-0.005*(Tp^2)*K1*K2/(J*L)) (0.1*Tp*K2/J-0.005*(Tp^2)*K2*R/(J*L));
0 (0.005*(Tp^2)*K1*R/(L^2)-0.1*Tp*K1/L) (1-0.1*Tp*R/L+0.005*(Tp^2)*(R^2/L^2-K1*K2/(J*L)))]
Bd=[0 ;(0.005*(Tp^2)*K2/(J*L)) ;(0.1*Tp-0.005*(Tp^2)*R/L)/L]
----------------------------------------------------------------------------------
AdBd=Ad*Bd;
Ad2Bd=Ad*Ad*Bd;
Sco=[Bd AdBd Ad2Bd]
R=[0 0 1]*inv(Sco)*(Ad*Ad*Ad)
Sob=[C;C*Ad;C*(Ad*Ad)]
H=(A*A*A)*inv(Sob)*[0;0;1]
-----------------------------------------------------------------------------------------
Ad-Bd*R
H*C-Bd*R
Ad-H*C