cal_rms.m

时间域有限差分法（TD-FDM）源代码

% to calculate normalized root mean square value

function rms = cal_rms(a,b,N)

kk = max(size(a));

sum = 0;
for i=1:kk
    if(b(i) ~= 0)
        sum = ((a(i)-b(i))/b(i))^2 + sum;
    end
end

rms = sqrt(sum/N);
end