Litres Baner
Название книги:

Математические модели в естественнонаучном образовании. Том I

Автор:
Денис Владимирович Соломатин
Математические модели в естественнонаучном образовании. Том I

000

ОтложитьЧитал

Шрифт:
-100%+

Доброй памяти Л.М.Мартынова, величайшего алгебраиста и криптографа современности, посвящается.

Предисловие

Связь между техническими и гуманитарными науками становится всё теснее. Классические задачи, такие как моделирование популяций и заболеваний, сменяются новыми проблемами моделирования машинного обучения, усложняющимися по мере накопления эмпирических данных, и вновь делают математику перспективной сферой человеческой деятельности. Не секрет, что естественный интеллект на сегодняшний день превосходит искусственного лишь эмоциональной составляющей и математической интуицией. Именно поэтому предполагается, что данная сфера будет по-прежнему одной из самых быстрорастущих.

Мы считаем, что межпредметные связи должны проявляться на всех уровнях математического образования. Студенты-математики получают определённый опыт и несомненную пользу, наблюдая приложения математики из неожиданных областей, открытых для них, извлекают выгоду из изучения того, как математические инструменты могут помочь им реализовать свои собственные проекты. Образ дидактики как нематематической науки, который сохраняется среди многих студентов педагогических колледжей, оказывает медвежью услугу тем, кто придерживается такой точки зрения. Настоящая монография является попыткой представить некоторые существенные темы математического моделирования на принципиально новом языке, адаптированном к математическому образованию. Надеемся, что это может мотивировать некоторых студентов-математиков педагогических специальностей продолжать свои математические исследования за пределами традиционного уровня. Такие студенты, как правило, имеют сильный интерес к математике и базовую математическую подготовку, достаточную для самостоятельного изучения её приложений. Таким образом, мы не предполагаем никакой дополнительной подготовки в области высшей математики за пределами классического курса; моделирование с помощью разностных уравнений позволяет свести к минимуму предварительную подготовку. Математические темы, обычно изучаемые на факультативных курсах, вводятся по мере необходимости для моделирования и последующего анализа полученных моделей. Несмотря на такой план изложения, мы знаем, что многие студенты изучают курс высшей математики и, возможно, другие специализированные курсы. Поэтому без колебаний включим вопросы для самопроверки и задачи для самостоятельного решения (они выделены курсивом отдельно), которые могут пригодиться тем, у кого есть дополнительная математическая подготовка. По нашему опыту работы в группах физико-математического профиля, студенты с фундаментальным образованием нашли для себя здесь много нового. Большая часть материала также апеллирует к знаниям студентов и по другим дисциплинам, которые просто интересны и занимательны сами по себе. Следовательно, монография может быть продуктивно использована как для проведения аудиторных занятий, так и для самостоятельного изучения на разных уровнях.

Нет уверенности в том, что к моменту издания будут охвачены все передовые направления столь бурно развивающейся отрасли, но стоит предположить, как поверхностное освещение какой-либо случайной темы вдохновит студентов на самостоятельное погружение в заинтересовавшую их область. В частности, когда вводятся некоторые модели, оставляем их анализ в качестве упражнений. Хотя подобное нельзя считать эффективным способом энциклопедичного изложения материала, мы надеемся, что это приведет к более глубокому пониманию изучаемого.

Поскольку компьютерные эксперименты с моделями бывают достаточно информативными, текст насыщен сериями программ. Простой интерфейс MATLAB, его широкая доступность как в профессиональной, так и в студенческой версиях, поддержка числовых и символьных вычислений, позволяет минимизировать усилия, затрачиваемые на решение типовых задач математического моделирования. Хотя компьютер и является инструментом, которым студенты должны уметь пользоваться, его использование ни в коем случае не является панацеей.

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

По всему тексту разбросаны контрольные вопросы, помеченные знаком «?». Они предназначены для дополнительного стимулирования сознания и предотвращения пассивного чтения. Ответы окажутся вполне очевидными после небольшого самостоятельного размышления, либо соответствующий вопрос будет обсуждаться в тексте позднее. Если находите такие вопросы раздражающими, пожалуйста, не стесняйтесь игнорировать их.

