Операторное уравнение по передаточной функции

Операторные передаточные функции

Содержание:

Операторные передаточные функции:

Практический смысл и назначение операторного метода в теории электрических цепей состоит, прежде всего, в представлении соотношения вход/выход в операторной форме, что даёт возможность существенно упростить процедуры анализа и синтеза электрических цепей и обеспечить связь между временным и частотным описаниями как колебаний, действующих в цепи, так и самой цепи.

Законы Ома и Кирхгофа в операторной форме

Покажем, что решение задач анализа колебаний в электрической цепи существенно упрощается при использовании операторного метода.

Законы Кирхгофа в операторной форме

Пусть токи

(17.1)

что говорит о формальной справедливости законов Кирхгофа для токов и напряжений, выраженных в операторной форме.

Операторные сопротивления и проводимости элементов электрических цепей

Убедимся в справедливости закона Ома для L-изображений колебаний на зажимах элементов R, L, С при нулевых начальных условиях (см.разд. 15.2) и найдём операторные изображения активного сопротивления, реактивного сопротивления индуктивности и ёмкости, а также их операторные проводимости

Для элемента активного сопротивления

(17.2)

т. е. операторное активное сопротивление равно самому активному сопротивлению, поэтому операторная активная проводимость равна самой активной проводимости

(17.3)

Для элемента индуктивности

правило дифференцирования даёт:

откуда операторные сопротивление и проводимость индуктивности равны:

(17.4)

(17.5)

Для элемента ёмкости

правило интегрирования даёт:

откуда операторные сопротивление и проводимость ёмкости равны:

(17.6)

(17.7)

Заметим, что поскольку оператор р согласно (16.2) определён как комплексное переменное

операторные сопротивления и проводимости элементов L и С получаются заменой оператора на оператор р при

Операторные сопротивление и проводимость последовательного и параллельного двухполюсников

Закон Ома при нулевых начальных условиях формально верен и для сложных двухполюсников, если в числе их элементов не содержатся независимые источники.

Определение:

