Задача на принятие решений в условиях риска! Delphi

[boss89]

Здраствуйте недавно столкнулся с задачей которую никак не получается реализовать. Обращаюсь к вам за помощью.
Задача.
Требуется принять решение о том, когда необходимо проводить профилактический ремонт ПЭВМ, чтобы минимизировать потери из-за неисправности. В случае если ремонт будет производится слишком часто, затраты на обслуживание будут большими при малых потерях из-за случайных поломок.

Так как невозможно предсказать заранее, когда возникнет неисправность, необходимо найти вероятность того, что ПЭВМ выйдет из строя в период времени t. В этом и состоит элемент ”риска”.

т.е. Математически это выглядит так: ПЭВМ ремонтируется индивидуально, если она остановилась из-за поломки. Через Tинтервалов времени выполняется профилактический ремонт всех n ПЭВМ. Необходимо определить оптимальное значение Т, при котором минимизируются общие затраты на ремонт неисправных ПЭВМ и проведение профилактического ремонта в расчёте на один интервал времени.

Пусть рt– вероятность выхода из строя одной ПЭВМ в момент t, а nt– случайная величина, равная числу всех вышедших из строя ПЭВМ в тот же момент. Пусть далее С1 – затраты на ремонт неисправной ПЭВМ и С2– затраты на профилактический ремонт одной машины.

Применение критерия ожидаемого значения в данном случае оправдано, если ПЭВМ работают в течение большого периода времени. При этом ожидаемые затраты на один интервал составят


где M(nt)– математическое ожидание числа вышедших из строя ПЭВМ в моментt. Так как nt имеет биномиальное распределение с параметрами (n, pt), то M(nt) = npt. Таким образом


Необходимые условия оптимальности T* имеют вид:

ОЗ (T*-1)? ОЗ (T*),

ОЗ (T*+1) ? ОЗ (T*).

Следовательно, начиная с малого значения T, вычисляют ОЗ(T), пока не будут удовлетворены необходимые условия оптимальности.

Выглядит довольно пугающе
Но всёже если ктото может помогите пожалуйста решить!

[Z1000000]

Теоретически так. Но критерий оптимальности никогда не выполняется. Либо он неправильный, либо я неправильно его вычисляю.

var
 i,N,T : integer;
 Pt,C1,C2 : Real;

function Zatrati ( t : Integer ) : Real;
begin
if t = 0 then begin Zatrati := 0; Exit; end;
Zatrati := N * ( C1*Pt*(t) + C2 ) / T;
end;

begin
//readln(N,Pt,C1,C2,T);
N := 100;
Pt := 0.1;
C1 := 10;
C2 := 50;
T := 100;
for i := 1 to T-1 do
 begin
 writeln ( Zatrati(i) );
 if ( Zatrati(i)<Zatrati(i+1) ) and ( Zatrati(i)< Zatrati(i-1) ) then
  begin
  writeln ('Искомый период найден : ', i); Exit;
  end;
 end;
writeln ('Критерий оптимальности не выполняется на заданном числе периодов T');
end.

[boss89]

Я вот не понимаю пользователь какие данные должен вводить?
Подскажите пожалуйста

А всё с этим разобрался))
Только разве мы Т сами задаем?