Radar Pulse Compression with Optimized Weighting Window for SAR Receivers

This paper proposes a novel design of a software-defined matched filter (MF) for digital receivers of synthetic aperture radar (SAR). The block diagram of the proposed receiver is described in detail. The purpose of this filter is to produce a SAR pulse with higher compression ratio (CR) and lower side lobe level (SLL) than that produced by the conventional MF. The proposed design is based on the idea of time windowing of the SAR pulse to construct the transfer function of the receiver filter. The shape of the proposed time-domain window is optimized to achieve the filter design goals including the minimization of the SLL and the realization of the target value of the CR. The transmitted SAR pulse is, first, subjected to linear frequency modulation and then subjected to the optimized window. The width (time duration) of the proposed window is divided into equal time intervals. The proposed time-domain window is constructed as a sequential continuous piecewise linear segments. The instantaneous value of the time-domain window at the start of each time interval is optimized so as to achieve the optimization goals. The width of the time-domain window is shown to be proportional to the width of the compressed pulse after optimization. The number of the time intervals into which the time duration of the window is divided is shown to have a significant effect on the optimization results. The particle swarm optimization (PSO) technique is then applied to get the window shape that minimizes the SLL for a specific predetermined value of the pulse CR. It is shown that the iterations of the PSO are fastly convergent and that the applied algorithm is computationally efficient. Also, it is shown that the desired value of the pulse CR is achieved with accuracy of 100%. Moreover, the achieved SLLs are about − 65dB\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$65\,\mathrm{dB}$$\end{document}, -90dB\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$- 90\,\mathrm{dB}$$\end{document}, -114dB\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$- 114\,\mathrm{dB}$$\end{document}, and -133dB\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$- 133\,\mathrm{dB}$$\end{document} for pulse CR of 5, 3, 2, and 1.5, respectively. Finally, for practical implementation of the introduced SAR pulse processing technique, the proposed optimized window is placed as a building block in a software-defined receiver of the SAR system.