Операторным сопротивлением (проводимостью ) двухполюсника называется отношение операторного напряжения на входе (операторного входного тока к операторному току на выходе (операторному напряжению на выходе)

(соответственно при нулевых начальных условиях.

Найти операторное сопротивление двухполюсника (рис. 17.1), состоящего из последовательно соединённых элементов R, L, С при нулевых начальных условиях.

Решение. Напряжение на зажимах двухполюсника при нулевых начальных условиях равно

Применим к полученному уравнению преобразование Лапласа:

откуда следует, что при последовательном соединении элементов их операторные сопротивления складываются, как и для комплексных сопротивлений, но оператор заменяется на оператор р (см. разд. 17.1.2):

(17.8)

Пример 17.2.

Найти операторную проводимость двухполюсника (рис. 17.2), состоящего из параллельно соединённых элементов R, L, С при нулевых начальных условиях.

Решение. Для тока согласно первому закону Кирхгофа имеем:

поэтому операторную проводимость заданного двухполюсника можно записать сразу:

(17.9)

В силу дуальности последовательного и параллельного контуров выражение (17.9) можно было записать сразу на основании формулы (17.8).

Выражения (17.8) и (17.9) представляют собой входные операторные функции цепи. Они дают основания определению операторного сопротивления и проводимости двухполюсника общего вида.

Операторные сопротивление и проводимость двухполюсника общего вида

Закон Ома, при нулевых начальных условиях, формально можно применить и для сколь угодно сложных двухполюсников. Ранее
(см. лекцию 5) было установлено, что если на входе двухполюсника действует источник напряжения с ЭДС то для контура (например, первого), замыкающегося через этот источник, по формуле Крамера можно записать:

Переходя к L-изображениям напряжений, токов и сопротивлений элементов цепи, получим представление двухполюсника в операторной форме (рис. 17.3), что позволяет записать L-изображение входного тока:

Теперь согласно определению операторной проводимости и операторного сопротивления имеем:

(17.10)

(17.11)

При этом нужно помнить, что определители и алгебраические дополнения в таких формулах записываются с учётом свойств преобразования Лапласа, как это сделано в примерах 17.1 и 17.2.

Определение операторной передаточной функции. Связь с импульсной и переходной характеристиками

В лекции 15 было показано, что во временной области соотношение вход/выход линейной электрической цепи при произвольном воздействии описывается уравнением свёртки:

где h(t) — импульсная характеристика, x(t) — воздействие, y(t) — реакция. При этом воздействие и реакция могут быть напряжениями или токами.

Для описания соотношения вход/выход в операторной форме воспользуемся L-изображением свёртки

(17.12)

откуда получаем соотношения вход/выход в операторной форме

(17.3)

которое называют передаточной функцией.

Определение:

Передаточной функцией линейной электрической цепи называется отношение L-изображения реакции к L-изображению воздействия при нулевых начальных условиях.

Выражение (17.13) говорит о том, что передаточная функция является L-изображением импульсной характеристики, т. е. импульсная характеристика является обратным преобразованием Лапласа передаточной функции:

(17.14)

(17.15)

Именно этими зависимостями объясняется содержащееся в определении передаточной функции требование нулевых начальных условий.

Связь между передаточной функцией и переходной характеристикой можно установить, если воспользоваться интегралом Дюамеля (15.20а) при нулевых начальных условиях:

когда Здесь, как и в случае импульсной характеристики, имеет место свёртка двух функций, которой в операторной области соответствует произведение L-изображений свёртываемых функций:

Первый сомножитель правой части полученного уравнения содержит L-изображение производной, поэтому окончательно можно записать:

(7.16)

(7.17)

что полностью соответствует связи импульсной и переходной характеристик (15.16).

Обратим внимание на то, что передаточная функция может быть получена из комплексных частотных характеристик формальным образом, а именно — простой заменой в КЧХ оператора на и наоборот: КЧХ может быть получена из передаточной функции заменой оператора на

В зависимости от того, какая величина выступает в качестве внешнего воздействия, а какая в качестве реакции цепи, различают четыре вида передаточных функций:

операторное передаточное сопротивление

(17.18)

операторную передаточную проводимость

(17.19)

передаточную функцию по току

(17.20)

передаточную функцию по напряжению

(7.21)

Последние две функции иногда называют операторными передаточными коэффициентами по току и по напряжению соответственно.

По любой из передаточных функций (17.18)—(17.21) нетрудно найти L-изображение реакции цепи, а затем и саму реакцию на заданное воздействие, поскольку любая передаточная функция Н(р) согласно (17.12) может рассматриваться как связующий коэффициент между L-изображения ми воздействия Х(р) и реакции Y(p).

Пример 17.3.

Записать передаточную функцию для последовательного колебательного контура (рис. 17.1, б) относительно напряжения на индуктивности.

Решение. По определению передаточной функции для индуктивности имеем

Но операторное напряжение на индуктивности равно:

Подставляя сюда операторное сопротивление (17.8), получаем искомую передаточную функцию:

(17.22)

Аналогично можно получить и другие передаточные функции для последовательной, параллельной или более сложной цепи. В последнем случае потребуется составить систему уравнений для L-изображений колебаний, воспользовавшись методом контурных токов или узловых напряжений.

Понятие о нулях и полюсах передаточной функции. Устойчивость передаточной функции

Задача 17.1.

Получить и исследовать общее выражение для передаточной функции цепи, когда воздействие представляет собой ЭДС источника напряжения, а реакцией является ток в выделенной ветви анализируемой цепи (рис. 17.4).

Решение. Выберем независимые контуры в цепи так, чтобы через источник напряжения замыкался ток только одного входного контура, а через интересующую нас ветвь — ток только одного выходного контура. На рис. 17.4 они обозначены индексами 1 и 2 соответственно.

Теперь, как и в задаче 5.2, необходимо положить При этих условия соответствующие операторные напряжения также оказываются равными нулю. Тогда операторный ток выходного контура получает вид:

откуда по определению передаточной функции имеем операторную передаточную проводимость

(17.23)

где — определитель системы операторных уравнений

a — операторный минор этого определителя относительно первой строки и второго столбца:

Заметим, что определитель и все его миноры представляют собой рациональные функции оператора р, все коэффициенты которых являются вещественными числами. Это объясняется тем, что при раскрытии определителя над его элементами совершаются только операции умножения, сложения и вычитания, а сами элементы представляют собой простейшие рациональные

функции с вещественными коэффициентами вида (17.11). Раскрывая определитель и минор и подставляя результаты в (17.23), получаем:

(17.24)

Полиномы числителя и знаменателя как и всякий полином, согласно основной теореме алгебры, могут быть представлены через их нули и соответственно следующим образом:

(17.25)

(17.26)

Отсюда передаточная функция (17.24) приобретает вид:

(П27)

—постоянный множитель;

— являются нулями числителя (корнями уравнения ) и называются нулями передаточной функции;

— называется характеристическим полиномом;

— являются нулями характеристического полинома (корнями уравнения ) и называются полюсами передаточной функции.

Названия корней уравнения нулями и корней уравнения полюсами связаны с тем, что при передаточная функция обращается в нуль, а при в бесконечность. Поскольку коэффициенты передаточной функции вещественны, то нули и полюсы могут быть или вещественными или составлять комплексно-сопряжённые пары:

Нули и полюсы наглядно отображаются на комплексной p -плоскости (рис. 17.5) значками ( ° ) и ( * ) соответственно.

На рис. 17.5 показаны:

  • вещественный положительный нуль и отрицательный полюс у которых частота
  • пара комплексно-сопряжённых нулей и и пара комплексно-сопряжённых полюсов

и у которых вещественные части отрицательны, а знаки соответствующих частот противоположны.

Отображение нулей и полюсов на p-плоскости называют картой нулей и полюсов. Различают левую и правую р-полуплоскости.

Карта нулей и полюсов позволяет оценить ряд свойств электрической цепи и, в частности, определить её устойчивость с точки зрения устойчивости передаточной функции.

Утверждение:

цепь является строго устойчивой тогда и только тогда, когда её передаточная функция имеет, полюсы только в левой р-полуплоскости, исключая мнимую ось.

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

вызывает ограниченную по величине реакцию

Но реакцию y(t) при нулевых условиях можно найти с помощью уравнения свёртки

Отсюда при заданных ограничениях имеем соотношение

из которого следует, что для получения равномерно ограниченной для всех t реакции, т. е. для обеспечения строгой устойчивости цепи должно выполняться условие абсолютной сходимости интеграла от импульсной характеристики:

(17.28)

Найдём расположение полюсов, которое соответствует полученному условию. Для этого представим импульсную характеристику h(t) цепи как обратное L-изображение передаточной функции (17.15) путём разложения последней на сумму простых дробей (16.28):

(17.29)

Подставим в интеграл (17.28) правую сумму (17.29)

(17.30)

и проведём ряд несложных преобразований.

Поскольку модуль суммы не превосходит суммы модулей, справедливо следующее неравенство:

В правой части полученного неравенства поменяем местами знаки суммирования и интегрирования и оставим только знак равенства:

Рассмотрим интеграл в правой части равенства, содержащий модуль экспоненты, при

во-вторых, первая экспонента под интегралом всегда неотрицательна, поэтому знак модуля можно опустить:

Остаётся исследовать сходимость интеграла при положительном и отрицательном показателе

Сходимость интеграла при означает, что для устойчивости передаточной функции (а потому и цепи), все полюсы должны иметь отрицательные действительные части т. е. лежать в левой р-полуплоскости, что и требовалось доказать.

Связь передаточной функции с частотными и временными характеристиками цепи

Как было показано в лекции 10, для определения частотных характеристик АЧХ и ФЧХ цепи необходимо знать комплексную частотную характеристику Получить КЧХ из передаточной функции несложно: необходимо лишь в (17.10) заменить оператор на поскольку частотные характеристики являются непрерывными функциями только частоты:

(17.31)

(17.32)

(17.33)

Эквивалентное выражение для КЧХ получается из (17.31), если воспользоваться комплексными функциями числителя и знаменателя:

(17.34)

(17.35)

(17.36)

Вследствие того, что функция

является иррациональной, обычно при анализе и синтезе цепей используют квадрат АЧХ:

(17.37)

Перечислим основные свойства передаточных функций и квадрата АЧХ пассивных цепей.

  1. Передаточная функция является L-изображением импульсной характеристики.
  2. Передаточная функция является дробно-рациональной функцией с вещественными коэффициентами.
  3. Полюсы устойчивой передаточной функции лежат в левой р-полуплоскости.
  4. Степени полиномов числителей передаточной функции и квадрата АЧХ не превышают степеней полиномов знаменателей; приневыполнении этого свойства АЧХ на бесконечно больших частотах должна принимать бесконечно большое значение, поскольку числитель в этом случае растёт быстрее знаменателя.
  5. Частотные характеристики цепи вычисляются по передаточной функции при подстановке
  6. Квадрат АЧХ является чётной рациональной функцией переменной с вещественными коэффициентами:
  7. По передаточной функции можно изобразить схему цепи

Обобщённая схема связи передаточной функции с характеристиками и свойствами цепи представлена на рис. 17.6.

Рекомендую подробно изучить предметы:
  1. Электротехника
  2. Основы теории цепей
Ещё лекции с примерами решения и объяснением:
  • Свободные колебания в пассивных электрических цепях
  • Цепи с распределёнными параметрами
  • Волновые параметры длинной линии
  • Колебания в линиях без потерь
  • Комплексные функции электрических цепей
  • Гармонические колебания в колебательном контуре
  • Частотные характеристики линейных электрических цепей
  • Частотные методы анализа и расчёта электрических цепей

При копировании любых материалов с сайта evkova.org обязательна активная ссылка на сайт www.evkova.org

Сайт создан коллективом преподавателей на некоммерческой основе для дополнительного образования молодежи

Сайт пишется, поддерживается и управляется коллективом преподавателей

Whatsapp и логотип whatsapp являются товарными знаками корпорации WhatsApp LLC.

Cайт носит информационный характер и ни при каких условиях не является публичной офертой, которая определяется положениями статьи 437 Гражданского кодекса РФ. Анна Евкова не оказывает никаких услуг.

Контрольная работа: Передаточные функции одноконтурной системы

Практическая работа № 1

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

Оценить устойчивость каждого из звеньев.

а) ; б).

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

