Дифференциация и интеграция подгонки
Этот пример показывает, как найти первые и вторые производные подгонки и интеграл подгонки, в значениях предиктора.
Создайте базовый синусоидальный сигнал:
xdata = (0:.1:2*pi)'; y0 = sin(xdata);
Добавьте шум в сигнал:
noise = 2*y0.*randn(size(y0)); % Response-dependent noise
ydata = y0 + noise;
Соответствуйте шумным данным пользовательской синусоидальной моделью:
f = fittype('a*sin(b*x)'); fit1 = fit(xdata,ydata,f,'StartPoint',[1 1]);
Найдите производные сглаженной функции в предикторах:
[d1,d2] = differentiate(fit1,xdata);
Отобразите на графике данные, подгонку и производные:
subplot(3,1,1) plot(fit1,xdata,ydata) % cfit plot method subplot(3,1,2) plot(xdata,d1,'m') % double plot method grid on legend('1st derivative') subplot(3,1,3) plot(xdata,d2,'c') % double plot method grid on legend('2nd derivative')
Обратите внимание на то, что производные могут также быть вычислены и построены непосредственно с методом графика cfit, можно следующим образом. Метод графика, однако, не возвращает данные по производным.
plot(fit1,xdata,ydata,{'fit','deriv1','deriv2'})
Найдите интеграл подгонки в предикторах:
int = integrate(fit1,xdata,0);
Отобразите на графике данные, подгонку и интеграл:
subplot(2,1,1) plot(fit1,xdata,ydata) % cfit plot method subplot(2,1,2) plot(xdata,int,'m') % double plot method grid on legend('integral')
Обратите внимание на то, что интегралы могут также быть вычислены и построены непосредственно с методом графика cfit, можно следующим образом. Метод графика, однако, не возвращает данные по интегралу.
plot(fit1,xdata,ydata,{'fit','integral'})
