hz_average3.m

Accelerating execution speed of fixed-point algorithms in MATLAB and Simulink

function y=hz_average3(u,acc)
% returns the scalar value for 
% avaergae of matrix u
% special case of avergaing 2x2 matrices using righ-shift
N=size(u,1)*size(u,2);
for i=1:N
acc(1) = acc(1) + u(i);
end;
y=bitshift(acc,-2);