Explicado como: 'la varianza de la señal es la media de sus cuadrados menos el cuadrado de su media'
La media de sus cuadrados (promedio de la potencia instantánea de todas las muestras) menos el cuadrado de su media (todas las muestras sumadas y cuadradas).
$$ \ sigma (x) \ equiv E [x ^ 2] - (E [x]) ^ 2 $$
¿No estoy seguro de si mi script es correcto? ¿Y qué es E en las ecuaciones de varianza?
%%
clc
clear
%% Signal
x = [ 0 5 5 0 -5 -5 ];
fs = 6;
%% Power Calculations
% Instantaneous power of each sample:
ins_pwr = x .^ 2;
disp(['Instantaneous Power (w) = ' mat2str(ins_pwr)]);
% Average power of signal:
avg_pwr = sum(ins_pwr)/length(x);
disp(['Average power (w) = ' num2str(avg_pwr)]);
% Variance of the signal:
var_pwr = (sum(ins_pwr)/length(x)) - (sum(x)^2);
disp(['Variance (w) = ' num2str(var_pwr)]);
Salida:
- Potencia instantánea (w) = [0 25 25 0 25 25]
- Potencia media (w) = 16.6667
- Varianza (w) = 16.6667