Вычитание тренда из данных

Введение

Функция MATLAB® detrend вычитает среднее значение или наилучшую эмпирическую кривую (в смысле наименьших квадратов) из ваших данных. Если ваши данные содержат несколько столбцов данных, detrend обработки каждый столбец данных отдельно.

Удаление тренда из данных позволяет вам сфокусировать свой анализ на колебаниях данных о тренде. Линейный тренд обычно указывает на систематическое увеличение или уменьшение в данных. Систематический сдвиг может быть результатом дрейфа датчика, например. Хотя тенденции могут быть значимыми, некоторые виды анализа дают лучшее представление после удаления тенденций.

Целесообразно ли это удалять эффекты тренда в данных, часто зависит от целей вашего анализа.

Удаление линейного тренда из данных

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

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

t = 0:300;
dailyFluct = gallery('normaldata',size(t),2); 
sdata = cumsum(dailyFluct) + 20 + t/100;

Найдите среднее значение данных.

mean(sdata)
ans = 39.4851

Постройте и пометьте данные. Заметьте систематическое увеличение курсов акций, которые отображают данные.

figure
plot(t,sdata);
legend('Original Data','Location','northwest');
xlabel('Time (days)');
ylabel('Stock Price (dollars)');

Примените detrend, который выполняет линейную подгонку к sdata и затем удаляет тренд из него. Вычитание выхода от входа дает к вычисленной линии тренда.

detrend_sdata = detrend(sdata);
trend = sdata - detrend_sdata;

Найдите среднее значение детрендированных данных.

mean(detrend_sdata)
ans = 1.8484e-14

Как ожидалось, данные без тренда имеет среднее значение очень близкое к 0.

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

hold on
plot(t,trend,':r')
plot(t,detrend_sdata,'m')
plot(t,zeros(size(t)),':k')
legend('Original Data','Trend','Detrended Data',...
       'Mean of Detrended Data','Location','northwest')
xlabel('Time (days)'); 
ylabel('Stock Price (dollars)');

Смотрите также

| | |