M: lab 4 5 6 7

Author: lnq99

Modeling/4/main.py

@@ -0,0 +1,81 @@
+from system.model import Model
+from system.generator import *
+from system.queue import Queue
+from system.service import Service
+
+from PyQt6 import uic
+from PyQt6.QtCore import pyqtSlot
+from PyQt6.QtWidgets import QApplication, QWidget, QMessageBox
+
+
+class MainWindow(QWidget):
+    def __init__(self):
+        super().__init__()
+        self.ui = uic.loadUi('main.ui', self)
+        self.init_ui()
+
+    def init_ui(self):
+        ui = self.ui
+        ui.a.setValue(2)
+        ui.b.setValue(8)
+        ui.m.setValue(5)
+        ui.sigma.setValue(1)
+        ui.return_factor.setValue(0.1)
+        ui.dt.setValue(10e-3)
+
+        ui.n_tasks.setRange(5, 5000)
+        ui.n_tasks.setValue(100)
+        ui.run.setStyleSheet('background-color:#388e3c')
+
+    @pyqtSlot(name='on_run_clicked')
+    def on_run(self):
+        ui = self.ui
+        try:
+            a = ui.a.value()
+            b = ui.b.value()
+            m = ui.m.value()
+            sigma = ui.sigma.value()
+            n_tasks = ui.n_tasks.value()
+            return_factor = ui.return_factor.value()
+            dt = ui.dt.value()
+
+            g_uniform = UniformGenerator(a, b)
+            g_normal = NormalGenerator(m, sigma)
+            queue = Queue()
+
+            print('Принцип Δt')
+            generator = RequestGenerator(g_uniform, n_tasks)
+            service = Service(g_normal, return_factor)
+            model = Model(generator, queue, service)
+            t, completed_tasks, queue_len_max = model.time_based(n_tasks, dt)
+            self.show_result(self.ui.result_time, t, completed_tasks, completed_tasks-n_tasks, queue_len_max)
+
+            print()
+
+            print('Событийный принцип')
+            generator = RequestGenerator(g_uniform, n_tasks)
+            service = Service(g_normal, return_factor)
+            model = Model(generator, queue, service)
+            t, completed_tasks, queue_len_max = model.event_based(n_tasks)
+            self.show_result(self.ui.result_event, t, completed_tasks, completed_tasks-n_tasks, queue_len_max)
+
+        except Exception as e:
+            QMessageBox.critical(self, 'Error', str(e))
+
+    @staticmethod
+    def show_result(text_edit, t, completed_tasks, returned_tasks, queue_len_max):
+        result_str = f'Обработанные заявки: {completed_tasks}\n\n'\
+            f'Повторно обработанные заявки: {returned_tasks}\n\n'\
+            f'Максимальная длина очереди: {queue_len_max}\n\n'\
+            f'Время работы: {t:.2f}'
+
+        text_edit.setPlainText(result_str)
+
+
+if __name__ == '__main__':
+    import sys
+
+    app = QApplication(sys.argv)
+    window = MainWindow()
+    window.show()
+    app.exec()

Modeling/4/main.ui