Introduction
Synthetic aperture radar (SAR) can be considered as the most effective land imaging system for earth remote sensing irrespective of the daytime and environmental conditions [1][2][3][4]. High imaging resolution and high performance of land target detection are essential requirements for efficient SAR systems [5]. Thus, in the design of SAR system transceiver, the main concern is the imaging resolution and the performance of land target detection. In pulsed SAR systems, a short pulse is required for high image resolution whereas a long pulse is required for high detection performance [5]. Consequently, there is always a trade-off between imaging resolution and detection performance. In this context, SAR pulse compression techniques are adapted to mitigate this trade-off problem by gathering some benefits of both high resolution and detection-performance. In pulse compression techniques, the transmitted SAR pulse has longer duration for high detection performance whereas its bandwidth corresponds to a shorter pulse to achieve high imaging resolution; see Fig. 1.
One of the most important issues in the SAR system design is the software-defined SAR transceiver design [6][7][8][9][10]. In a typical SAR receiver that employs frequency chirping for pulse compression, the received SAR echo is processed using a matched filter (MF) to get its time-domain form at the MF output as a pulse that has a main lobe and multiple sidelobes. The matched filter of the SAR receiver maximizes the signal-to-noise ratio at the receiver output leading to increase the receiver sensitivity and, thus, enhancing the radar system performance. Moreover, for a compressed SAR pulse by applying frequency chirping, the MF leads to suppress the sidelobe level of the received SAR echo at the receiver output and, thus, improves the SAR imaging resolution.
Satisfactory design of the frequency chirping of the transmitted pulse results in a received echo pulse with a main lobe of much higher level than that of the side lobes. Also, the width of the main lobe is much narrower than that of the transmitted pulse. The width of the received echo pulse at the MF output can be calculated as the time difference between the first two nulls on the sides of the main lobe of the signal at the MF output. For satisfactory performance of the operation of SAR pulse compression, it is required to get a pulse compression ratio (CR) that is equal to a predetermined (desired) value. The CR can Fig. 1 Scheme of SAR pulse compression for enhancing the imaging resolution be defined as the ratio between the width of the main lobe of the received signal at the MF output and the width of the transmitted pulse before chirping.
Over the last decade, a lot of research works have been concerned with issue of SAR pulse compression. For example, the work of [9] proposes a radar pulse compression scheme through an approach that is related to the design of efficient non-linear frequency modulation (NLFM) waveforms namely, a temporal predistortioning method of LFM signals by some nonlinear frequency laws. This method produces a compressed radar pulse with SLL of about−40dB . In [10], standard tapered windows are modified and applied to side lobe suppression in compressed pulses with LFM chirping to produce a compressed radar pulse with maximum SLL of about −43dB . In [11], side lobe suppression is performed by some NLFM laws achieved by applying of a stationary phase technique (SPT) synthesis algorithm, where a SLL of about. In [12], different types of weighting window are applied for SLL reduction. This work demonstrates some important numerical results for the achieved SLL for the applied weighting windows of different shapes to show that triangular, Gaussian and Hann give the highest performance. In [13], stretch and correlation processing techniques are applied for pulse compression to achieve a highly compressed pulse, however, the achieved maximum SLL is about −23.8dB . In [14], a variety of adaptive processing techniques are applied to range compression in SAR image focusing and evaluated under real and ideal conditions. Some research work present combined techniques for radar pulse compression. For example, in [15], the SAR pulse compression is performed using spectrum modification and window weighting technique. The application of such a combined technique is shown to be capable of reducing the maximum SLL to about −49dB . The work of [16] introduces combination of phase predistortion and spectrum modification for SAR pulse compression to obtain a maximum SLL of about − 62 dB.
In the present work, a novel design of the MF is proposed for the receiver of SAR to produce a SAR pulse with high CR and low SLL. The proposed MF design employs time windowing of the SAR pulse to construct the transfer function of the receiver filter. The time waveform of the proposed compression window is optimized to achieve the filter design goals including minimization of the SLL and realization of the desired CR. The proposed compression window is placed as a building block in a software defined receiver of the SAR system. The proposed pulse compression process requires the transmitted SAR pulse to be subjected to LFM. The time-domain window is constructed as a sequential continuous piecewise linear segments. The particle swarm optimization (PSO) [17] technique is then applied to get the window shape that minimizes the SLL for a specific target value of the CR.
The PSO is an evolutionary multi-objective optimization technique that can arrive at the best shape of the compression window so as to efficiently satisfy the multiple optimization goals including the minimization of the SLL and the achievement of a specific SAR pulse compression ratio. It can operate with small or large number control parameters even when they have inhomogeneous types. Moreover, the PSO is a simple, computationally efficient, and robust iterative technique that can arrive at the design goals even when the initial values of the control parameters are far from the optimum values. Owing to its fast convergence, the PSO algorithm can arrive at the design goals in a small number of iterations. Among the other evolutionary optimization technique, the PSO can be considered as the most effective and efficient technique for optimizing the shape of the compression window.
The remaining of this paper is organized as follows. The next section gives a mathematical formulation of the frequency chirping of SAR pulse. Section 3 provides a description of the conventional MF receiver used in the SAR system receiver for pulse compression. Section 4 gives a detailed description of the novel SAR system receiver proposed in the present work with the optimized compression window. Section 5 explains how the PSO algorithm is applied to construct the time-domain compression window. Section 6 presents the numerical results with interesting discussions and conclusions. Finally, Sect. 7 summarizes the most important conclusions of the present work.

Frequency Chirping of the SAR Pulse
A conventional method for radar pulse compression is by frequency chirping using LFM. In LFM [5,6], the radar pulse is constructed as a sinusoidal signal whose amplitude is constant over the pulse duration and zero otherwise. The frequency of the sinusoidal signal is f b at the start of radar pulse and increases linearly with the time until it reaches f e at the end time of the pulse. If the pulse duration is T then the slope of increase of the instantaneous frequency is f e − f b ∕T.

