Банк рефератов содержит более 364 тысяч рефератов, курсовых и дипломных работ, шпаргалок и докладов по различным дисциплинам: истории, психологии, экономике, менеджменту, философии, праву, экологии. А также изложения, сочинения по литературе, отчеты по практике, топики по английскому.
Полнотекстовый поиск
Всего работ:
364139
Теги названий
Разделы
Авиация и космонавтика (304)
Административное право (123)
Арбитражный процесс (23)
Архитектура (113)
Астрология (4)
Астрономия (4814)
Банковское дело (5227)
Безопасность жизнедеятельности (2616)
Биографии (3423)
Биология (4214)
Биология и химия (1518)
Биржевое дело (68)
Ботаника и сельское хоз-во (2836)
Бухгалтерский учет и аудит (8269)
Валютные отношения (50)
Ветеринария (50)
Военная кафедра (762)
ГДЗ (2)
География (5275)
Геодезия (30)
Геология (1222)
Геополитика (43)
Государство и право (20403)
Гражданское право и процесс (465)
Делопроизводство (19)
Деньги и кредит (108)
ЕГЭ (173)
Естествознание (96)
Журналистика (899)
ЗНО (54)
Зоология (34)
Издательское дело и полиграфия (476)
Инвестиции (106)
Иностранный язык (62791)
Информатика (3562)
Информатика, программирование (6444)
Исторические личности (2165)
История (21319)
История техники (766)
Кибернетика (64)
Коммуникации и связь (3145)
Компьютерные науки (60)
Косметология (17)
Краеведение и этнография (588)
Краткое содержание произведений (1000)
Криминалистика (106)
Криминология (48)
Криптология (3)
Кулинария (1167)
Культура и искусство (8485)
Культурология (537)
Литература : зарубежная (2044)
Литература и русский язык (11657)
Логика (532)
Логистика (21)
Маркетинг (7985)
Математика (3721)
Медицина, здоровье (10549)
Медицинские науки (88)
Международное публичное право (58)
Международное частное право (36)
Международные отношения (2257)
Менеджмент (12491)
Металлургия (91)
Москвоведение (797)
Музыка (1338)
Муниципальное право (24)
Налоги, налогообложение (214)
Наука и техника (1141)
Начертательная геометрия (3)
Оккультизм и уфология (8)
Остальные рефераты (21692)
Педагогика (7850)
Политология (3801)
Право (682)
Право, юриспруденция (2881)
Предпринимательство (475)
Прикладные науки (1)
Промышленность, производство (7100)
Психология (8692)
психология, педагогика (4121)
Радиоэлектроника (443)
Реклама (952)
Религия и мифология (2967)
Риторика (23)
Сексология (748)
Социология (4876)
Статистика (95)
Страхование (107)
Строительные науки (7)
Строительство (2004)
Схемотехника (15)
Таможенная система (663)
Теория государства и права (240)
Теория организации (39)
Теплотехника (25)
Технология (624)
Товароведение (16)
Транспорт (2652)
Трудовое право (136)
Туризм (90)
Уголовное право и процесс (406)
Управление (95)
Управленческие науки (24)
Физика (3462)
Физкультура и спорт (4482)
Философия (7216)
Финансовые науки (4592)
Финансы (5386)
Фотография (3)
Химия (2244)
Хозяйственное право (23)
Цифровые устройства (29)
Экологическое право (35)
Экология (4517)
Экономика (20644)
Экономико-математическое моделирование (666)
Экономическая география (119)
Экономическая теория (2573)
Этика (889)
Юриспруденция (288)
Языковедение (148)
Языкознание, филология (1140)

Реферат: Факультет: Кибернетики Кафедра: Биомедицинская Электроника работа “Расчёт асцитной гепатомы Зайдела” Дисциплина: “Моделирование в медицине”

Название: Факультет: Кибернетики Кафедра: Биомедицинская Электроника работа “Расчёт асцитной гепатомы Зайдела” Дисциплина: “Моделирование в медицине”
Раздел: Остальные рефераты
Тип: реферат Добавлен 07:12:31 29 сентября 2011 Похожие работы
Просмотров: 53 Комментариев: 6 Оценило: 0 человек Средний балл: 0 Оценка: неизвестно     Скачать

МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ ИНСТИТУТ

РАДИОТЕХНИКИ, ЭЛЕКТРОНИКИ И АВТОМАТИКИ (ТУ)

Факультет: Кибернетики

Кафедра: Биомедицинская Электроника

КУРСОВАЯ РАБОТА

“Расчёт асцитной гепатомы Зайдела”

Дисциплина: “Моделирование в медицине”

Работу выполнил студент группы КМ-1-хх:

ФИО _______________________

Проверил преподаватель:

Бабушкина Нина Александровна ______________________

МОСКВА 200х

Содержание

Титульный лист…………………………………………………………………………………1

Содержание……………………………………………………………………………………...2

Введение…………………………………………………………………………………………3

1. Заболевание…………………………………………………………………………………..3

2. Аппроксимация исходных данных…………………………………………………………..3

2.1. Подбор аппроксимирующей функции…………………………………………………….3

2.2. Подбор коэффициентов функции ПК……………………………………………………..3

2.3. Выбор функции, обеспечивающей наилучшую аппроксимацию исходных данных…..4

2.4. Аналитический расчёт коэффициентов методом наименьших квадратов……………...5

2.4.1. Метод наименьших квадратов…………………………………………………………...5

2.4.2. Расчёт коэффициентов степенной функции……………………………………………6

2.5. Расчёт коэффициентов итерационным методом Ньютона……………………………….7

2.5.1 Алгоритм метода Ньютона………………………………………………………………..7

2.5.2 Расчёт коэффициентов степенной функции методом Ньютона………………………..8

2.6. Результаты………………………………………………………………………………….10

3. Расчёт биологических параметров………………………………………………………….10

3.1. Время жизни организма без лечения и запас жизненных сил…………………………..10

3.2. Дозовая зависимость…………………………………...…………………………………..11

3.3. Расчёт времени жизни организма после курса лечения…………………………………12

Результаты………………………………………………………………………………………13

Список литературы……………………………………………………………………………..14

Введение

Необходимо рассчитать время жизни организма с асцитной гепатомой Зайдела при следующих исходных данных:

Количество введения доз препарата: 2 раза

Вводимая доза: D=0.8 МПД

Задержка введения доз:суток

1. Заболевание

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

Асцитная гепатома Зайдела (ascite Zajdel hepatoma) – заболевание животных. Обычно её прививают специально в исследовательских целях.

2. Аппроксимация исходных данных.

2.1. Подбор аппроксимирующей функции

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

Экспоненциальная функция: (2.1)

Степенная функция: (2.2)

2.2. Подбор коэффициентов функции ПК

Подбор коэффициентов производился в пакете MathCAD 2001 с помощью функции genfit(), которая работает по итерационному методу градиентного спуска.

Экспоненциальная:

(2.3)

Степенная функция:

(2.4)

Рис. 2.1. Аппроксимация исходных данных экспоненциальной и степенной функциями

2.3. Выбор функции, обеспечивающей наилучшую аппроксимацию исходных данных.

Чтобы выбрать функцию, обеспечивающую наилучшую аппроксимацию исходных данных, необходимо просчитать коэффициент корреляции Пирсона R.

, где (2.5)

N(t) – исходные данные
G(t) – аппроксимирующая функция

(2.6) - средние значения популяции опухоли, рассчитанное по исходным данным.

(2.7) - средние значения популяции опухоли, рассчитанное по аппроксимированным данным..

Коэффициент корреляции показывает статистическую связь между выборками. Он может принимать значения между -1 и +1, причём если значение находится ближе к 1, то это означает наличие сильной связи, а если ближе к 0, то слабой. Если коэффициент корреляции отрицательный, это означает наличие противоположной связи: чем выше значение одной переменной, тем ниже значение другой.

Для полученной экспоненциальной функции (2.3) коэффициент корреляции R1 по (2.5) равен:

R1 = 0.99654

Для полученной степенной функции (2.4) коэффициент корреляции R2 по (2.5) равен:

R2 = 0.99907

R2 > R1 , следовательно, степенная функция лучше аппроксимирует исходные данные.

2.4. Аналитический расчёт коэффициентов методом наименьших квадратов

2.4.1. Метод наименьших квадратов