@@ -0,0 +1,132 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<ui version="4.0">
+ <class>Form</class>
+ <widget class="QWidget" name="Form">
+  <property name="geometry">
+   <rect>
+    <x>0</x>
+    <y>0</y>
+    <width>451</width>
+    <height>490</height>
+   </rect>
+  </property>
+  <property name="windowTitle">
+   <string>Form</string>
+  </property>
+  <widget class="QWidget" name="gridLayoutWidget">
+   <property name="geometry">
+    <rect>
+     <x>10</x>
+     <y>10</y>
+     <width>431</width>
+     <height>151</height>
+    </rect>
+   </property>
+   <layout class="QGridLayout" name="gridLayout">
+    <item row="3" column="0">
+     <widget class="QLabel" name="label_4">
+      <property name="text">
+       <string>Вероятность возврата заявки</string>
+      </property>
+     </widget>
+    </item>
+    <item row="0" column="0">
+     <widget class="QLabel" name="label_2">
+      <property name="text">
+       <string>Генератор U(a,b)</string>
+      </property>
+     </widget>
+    </item>
+    <item row="4" column="0">
+     <widget class="QLabel" name="label_5">
+      <property name="text">
+       <string>Δt</string>
+      </property>
+     </widget>
+    </item>
+    <item row="1" column="0">
+     <widget class="QLabel" name="label_3">
+      <property name="text">
+       <string>Обслуживающий аппарат N(m,σ^2)</string>
+      </property>
+     </widget>
+    </item>
+    <item row="2" column="0">
+     <widget class="QLabel" name="label">
+      <property name="text">
+       <string>Количество заявок</string>
+      </property>
+     </widget>
+    </item>
+    <item row="0" column="1">
+     <widget class="QDoubleSpinBox" name="a"/>
+    </item>
+    <item row="0" column="2">
+     <widget class="QDoubleSpinBox" name="b"/>
+    </item>
+    <item row="1" column="1">
+     <widget class="QDoubleSpinBox" name="m"/>
+    </item>
+    <item row="1" column="2">
+     <widget class="QDoubleSpinBox" name="sigma"/>
+    </item>
+    <item row="3" column="1" colspan="2">
+     <widget class="QDoubleSpinBox" name="return_factor"/>
+    </item>
+    <item row="4" column="1" colspan="2">
+     <widget class="QDoubleSpinBox" name="dt"/>
+    </item>
+    <item row="2" column="1" colspan="2">
+     <widget class="QSpinBox" name="n_tasks"/>
+    </item>
+   </layout>
+  </widget>
+  <widget class="QPushButton" name="run">
+   <property name="geometry">
+    <rect>
+     <x>120</x>
+     <y>190</y>
+     <width>211</width>
+     <height>31</height>
+    </rect>
+   </property>
+   <property name="text">
+    <string>Моделировать</string>
+   </property>
+  </widget>
+  <widget class="QWidget" name="gridLayoutWidget_2">
+   <property name="geometry">
+    <rect>
+     <x>10</x>
+     <y>250</y>
+     <width>431</width>
+     <height>221</height>
+    </rect>
+   </property>
+   <layout class="QGridLayout" name="gridLayout_2">
+    <item row="0" column="0">
+     <widget class="QLabel" name="label_6">
+      <property name="text">
+       <string>Принцип Δt</string>
+      </property>
+     </widget>
+    </item>
+    <item row="0" column="1">
+     <widget class="QLabel" name="label_7">
+      <property name="text">
+       <string>Событийный принцип</string>
+      </property>
+     </widget>
+    </item>
+    <item row="1" column="0">
+     <widget class="QPlainTextEdit" name="result_time"/>
+    </item>
+    <item row="1" column="1">
+     <widget class="QPlainTextEdit" name="result_event"/>
+    </item>
+   </layout>
+  </widget>
+ </widget>
+ <resources/>
+ <connections/>
+</ui>

Modeling/4/system/__init__.py

@@ -0,0 +1,69 @@
+from numpy.random import default_rng
+
+
+class IGenerator:
+    generator = default_rng()
+
+    def generate(self) -> float:
+        raise NotImplementedError
+
+
+class UniformGenerator(IGenerator):
+    def __init__(self, low, high):
+        if low > high: low, high = high, low
+        self.low = low
+        self.high = high
+
+    def generate(self):
+        return self.generator.uniform(self.low, self.high)
+
+
+class NormalGenerator(IGenerator):
+    def __init__(self, loc, scale):
+        self.loc = loc
+        self.scale = scale
+
+    def generate(self):
+        return self.generator.normal(self.loc, self.scale)
+
+
+class RequestGenerator:
+    def __init__(self, generator: IGenerator, n_task: int):
+        self._generator = generator
+        self._remaining_time = 0
+        self._is_ready = False
+        self._generated_tasks = 0
+        self._n_task = n_task
+
+    def generate(self):
+        r = -1
+        while r <= 0:
+            r = self._generator.generate()
+        self._remaining_time = r
+        self._is_ready = True
+        self._generated_tasks += 1
+        return self._remaining_time
+
+    def elapse(self, time):
+        if not self._is_ready:
+            if self._generated_tasks == self._n_task:
+                return 0
+            r = self._remaining_time - time
+            if r <= 10e-4:
+                r = self.generate() - r  # suppose >0
+            self._remaining_time = r
+            return self._remaining_time
+
+    def pop(self):
+        if self._is_ready:
+            self._is_ready = False
+            return self._generated_tasks    # task id
+        return None
+
+    @property
+    def is_ready(self):
+        return self._is_ready
+
+    @property
+    def generated_tasks(self):
+        return self._generated_tasks