.

1. а). Дифференциальное уравнение можно записать в виде:

.

Обозначим Y(s) и F(s) как изображения сигналов соответственно y и f , тогда операторное уравнение (при нулевых начальных условиях) примет вид:

1,25s3Y(s) – 4s2Y(s) + 5sY(s) = 3F(s) – sF(s).

Данное уравнение можно преобразовать, вынеся Y(s) и F(s) за скобки:

Y(s). (1,25s3 – 4s2 + 5s) = F(s). (3 – s).

.

Очевидно, что входной сигнал x отсутствует, и выходной сигнал у определяется только внешним воздействием f (система, действующая по возмущению): , то получается уравнение Y(s) = WF(s).F(s). Структурная схема объекта приведена на рис. 1.

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

A(s) =.

Если приравнять данное выражение к нулю, то образуется характеристическое уравнение , корни которого:

, и .

Распределение корней на комплексной плоскости показано на рис. 2. По рисунку видно, что корни лежат в правой полуплоскости, следовательно, объект неустойчив.

б) Дифференциальное уравнение можно записать в виде:

.

Обозначим Y(s), X(s) и F(s) как изображения сигналов соответственно y , x и f , тогда операторное уравнение (при нулевых начальных условиях) примет вид:

2s2Y(s) + 4sY(s) + 10Y(s) = 3X(s) + 4sF(s).

