diff --git a/applications/lofar2/model/run_pfir_coeff.m b/applications/lofar2/model/run_pfir_coeff.m
index a16ef12190fa45322ed43d7a94c957b1a5ea4f33..da526ae98a499551c2da7d7118a905e98da34fb8 100644
--- a/applications/lofar2/model/run_pfir_coeff.m
+++ b/applications/lofar2/model/run_pfir_coeff.m
@@ -379,6 +379,8 @@ end
 title(['FIR filter coefficients for ', config.design, strNL]);
 xlabel(['Taps 1:', num2str(L)]);
 ylabel(['Channels 1:', num2str(N)]);
+file_name = ['plots/', file_name_prefix, '_coefficients.jpg'];
+print(file_name, '-djpeg')
 
 %% Plot FIR-filter DC response per polyphase
 if coeff_w>0
@@ -402,6 +404,8 @@ if coeff_w>0
     xlabel(['Polyphase 1:', num2str(N)]);
     ylabel(['Sum of taps (with median at zero is ', num2str(dc_polyphases_median), ')']);
     grid on;
+    file_name = ['plots/', file_name_prefix, '_dc_response_per_polyphase.jpg'];
+    print(file_name, '-djpeg')
 end
 
 % Plot DC adjustment for LOFAR subband FIR-filter coefficients
@@ -433,7 +437,8 @@ grid on;
 title(['Full FIR filter transfer function for ', config.design, strNL]);
 xlabel('Frequency [channels]');
 ylabel('Power [dB]');
-file_name = ['plots/', file_name_prefix, '_fir_transfer_function.jpg'];
+file_name = ['plots/', file_name_prefix, '_transfer_function_full.jpg'];
+print(file_name, '-djpeg')
 
 % . zoomed plot in power
 fig=fig+1;
@@ -447,7 +452,7 @@ grid on;
 title(['Zoomed FIR filter transfer function for ', config.design, strNL]);
 xlabel('Frequency [channels]');
 ylabel('Power');
-file_name = ['plots/', file_name_prefix, '_fir_transfer_function.jpg'];
+file_name = ['plots/', file_name_prefix, '_transfer_function_zoom.jpg'];
 print(file_name, '-djpeg')
 
 % Plot FIR-filter amplitude and phase characteristic using freqz()
@@ -491,3 +496,6 @@ title(['FIR filter sum of power responses between channels for ', config.design,
 xlabel('Frequency [channels]');
 ylabel('Power');
 
+file_name = ['plots/', file_name_prefix, '_transfer_function_two_channels.jpg'];
+print(file_name, '-djpeg')
+