Поможем написать учебную работу
Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
Если у вас возникли сложности с курсовой, контрольной, дипломной, рефератом, отчетом по практике, научно-исследовательской и любой другой работой - мы готовы помочь.
ФЕДЕРАЛЬНОЕ АГЕНСТВО ПО ОБРАЗОВАНИЮ
АРЗАМАССКИЙ ПОЛИТЕХНИЧЕСКИЙ ИНСТИТУТ (филиал)
ГОСУДАРСТВЕННОГО ОБРАЗОВАТЕЛЬНОГО УЧРЕЖДЕНИЯ ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ
«НИЖЕГОРОДСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХГИЧЕСКИЙ УНИВЕРСИТЕТ
ИМЕНИ Р.Е. АЛЕКСЕЕВА»
Лабораторная работа
по «Моделированию процессов и систем»
«Формирование и исследование матричных моделей динамической системы»
Вариант №7
Выполнил:
студент группы АСП 11-1
Кошелев С.И.
Проверил:
Волков В.Л.
Арзамас, 2013 г.
Исходные данные
Mu=56, R=1
1)Непрерывная матричная модель
Код программы в Матлабе
clc
Sg=tf([16 0 -0.176 0 0.00016],[0.035 0 -7.94 0 87.61 0 -27.72 0 2.18]);
Wg=fact(Sg)
R=1;
Mu=56;
[b,a]=tfdata(Wg,'v');
sys=ss(Wg);
[A,B,C,D]=ssdata(sys)
W0=max(damp(Wg))
Wm=W0*10^0.75
dt=pi/Wm
sysd=c2d(sys,dt)
[F,G,C,T,dt]=ssdata(sysd)
Результат работы программы
Transfer function:
16 s^4 - 0.176 s^2 + 0.00016
---------------------------------------------------
0.035 s^8 - 7.94 s^6 + 87.61 s^4 - 27.72 s^2 + 2.18
Transfer function:
4 s^2 + 0.5265 s + 0.01265
----------------------------------------------------
0.1871 s^4 + 3.523 s^3 + 11.96 s^2 + 7.939 s + 1.476
A =
-18.8334 -7.9901 -2.6523 -0.9865
8.0000 0 0 0
0 2.0000 0 0
0 0 0.5000 0
B =
2
0
0
0
C =
0 1.3363 0.0879 0.0042
D =
0
Схема испытаний непрерывной матричной модели
График выходного процесса матричной модели
2) Дискретная матричная модель
Код программы в Матлабе
W0=max(damp(Wg))
Wm=W0*10^0.75
dt=pi/Wm
sysd=c2d(sys,dt)
[F,G,C,T,dt]=ssdata(sysd)
W0 =
14.6712
Wm =
1.5080e+006
dt =
2.0832e-006
a =
x1 x2 x3 x4
x1 1 -1.664e-005 -5.525e-006 -2.055e-006
x2 1.667e-005 1 -4.604e-011 -1.713e-011
x3 3.472e-011 4.166e-006 1 -2.378e-017
x4 1.205e-017 2.17e-012 1.042e-006 1
b =
u1
x1 4.166e-006
x2 3.472e-011
x3 4.822e-017
x4 1.256e-023
c =
x1 x2 x3 x4
y1 0 1.336 0.08794 0.004226
d =
u1
y1 0
Sampling time: 2.0832e-006
Discrete-time model.
F =
1.0000 -0.0000 -0.0000 -0.0000
0.0000 1.0000 -0.0000 -0.0000
0.0000 0.0000 1.0000 -0.0000
0.0000 0.0000 0.0000 1.0000
G =
1.0e-005 *
0.4166
0.0000
0.0000
0.0000
C =
0 1.3363 0.0879 0.0042
T =
0
dt =
2.0832e-006
Схема испытаний дискретной непрерывной матричной модели
График выходного процесса матричной модели
Контрольный вопрос: Как проверить правильность вычисленного значения максимальной частоты спектра процесса?
Ответ: Подставить полученное значение максимальной частоты спектра процесса в формулу :
Если значение сойдутся, то максимальная частота спектра процесса вычислена правильно.