Задача по мат моделированию

@DWfrya · Регистрация: 18.03.2014

Студворк — интернет-сервис помощи студентам

Нужно написать программу на Java(либо любом ООП) и любом языке типа Simula.
Задачка такая:
В школе имеется 100 компьютеров. Компьютеры периодически выходят из строя. Время нормальной работы распределено нормально с матожиданием 30 дней и дисперсией 15 дней. Время ремонта распределено равномерно от 1 до 3 дней. Сколько нужно ремонтников, для того что бы постоянно работало 90 компьютеров в будни дни.

Буду очень признателен за помощь, я просто даже не знаю с чего начать
Написал только методы генераторов чисел необходимые.

Нашел что-то на питоне только

Python

# -*- coding: cp1251 -*-
import math
import random
 
compsPos = []
queue = {}
maxComp = 100
 
def getWorkTime():
    cv = math.floor(random.gauss(30, 15))
    if (cv < 0):
        cv = 0
    return cv
 
def addqueue(item):
    l = math.floor(random.gauss(2, 1))
    if l < 1:
        l = 1
    if l > 3:
        l = 3
    queue[item] = l
    print 'comp ',item,' added to queue for ',l,' days'
    
 
# первично заполняем массив компов их временами работы
for i in range(0, maxComp, 1):
    cv = getWorkTime()
    compsPos.append(cv)
 
print compsPos
rl = []
# дни
for i in range(365):
    print 'day ',i
    # цикл выходов из строя
    for j in range(0, maxComp, 1):
        # этот комп сейчас в ремонте
        if (j in queue.iterkeys()):
            continue
        # этот комп сломается 100%, так как ему работать 0 дней
        if (0 == compsPos[j]):
            addqueue(j)
            continue
        # иначе оцениваем вероятность поломки
        pos = 1-(1 / compsPos[j])
        # бросаем случайное числе
        # если оно больше вероятности поломки, то комп сломался
        # иначе работает
        cur = random.random()
        if (cur > pos):
            addqueue(j)
            continue
        # уменьшаем на 1 время работы компа        
        compsPos[j] = compsPos[j] - 1
    rl.append(len(queue))
    # цикл возвратов в строй
    for j in queue.keys():
        queue[j] = queue[j]-1
        # комп починен, возвращаем в строй
        if (0 == queue[j]):
            # задаём новое время работы
            cv = getWorkTime()
            compsPos[j] = cv
            print 'comp ',j,' returned with new time ',cv
            # удаляем из очереди
            del queue[j]
 
print rl
print max(rl)

@xoraxax · 06.12.2016, 13:55

с теории начни, опиши русским языком как эта хрень считается, кто-нибудь с программой поможет

@DWfrya · 07.12.2016, 11:47 **[ТС]**

В общем на java написал. Кому интересно, выкладываю код.
Теперь нужно это как то на имитационный язык перевести.
Тогда в другом разделе открою тему. Задача по мат моделирвоанию

Java

package com.company;
 
import java.util.ArrayList;
 
public class Main {
    private long count = 0;
    private static long a = 4096;
    private static long c = 150889;
    private static long m = 714025;
    private static long seed = 23;
    static int maxComp = 100;
    static ArrayList<Comp> compList = new ArrayList<>();
    static ArrayList<Integer> repairList = new ArrayList<>();
    public static void main(String[] args) {
        for (int i = 0; i < maxComp; i++) {
            compList.add(new Comp(i, getWorkTime(), -1));
        }
        int kol[] = new int[365];
        int day=0;
        for (int i = 0; i < 365; i++) {
            System.out.println("day "+i);
            for (int j = 0; j < maxComp; j++) {
                if (compList.get(j).repairTime>0) {continue;}
                if (compList.get(j).workTime == 0 && compList.get(j).repairTime==-1) {
                    addqueue(j, compList.get(j));
                    repairList.add(compList.get(j).id);
                    continue;
                }
                if (day<5) compList.get(j).workTime--;
            }
            day++;
            kol[i]=repairList.size();
 
            for (int j = 0; j < repairList.size(); j++) {
                int idComp = repairList.get(j);
                if (compList.get(idComp).repairTime>0) compList.get(idComp).repairTime--;
                else if (compList.get(idComp).repairTime==0) {
                    compList.get(idComp).workTime=getWorkTime();
                    compList.get(idComp).repairTime=-1;
                    System.out.println("comp "+idComp+" returned with new time "+compList.get(idComp).workTime);
                    repairList.remove(j);
                }
            }
            if (day==7) day=0;
 
        }
        int max = kol[0];
        for (int i = 1; i < kol.length; i++) {
            if (kol[i]>max) max=kol[i];
        }
        System.out.println(max+" computers are broken max");
        System.out.println("need "+(max-10)+" workers");
 
    }
 
    public static int getWorkTime() {
        int cv = Normal(30, 15);
        if (cv<0) cv=30;
        return cv;
    }
 
    public static void addqueue(int item, Comp comp) {
        double l = Ravnomern(1,4);
        if (l<1) l =1;
        else if (l>3) l=3;
        int il = (int)(Math.round(l));
        comp.repairTime=il;
        System.out.println("comp "+item+" added to queue for "+il+" days");
    }
 
    public static int Normal(double mat, double dis)
    {
        double U = Math.random();
        double V = Math.random();
 
        double R1 = Math.sqrt(-2 * Math.log(U)) * Math.sin(2 * Math.PI * V);
        double R2 = Math.sqrt(-2 * Math.log(U)) * Math.cos(2 * Math.PI * V);
        double R = (R1 + R2) / 2;
        return (int)(R * dis + mat);
    }
 
    public static double getRand() {
        seed = (a * seed + c) % m;
        return seed / (double) m;
    }
 
    public static double Ravnomern (int a, int b){
        double ret = 0;
        ret = (b - a)*getRand() + a;
        return ret;
    }
}
 
package com.company;
 
/**
 * Created by Денис on 06.12.2016.
 */
public class Comp {
    int workTime, repairTime, id;
    public Comp(int id, int workTime, int repairTime) {
        this.id=id;
        this.workTime=workTime;
        this.repairTime=repairTime;
    }
}

@xoraxax 3638 / 2970 / 918 Регистрация: 05.07.2013 Сообщений: 14,220
	06.12.2016, 13:55	2
	с теории начни, опиши русским языком как эта хрень считается, кто-нибудь с программой поможет 0

	07.12.2016, 11:47

Опции темы

Задача по мат моделированию

Решение