Данное уравнение можно преобразовать, вынеся Y(s) и X(s) за скобки:

Y(s). (5s2 + 4s + 10) = 3X(s) + 4sF(s).

.

Если обозначить передаточные функции объекта как

и ,

то получается уравнение Y(s) = Wx(s).X(s) + WF(s).F(s). Структурная схема объекта приведена на рис. 3.

Характеристическая функция имеет вид:

,

а характеристическое уравнение:

.

Корни этого уравнения равны:

и .

Распределение корней на комплексной плоскости показано на рис. 4:

Все корни характеристического уравнения лежат в левой полуплоскости, очевидно, что объект устойчив.

2. Дана передаточная функция вида:

Зная, что по определению, , получим:

, тогда:

.

Применяя к полученному выражению обратное преобразование Лапласа, находим искомое дифференциальное уравнение:

.

Практическая работа № 2

Дана одноконтурная АСР, для которой определена передаточная функция регулятора (Р) с настройками и дифференциальное уравнение объекта управления (ОУ). Требуется определить:

— передаточную функцию разомкнутой системы W∞(s),

— характеристическое выражение замкнутой системы (ХВЗС),

— передаточные функции замкнутой системы Фз(s) – по заданию, Фв(s) – по возмущению, ФЕ(s) – по ошибке,

