Арзамас, 2013 г

Волков В.Л.

Кошелев С.И.

Лабораторная работа

Форма одежды для женщин

Для лиц высшего начальствующего состава

Демисезонная повседневная форма одежды

Летняя повседневная форма одежды

Примечание: образцы форменной одежды, отмеченные желтым цветом находятся на стадии утверждения.

по «Моделированию процессов и систем»

«Формирование и исследование матричных моделей динамической системы»

Вариант №7

Выполнил:

студент группы АСП 11-1

Проверил:

Исходные данные

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 =

C =

0 1.3363 0.0879 0.0042

D =

Схема испытаний непрерывной матричной модели

График выходного процесса матричной модели

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 =

dt =

2.0832e-006

Схема испытаний дискретной непрерывной матричной модели

График выходного процесса матричной модели

Контрольный вопрос: Как проверить правильность вычисленного значения максимальной частоты спектра процесса?

Ответ: Подставить полученное значение максимальной частоты спектра процесса в формулу:

Если значение сойдутся, то максимальная частота спектра процесса вычислена правильно.

Дата добавления: 2015-05-09