В тексте больше фактического материала, чем можно было бы охватить за полноценный семестр, предлагая преподавателям множество вариантов изучения тем. Например, темы глав 1, 2, 3 и 7, пожалуй, являются наиболее привычными для академического курса математического моделирования и уравнения математической физики, охватывающие классические модели, как линейные, так и нелинейные. Главы 4 и 5 предлагают студентам введение в новые темы эволюционного моделирования и построения дерева решений, которые являются перспективными и полезными при построениях различных классификаторов. Глава 6, посвященная генетическим алгоритмам, дает представление о другой области, в которой математика и биология оказались тесно переплетены, вдохновляясь идеями друг от друга. Глава 8 и Приложение содержат краткое введение в основные инструменты подгонки кривых и статистики.

В главе 1 вводится понятия динамического моделирования с помощью разностных уравнений от одной переменной, включая ключевые понятия равновесия, линеаризации и стабильности. В главе 2 эксплуатируем матричную алгебру, вычисление собственных чисел и векторов, с помощью линейных моделей от двух переменных. Эти главы являются основой для изложения всего последующего материала.

Экскурс в теорию вероятностей появляется в двух разделах главы 4, чтобы смоделировать эволюционное развитие математического образования, а затем расширяется в главе 6 на применение генетических алгоритмов. Глава 5, которая имеет алгоритмический привкус, отличный от остальной части текста, частично зависит от формул расстояния, полученных в главе 4. Трактовка моделей в главе 8, естественным образом зависит от введённой в главе 3 модели взаимодействующих факторов.

Разработка этой монографии началась давно, в течение нескольких лет появлялись краткие версии нескольких глав. После многих дополнительных пересмотров конспекты курса достигли критической массы, и публикация их для использования другими преподавателями оказалась лишь вопросом времени.

Многие студенты помогали, как в качестве критиков, так и в качестве экспериментаторов, пытаясь решать предлагаемые задачи, вычитывая текст, задавая много вопросов. Несмотря на все старания, некоторые ошибки наверняка просочились, пожалуйста, сообщите автору обо всем, что сочтёте неверным.

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

MATLAB по сути является полноценным языком программирования с отличными графическими возможностями. Мы воспользовались этими функциями, чтобы предоставить несколько программ, упрощая изучение математических моделей для начинающих. В упражнениях и проектах используются некоторые программы и файлы данных, приводимые ниже.

Готовые тексты программ приводятся в целях минимизации необходимого багажа знаний синтаксиса MATLAB. Чтобы запустить большинство программ ниже, просто скопируйте их из электронного варианта книги или введите в окне редактора скриптов MATLAB.

Можно просто сохранить текст программы в файл с расширением m. Для запуска такого скрипта в окне водится имя запускаемого файла. После запуска будет задан ряд вопросов о конструируемых моделях и их параметрах. Команда help имя_файла также предоставляет краткое описание функций программы, получаемое из комментариев в первых строках программы. Поскольку m-файлы являются текстовыми файлами, они могут быть прочитаны и изменены любым заинтересованным пользователем.

Некоторые m-файлы определяют функции, которые принимают аргументы. Например, такая команда, как compseq(seq1,seq2), запускает программу compseq.m для сравнения двух последовательностей ДНК seq1 и seq2. Набрав help compseq, можно увидеть объяснение синтаксиса такой функции. A mat-файл содержит данные, доступ к которым возможен только из MATLAB. Чтобы загрузить такой файл, скажем, seqdata.mat, введите load seqdata. Имена всех новых переменных, которые будут созданы, можно увидеть вводя команду who, в то время как значения, хранящиеся в этих переменных, можно увидеть просто вводя имя переменной.

Некоторые файлы данных могут быть предоставлены в виде m-файлов, тогда вспомогательные комментарии и пояснения сохраняются вместе с данными. Для них запуск m-файла создает переменные, так же, как и загрузка mat-файла. Комментарии можно прочитать с помощью любого текстового редактора.

В ходе выполнения задач для самостоятельного решения предлагается использовать следующие файлы скриптов MATLAB, доступных из открытых источников:

 aidsdata.m – содержит данные числа случаев синдрома приобретенного иммунодефицита;

 cobweb.m и cobweb2.m – рисуют графики с паутинной диаграммой для итераций модели с одной популяцией; первая программа оставляет все рисуемые линии, а вторая программа постепенно стирает их;

 

 compseq.m – функция сравнивает две последовательности ДНК, получая частотную таблицу количества фрагментов с каждой из возможных базовых комбинаций;

