Архив за месяц: Февраль 2013

На некоторых российских колбасных заводах остро стоит вопрос с грызунами (rodental problem). В тех местах, которые редко посещает санитарно-эпидемиологическая служба, всегда возникает стимул сэкономить на гигиене. При этом одна проблема может быть решена совершенно экзотическим способом.

Никто не в состоянии идеально вычистить куттеры от кусочков мяса. Поэтому по ночам крысы любят залезать в механизмы и лакомиться свежей (или не очень) говядиной, свининой и бараниной. Бороться с крысами почти невозможно в силу ветхости колбасных заводов и наличия крысиных ходов под полом, под землёй и в вентиляционных системах. Что делать с грызунами, которые всю ночь пачкали лапами и хвостами агрегаты и гадили в мясо?

На некоторых заводах с проблемой борются очень оригинально. Ранним утром, вместо того чтобы прогнать крыс, работник измельчительного цеха включает куттеры. В ту же секунду в цехе раздаётся визг: «Уи-и-и-и-и!» — а затем хруст. Крысы погибают под ножами. Их кости достаточно тонкие, чтобы измельчиться без предварительной обвалки. Мне неизвестно, что происходит с этими покрошенными крысами далее. Скорее всего, из соображений гигиены эту массу, содержащую кишечник, кишечные массы, острые зубы, шерсть и мелкие косточки, выбрасывают. Известно, что труп становится непригодным в пищу именно из-за образования токсинов (в т. ч. кадеверина) в кишечнике. Колбаса с измельчённой целой крысой обладала бы слишком малым сроком годности. Поэтому разумно предположить, что крысиные массы изымаются из технологического процесса, а потребителю крыса может попасться в совершенно ничтожной концентрации в виде мелких остатков, не счищенных с ножей.

Естественно, я не перестану кушать колбасу. Естественно, при каждом употреблении колбасы я держу эту картину в голове. Более того, я пытаюсь определить состав мясного продукта и распробовать какие-нибудь необычные кусочки в надежде, что это окажется именно мышь или крыса (в случае с шаурмой — кошка, собака или бомж). Брезгливость — это одно из худших качеств, которое, к сожалению, регулярно культивируется человечеством в рекламе дезинфицирующих средств или стерильных предметов потребления. Господа, что бы вы кушали во время кораблекрушения? Анчоусы? Или дрожащими руками бы разделали чью-нибудь человечину и прибыли на берег совершенно невменяемым маньяком или, наоборот, недотрогой, утратившим человеческий язык и человеческий разум? Провели бы остаток жизни в комнате с мягкими стенами? Необходимо быть заранее готовым к любой продовольственной ситуации. В одной из следующих публикаций я расскажу, как правильно приготовить падаль, чтобы не потравиться.

Рис. 1. Тушка здорового деревенского петуха

Любое мясо может попасться в пищу. Даже трупы бездомных людей могут быть пущены в ход. Однако не стоит из-за этого отказываться от потребления колбасы или сосисок: доля незадекларированного мяса в продукте статистически незначима, поэтому можно ею пренебречь. (Но она неизбежно присутствует!)

Отличная задачка для тех, кто любит боль.

Задача. Агафья Степанна стала ипохондриком и теперь очень боится заболеть — она каждый час измеряет свою температуру градусником. Она бывает больше либо равной 36,6° (т. е. 36,6° плюс «надбавка».), причём «надбавка» распределена экспоненциально с параметром \(\lambda=4\). Если результат измерения превысил 37°, то Степанна начинает волноваться и оттого температурит так, что каждый следующий час результат измерения увеличивается по сравнению с предыдущим на случайную величину, имеющую то же самое распределение, что и надбавка (\(\mathop{\mathrm{Exp}}(4)\)). Как только достигнут или превышен потолок в 40°, Степанна теряет сознание, холодеет, и через час её температура возвращается на отметку 36,6°, а сама Степанна приходит в сознание и радуется тому, что абсолютно здорова. После этого она снова каждый час мерит температуру по тому же принципу.

• Каково матожидание её температуры за всё время?
• Какой процент всего времени Степанна температурит (\(t\geqslant37^{\circ}\))? Какова средняя продолжительность её здорового состояния? Какова средняя продолжительность её болезни? (Длина блуждания.)
• Если ипохондрик потеряла сознание, то какова вероятность, что она увидела на градуснике отметку выше, чем 40,5°? (Условная вероятность.)
• Новое условие: предположим, что если температура подскакивает до 41° или выше, то старушка умирает. (Это вполне возможное событие: если было 39,9°, а потом раз — и резко подскочила на 1,15° или выше; вероятность скачка на 1,15° или более равна ≈0,01.) Сколько часов в среднем живут старушки с того момента, как становятся ипохондриками? Каково матожидание температуры живых старушек? Если температура в пределах нормы, то какова вероятность умереть при следующем взлёте температуры? А в один из десяти следующих взлётов? (Длина блуждания. Теперь из матожидания будут выброшены аутлайеры, оно чуть понизится, но совсем на кроху.)
• Как изменятся все ответы в предыдущих пунктах, если после потери сознания старушкина температура падает не до 36,6°, а до 36,6° плюс надбавка (т. е. сразу разыгрывается случайным образом)?