Time Waveform of the LFM Chirped Pulse
In LFM, the transmitted signal is a chirped pulse that can be expressed as follows.
where (t) is the instantaneous value of the angle.
The instantaneous frequency, f i (t) can be obtained by differentiating (t) with respect to the time To obtain LFM chirping the instantaneous frequency, f i (t) , should take the following form where t b is the start time of the pulse and t e = T + t b is the end time of the pulse; f b and f e are the start and stop frequencies.
Thus, the angle (t) can be expressed as (t) = ∫ f i (t)dt ; this gives, It should be noticed that constant of integration in (4) is set to get the phase angle of the sinusoidal signal equal to zero at the start of the pulse.

Time Discretization for Simulation
For simulation of the SAR pulse transmission, reception, and processing, the time should be discretized so as to apply fast Fourier transform (FFT) and inverse fast Fourier transform (IFFT) operations. Let N s be the number of time samples of the transmitted SAR pulse, s(t) . Let the time interval between the successive samples be Δt . Thus, the start time of the nth period is t n = (n − 1)Δt + t b and t b = n b − 1 Δt , where n b is the number of the time sample at which the transmitted pulse starts. Thus,. The total time for simulation is T T = LΔt where L is the total number of sampling periods over which the simulation is performed. The center frequency, f c = f b + f e ∕2 , is the operating frequency of the SAR. The bandwidth of operation is B = f e − f b , then the start and stop frequencies f b , f e can be, respectively, expressed as follows.
Thus, the nth frequency components of the transmitted pulse can be expressed as, For accurate simulation, the sampling frequency, f s , should be much greater than the center frequency f c , i.e. f s ≫ f c . The minimum sampling frequency f smin = 2f e .

Conventional Design of the MF Receiver for SAR Pulse Compression
When the transmitted pulse travels from the SAR antenna to the target and is reflected back to the antenna it gets affected by noise. The main role of the MF [5,6] is to efficiently retrieve the known transmitted signal from the received noisy signal, so MF-based receiver, the system is designed to search for the signal which has the similar characteristics as the transmitted pulse, in the received signal. The MF is an optimal filter which maximizes the signal to noise ratio (SNR) at the receiver output.

Conventional Design of MF Receiver for SAR Systems Employing FM Chirping
The transfer function of the MF can be described by the frequency response H(f ) . A software-defined SAR receiver has its all functionalities applicable as software modules except for the low-noise amplifier (LNA). Such a receiver is shown in Fig. 2; the received echo is fed into a MF whose transfer function is the conjugate of the transmitted signal, H(f ) = S * (f ) . The signal, P(f ) , at the output of the MF is the compressed pulse which is the IFFT of the product of the received signal spectrum, R(f ) , and the transfer function of the MF, H(f ) . Thus, the signal spectrum at the MF output can be expressed as, Under the assumption that the time-domain waveform of received pulse r(t) is identical to that of the transmitted signal,s(t) , one can put R(f ) = Ae j S(f ) in (6), where A and are the magnitude and phase of the coefficient of backscattering on the radar target. Thus, the frequency-domain expression at the matched filter output can be expressed as, From (7), it is shown that the bandwidth of the processed echo pulse, p(t) , at the MF output is the same as that of the transmitte chirped pulse, s(t).