distances.m – функция вычисляет расстояния Джукса-Кантора, 2-параметрическое расстояние Кимуры и логарифмическое расстояния между всеми парами в коллекции последовательностей ДНК;

 distJC.m, distK2.m и distLD.m – функции вычисляют расстояние Джукса-Кантора, 2-параметрическое расстояние Кимуры и логарифмическое расстояние для одной пары последовательностей, описываемых частотным массивом сайтов ДНК с каждой комбинацией оснований;

 flhivdata.m – содержит последовательности ДНК гена оболочки вируса иммунодефицита человека из «случая стоматолога во Флориде»;

 genemap.m – моделирует данные тестового скрещивания для проекта генетического картирования, используя гены мухи или мыши;

 genesim.m – производит временной график частоты аллелей гена в популяции фиксированного размера; относительные значения приспособленности для генотипов могут быть установлены для моделирования естественного отбора;

 informative.m – функция находит участки в выровненных последовательностях ДНК, которые информативны для метода максимальной экономии;

 longterm.m – рисует диаграмму бифуркации для модели с одной популяцией, показывая долгосрочное поведение по мере изменения значения одного параметра;

 markovJC.m и markovK2.m – эти функции осуществляют получение марковской матрицы Джукса-Кантора или 2-параметрической модели Кимуры с заданными значениями параметров;

 mutate.m и mutatef.m – моделирует мутации последовательности ДНК по марковской модели замещения оснований; вторая программа является функциональной версией первой;

 nj.m – функция реализует алгоритм присоединения соседей для построения дерева из массива расстояний;

 onepop.m – отображает графики итераций модели с одной популяцией;

 primatedata.m – содержит последовательности митохондриальной ДНК из 12 приматов, а также вычисленные расстояния между ними;

 seqdata.mat – содержит смоделированные данные последовательности ДНК;

 seqgen.m – функция генерирует последовательности ДНК с заданной длиной и распределением оснований;

 sir.m – отображает итерации эпидемиологической модели SIR, включая графики временной и фазовой плоскости;

 twopop.m – отображает итерации 2-популяционной модели, включая графики временной и фазовой плоскости.

Глава 1. Динамическое моделирование разностными уравнениями

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

Как мы можем понять такие сложные системы, как те, которые возникают в социальных науках? Как мы можем проверить, достаточно ли нашего предполагаемого понимания ключевых процессов, чтобы описать, как ведет себя система? Математический язык предназначен для точного описания, и поэтому описание сложных систем часто требует математической модели.

В этой главе мы рассмотрим некоторые способы, которыми математика используется для моделирования динамических процессов в обучении математике. Простые формулы связывают, например, количество абитуриентов в определенном году с выпускниками последующих лет. Мы учимся понимать последствия, которые можно прогнозировать, составляя уравнение, средствами математического анализа, при этом наша формализация может быть проверена эмпирическими наблюдениями. Хотя многие из моделей, которые мы рассматриваем, могут на первый взгляд показаться грубыми упрощениями, их сила в простоте. Чем проще модель, тем яснее становятся предсказываемые её последствия исходя из самых базовых предположений.

Начнем с того, что сосредоточимся на моделировании того, как количество выпускников физико-математических классов растёт или сокращается с течением времени. Поскольку математические модели должны основываться на вопросах, вот несколько вопросов, которые следует учитывать: почему число выпускников иногда растёт, а иногда сокращается? Должны ли объемы выпусков вырасти до такой степени, что они станут неустойчиво большими, а затем сойдут до нуля? Если нет, то должно ли количество выпускников достичь некоторого равновесия? Если равновесие существует, какие факторы ответственны за него? Является ли такое равновесие настолько тонким, что любое нарушение может положить ему конец? Что определяет, следует ли данная тенденция одному из этих курсов или другому?

Начнём разбирать перечисленные вопросы с помощью самой простой математической модели изменяющейся численности населения.

1.1. Мальтузианская модель

