blob: dde2b25d60a55bb6bfba7ad7972b076e45c7a7f6 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
function [skew, kurtosis] = pdf_stats(x,p)
% Define some useful bits
dx = x(2) - x(1);
% The expected value (From LibreTexts 4.1)
E = @(x) sum(x .* p) * dx;
% First moment (mean) (From LibreTexts 4.3)
mu = E(x);
% Second moment (variance) (From LibreTexts 4.3)
sigma = sqrt(E((x - mu).^2));
% Third moment (skewness) (From LibreTexts 4.4)
skew = E(((x - mu)/sigma).^3);
% Fourth moment (excess kurtosis) (From LibreTexts 4.4)
kurtosis = E(((x - mu)/sigma).^4) - 3;
end
|
