Расчет запаса жизненных сил. Время жизни после лечения

 

План работы

1. Введение

2. Задание

3. Аппроксимация данных. Построение модели

3.1 Выбор функции

3.2 Аппроксимация в MATLAB

3.3 Выбор оптимальной функции

3.4 Регрессионный анализ экспоненциальной функции

3.5 Регрессионный анализ степенной функции

4. Расчет запаса жизненных сил. Время жизни после лечения

5. Заключение

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

Приложение

1. Введение

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

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

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

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

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

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

2. Задание

Цель работы: построить модель лечения солидной саркомы в программе MATLAB7; найти время жизни существа после лечения, ориентируясь на время жизни объекта до лечения.

Основные задачи: найти функцию с минимальным СКО исходя из заданного графика (Рис.1); найти минимальные значения параметров a и b итоговой функции; вычислить запас жизненых сил объекта; включить в программу дозовую зависимость жизни объекта.

Рис. 1. Экспериментальная кривая роста солидной саркомы

Заданные экспериментальные данные:

V (t), см31.523.53.74.85788.29.513.51515.518t, сутки6.57910111212.81415.21717.8202122

Дополнительные данные:

D=0.8, где D - водимая доза.

T (ж) =31, где T (ж) - время жизни или количество прожитых дней.

t1=7, где t1 - начало лечения или введение первой дозы.=3, где n - общее количество вводимых доз/инъекций.=6, где td - промежуток или количество дней между вводимыми дозами.

E=5, где E - эффект при максимально переносимой дозе.

3. Аппроксимация данных. Построение модели

3.1 Выбор функции

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

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

(2.1)

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

(2.2)

3.2 Аппроксимация в MATLAB

Нахождение коэффициентов производился в пакете MATLAB 7.0 с помощью функций fminsearch (). (см. приложение) x = fminsearch (fun,x0,options,P1,P2,.) - передает зависимые от задачи параметры P1, P2 и т.д. непосредственно в функцию fun. Осуществляет безусловную минимизацию функций. Если опции не определены, используется опция = [] как структурный ноль.

Расчеты в MatLab дали следующие результаты:

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

;

0.820

? = 0.1185

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

;

0.075

? =0.029

? =2.065

Рис. 2. Динамические кривые аппроксимации экспериментальных данных роста опухоли двумя функциями.

3.3 Выбор оптимальной функции

Для выбора функции с большей корреляцией с исходными данными, необходимо рассчитать сумму квадратов отклонений (СКО) между исходными точками и аппроксимированными аналитической функцией.

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

СКО1=21.35

при 0.820

? =0.1185

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

СКО2=0.922

при 0.0755

? =0.028

? = 2.065

Т.к. СКО2 < СКО1, степенная функция лучше аппроксимирует исходные данные. Следовательно, для дальнейших расчетов используем степенную функцию.

Рис. 3. График изменения СКО при различных значениях объема.

Рис. 4. График изменения СКО при различных значениях объема.

3.4 Регрессионный анализ экспоненциальной функции

Формула экспоненциальной функции:

.

Приводим ее к линейному виду y=a+b*x, где вместо x выступает t, логарифмируем:

Пусть , тогда:

Продифференцируем по a и приравняем к 0, т.е.:

Пренебрегаем множителем, раскрываем скобки и выносим а за знак суммирования:

Выведем формулу для a:

(1)

Произведем численный расчет a. Для этого составим таблицу необходимых значений, учитывая, что - значения размеров опухоли, равных Взятое число точек m=14.

Таблица

V (t) 1.351.82.252.74.14.457.37.78.29.111.81515.418.2t6791011121314151718192122-0.020.270.490.671.171.661.721.781.882.142.382.322.412.58-0.121.894.416.711.9914.0421.5824.0826.731.9638.5245.2250.656.76

Из таблицы получаем:

;

;

Подставим полученные значения в формулу (1) и получим:

3.5 Регрессионный анализ степенной функции

Формула степенной функции:

;

Действуя аналогичными методами получаем:

