Диссертация (1173106), страница 22
Текст из файла (страница 22)
- ЛиАЗ 5256.53№Гос номерГодвыпускаLгод пробег за 2016 год10463тт05201328411 км234567891011120464тт05Lг= с нач. экс. Поданным тех. отдел66092 кмПримечание4 шины заменяли280/70 R 57230515 км72519 кмНе заменяли29208 км70078 кмНе заменяли0475тт0522857 км55237 кмНе заменяли0476тт0521512 км52257 кмНе заменяли0477тт0525129 км71148 кмНе заменяли0481тт0529685 км77514 кмНе заменяли0482тт0523192 км61695 кмНе заменяли0483тт0531700 км69876 кмНе заменяли0485тт0532730 км75000 км3 шины заменяли0486тт0537296 км83278 кмНе заменяли0489тт0530943 км77355 км2 шины заменялиСтоимость шин от 2200 до 2800 сомони или 15714 руб.
до 20000 руб.0473тт0520132013201320132013201320132013201320132013Сделать измерения для Шины:Вторая группа - автобусы год поставки в АТП 2007 г. - ЛиАЗ 5256№Гос номер12345678910111213140176тт010180тт010182тт010204тт010205тт010212тт010238тт010240тт010259тт010200тт010218тт010252тт010185тт010211тт01Годвыпуска20062006200620062006200620072007200720062006200720062006Lгод пробег за 2016 год16708 км20629 км16226 км13 881 км16240 км34564 км18377 км7155 км32139 км18241 км21564 км36745 км29284 км9180 кмLг= с нач. экс.
Поданным тех. отдел134250 км238437 км202264 км191999 км181360 км332863 км253489 км294406 км368393 км280190 км241281 км260081 км217850 км217042 кмПримеч15 шин заменял18 шин заменял16 шин заменял13 шин заменял15 шин заменял13 шин заменял15 шин заменял19 шин заменял22 шин заменял21 шин заменял17 шин заменял16 шин заменял20 шин заменял15 шин заменял178Приложение ВВыдержки из программы по расчету математической моделиlm0 = 74; W = 35000; Tp = 8;a1 = 726.3; s1 = 261.5; a2 = 942; s2 = 342; lam1 = 0.00114; lam2 = 0.00287;lmfl = 0.0006; Czf = 23.63; k1fl = 0;ltor = 0.00264; Ctor = 6.846; k1tr = 0;Cavt = 6400; Cbnk = Log [1 + 0.07]; k1f1 = 0;comment = "здесь lm0 - макс.годовой пробег в точке t=0; W - объем работ, \которые парк должен выполнить за год; Tp - продолжительность \эксплуатации в годах (основной параметр для оптимизации); a1 - \мат.ожидание закона распред.
вероятностей пробега до естественного \выбывания автобуса (списания по невозможности продолжать работу), s1 \- средн.квадр.отклонения этого ЗРВ; a2 - мат.ождн. ЗРВ нормалн. \составляющей функции, описывающей поведение бетта - т.е. поддавления \альфа(т) до коэффициента выпуска, s2 - сигма для того же ЗРВ, где а2; \lam1 - показатель степени для экспонентц. зависимости для \ктг=альфа-те; lam2 - показатель степени для экспонентц.
зависимости \для lmx[.,.]- максимальный годовой пробег в зависимости от \сумм. пробега; lmf1 - показатель степени для экспонентц. зависимости \для удельной стоимости топлива; Czf - удельная стоимость расхода \топлива в тыс. руб. на тыс. км., в нулевой точке; k1f1 - коэффициент \смещения для функции уд. стоимости топлива от сумм. пробега; ltor - \показатель степени для экспонент.
зависимости для удельной стоимости \ТО и ТР в тыс. руб. на тыс. км.; Ctor - то же, что и Czf, но для ТО и ТР; \k1tr - то же, что и k1f1, но для ТО и ТР; Cavt - стоимость нового \автобуса в тыс.руб.; Cbnk - показатель банковского роста";comment = "функции: R[] - функция безотказности парка, доля оставшихся вживых автобусов от сум. пробега; alft[] - коэф. техн. готовности ; \bet[] - коэфф. использования ; lmx[] - годовой пробег максимальный от \сумм.
пробега.";179R [l1_, a1_, s1_] := (1/2 (1 - Erf[(-a1 + l1)/(Sqrt[2] s1)]));alft[l1_, lm1_] := (Exp[-lm1*l1]);bet [l1_, a2_, s2_] := (1/2 (1 - Erf[(-a2 + l1)/(Sqrt[2] s2)]));lmx[l1_, lm2_] := lm0*(Exp[-lm2*l1]);comment ="значение Ai[y] - начальная закупка, стартовое значение количестваавтомобилей в момент 0. В зависимости от 'y' - предельного пробега списанияавтомобиля ";Ai[y_] := (W/(lmx[y, lam2]*R[y, a1, s1]));U[l11_] := NIntegrate[R[x1, a1, s1], {x1, 0, l11}];U2 = Interpolation [Table[{x, U[x]}, {x, 0, 2000, 10}]];l0[l11_] := (R[l11, a1, s1]*lmx[l11, lam2]);…………………………..For [y2 = 1, y2 < 26, y2++,ghz1[y2] =FindRoot[(IZflf2[x1, y2] + IZtrf2[x1, y2]) - x1*(1 + y2/4)/70 +Zaf2[x1, y2] - tbvlzz[y2] - 0.5 == 0, {x1, 1, tbvlyz[y2]}];ghz2[y2] =FindRoot[(IZflf2[x1, y2] + IZtrf2[x1, y2]) - x1*(1 + y2/4)/70 +Zaf2[x1, y2] - tbvlzz[y2] - 0.3 == 0, {x1, tbvlyz[y2], 80}];];tbvlxz1 =Interpolation [Table[{y3, ghz1[y3][[1]][[2]]}, {y3, 1, 25}],InterpolationOrder -> 2];tbvlxz2 =Interpolation [Table[{y3, ghz2[y3][[1]][[2]]}, {y3, 1, 25}],InterpolationOrder -> 2];zonaz = ParametricPlot3D[{x1,y2, (IZflf2[(x1 - 1.75)*4, y2] +180IZtrf2[(x1 - 1.75)*4, y2]) - ((x1 - 1.75)*4)*(1 + y2/4)/70 +Zaf2[(x1 - 1.75)*4, y2] + 0.01}, {y2, 1, 16}, {x1, 2, 25},RegionFunction ->Function [{x, y, z, y2, x1},And [x1 > tbvlxz1[y2]*0.25 + 1.75, x1 < tbvlxz2[y2]*0.25 + 1.75]],PlotPoints -> 50,PlotStyle -> {RGBColor[0.05, 0.9, 0.7], Opacity -> 0.7}]…………………………..For [y2 = 1, y2 < 15,y2++, (Print["N возраст.групп=", y2, " лев.граница=",tbvlxz1w[y2]*0.25 + 1.75, " центр.решение=",tbvlyzw2[y2]*0.25 + 1.75, " прав.граница=",tbvlxz2w[y2]*0.25 + 1.75, " сумм.издержки=", tbvlzzw2[y2]]);].