Решение задачи СМО с использованием системы MathCAD
Решим задачу СМО, приведенную в примере 4, с использованием системы MathCAD.
Ввод текста на всех этапах решения задачи осуществляют с помощью комбинаций клавиш «Shift + "» , что позволяет создать текстовую область. Введем на рабочем листе первый пункт расчета:
1. Задание исходных данных.
Затем последовательно введем исходные данные в поле экрана (рис. 6.18):
X -0.25 ц :=0.33 п ~2
Для решения задачи воспользуемся блоком функций Given…Find, Для этого необходимо задать начальные приближения. Введем на рабочем листе второй пункт расчета:
Y. Задание исходных данных: X := 0.25 Ji :=0 ВВзЗГ"
2. Начальные приближения. РО :=0 2Ь PI :Ч).25 Р2 :=0.25 РЗ :=0 25
3. Запись системы уравнений, описывающей функционирование многоканальной СМО
Given
—X-?Q + I Р1 = 0 А-Р0-(А + ц)-Р1+2-й Р2-0 А-Р1-(Я. + 2 [х) Р2н ЗцРЗ = 0 Р0+Р1+Р2+Р3 = 1
Л
4. Результаты решения — РО = 0.476 Р1 и 0.357 Р2 = 0.134 РЗ = 0.033
Рис.. 6. IS Определение параметров функционирования многоканальной замкнутой
СМО с отказами в системе MaihCAD
2. Начальные приближения
Затем последовательно введем начальные приближенные значения искомых парамет-ров:
Р0:=0.4 Р1-0.2 Р2:=0.2 РЗ-0.2
Введем на рабочем листе третий пункт расчета.
3. Запись системы уравнений, описывающей функционированиемногока нальнои СМО.
:= lrind(P0,P],P2P3) |
Вначале вводим ключевое слово С Yen (Дано), которое может быть напечатано прописными, строчными буквами или начинаться с прописной, затем исходную систему уравнений с жлрныч «=», а в заключение — вектор искомых величин. Для этого в ноле рабочего листа определяем местоположение вектора. Если в окне выведена панель инструментов Math (Математика), нужно щелкнуть но кнопке с изображением матрицы. Появится аналоговое окно Matrix (Матрица). Здесь щелкаем по кнопке с аналогичным изображением или нажимаем комбинацию клавиш «Ctrl + М». В обоих случаях появится диалоговое окно Insert Matrix (Вставить матрицу). В его текстовых полях Rows (Строки) и Columns (Столбцы) вставим нужное число ст рок и столбцов в нашей задаче 4 и 1 соответственно. После щелчка по кнопке ОК появится шаблон с метками для ввода искомых данных, Подводя курсор или указатель мыши к каждой из
Них, зададим искомые параметры, затем знак присваивания и имя встроенной функции Friid(P0,Pl, P2,P3,P4,P5). Выведем на рабочий лист четвертый пункт расчета: 4. Результаты решения.
Для получения результатов расчета искомых величин достаточно набрать имя H)Oki Юго параметра и знак J >авенсгва, па жав соо гветс гвук>- щук > клавишу или щелкнув по кнопке со знаком равенства, расположенной в верхнем левом уг. ту панели инструментов Evalu… (Вычисления;.
Рассмотрим неустановившийся режим работы СМО когда ее ос — н овные характеристики завися г от времени В этом случае функционирование СМО будем описывать системой обыкновенных дифференциальных уравнений:
#h=-(A] *-ix)P](t) + }.PQ(t) + 2iiP2(t), Щ = -(к+2хЩ (0 + XPl(t) + ЗцР, (/),
P‘{t)=-здею(6.113)
На рис. 6.19 представлены исходные данные и система дифференциальных уравнений, описывающая функционирование многоканальной замкнутой СМО при неустановившемся режиме работы.
1 Задание исходных данных многоканальной замкнутой СМО с отказами при неустановившемся режиме.
X := 0 25 ц := 0 3333
2 Система обыкновенных дифференциальных уравнений, описывающих функционирование многоканальной замкнутой СМО с отказами
Ip0(t)=-x-po(t)+li-p1(t) dt
TlltJHP P0(t)-(X +м) Pi(t) + 2-M. P2(t) dt
Dt
^{t)=A. p2(t)-3.M. p3(0
Pile. 6.19. Олисание функционирования многоканальной замкнутой СМО
При неустановившемся режиме
На рис. 6 20 показаны правые части системы уравнении в виде вектор-столбца. Каждый его элемент определяет значение правой части соответствующего дифференциального уравнения на любом шаге ин гегрирования (решения). Здесь же даны начальные значения искомых параметров в виде вектор — столбца. В нижнеи части рисунка определены начальное и конечное время интегрирования и число шаг ов решения системы дифференциальных уравнений.
3 Формирование функции, которая определяет вектор отклонений искомых величин D(t, Р) и векгор начальных значений Р.
-х-РоСО+цРИО |
‘п |
||
XP0-(0-(^) iY(0*-2 MY(t) |
Р— |
0 |
|
Х-Р1-(0-(Х4 2.ц)Р2(1) + 3.цР3(1) XP2(t)-3^P3(t) |
Г.— |
0 0 V / |
4. Определение дополнительных параметров для решения системы обыкновенных дифференциальных уравнений" начальное и конечное время исследования системы и число тагов интегрирования t0 := 0 tl := 6 N := 1000
Рис, 6 20. Представление системы дифференциальных уравнений в виде, доступном
Для решения ее в системе MathCAD
На рис. 6.21 приводится решение системы дифференциальных уравнений многоканачьнои замкнутой СМО с использованием встроенной функции Rkfixed(P,TO,Tl,N,D), реализующей метод Рунге — Кут — та с фиксированным шагом. Для вызова этой функции необходимо: щелкнуть по пункту Insert (Вставка) главного меню, а затем по пункту Function (Функция) падающего меню или нажать комбинацию клавиш «Ctrl + Е» Появится диалоговое окно Insert Function (Вставить функцию):
Паити в разделе Function Category (Категория функции) строку Differential Equation Solving (Решение дифференциального уравнения) и щелкнуть по пей левой кнопкой мыши. В правом поле в разделе Function Name (Имя функции) появится список функций для решения дифференциальных уравнений;
Найти в списке функцию Rkfixed и щелкнуть по ней мышью. После этого щелкнуть по кнопке ОК. В месге установки визира появится шаблон функции для решения системы уравнений методом Рунге-Кутга.
5. Решение системы обыкновенных дифференциальных уравнений с использованием встроенной функции rkhxed(P, tO, tl, N,D).
S:=rkfixcd(P, tO, tl, N,D) t:=S<0) PO:=S(1) PL=S<2> P2:=S<3> P3:=S<4>
6. Представление результатов решения в графическом виде.
Рис. 6.21 Решение системы дифференциальных уравнений, описывающих функционирование многоканальной замкнутой СМО с отказами
Можно сразу найти функцию rkfixed в правом поле раздела Function Name, после чего щелкнуть мышью по ней, а затем по кнопке ОК, но это займет больше времени. В обоих случаях в нижних полях диалогового окна будет дано правильное обозначение выбранной функции со всеми аргументами, а также краткое описание ее действий. На рис. 6.21 приведено графическое решение системы дифференциальных уравнений для первых двух искомых параметров, т. е. представлено поведение параметров РО, PI, Р2 и РЗ — вероятности отсутствия требований и наличия в системе одного, двух и трех требований соответственно в зависимости от длительности процесса.
Анализируя графическое решение системы обыкновенных диф — фс ренциальныхуравнений, описывающей функционирование заданной многоканальной замкнутой СМО, можно заметить, что примерно через 0,2 ч она переходит в установившийся режим работы. При этом значения вероятностей состояний режима работы системы при решении совокупности обыкновенных дифференциальных уравне
ний практически полностью соответствуют решению системы алгебраических уравнений для установившегося режима работы.