Moving average of semi-deviation
2 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Hello,
I am trying to calculate the moving average of the semi-deviation of 1 column of data. I have been able to find a semi-deviation add-in but cannot find anything about making a moving average of the semi-deviation.
Any help appreciated.
0 commentaires
Réponses (4)
Image Analyst
le 15 Juin 2016
What is the definition of " semi-deviation"?
If you mean "Standard Deviation" then try stdfilt() or movstd().
1 commentaire
Image Analyst
le 16 Juin 2016
If you have a column vector of data and want the sd of data below the mean, just to
theMean = mean(data);
semiDeviation = std(data(data<theMean));
Star Strider
le 16 Juin 2016
I wrote my own function from your link and implemented an unfortunately inefficient filter to get the moving semideviation (using my own semideviation function):
semidev = @(x) sqrt(sum((mean(x(:)) - x(x(:) < mean(x(:)))).^2) / length(x(:)));
intvl = 10; % Index Interval For Moving Average
t = 1:100; % Time Vector
data = randn(1, 100); % Create Data
datav = [data, data(end-intvl+1:end)]; % ‘Pad’ ‘data’ With Repeat End Data
for k1 = 1:length(datav)-intvl;
wndw = k1:k1+intvl; % Subscript ‘Window’
sv(wndw) = semidev(datav(wndw)); % Calculate Semideviation Over ‘Window’
end
sv = sv(1:length(data)); % Trim To Correct Length
figure(1)
plot(t, data, 'bp')
hold on
plot(t, sv, '-r')
hold off
grid
It runs. I’ll let you determine if it produces the correct results.
0 commentaires
N/A
le 17 Juin 2016
2 commentaires
Star Strider
le 17 Juin 2016
If it gives you the result you want, go for it!
I based my anonymous function code on the expression you linked to. It should give you the result you want, since all the x-values meeting the criteria are set to their appropriate values, and those that do not meet the criteria are set to zero by default anyway.
Chris Turnes
le 17 Juin 2016
The movstd call here will include the zeros in the calculation. From the definition, it sounds like you want to ignore those points, rather than treat them as zero. I think you could do this by:
data_ds( data > data_movmean ) = nan;
data_movdsstd = movstd( data_ds, [intvl, 0], 'omitnan' );
This should actually ignore those points rather than just consider them to be zero.
Voir également
Catégories
En savoir plus sur Axis Labels dans Help Center et File Exchange
Produits
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!