Метод наименьших квадратов — один из методов теории ошибок для оценки неизвестных величин по результатам измерений, содержащим случайные ошибки. Метод наименьших квадратов применяется также для приближенного представления заданной функции другими (более простыми) функциями и часто оказывается полезным при обработке наблюдений.[1]

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

Пусть дано решить систему уравнений, число которых более числа неизвестных x, у, z…

ax + by + cz… + n = 0

a1 x + b1 y + c1 z… + n1 = 0

a2 x + b2 y + c2 z… + n2 = 0 (2.8)

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

[aa] = a1 a1 + a2 a2 +…

[ab] = a1 b1 + a2 b2 +…

[ac] = a1 c1 + a2 c2 +… (2.9)

[bb] = b1 b1 + b2 b2 +…

[bc] = b1 c1 + b2 c2 +…

то нормальные уравнения представятся в следующем простом виде:

[aa]x + [ab]y + [ac]z +… [an] = 0

[ab]x + [bb]y + [bc]z +… [bn] = 0 (2.10)

[ac]x + [bc]y + [cc]z +… [cn] = 0

Число этих уравнений равно числу неизвестных и легко решается.

2.4.2. Расчёт коэффициентов степенной функции

Линеаризуем степенную функцию:

(2.11)

Линеаризуем исходные данные:

В соответствии с п.2.4.1 переназначим переменные:

Тогда исходные данные можно записать в форме (2.8):


a1 x + y + n1 = 0

a2 x + y + n2 = 0

a3 x + y + n3 = 0

a4 x + y + n4 = 0

a5 x + y + n5 = 0

0.69x + y - 1.1 = 0

1.1x + y - 2.83 = 0

1.39x + y – 3.9 = 0

1.61x + y - 5 = 0

1.79x + y – 5.7 = 0

[aa] = a1 a1 + a2 a2 +…+ a5 a5 = 9.41

[ab] = a1 b1 + a2 b2 +…+ a5 b5 = 6.58

[an] = a1 n1 + a2 n2 +… + a5 n5 = -27.62

[bb] = b1 b1 + b2 b2 +…+ b5 b5 = 5

[bn] = b1 n1 + b2 n2 +…+ b5 n5 = -18.63

Следовательно, нормальные уравнения записываются по форме (2.10) так:

[aa]x + [ab]y + [an] = 0

[ab]x + [bb]y + [bn] = 0

9.41x + 6.58y – 27.62 = 0

6.58x + 5y – 18.63 = 0

Решая эту систему, получаем:

Следовательно, аппроксимирующая функция будет иметь вид:

2.5. Расчёт коэффициентов итерационным методом Ньютона

2.5.1 Алгоритм метода Ньютона.

Метод Ньютона представляет собой метод наискорейшего спуска с шагом, длина которого зависит от свойств минимизируемой функции. Метод Ньютона основан на квадратической аппроксимации минимизируемой функции в окрестности точки x(k) , где (k) – номер итерации Минимум квадратической функции легко найти, приравнивая ее градиент нулю. Можно сразу же вычислить положение экстремума и выбрать его в качестве следующего приближения к точке минимума. Новая итерация вычисляется по формуле:

(2.12)

Пусть f(x) - минимизируемая фукнция с векторным аргументом . Алгоритм наискорейшего спуска реализует итерационную процедуру движения к минимуму из произвольно выбранной точки начального приближения в направлении наиболее сильного уменьшения функции, определенном в окрестности текущего значения аргумента минимизируемой функции. Такое направление противоположно направлению, задаваемому вектором градиента минимизируемой функции f(x):

(2.13)

Вычисляя точку нового приближения по формуле (2.12) и разлагая f(x(k+1) ) в ряд Тейлора, получим формулу квадратической аппроксимации fкв (x(k+1) ):

, где

(2.14)

- матрица вторых производных:

(2.15)

Условие минимума fкв (x(k+1) ) по . Вычислим градиент из (2.14):

(2.16)

Для учета фактических особенностей минимизируемой функции будем использовать в (2.16) значения градиента и матрицы вторых производных, вычисленных не по аппроксимирующей fкв (x), а непосредственно по минимизируемой функции f(x). Заменяя fкв (x) в (2.16), найдем длину шага

(2.17)

