Fs = 200000000; % Sampling Frequency
Fstop1 = 78500000; % First Stopband Frequency
Fpass1 = 78800000; % First Passband Frequency
Fpass2 = 79200000; % Second Passband Frequency
Fstop2 = 79500000; % Second Stopband Frequency
Dstop1 = 0.001; % First Stopband Attenuation
Dpass = 0.057501127785; % Passband Ripple
Dstop2 = 0.0001; % Second Stopband Attenuation
dens = 20; % Density Factor
% Calculate the order from the parameters using FIRPMORD.
[N, Fo, Ao, W] = firpmord([Fstop1 Fpass1 Fpass2 Fstop2]/(Fs/2), [0 1 ...
0], [Dstop1 Dpass Dstop2]);
% Calculate the coefficients using the FIRPM function.
b = firpm(N, Fo, Ao, W, {dens});
Hd = dfilt.dffir(b);
[b,a]=tf(Hd);
yt75=filter(b,a,out1 + out2*i);