ep" FOR OUTPUT AS # 2
PRINT "IOY"
INPUT "BB W1, a1, f1, W2, A2, F2", W1, A1, F1, W2, A2, F2 = F1/57.3: F2 = F2/57.3
W = W1: A = A1: F = F1: B = 0: T1 = 1000: N = 2: Z = .0000001
IF ABS (T1 * W)
T1 = T1/2: GOTO 40
IF B = 1 THEN 110
T = (TAN ((-F - T1 * W)/N))/W
K = A * (((T * W) ^ 2 + 1) ^ (N/2))
'PRINT "K ="; K; "T ="; T
W = W2: A = A2: F = F2: B = 1
GOTO 140S9
N = 2 * (LOG (K/A)/LOG ((T * W) ^ 2 + 1))
T1 = (-F - N * ATN (T * W))/W
'PRINT "N ="; N; "T1 ="; T1
W = W1: A = A1: F = F1: B = 0
A0 = K/(((T * W) ^ 2 + 1) ^ (N/2))
F0 =-N * ATN (T * W) - T1 * W
IF ABS ((A - A0)/A)> Z THEN 60
IF ABS ((F - F0)/F)> Z THEN 60
PRINT # 1, "K ="; K; "T ="; T; "T1 ="; T1; "N ="; N = K/ (((T * W1) ^ 2 + 1) ^ (N/2)): Fm1 =-N * ATN (T * W1) - T1 * W1 # 1, "Am1 ="; Am1; "Fm1 = "; Fm1 * 57.3;" A1 = "; A1;" F1 = "; F1 * 57.3 = K/(((T * W2) ^ 2 + 1) ^ (N/2)): Fm2 =-N * ATN (T * W2) - T1 * W2 # 1, "Am2 ="; Am2; "Fm2 ="; Fm2 * 57.3; "A2 ="; A2; "F2 = "; F2 * 57.3
W = 0 = 2 * W2
H = (W1 - W)/15
A0 = K/(((T * W) ^ 2 + 1) ^ (N/2))
F0 =-N * ATN (T * W) - T1 * W
'240 PRINT "W ="; W; "P ="; A0 * COS (F0); "Q ="; A0 * SIN (F0) = A0 * COS (F0) = A0 * SIN (F0) # 2, P, Q
IF W> W1 THEN 280
260 W = W + H
GOTO 220
END
За допомогою програми, були отримані наступні параметри:
K = .9761574, T = 28.74019, T1 = 4.299038, N = 2.073122, Am1 = .57
Fm1 = -88.99999, A1 = .57, F1 = -89, Am2 = .25, Fm2 = -136, A2 = .25, F2 = -136
Графік КЧХ представлений на рісунке1.3.
В
Малюнок 1.3-Комплексно-частотна характеристика об'єкта
ВИСНОВКИ
У процесі виконання курсової роботи проведена ідентифікація об'єкта управління по двох точках експериментальної амплітудно - фазової характеристиці об'єкта.
Отримано, що математична модель об'єкта у вигляді передавальної функції, визначені значення невідомих коефіцієнтів
K = 0.9761574, T = 28.74019, ? = 4.299038, N = 2.073122
і проведена перевірка адекватності моделі та об'єкта шля...
