Порядок выполнения лабораторной работы. · Вызвать MATLAB

Функцией дожития (надёжности) называется вероятность того, что изучаемый процесс продлится не меньше заданного времени t

,

Т обозначает случайную величину, имеющую смысл длительности процесса (длительности безотказной работы прибора, длительность болезни и т.д.).

Для модели Гомпертца функция дожития имеет вид

. (2.1)

1.1 Получите у преподавателя значения параметров a и b. Сформируйте вектор значений переменной t в диапазоне от 0 до 30 с шагом 0.1, для которых будут вычисляться значения функция дожития, оператором

t=0:0.1:30;

1.2 Вычислите значения функции дожития для модели Гомпертца по формуле (2.1) с помощью оператора

S=exp(-a/b*(exp(b*t)-1));

1.3 Запишите оператор построения кривой дожития

plot(t,S); xlabel(‘t’); ylabel(‘S(t)’); title(‘Survival by Gompertz using plot’);

1.4 Откройте в новом окне редактора пустой файл, запишите в него оператор вычисления функции дожития для модели Гомпертца по формуле (2.1). Оформите файл как файл-функцию GMP с входными параметрами а, b и t и сохраните его в рабочей директории в файле с именем GMP.m

function out=GMP(a,b,t)

out=exp(-a/b*(exp(b*t)-1));

1.5 В ранее открытом файле Lab2.m запишите операторы открытия нового окна, построения графика функции дожития для модели Гомпертца на интервале [0,30] с помощью оператора fplot и функции GMP

figure;

fplot(@GMP,[0,30]);

xlabel(‘t’); ylabel(‘S(t)’); title(‘Survival by Gompertz using fplot’);

1.6 С помощью команды run Lab2 выполните записанные операторы. Значения параметров, построенные графики внесите в отчёт о выполнении лабораторной работы. Сделайте вывод о близости графиков, построенных с помощью операторов plot и fplot.

Для модели Гомпертца-Майкхема функция дожития имеет вид

. (2.2)

2.1 Получите у преподавателя значения параметров a, b и c. Сформируйте вектор значений переменной t в диапазоне от 0 до 30 с шагом 0.1, для которых будут вычисляться значения функция дожития, оператором

t=0:0.1:30;

2.2 Вычислите значения функции дожития для модели Гомпертца по формуле (2.1) и для модели Гомпертца-Майкхема по формуле (2.2) с помощью операторов

SG=exp(-a/b*(exp(b*t)-1));

SGM=exp(-a/b*(exp(b*t)-1)-c*t);

2.3 Запишите оператор построения кривых дожития на одном графике с выводом подписей к кривым

plot(t,SG,t,SGM);

xlabel(‘t’); ylabel(‘S(t)’); title(‘Survival by Gompertz and Gompertz-Mekham’);

legend(‘Gompertz’,’Gompertz-Mekham’);

2.5 С помощью команды run Lab2 выполните записанные операторы. Значения параметров и построенный график внесите в отчёт о выполнении лабораторной работы. Сделайте вывод о влиянии величины параметра с на крутизну падения кривой дожития.

Задание 3. Построить в одном окне графики функции надёжности для модели Гомпертца, график логарифма интенсивности отказа, график кумулятивного риска