Итак, последовательность вычислений для реализации алгоритма метода Ньютона:

  1. Произвольно задать точку начального приближения x(0)
  2. В цикле по номеру итерации k=0,1… вычислить:
    1. Значение вектора градиента по формуле (2.13)
    2. Значение матрицы вторых производных по формуле (2.15)
    3. Значение матрицы, обратной матрице вторых производных
    4. Значение шага по формуле (2.17)
    5. Новое значение приближения x(0) по формуле (2.12)
  3. Закончить итерационный процесс при достижении нужного приближения.[2]

2.5.2 Расчёт коэффициентов степенной функции методом Ньютона

За минимизируемую функцию возьмём среднеквадратичное отклонение (СКО) между исходными точками и аппроксимированными аналитической функцией:

(2.18)

За начальные приближения выберем:

Для случая со степенной функцией формулы (2.13) и (2.15) имеют следующий вид:

Итерации:

  1. Первая итерация.

  1. Вторая итерация.

  1. Третяя итерация.

  1. Четвёртая итерация.

  1. Пятая итерация.

  1. Шестая итерация.

Итерационный цикл закончен, т.к. результат вычисления коэффициентов совпал с рассчитанными на ПК (п.2.2).

2.6. Результаты

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

3. Расчёт биологических параметров.

3.1. Время жизни организма без лечения и запас жизненных сил.

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

Запас жизненных сил определяют как площадь под аналитической кривой от начала заболевания до летального исхода.

(3.1)

Рис. 3.1. Иллюстрация по запасу жизненных сил.

3.2. Дозовая зависимость.

Рис. 3.2. График дозовой зависимости

Задержку роста опухоли определяют по данным дозовой зависимости.

Вводимая доза: D=0.8 МПД

Задержка введения доз:суток

Задержка роста опухоли:

(3.2)

3.3. Расчёт времени жизни организма после курса лечения.

Величина запаса жизненных сил не меняется со временем и является величиной постоянной. На основании этого факта можно произвести расчёт времени жизни организма после курса лечения. Поскольку при введении дозы препарата происходит задержка роста опухоли, представим этот процесс в виде аналитической кусочно-прерывной функции из трёх интервалов: первый – до первого введения дозы, второй и третий - после соответствующего введения дозы с задержкой в 5 суток (Рис.3.3).

(3.3)

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

(3.4)

Рис. 3.3. График роста опухоли после лечения.

Необходимо проверить, не умрёт ли организм до окончания цикла лечения. Для этого просчитаем расход запаса жизненных сил до второго введения:

,

следовательно организм не умрёт до завершения цикла лечения.

По формуле (3.4) рассчитываем время жизни после цикла лечения:

Результаты

Развитие опухоли лучше всего аппроксимирует степенная функция. При двукратном введении дозы уровнем 0.8 МПД, организм сможет прожить 10.5 суток.

Рис. 4.1. Смоделированный график развития опухоли до летального исхода

Список литературы

  1. Пыльнов Ю.В. Регрессионный анализ полиномиальных моделей. – М.: МИРЭА, 1994, 56 с.
  2. Ван дер Варден Б.Л. Математическая статистика. – М: Эдиториал УРСС, 2006, 435 c.
Оценить/Добавить комментарий
Имя
Оценка
Комментарии:
Привет студентам) если возникают трудности с любой работой (от реферата и контрольных до диплома), можете обратиться на FAST-REFERAT.RU , я там обычно заказываю, все качественно и в срок) в любом случае попробуйте, за спрос денег не берут)
Olya17:55:50 01 сентября 2019
.
.17:55:49 01 сентября 2019
.
.17:55:48 01 сентября 2019
.
.17:55:48 01 сентября 2019
.
.17:55:47 01 сентября 2019

Смотреть все комментарии (6)
Работы, похожие на Реферат: Факультет: Кибернетики Кафедра: Биомедицинская Электроника работа “Расчёт асцитной гепатомы Зайдела” Дисциплина: “Моделирование в медицине”

Назад
Меню
Главная
Рефераты
Благодарности
Опрос
Станете ли вы заказывать работу за деньги, если не найдете ее в Интернете?

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



Результаты(258798)
Комментарии (3487)
Copyright © 2005-2020 BestReferat.ru support@bestreferat.ru реклама на сайте

Рейтинг@Mail.ru