Save cl.

applications/lofar2/model/pfb_bunton_annotated/Gen_filter12.m

@@ -30,6 +30,8 @@ Norder = Kcoefs - 1;
 W0 = 1 / 15.3;
 %W0 = 1 / 15.19;
 c = fircls1(Norder, W0, .02, .0003);
+dc_gain_c = sum(c)
+% c = c / dc_gain_c;  % force DC gain = 1
 
 figure(1)
 subplot(3, 1, 1)
@@ -79,8 +81,7 @@ Q = 12;
 % Y = interpft(X, N) returns a vector Y of length N obtained by interpolation
 % in the Fourier transform of X.
 cl = interpft(c, length(c) * Q) / Q;
-dc_gain_c = sum(c);
-dc_gain_cl = sum(cl);
+dc_gain_cl = sum(cl)
 Ncoefs = Kcoefs * Q;
 Nfft = Kfft * Q;
 
@@ -98,3 +99,6 @@ xlabel('Frequency (FFT bins)','FontSize',10)
 ylabel('Magnitude','FontSize',10)
 title ('Response of fircls1 filter c and interpolated cl')
 grid
+
+clT = cl';
+save('cl.txt', 'clT', '-ascii', '-double')