Conventional Design of MF Receiver Employing Windowing Technique for SAR Pulse Compression
The MF applies the process of FFT to the received SAR signal as already shown in Fig. 2. The application of this process is performed under the assumption that the received signal is periodic, which is not the general case. When the received signal is non-periodic, then the FFT results in a spread of the received signal spectrum (in the frequency domain) that makes the frequency content of the received signal difficult to identify. To deal with this problem, the proposed receiver design employs a time-domain window function that results in a periodic signal when multiplied by the received signal. As a consequence of employing this window, the resulting sidelobes of the signal at the receiver output are suppressed. Matched filtering can be done in the frequency domain, by combining the MF operation with the windowing so that required computation is reduced. Figure 3 shows the block diagram for frequency domain windowing for LFM side lobe suppression. In this diagram it is shown how the window is applied to the filter frequency response before the true filtering step.

Proposed SAR System Receiver with Optimized Compression Window: Novel Design of the MF Receiver
A time-domain window of a standard shape such as Hamming or Hanning or any other standard pulse can be employed in the proposed receiver design. However, it is proposed to construct an optimized window shape by application of the PSO for the purpose of minimizing the SLL of the signal at the receiver output and achieving a specific pulse CR. The proposed windowing method for SAR pulse compression can be explained in view of the block diagram shown in Fig. 4.
In the proposed software-defined transceiver of the SAR whose block diagram is shown in Fig. 3, the transmitted pulse s(t) is obtained by using phase-locked-loop (PLL) that employs a voltage-controlled oscillator (VCO) with saw tooth control input to get the desired frequency ramp form in the time domain with feedback for frequency controlled. The frequency-domain form of the transmitted signal s(t) is obtained by The transfer function (f ) is constructed by multiplying the FFT of the reference (LFM) pulse, s(t) , which is the pulse to be transmitted, by t W(f ) as follows.
The function R(f ) is obtained as the FFT of the received echo, r(t) and is, then, multiplied by the FFT, W(f ) , of the window function to get the frequency domain function Y(f ).
Finally, the function Y(f ) is multiplied by the conjugate of the transfer function (f ) to get the frequency-domain form the pulse at the output of the SAR receiver, P(f ).
The function P(f ) is then subjected to IFFT to get the received pulse, p(t) , at the output of the MF. According to (9) through (11), the spectrum of receiver output can be expressed as follows.
where H(f ) is the overall transfer function of the SAR receiver.
Making use of (10) the expression (11) can be reformulated as follows.
Making use of (9), the conjugate of (f ) can be substituted into the expression (13) to get.
Finally, P(f ) can be expressed as follows.
For simulation, it is assumed that the received signal r(t) is identical to the transmit signal, s(t) , and hence one can consider R(f ) = Ae j S(f ) as mentioned above. Under this assumption, the substitution for R(f ) into (15) gives following expression for the signal at the MF output.
As the window w(t) is a narrow pulse in time domain, it has a wide spectrum W(f ) in the frequency domain. Thus, according to (16), the pulse p(t) = IFFT{P(f )} is much narrower than s(t) in the time domain.

Optimized Time-Domain Window for SAR Pulse Compression
In this section, the method of compression window optimization to produce minimum SLL for a specific desired value of the pulse CR is described in detail. First, the time-domain window is constructed as piecewise linear curve and then the PSO algorithm is applied after the formulation of the cost function. The following subsections are dedicated for this purpose.

Piecewise-Linear Segmentation for Construction of the Time Waveform of the Compression Window
To get the optimum shape of the windowing function that minimizes the SLL and achieves a desired value of the pulse CR, this window constructed as successive time samples of arbitrary values as shown in Fig. 5, the compression window can be described as a vector w, which gives the values of the window at the corresponding times samples given by the time vector .
where, K, is the number of time samples. Also, it is assumed that the interconnections between the sample points of the compression window are straight-line segments. In other words, each pair of successive points t k−1 , w k−1 and t k , w k are connected by a linear segment. It should be noted that for simulation, the application of the time window starts with the start of the chirped pulse. This means that the duration of the application of the compression window w(t) is the same as the duration of the SAR pulse. However, the window width can be different from the SAR pulse width. This means that t K = T + t 0 and K = N s − N 0 . In general, the compression window can take arbitrary shape. However, the width of w(t) can be given as the 3 dB power width which is the time interval over which the window samples have one-half the total power of the compression window. The width,T w , of the compression window is defined as the central time duration over which the area under the squared curve is one-half the total area under the squared curve. This time interval is centered at t c = N 2 Δt + t 0 .