— коэффициенты усиления АСР,

Р — ПИ-регулятор с ПФ вида ;

дифференциальное уравнение объекта управления:

.

Определим передаточную функцию объекта:

W об( s ) .

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

Характеристическое выражение замкнутой системы:

;

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

— по заданию;

— по ошибке;

— по возмущению.

По передаточным функциям определим коэффициенты усиления путем подстановки в них s = 0:

К3 = Ф3(0) = 1 – по заданию;

КЕ = ФЕ(0) = 0 – по ошибке;

Кв = Фв(0) = 0 – по возмущению.

Определим устойчивость АСР по критерию Гурвица.

Так как коэффициенты ХВЗС а3 = 4, а2 = 6, а1 = 18, а0 = 4 (степень полинома n = 3), то матрица Гурвица имеет вид:

Диагональные миноры матрицы равны соответственно:

Поскольку все определители положительны, то АСР является устойчивой.

Практическая работа № 3

По табличным данным построить переходную кривую объекта, определить параметры передаточной функции объекта, рассчитать настройки ПИД-регулятора, обеспечивающие 20%-е перерегулирование.

DXвх = 5,5 кПа; DY = 0,149 %; tзап = 40 сек

Операторное уравнение по передаточной функции

Понятие линейного динамического звена

САУ удобно представлять для анализа и при синтезе в виде взаимосвязанной совокупности отдельных элементов – динамических звеньев.

Под динамическим звеном понимают в общем случае абстрактное устройство, имеющее вход и выход, и для которого задано уравнение, связывающее сигналы на входе и выходе, как это показано на рис. 1.

Подробное изучение свойств реальных объектов управления и систем автоматического управления приводит к описанию динамических звеньев в виде нелинейных дифференциальных уравнений. Но во многих случаях их можно линеаризовать, то есть заменить нелинейные уравнения линейными, приближенно описывающими процессы в системах. Тем самым осуществляется декомпозиция задач анализа и синтеза систем, то есть первоначально используют линейное представление, а затем осуществляют учет вносимых нелинейностями особенностей. Такому подходу способствует то, что, в большинстве случаев, нормально функционирующая система работает в режиме малых отклонений, при которых нелинейности не проявляются. В дальнейшем мы будем рассматривать преимущественно аппарат изучения линейных систем, а особенности систем других классов: нелинейных, импульсных, цифровых и стохастических, будут излагаться позднее в других учебных дисциплинах.

Если уравнение, связывающее сигналы и , линейно, то говорят о линейном динамическом звене

Уравнение линейного динамического звена имеет следующий общий вид:

где — постоянные коэффициенты, .

Использовать такое описание динамического звена в задачах анализа и синтеза систем и объектов управления не рационально, поэтому существуют и иные формы описания и представления динамических звеньев и систем в целом.

Подвергнем уравнение (1) преобразованию Лапласа, считая начальные условия нулевыми и заменяя оригиналы сигналов их изображениями по Лапласу

.

Используя теоремы преобразования Лапласа линейности и дифференцирования, получим операторное уравнение, связывающие изображения входного и выходного сигналов

Преобразуем уравнение (2) к следующему виду

Получим из (3) отношение изображений выходного и входного сигналов

Отношение (4) не зависит от изображений сигналов, определяется только параметрами самого динамического звена (), имеет вид дробно-рациональной функции.

Отношение изображений выходного и входного сигналов называют передаточной функцией динамического звена

.

,

называют характеристическим уравнением динамического звена, так как знаменатель передаточной функции – это характеристический полином дифференциального уравнения, описывающего динамическое звено.

Определим передаточную функцию динамического звена по его принципиальной электрической схеме

По второму закону Кирхгоффа запишем уравнения описывающие схему

С учетом того, что

