function H=mfrdadd(H1,H2,n)
% Add up two MFD freq response arrays
% Copyright (c) Dingyu Xue, Northeastern University, China
% Last modified 18 May, 2022
[H1,H2]=mfrdmat(H1,H2,n); H=H1+H2;
end