Application of the PSO to Minimize the SLL of the Received Echo Pulse
The objective of the optimization problem is to arrive at the optimal shape of time-domain compression window in its discrete form under the assumption that the discrete points of the window are interconnected using linear segments. The final goal is to minimize the SLL of the SAR pulse at the output of the MF of the SAR receiver and to realize a predetermined value of the CR. The swarm is composed of N S particles; each particle can be considered as a vector has N components (time samples) as given by (20).

Formulation of the Optimization Problem
Let the achieved time duration between the first two nulls on the sides of the main lobe of the signal r(t) at the output of the MF be T P and let the desired value (optimization goal) of this duration to be T D . Thus, the target value of CR D is T∕T D . The CR of the radar pulse is defined as the ratio between the chirped pulse width, T , and the time duration, T p , between the first nulls on the sides of the main lobe of the compressed pulse at the receiver output. Thus, the achieved CR can be defined as CR = CR P = T∕T P . As the objective of the PSO algorithm is to minimize the SLL and to get the compressed pulse width, T P , as close as possible to the desired width, T D , the cost function can be formulated as follows.
where max ( ) , is the maximum SLL achieved by a window shape that is given by the vector , T P ( ) is the achieved width, T D is the desired pulse width, F S and F P are weighting factors to determine the contributions of the maximum SLL and the deviation from the desired window width to the cost function C( ). According to (19), the PSO has two objectives; minimization of the maximum SLL and achieving a predetermined value of the pulse CR. It may be logic to intuit that relaxing the condition to achieve the desired pulse CR enables the PSO algorithm to achieve a lower value of the maximum SLL whereas aiming at a higher CR will cause the achievable SLL to increase.
The optimization problem can be formulated as follows. It is required to get the samples that achieve the following The constraint (20-b) means the window samples should be positive while the PSO iterations are running.

Implementation of the PSO
The optimization space (swarm of particles and their positions) contains N P particles. The nth particle has its position within the swarm in the th iteration given by the vector ( ) n . The implementation of the PSO algorithm can be divided into four stages: (i) initialization of positions, n (0) and velocities, n (0) , for the particles of the swarm. The initial values of the particles w n (0) where ∀n ∈ 1, 2, .., N s can be set by assigning the sample values of W n,k (0) to take the shape of Hamming window with the width, W H . The samples of where k is the time sample index, n is the particle number, and r n,k is a random ∈ [0, 1]. The subscripts n, k of the r n,k mean that this random number is generated for each pair of the indices(n, k) . (ii) calculation of the local best positions, n , for the particles, (iii) Calculation of the global best position, g, (iv) calculation of the particles velocities, n (1) , and positions, n (1) , for the next iteration. The superscript (H) refers to the samples of the Hanning window that is suggested as the initial values of the window samples.
The following equations are used to implement an iterative PSO algorithm: where is the iteration number (time index), u is the inertia weight parameter,c 1 , c 2 are acceleration factors and r 1 , r 2 are random numbers between 0 and 1. For each particle, the initial position, n (0) , is determined by assigning random values (varying around the unity) to the amplitudes of the excitation voltages of the array elements. The initial local best position of each particle, n (0) , is assigned the same initial value of the particle position, n (0) . The initial value of the velocity of each particle is set to zero. The local best position, ( ) n , for each particle is the position of this particle that results in the minimum value of the cost function over the progressive iterations during the run of the PSO algorithm. The global best position, ( ) , for the particles in the swarm is the position among the local best positions that results in the absolute minimum value of the cost function over the successive iterations during the run of the PSO algorithm. In each iteration, the velocity, n , and position, n ,for each particle in the swarm are updated, as given by (22-a) and (22-b).