,

Получим операторные уравнения

Из второго уравнения выразим значение изображения тока

Подставим полученное выражение в первое уравнение системы

.

В итоге получаем искомую передаточную функцию

.

Графически передаточные функции динамического звена представляют в следующем виде:

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

.(5)

В общем случае САУ состоит из множества динамических звеньев, сигналы с выходов звеньев могут суммироваться или вычитаться, суммироваться с внешними для САУ сигналами. Суммирование и вычитание изображений сигналов могут быть представлено графически с помощью суммирующих звеньев:

Показанная выше неоднозначность графического представления вычитания изображений на суммирующем элементе связана с различием в стандартах разных стран.

Используя графическое представление передаточных функций звеньев и суммирующие звенья, можно в графической форме представить операторные уравнения, описывающие САУ. Такое графическое представление операторных уравнений в ТАУ называют структурной схемой.

По математической модели объекта управления в форме системы дифференциальных уравнений определить структурную схему объекта.

Получим систему операторных уравнений, подвергнув исходную систему дифференциальных уравнений преобразованию Лапласа и заменив оригиналы изображениями,

Из первого уравнения системы операторных уравнений, которое описывает динамическое звено объекта управления, после преобразований получим

.

Тогда передаточная функция этого звена имеет вид

,

а выражение описывает суммирующее звено . Таким образом, получены два фрагмента структурной схемы

Из второго уравнения системы операторных уравнений, которое описывает динамическое звено объекта управления, после преобразований получим, вводя обозначение,

.

Тогда передаточная функция этого звена имеет вид

,

а выражение описывает суммирующее звено . Таким образом, получены еще два фрагмента структурной схемы

Соединим все фрагменты структурной схемы объекта управления, объединяя одноименные сигналы, либо разветвляя их с помощью точек ветвления , показанных на схеме. В результате получим

Временные характеристики динамического звена

Временной или импульсной характеристикой динамического звена называют реакцию звена на , обозначая ее как . При этом схема эксперимента имеет вид –

Выясним, что представляет собой временная характеристика, то есть почему ее называют характеристикой динамического звена?

Для этого рассмотрим динамическое звено с передаточной функцией

В этом случае, в соответствии с (5), имеем

.

Получаем, что передаточная функция звена – это изображение по Лапласу импульсной характеристики динамического звена. В свою очередь, импульсная характеристика может быть определена по передаточной функции

,

при использовании разложения в форму Хэвисайта и обратное преобразование Лапласа.

Знание импульсной характеристики позволяет определить реакцию динамического звена на сигнал любой формы.

Для динамического звена с передаточной функцией преобразуем (5), используя теорему об умножении изображений преобразования Лапласа,

,

а если легко получить , тогда

.

Переходной характеристикой или переходной функцией динамического звена называют реакцию динамического звена на , обозначая ее как . При этом схема эксперимента имеет вид –

Для анализа переходной характеристики рассмотрим динамическое звено с передаточной функцией

В этом случае, в соответствии с (5), имеем

.

По теореме об интегрировании оригинала имеем

Переходная функция является интегралом по времени от импульсной характеристике и наоборот

.

Переходная характеристика динамического звена может быть определена по передаточной функции

Контрольные вопросы и задачи

Что такое линейное динамическое звено?

Как определить передаточную функцию линейного динамического звена?

Перечислите основные элементы структурных схем систем управления.

Как определить по передаточной функции динамического звена его временные характеристики: импульсную и переходную?

Как по переходной характеристике определить импульсную характеристику динамического звена?

Определите передаточную функцию динамического звена по его принципиальной электрической схеме

.

Определите передаточную функцию динамического звена по его принципиальной электрической схеме

.

По математической модели объекта управления в форме системы дифференциальных уравнений определить структурную схему объекта.


источники:

http://www.bestreferat.ru/referat-292845.html

http://drive.ispu.ru/elib/lebedev/3.html

Название: Передаточные функции одноконтурной системы
Раздел: Рефераты по математике
Тип: контрольная работа Добавлен 21:52:15 24 декабря 2010 Похожие работы
Просмотров: 508 Комментариев: 14 Оценило: 3 человек Средний балл: 5 Оценка: неизвестно Скачать