Интервалы предсказания для 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')
