Интервалы прогнозирования для cfit или sfit объект
ci = predint(fitresult,x)cfit объект fitresult при новых значениях предиктора, заданных вектором x. fitresult должен быть выводом из fit функция, содержащая необходимую информацию для ci. ci является n-by-2 массив, где n = length(x). Левый столбец ci содержит нижнюю границу для каждого коэффициента; правый столбец содержит верхнюю границу.
ci = predint(fitresult,x,level,intopt,simopt)
Границы наблюдения шире, чем функциональные границы, потому что они измеряют неопределенность прогнозирования подогнанной кривой плюс случайное изменение в новом наблюдении. Несимметричные границы для отдельных элементов x; одновременные границы для всех элементов x.
Вычислите и запишите интервалы наблюдения и функционального прогнозирования для соответствия шумным данным.
Создание шумных данных с экспоненциальным трендом.
x = (0:0.2:5)'; y = 2*exp(-0.2*x) + 0.5*randn(size(x));
Поместите кривую в данные с использованием одной экспоненциальной точки.
fitresult = fit(x,y,'exp1');Вычислите 95% интервалы наблюдения и функционального прогнозирования, как одновременные, так и несимметричные. Несимметричные границы для отдельных элементов x; одновременные границы для всех элементов x.
p11 = predint(fitresult,x,0.95,'observation','off'); p12 = predint(fitresult,x,0.95,'observation','on'); p21 = predint(fitresult,x,0.95,'functional','off'); p22 = predint(fitresult,x,0.95,'functional','on');
Постройте график интервалов данных, аппроксимации и прогнозирования. Границы наблюдения шире, чем функциональные границы, потому что они измеряют неопределенность прогнозирования подогнанной кривой плюс случайное изменение в новом наблюдении.
subplot(2,2,1) plot(fitresult,x,y), hold on, plot(x,p11,'m--'), xlim([0 5]), ylim([-1 5]) title('Nonsimultaneous Observation Bounds','FontSize',9) legend off subplot(2,2,2) plot(fitresult,x,y), hold on, plot(x,p12,'m--'), xlim([0 5]), ylim([-1 5]) title('Simultaneous Observation Bounds','FontSize',9) legend off subplot(2,2,3) plot(fitresult,x,y), hold on, plot(x,p21,'m--'), xlim([0 5]), ylim([-1 5]) title('Nonsimultaneous Functional Bounds','FontSize',9) legend off subplot(2,2,4) plot(fitresult,x,y), hold on, plot(x,p22,'m--'), xlim([0 5]), ylim([-1 5]) title('Simultaneous Functional Bounds','FontSize',9) legend({'Data','Fitted curve', 'Prediction intervals'},... 'FontSize',8,'Location','northeast')