Results and Discussions
This section is concerned with the presentation and discussion of the numerical results of the achieved SLL and pulse CR using the conventional LFM chirping and the conventional weighting window methods. Also various the achievements of SAR pulse compression techniques proposed in some recent publications. Finally the SLL and pulse CR achieved by the optimized window technique introduced in the present work are presented, discussed, and compared to those achievement available in the most recently published work.

Chirped SAR Pulse using Conventional LFM
The LFM is a conventional method for improving the resolution of the radar, the unchirped radar pulse, shown in Fig. 6a is subjected to chirping using LFM. Thus, the square pulse becomes the envelope of a sinusoidal wave whose instantaneous frequency increases linearly with time as shown in Fig. 6b. It is shown that the phase of the sinusoidal signal is zero at the start time of the chirped pulse. A numerical example can be used for quantitative explanation; a square pulse of duration T = 50ns , Fig. 6a is subjected to chirping using LFM as described in Sect. 2 to get the chirped pulse whose time waveform is shown in Fig. 6b. For discretization, the sampling frequency is set as f s = 15f c , where f c = 1.27GHz is the central frequency, N s = 15 , and L = 12700 . The application of the FFT results in the frequency spectrum of the chirped pulse shown in Fig. 6c. The process of the conventional MF described in Sect. 3 is applied to receive and process the echo pulse, r(t) , reflected from the radar target due to the transmitted chirped pulse, s(t) , described above. The time waveform of the conventionally LFM-chirped pulse at the receiver MF output is shown in Fig. 7. It is shown that the processed echo pulse, p(t) , at the MF output has a main lobe and many side lobes of very low level relative to that of the main lobe. The main lobe is compressed, i.e., it has very short time duration relative to the duration T of the transmitted radar pulse. However, many side lobes exist, the highest of them is the second lobe. It shown that the maximum SLL is about − 13.1 dB . The width of the pulse p(t) can be calculated as the time difference between the first two nulls on the sides of the main lobe. It is shown that the compressed pulse width is 0.39ns . The pulse compression can be obtained by dividing the width of the pulse before compression, Fig. 6a by the width of the main lobe of p(t) ; this gives a CR of 128 . Thus, the frequency chirping using LFM results in very high CR but such SLL may not be adequate for accurate SAR imaging.

Compression Window Technique to Improve the Performance of LFM-Chirping
As mentioned in Sect. 4, the application of time or frequency windowing to LFM-chirped pulse leads to improve the performance of the pulse compression regarding the CR and the maximum SLL. One of the well-known and most effective time-domain windows is the Hanning window. The application of Hanning window of width 1ns , shown in Fig. 8a, to the LFM-chirped pulse, shown in Fig. 6, results in the time waveform shown in Fig. 8b for the received signal at the output of the MF. The width of the main lobe is 1.2ns , which corresponds to CR of about 42, whereas the maximum SLL is about − 22.6dB.
For investigating the impact of the compression window width on the effectiveness of the proposed pulse compression process, the following example is demonstrated. In this example, the application of Hanning window of width 6ns , whose time waveform is presented in Fig. 9a, results in the received signal at the output of the MF whose time waveform is shown in Fig. 9b. It is shown that the width of the main lobe is 6.6ns , which corresponds to CR of about 7.5 , whereas the maximum SLL is about −28.2dB.
From the previous two examples, it is clear that to obtain lower SLL of the received pulse at the MF output, the width of the Hanning window should be increased which, in turn, results in lower CR of the SAR pulse. Thus, a compromise should be done to choose the window width that is more appropriate for a specific application.

