rms.m

非常好的数字处理教程

function rms_value = rms(signal)

    sigsum = 0;
	N = length(signal);
    for i=1:N
		sigsum = sigsum + signal(i)^2;
    end

	rms_value = (1/N)*sqrt(sigsum);