Modeling/4/system/generator.py

@@ -0,0 +1,75 @@
+from .generator import RequestGenerator
+from .queue import Queue
+from .service import Service
+
+
+LOG_FREQ = 5
+
+
+class Model:
+    def __init__(self, generator: RequestGenerator, queue: Queue, service: Service):
+        self._generator = generator
+        self._queue = queue
+        self._service = service
+
+    def time_based(self, n_tasks, dt=10e-3):
+        generator, queue, service = self._generator, self._queue, self._service
+        t_current = 0
+        n_tasks_for_logging = n_tasks // LOG_FREQ
+
+        while generator.generated_tasks < n_tasks \
+                or not queue.is_empty or not service.is_ready:
+            generator.elapse(dt)
+            _, return_task = service.elapse(dt)
+
+            if return_task:
+                queue.enqueue(return_task)
+
+            if generator.is_ready:
+                task = generator.pop()
+                queue.enqueue(task)
+
+            if service.is_ready and not queue.is_empty:
+                if service.completed_tasks % n_tasks_for_logging == 0:
+                    print(f'Processed {service.completed_tasks:3d} tasks, t = {t_current:.4f}')
+                task = queue.dequeue()
+                service.process(task)
+
+            t_current += dt
+
+        return t_current, service.completed_tasks, queue.len_max
+
+    def event_based(self, n_tasks):
+        """
+        events: task_generated, task_completed
+        """
+        generator, queue, service = self._generator, self._queue, self._service
+        t_current = 0
+        n_tasks_for_logging = n_tasks // LOG_FREQ
+        dt = 0
+
+        while generator.generated_tasks < n_tasks \
+                or not queue.is_empty or not service.is_ready:
+            t_current += dt
+            t_remain_task_generated = generator.elapse(dt)
+            t_remain_task_completed, return_task = service.elapse(dt)
+
+            if return_task:
+                queue.enqueue(return_task)
+
+            if generator.is_ready:
+                task = generator.pop()
+                queue.enqueue(task)
+
+            if service.is_ready and not queue.is_empty:
+                if service.completed_tasks % n_tasks_for_logging == 0:
+                    print(f'Processed {service.completed_tasks:4d} tasks, t = {t_current:.4f}')
+
+                task = queue.dequeue()
+                t_remain_task_completed = service.process(task)
+
+            t_list = [t_remain_task_generated, t_remain_task_completed]
+            t_list = list(filter(lambda x: x > 0, t_list))
+            dt = min(t_list) if t_list else 0
+
+        return t_current, service.completed_tasks, queue.len_max

Modeling/4/system/model.py

@@ -0,0 +1,34 @@
+from collections import deque
+
+
+class Queue:
+
+    def __init__(self, maxlen=None):
+        self._queue = deque(maxlen=maxlen)
+        self._len_max = 0
+        self._len = 0
+
+    def enqueue(self, item):
+        self._queue.append(item)
+        self._len += 1
+        if self._len > self._len_max:
+            self._len_max = self._len
+
+    def dequeue(self):
+        if self._len > 0:
+            self._len -= 1
+            return self._queue.pop()
+        else:
+            return None
+
+    @property
+    def is_empty(self):
+        return self._len <= 0
+
+    @property
+    def len_max(self):
+        return self._len_max
+
+    @property
+    def len(self):
+        return self._len