generateHighPassFilter(Fpass, Fstop, Fs, N)

function Hd = generateHighPassFilter(Fpass, Fstop, Fs, N)


Dstop = 0.0001;          % Stopband Attenuation
Dpass = 0.057501127785;  % Passband Ripple
dens  = 20;              % Density Factor

% Calculate the order from the parameters using FIRPMORD.
[N, Fo, Ao, W] = firpmord([Fstop, Fpass]/(Fs/2), [0 1], [Dstop, Dpass]);

% Calculate the coefficients using the FIRPM function.
b  = firpm(N, Fo, Ao, W, {dens});

% Calculate the coefficients using the FIRPM function.
Hd = b;