The Optimized Window for SAR Pulse Compression
In this subsection, six numerical examples are presented to show the performance of SAR pulse compression using the optimized window technique proposed in the present work. The PSO algorithm runs to construct the optimum window with the optimization goal to achieve a predetermined value of the CR and to minimize the SLL. Let the initial window shape for the PSO algorithm be that of a Hanning window of width 1ns as that shown in Fig. 10a. In this case, the PSO algorithm takes about 300 iterations to arrive at the optimized window shape shown in Fig. 10b whose width is still 1ns like that of the initial window. The resulting time waveform of the compressed pulse at the MF output is shown  Fig. 10c where the main lobe width is about 2.5ns . The resulting CR is 20 , whereas the maximum SLL is about − 40.5dB . Thus, when compared with the case of applying the Hanning window of width 1ns as previously discussed; see Fig. 8, the application of an optimized window of the same width results in reducing the SLL from − 22.6dB to −40.5dB whereas the CR is decreased from 42 to 20 . As the SLL is more critical for the radar operation than the CR, the application of the optimized window can be considered much superior to the application of the Hanning window.
Other cases are for compressed pulse using the optimized window present in Figs. 11, 12, 13 and 14. In these examples, it is shown that the lower the predetermined (desired) pulse CR the lower the achieved SLL.
The last example of the application of PSO on Hanning window is shown in Fig. 15 with width 6ns . The optimized window shape is shown in Fig. 15a. The output of the MF is  Fig. 15b where the main lobe width is about 25ns . The resulting CR is 2 , whereas the maximum SLL is − 115 dB . Thus, when compared with the case of applying the Hanning window of width 6ns , the application of the optimized window of the same width results in reducing the SLL from − 28.2dB to − 114dB whereas the CR is decreased from 7.4 to 2.
For comparative assessment of the effectiveness of the proposed method, the well known Hanning window, as one of the conventional time-domain shapes of the radar pulse compression windows, is applied to compress the SAR pulae presented in Fig. 6. so as to compare its performance regarding the resulting pulse CR and maximum SLL with the performance obatined by the application of the optimizaed window proposed in the present work. A Matlab® is written specially for the purpose of comparizon. It should be noted that some results that have been obtained by applying the Hanning window are presented in Figs. 8 and 9. Table 1 shows a list of the achieved SLL and the corresponding CR using the Hanning window compared to those achieved by the optimized window technique proposed in the present work. It is shown that, for the same value of the pulse CR, the proposed window achieves much lower value of the SLL than that achieved by the Hanning window. It should be noted that, a narrower compression window gives higher pulse CR. For example, for a pulse CR of 20 , Hanning window of width 1ns achieves a SLL of − 23.5dB , whereas the proposed window (of the same width) achievs a SLL of − 42dB . For relatively low pulse CR of 5, Hanning window of width 5.2ns achieves a SLL of −36dB , whereas an optimized window of the same width achievs a SLL of −65dB . Thus, with increasing the compression window width the maximum SLL is improved whereas pulse CR is decreased. However, the rate of reduction (improvement) of the SLL with increasing the width of the compresion window using the proposed optimized window is much better than that obtained by the application of Hanning window provided that both windows have the same width.

Relations Among the Achieved CR, SLL and Width of the Compression Window
Two objectives are to be achieved Predetermined value of the CR and low SLL of the received SAR pulse at the MF output. Figure 16 presents the dependencies of the max SLL of the received pulse at the MF output and the corresponding CR on the width of  the compression window that has been optimized using the PSO. It is shown that to obtain lower SLL of the received pulse at the MF output the width of compression window should be increased which results in lower CR of the SAR pulse. Thus, the objectives of high CR and low SLL seem to be contradicting with each other. A compromise should be carried out to select the optimum width of the compression window.