Я боль не очень люблю, поэтому приведу не решение, а собственные размышления на тему.

Мне кажется, имеет смысл рассмотреть такой случайный процесс: есть нижняя граница коридора (36,6), средняя (37), верхняя (40) и смертельная (41). В коридоре в каждый момент времени с пола подпрыгивает частица. Как только в момент времени величина выскакивает из средней границы коридора, то начинается случайное блуждание вверх. Как только в результате случайного однонаправленного блуждания частица пробивает потолок, то частица падает на пол, и снова начинаются случайные прыжки с пола. Если же, вылетая, частица пересекла смертельную границу, то игра заканчивается. Требуется найти среднюю вертикальную координату этой частицы, а также долю наблюдений, когда частица вне меньшего коридора... Вроде всё понятно.
Кстати, процент времени болезни не может быть равен 100, так как после каждого «рекорда» следующий период она здорова (36,6), а длина температуренья конечна (кстати, это потребуется доказать в процессе решения).

Реализация Агафьи Степанны в программе Mathematica.

t = 1;    (* Начальный момент времени... 1 или 0 на вкус*)
x = 36.6; (* Пол *)
i = 1;    (* Для того чтобы делать циклы *)
cycles = 10000; (* Как бы количество циклов *)
coord = { }; (* Резервируем место под запись наблюдений *)
While[i < cycles, (* Главный цикл *)
 x += RandomVariate[ExponentialDistribution[4]]; (* Разыгрываем прыжок темп. с пола *)
 If[ (* НАЧАЛО ЕСЛИ *)
  x < 37, (* ЕСЛИ не выскочили, то *)
  AppendTo[coord, {t, x}]; (* ТОГДА Пишем в хронику номер наблюдения и темп., *)
  x = 36.6; t++, (* Обнуляем частицу, разыгрываем следующее наблюдение *)
      While[x < 40, (* ИНАЧЕ ЦИКЛ2: будем расти, пока меньше 40 *)
     AppendTo[coord, {t, x}];    (* Пишем в хронику номер наблюдения и темп., *)
        x += RandomVariate[ExponentialDistribution[4]]; (* Снова растём, ибо больны *)
        t++ (* Разыгрываем следующее наблюдение *)
       ]; (*КОНЕЦ ЦИКЛ2, когда стало выше 40 *)
      AppendTo[coord, {t, x}]; (* Пишем точку, которая выше 40 *)
            t++; (* Ждём 1 час, пока без сознания *)
      x = 36.6;  (* Охлаждаем тело *)
      AppendTo[coord, {t, x}]; (* Записываем здоровое наблюдение *)
(* Закомментировать предыдущую строку для последнего пункта задачи *)
      t++; (* Разыгрываем следующее наблюдение *)
(* Закомментировать предыдущую строку для последнего пункта задачи *)
  ]; (* КОНЕЦ ЕСЛИ
 i++ (* Перейти к следующей процедуре проверки
 ]; (* Конец главного цикла
(*Print[TraditionalForm[coord]];*)
(* Предыдущая строка, если снять комментарий, выводит таблицу наблюдений *)
(* Я по ней проверял, не потерялось ли чего *)
 
(* А теперь рисуем все наблюдения *)
Show[
 ListPlot[coord, PlotStyle -> Directive[PointSize[Medium], Black], (* Наблюдения чёрными точками *)
  PlotRange -> All],
 ListLinePlot[{{0, 37}, {t, 37}}, (* Горизонтальные линии пороговых уровней *)
  PlotStyle -> {Thickness[0.002], Red, Dashed}],
  ListLinePlot[{{0, 40}, {t, 40}},
   PlotStyle -> {Thickness[0.002], Red, Dashed}],
 ListLinePlot[{{0, 41}, {t, 41}},
  PlotStyle -> {Thickness[0.002], Red, Dashed}]
 ]

20 наблюдений:

200 наблюдений:

2 000 наблюдений:

20 000 наблюдений:

Посмотрев на скриншоты, можно прийти к выводу (правда, это надо ещё и обосновать), что если температура подскочит выше 42°, то старушка взрывается... шутка! Можно прийти к выводу, что если принять условие пункта 4, то дни старушки сочтены, длина жизни конечна; ведь так много аутлайеров...

Особенность программы: количество циклов становится количеством наблюдений, когда частица в коридоре. В итоге реально всего шагов в симуляторе становится не столько, сколько циклов, а раза в 3,3 больше. Можно было бы реализовать так, чтобы сколько задали, столько и было, но это нужно подумать и что-то основное поменять (наложить проверку не на номер цикла, на t).

P.S. MathJax поддерживает знаки нестрогого неравенства с наклонными чёрточками, кои так характерны для русской типографики: \leqslant и \geqslant. Не забудьте включить AMS в конфигурации MathJax, и будет вам счастье: \(\leqslant\), \(\geqslant\).