Предположим, мы выращиваем не будущих математиков, а популяцию какого-то организма, скажем, мух, в лаборатории. Представляется разумным, что в любой данный день численность населения будет меняться из-за новых рождений, так что оно увеличивается за счет добавления определенной доли f от имеющегося населения. При этом часть d от имеющегося населения погибнет, условно, как бы цинично это не звучало, но многие профессиональные математики после выпуска вынуждены работать не по специальности, что смерти подобно.

Рассмотрим простейшую прикладную модель, которую предложил Томас Мальтус в своём очерке 1798 года о принципе народонаселения, неоднократно подвергавшемся всесторонней критике. Если люди живут в течение 70 лет, то мы ожидаем, что из большой популяции примерно 1/70 населения будет умирать каждый год; таким образом, . Если, с другой стороны, мы предположим, что на каждые сто человек приходится около четырех рождений в год, мы имеем . Обратите внимание, что в этом случае мы выбрали год в качестве единиц времени.

Вопросы для самопроверки:

– Объясните, почему для любой популяции  должно быть в диапазоне от 0 до 1.  Что будет означать ?  Что будет означать ?

– Объясните, почему  должно быть не менее 0, но может быть больше 1. Можете ли вы назвать реальные популяции (при должном выборе единицы времени), для которых  будет больше 1?

– Используя годы в качестве единицы времени, какие значения f и d будут уместны для моделирования числа выпускников естественно-научного профиля? Гуманитарного? Социально-экономического? Технологического и универсального?

Чтобы смоделировать значения P сфокусируемся на следующем за P изменении численности. Формально . Это означает, что, учитывая текущее значение , скажем, , а также  и , например,  и , можно предсказать изменение . Таким образом, в начале следующего временного периода суммарная численность составляет .

Введём несколько вспомогательных обозначений для упрощения восприятия математической модели. Пусть  – размер популяции, измеренный в момент времени , тогда  это приращение или изменение численности между последовательными моментами времени.

Ясно, что  зависит от , поэтому можно встретить подстрочный индекс  рядом с , так как для разных значений  приращение  оказывается разным. Тем не менее, этот индекс не редко пропускают.

Теперь то, что нас в конечном итоге волнует, это понимание динамики популяции , а не только приращения . Но . Объединив константы вместе, обозначив за , модель стала гораздо проще: .

Популяризаторы науки часто называют константу  конечной скоростью роста населения. (Слово «конечный» используется, чтобы отличить это число от любого вида мгновенной скорости, которая включала бы производную, как вы знаете из курса дифференциального исчисления. Для значений , , и  использованных ранее, вся модель теперь имеет вид , где . Первое уравнение, выражающее  через , называется разностным уравнением, а второе, задающее , является его начальным условием.  С этими двумя уравнениями легко составить таблицу значений численности  с течением времени, как в таблице 1.1.

Таблица 1.1. Рост популяции по простой модели

Момент времени         Численность

0                                          500

1                                          (1. 07)500 = 535

2                                          (1. 07)2500 = 572.45

3                                          (1. 07)3500 ≈ 612.52

…                                         …

По закономерностям в таблице 1.1 легко перейти от рекуррентного соотношения для  к замкнутой форме записи, чтобы осталась только зависимость от  в явном виде: . На этой модели теперь легко предсказать численность популяции в любое время.

Может показаться странным называть  разностным уравнением, когда разность  там не появляется. Однако уравнения  и  эквивалентны, поэтому любое из них разумно определять одним и тем же термином.

Пример. Предположим, что система математического образования имеет очень жесткие ограничения на целевые цифры приёма в ВУЗы (что вполне реалистично на просторах СНГ), по которым каждый год выпускается 200 молодых специалистов и все сотрудники пенсионного возраста уходят на заслуженный отдых. После того, как состоялся очередной выпуск, только 3% остаются работать по специальности, чтобы связать свою профессиональную деятельность с математикой, остальные либо эмигрируют, либо находят выше оплачиваемую работу. Чтобы написать разностное уравнение в этой системе, где будем измерять  в поколениях, нужно просто заметить, что уровень «смертности» равен , в то время как эффективная «плодовитость» системы равна . Следовательно, .

Вопросы для самопроверки:

– Будет ли общая численность математиков расти, а не уменьшаться при таких условиях?

– Предположим, вы не знаете эффективной «плодовитости», но знаете, что численность  стабильна (неизменна) с течением времени. Какой должна быть ? (Подсказка: поймите, что такое , если численность стабильна?) Если каждый год выпускается 200 молодых специалистов, какая их часть должна оставаться в системе и обучать математиков следующего поколения?

 

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

Вопросы для самопроверки:

– Можете ли вы представить себе обстоятельства, при которых игнорирование уменьшения числа профессионалов той или оной области было бы хорошей идеей?

Так что же такое разностное уравнение? Теперь, когда увидели разностное уравнение на примере, можно попытаться дать строгое определение: разностное уравнение – это формула, выражающая значения некоторой величины  в терминах предыдущих значений . Таким образом, если  является какой-либо функцией, то  называется разностным уравнением. В предыдущем примере использовалась , но часто  будет более сложным.

Изучая разностные уравнения и их приложения, рассмотрим два основных вопроса: 1) Как найти подходящее разностное уравнение для моделирования ситуации? 2) Как понять поведение модели разностных уравнений после того, как её нашли?

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

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