Пусть , тогда:

Продифференцируем по a и по b и приравняем к 0, т.е.:

Пренебрегаем множителями и раскрываем скобки:

Преобразуем эти выражения:

Из полученного выразим a и b:

Произведем численный расчет a и b. Для этого составим таблицу необходимых значений, учитывая, что - значения размеров опухоли, равных Взятое число точек m=14.

Таблица

V (t) 1.351.82.252.74.14.557.37.78.29.111.81515.418.2t67910111213141517181921220.30.5870.810.991.411.4921.9452.0792.1042.22.4682.7082.742.89ln 1.7911.9452.1972.3022.3942.4842.5492.6392.7212.8332.8792.9953.0443.091ln 3.23.7834.8265.2995.7456.176.4976.9647.4038.0258.2888.979.2659.554

Из таблицы получаем:

;

;

;

;

Подставим полученные значения в (2) и (3):

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

4. Расчет запаса жизненных сил. Время жизни после лечения

Геометрический смысл переменой - площадь под выбранной кривой. Следовательно, расчет запаса жизненных сил производится по формуле:

При подстановке значений в нашу программу получаем:

Далее находим эффект вводимой инъекции в зависимости от максимально допустимой дозы (МДП).

, где

Отсюда:

E (D) = D*E=4

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

Учитывая эффект от дозы, можно найти время жизни организма после лечения:

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

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

После подстановки данных получен следующий результат: T (ж) = 44,6.

Рис. 6. Кривые роста опухоли до и после лечения (при многократном введении дозы n=3)

5. Заключение

Проанализировав исходные данные и сделав определенные расчеты времени жизни при однократном введении максимально переносимой дозы и при многократном введении (n=5) дозы D=0.8 МПД с интервалом t=6 суток, я пришел к выводу, что проведение длительного лечения эффективнее, следовательно, продолжительность жизни организма увеличивается в данном случае на 10 суток.

На основании исходных данных и расчетов СКО можно сделать вывод, что развитие опухоли лучше всего аппроксимирует степенная функция.

онкологической заболевание время жизнь

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

1. Бабушкина Н.А. Математическое моделирование как метод изучение функционирования организма. М.: МИРЭА, 1995. - 63с.

. Ван дер Варден Б.Л. Математическая статистика. - М: Эдиториал УРСС, 2006, 435 c.

. Дрейпер Н., Смит Г. Прикладной регрессионный анализ. Множественная регрессия. Изд.2, перераб. и доп. - М.: Диалектика, 2007, 212 с.

. Потемкин В. Введение в MATLAB CHM. М.: Диалог-МИФИ, 2000. - 256 с.

. Пыльнов Ю.В. Регрессионный анализ полиномиальных моделей. - М.: МИРЭА, 1994, 56 с.

. Михальский А.И. Лекционные материалы по курсу КТ в МБС, 2013 год.

Приложение

Текст программы в MATLAB:

= [1.35,1.8,2.25,2.7,4.1,4.55,7.3,7.75,8.2,9.1,11.8,15,15.45,18.2]; % Значения объёма по точкам.= [6,7,9,10,11,12,13,14,15,17,18, 19,21,22]; % Значения времени по точкам.=1.35;=14;=fminsearch (@G1,1, [],Y1,T1,V01) % Экспоненциальная.

Y2=V01. *exp (T1*P (1));= (1/n) * sum ( (Y1 - Y2). ^2);=fminsearch (@G2, [1,1], [],Y1,T1,V01) % Степенная.=V01. * (1+T1/ (P1 (1))). ^ (P1 (2));= (1/n) * sum ( (Y1 - Y3). ^2);(T1,Y1,'ko-',T1,Y2,'k+: ',T1,Y3,'k-. '),grid;('Исходные данные','Экспоненциальная','Степенная');('Время');('Объём');= [];

m2= [];V01=0.0001: 0.01: 3;