Computational Efficiency and Rate of Convergence of the PSO Algorithm
The cost function C( ) given by (19) decays with progressive iteration of the PSO algorithm as shown in Fig. 17 for three different cases. As shown in this figure, the PSO consumes from 200 -400 iterations to reach the steady state at which the best achievable performance is obtained. The stair case in the curve is attributed to the application of the constraint given by (20-b). In each iteration, if a time sample of the compression window is negative it is overridden to be zero. Of course, the application of this constraint causes the decay curve to be locally horizontal which results in the staircase appearance of the curve and also, slows down the decay of the cost function with the PSO iterations. The fast convergence of the PSO algorithm and the low value of the cost function at the steady state of the PSO algorithm reflects the computational efficiency of the proposed method including the suitable formation of the optimization problem together with the cost function and the robustness of the developed algorithm for the application of the PSO.

Comparative Performance Summary
Comparisons among the achieved SLL using different pulse compression techniques proposed in some recent publications are demonstrated in Table 2. It is shown that the optimized compression window method proposed in the present work gives much lower SLL than those obtained in the other published work. However, the achieved SLL in the present work is variable and depends on the desired value of the pulse CR as listed in Table 2. It should be noted that the methods listed in Table 2 are briefly mentioned in the "Introduction" section of the present paper. Also, it is noticed that the best results considering the suppression of the sidelobes result in SLLs of those obtained using the methods that employ combined techniques like those presented in [15] and [16] to obtain SLLs of about −49dB and −62dB , respectively. Recall that the method presented in [15], applies the combined spectrum modification and window weighting technique for SAR pulse compression, whereas the work of [16] introduces the combined technique of phase predistortion and spectrum modification for SAR pulse compression. However, the maximum SLLs achieved in the present work depend on the pulse CR. To achieve better SLL, the pulse CR should be reduced. As shown in Table 1, a maximum SLL of −65dB is achieved for pulse CR of 5, which is better than those obtained in [15] and [16].

Conclusion
A novel design of a software-defined filter for digital receivers of SAR systems has been proposed. The proposed filter produces a SAR pulse with higher CR and lower SLL than that produced by the conventional MF. The novel filter design is based on the time windowing of the SAR pulse, where the shape of the proposed time-domain window is optimized to minimize the SLL and to realize the desired pulse CR. The transmitted SAR pulse is, first, subjected to LFM and then subjected to the optimized window. The proposed window is constructed as a sequential continuous piecewise linear segments. The PSO algorithm is then applied to get the window shape that minimizes the SLL for a specific predetermined value of the pulse CR. The iterative PSO is shown to be fastly convergent and computationally efficient. Also, it is shown that the desired value of the pulse CR is exactly achieved. On the other hand, the achieved SLLs are about − 65 dB , − 90 dB , − 114 dB , and − 133 dB for pulse CRs of 5, 3, 2, and 1.5, respectively. By comparsison, the method proposed in the present work has been proved to be superior to those conventional windowing techniques used for radar pulse compression. For example, a pulse CR of 20 , Hanning window of width 1ns achieves a SLL of − 23.5dB , whereas the proposed window achievs a SLL of −42dB . For relatively low pulse CR of 5, Hanning window of width 5.2ns achieves a SLL of −36dB , whereas an optimized window of the same width achievs a SLL of −65dB . Als, when compared with the numerical results presented in some recently published work, the SLL achieved in the present work for specific pulse CR is shown to be better than those than those obtained by the other techniques. As the proposed technique has been applied in the present paper for only LFM chirped pulse, it is suggested, as a future extension, to apply the optimized window technique, for compression of a SAR pulse that has been chirped using NLFM. This is expected to give better performance of the SAR pulse compression than that obtained using LFM chirping.
Funding Open access funding provided by The Science, Technology & Innovation Funding Authority (STDF) in cooperation with The Egyptian Knowledge Bank (EKB). Funding not applicable to this article (no funds for this research).
Data Availability Data sharing not applicable to this article as no datasets were generated or analyzed during the current study. Availability of data and material not applicable as there is no data sets used.
Code Availability Code availability not applicable, no codes developed in this research.

Declarations
Conflict of interest Authors declare that there are no Conflicts of interest for this research.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.