Задачи для самостоятельного решения:

1.1.1. Популяция изначально составляла 100 особей, но из-за комбинированного воздействия рождений и смертей она утраивается каждый час.

а. Составьте таблицу численности популяции для  пробегающего значения от 0 до 5, где  измеряется в часах.

б. Приведите два уравнения, моделирующих рост популяции, сначала путем выражения  через , а затем выразив  через .

в. Что можно сказать об уровнях рождаемости и смертности среди населения вашей страны? Земного шара?

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

а. Запишите уравнение, моделирующее эту ситуацию. Уточнив, сколько реального времени представлено шагом 1 в параметре  и каково было начальное количество школ в период новейшей истории.

б. Заполните таблицу и нарисуйте график числа школ в зависимости от .

в. Сопоставьте полученные результаты с официальными данными Росстата. Это соответствует вашей модели? Какие выводы и/или вопросы это вызывает?

1.1.3. С помощью ручного калькулятора составьте таблицу значений численности населения выбирая  в диапазоне от 0 до 6 для следующих моделей. Затем отобразите табличные значения на графике.

а. ,

б. ,

в. ,

1.1.4. Повторите решение задачи 1.1.3(а) с помощью MATLAB, введя последовательность команд, например:

p=1

x=p

p=1.3*p

x=[x p]

p=1.3*p

x=[x p]

Возврат к предыдущим командам для их повторения можно осуществлять нажатием клавиши "↑". Объясните, как это работает. Теперь повторите решение с использованием цикла, например:

p=1

x=1

for i=1:10

p=1.3*p

x=[x p]

end

Отступ не является обязательным, но помогает сделать цикл for-end понятнее для чтения. Объясните, как это работает. Визуализируйте полученные данные на графике с помощью команды:

plot([0:10],x)

1.1.5. Для модели, указанной в задаче 1.1.3 а), сколько времени должно пройти, прежде чем популяция превысит 10, превысит 100 и превысит 1 000? Используйте MATLAB, чтобы вычислить это экспериментальным путём, а затем вычислите аналитически, используя логарифмирование и тот факт, что . Обнаруживается ли закономерность в изменениях вычисленной продолжительности? Объясните, когда и почему значение стабилизируется.

1.1.6. Если бы данные в таблице 1.2 о численности докторов физико-математических наук были собраны по десятилетиям с момента основания института математики, соответствовали бы они геометрической модели? Будет ли численность соответствовать геометрической модели хотя бы в некотором временном интервале? Объясните наблюдаемое явление.

Таблица 1.2. Численность учёных в стране (сотни)


0 1             2             3             4             5             6             7             8             9             10



           1,94       3,04       4,62       6,72       9,26       11,88     14,08     15,52     16,26     16,60     16,72

1.1.7. Заполните пропуски:

а. Модели  и  представляют растущие значения, когда  – любое число в диапазоне _______, а  – любое число в диапазоне _______.

б. Модели  и  представляют уменьшающиеся значения, когда  – любое число в диапазоне _______, а  – любое число в диапазоне _______.