[P3,f] =fminsearch (@G1,1, [],Y1,T1,V01);= [m1; V01];= [m2; f];(m1,m2,'k-. ','LineWidth',2),grid;('Мин. Объём');('СКО');= [];= [];V01=0.001: 0.01: 2;

[P4,f2] =fminsearch (@G2, [1,1], [],Y1,T1,V01);= [n1; V01];= [n2; f2];(n1,n2,'k-. ','LineWidth',2),grid;('Мин. Объём');

ylabel ('СКО');

% До лечения.

V01=0.8200;=0.1440;=V01*exp (T1. *a1);=0.0755;2=0.029;=2.065;=0.1; %Промежуток развития=V02*diag (1+ (T1/a2)) ^b; % Создадим массив объёма созданного по формуле с учётом оценки.

Y3=diag (Y3);=Y3';(T1,Y1,'ko',T1,Y2,'k-. ',T1,Y3,'k'); % Сравниваем получившиеся данные, на графике без лечения. В качестве проверки строим экспоненциальную оценку.on('Время');('Объём');('Исходные данные', 'Экспоненциальная функция','Степенная функция');

% Без лечения.m= [0: h: 32]; % Время жизни без лечения.=V02*diag (1+ (T2m/a2)) ^b; % Создадим массив объёма созданного по формуле с учётом оценки.

Y8=diag (Y8);=Y8';=sum (Y8. *h); % Считаем примерный запас жизненных сил, на примере без лечения.

% С лечением.= [0: 1: 5]; % Зависимость эффекта, от максимально переносимой дозы.

MPD= [0: 0.2: 1];(MPD,E,'k-. ');on('Переносимая доза');

ylabel ('Эффект');

%Найдем эффект, при дозе n от переносимой.=0.8;=E/MPD; %Находим тангенс при макс. дозе.=1;=D*MPD*tgEMPD; %Находим эффект при дозе n.=7; %Время начала лечения.=3; %Количество инъекций.=6; %Промежуток между инъекциями.=t1+ (n1-1) *td; %День последней инъекции.

Y41= [];= [];= [];= [];= [];= [];= [];

Tx4= [];=0; % Введём время развития, максимальное значение которого, на выходе из последнего цикла и будет временем жизни с лечением.=t1; % Время 1-ой инъекции.=t1+ (n1-2) *td; % Время 2-ой инъекции.=tj1; % Время последней инъекции.=31; % Время жизни.Tj<Taux1

Y4=V02*diag (1+ (Tj/a2)) ^b;= [Y41; Y4];= [Tx1; Tj];=Tj+h;;=Tj-ED; % Ввод дозы 1.Tj<Taux2=V02*diag (1+ (Tj1/a2)) ^b;= [Y41; Y4];= [Tx1; Tj];=Tj1+h;=Tj+h;;=Tj1-ED; % Ввод дозы 2.Tj<Taux3=V02*diag (1+ (Tj2/a2)) ^b;= [Y41; Y4];= [Tx1; Tj];=Tj2+h;=Tj+h;;=Tj2-ED; % Ввод дозы 3.

Zp=sum (Y41. *h);Zp<zp3; % Рассчитаем график, до летального объёма.

Zp=Zp+Y4*h;=V02*diag (1+ (Tj3/a2)) ^b;= [Y41; Y4];= [Tx1; Tj];=Tj3+h;=Tj+h;;(T2m,Y8,'k-. ',Tx1,Y41,'k');on('Время');('Объём');('Без лечения','С лечением');-файл G1 для расчета экспоненциальной функции:

function out=G1 (P,V1,T1,V01)=P (1);=V01. *exp (T1*a);=V1-Y;

out=sum (out. *out);

M-файл G2 для расчета степенной функции:

function out=G2 (P1,V1,T1,V01)=P1 (1);=P1 (2);=V01. * (1+T1/a). ^b;=V1-Y;=sum (out. *out);

План работы 1. Введение 2. Задание 3. Аппроксимация данных. Построение модели 3.1 Выбор функции 3.2 Аппроксимация в MATLAB 3.3 Выбор оптим

Больше работ по теме: