samplingRate = 48000;
width = 550;
cutoff = 1500;
ripple = -80;
A = -ripple;

transitionWidthNorm = 2 * pi * width / samplingRate;
cutoffNorm = cutoff / samplingRate;



nTemp = (A - 7.95) / (2.285 * transitionWidthNorm);
N = ceil(nTemp);

if(mod(N, 2) == 0)
    N = N + 1;
end


beta = 0.1102 * (A - 8.7);



%Working out window starts here.
hN = (N - 1) / 2;
n= -hN : hN;


wd = kaiser(N, beta)';
hd = 2 * cutoffNorm * sin(2 * pi * cutoffNorm * n)./(n * 2 * pi * cutoffNorm);
hd(hN+1) = 2 * cutoffNorm;


hdw = hd.*wd;

%Working out window ends here.

lowpassFile = fopen('LowPassCoef.asm', 'w');
for c = 1 : N
    fprintf(lowpassFile, 'lpcoef_%d\tdc\t%1.8f\n', c, hdw(c));
end
fclose(lowpassFile);

%Start working out high pass filter here.

transitionWidthNorm = 2 * pi * width / samplingRate;
cutoffNorm = cutoff / samplingRate;



nTemp = (A - 7.95) / (2.285 * transitionWidthNorm);
N = ceil(nTemp);

if(mod(N, 2) == 0)
    N = N + 1;
end


beta = 0.1102 * (A - 8.7);



%Working out window starts here.
hN = (N - 1) / 2;
n= -hN : hN;


wd = kaiser(N, beta)';
hd = -2 * cutoffNorm * sin(2 * pi * cutoffNorm * n)./(n * 2 * pi * cutoffNorm);
hd(hN+1) = 1 - (2 * cutoffNorm);


hdw = hd.*wd;

%Working out window ends here.

%Plot graphs of high pass filter.
figure(1)
subplot(3, 1, 1), plot(n, hd);
subplot(3, 1, 2), plot(n, wd);
subplot(3, 1, 3), plot(n, hdw);
[h1,w] = freqz(hd,1,4096,samplingRate);
[h2,w] = freqz(hdw,1,4096,samplingRate);
figure(2)
subplot(2,1,1)
plot(w,20*log10(abs(h1)),w,20*log10(abs(h2)))
subplot(2,1,2)
plot(w,unwrap(angle(h1)),w,unwrap(angle(h2)))

highPassFile = fopen('HighPassCoef.asm', 'w');
for c = 1 : N
    fprintf(highPassFile, 'hpcoef_%d\tdc\t%1.8f\n', c, hdw(c));
end
fclose(highPassFile);