в. Модели  и  представляют стабильные значения, когда  – любое число в диапазоне _______ и когда  – любое число в диапазоне _______.

1.1.8. Объясните, почему модель  не может иметь смысла для описания численности популяции, когда .

1.1.9. Предположим, что популяция описывается моделью  и . Найдите  для .

1.1.10. Говорят, что модель имеет устойчивое состояние или точку равновесия при  если всякий раз, когда , имеем .

а. Перефразируйте определение следующим образом: модель имеет устойчивое состояние при  если всякий раз, когда , имеем  .

б. Перефразируйте определение неформально: модель имеет устойчивое состояние , если ___.

в. Может ли модель, описываемая равенством  иметь устойчивое состояние? Объясните почему.

1.1.11. Объясните, почему модель  приводит к формуле .

1.1.12. Предположим, что на численность определенного населения влияют только рождение, смерть, иммиграция и эмиграция, каждая из которых происходит ежегодно в размере, прямо пропорциональном численности населения. То есть, если население составляет , то в течение периода времени в 1 год число рождений составляет , число смертей , число иммигрантов равно , а число эмигрантов равно , для некоторых , ,  и . Покажите, что популяция все еще может быть смоделирована равенством  и выведите формулу для вычисления .

1.1.13. Как хорошо известно лимнологам и океанографам, количество солнечного света, проникающего на различные глубины воды, может сильно повлиять на численность живущих там организмов. Предположим, что вода имеет равномерную мутность, а количество обитателей на каждом метре в глубину пропорционально количеству поступающего света.

а. Объясните, почему это приводит к модели вида , где  обозначает количество света, проникшего на глубину  метров.

б. В каком диапазоне должны находиться параметры этой модели, чтобы иметь физический смысл?

в. При  и  постройте график  для .

г. Применима ли аналогичная модель к фильтрации света через полог леса? Применимо ли там предположение о «равномерной мутности»?

1.1.14. В таблице 1.3 приведены данные о численности обучающихся физмат школ.

а. Изобразите данные на графике. Соответствуют ли эти данные геометрической модели роста? Объясните почему да или почему нет, используя графические и численные методы оценки. Можете ли придумать факторы, которые приведут к отклонению от геометрической модели?

б. Используя данные только за 1980 и 1985 годы для оценки скорости роста геометрической модели, посмотрите, насколько хорошо результаты модели согласуются с данными последующих лет.

в. Вместо того, чтобы просто использовать данные 1980 и 1985 годов для оценки показателя роста числа школьников, найдите способ использовать все данные, чтобы получить то, что (предположительно) должно быть лучшей геометрической моделью. Проявите творчество. Есть несколько разумных подходов. Соответствует ли ваша новая модель данным лучше, чем модель из части (б)?

Таблица 1.3. Оценки числа школьников

Год        Численность школьников (в 1 000 человек)

1980                     213,260

1985                     231,658

1990                     245,976

1995                     254,504

2000                     263,368

2005                     263,952

2010                     302,690

2015                     328,602

2020                     359,980

1.1.15. Предположим, что популяция моделируется уравнением , где  измеряется в единицах. Если решим измерить численность популяции в тысячах единиц, обозначив это число за , то уравнение, моделирующее популяцию, могло измениться. Объясните, почему модель по-прежнему будет простой . Подсказка: обратите внимание на то, что .

1.1.16. В данной задаче исследуем, как изменится модель, если изменить количество времени, представленное приращением переменной  на единицу. Важно отметить, что эта ситуация не всегда имеет биологический смысл. Например, для организмов, таких как многие насекомые, поколения не перекрываются. Дрозофилы не воспитывают себе преемников. Но время их размножения имеет регулярное распределение, поэтому использование приращения времени меньшее, чем промежуток между двумя последовательными временами рождения, было бы бессмысленным. Однако для более сложных организмов, таких как люди, с перекрывающимися поколениями и практически непрерывным размножением, нет естественного ограничения на выбор значения приращения времени. Таким образом, популяции иногда моделируются с «бесконечно малым» приращением времени (т.е. дифференциальными уравнениями, а не разностными). Эта ситуация иллюстрирует связь между двумя типами моделей: дискретная и континуальная.


Издательство:
Автор
Поделиться: