/ / Article
Special Issue

## Iterative Methods for Nonlinear Equations or Systems and Their Applications

View this Special Issue

Research Article | Open Access

Volume 2013 |Article ID 949084 | https://doi.org/10.1155/2013/949084

Bin Wang, Jinkuan Wang, Xin Song, Fengming Xin, "Transmitted Waveform Design Based on Iterative Methods", Journal of Applied Mathematics, vol. 2013, Article ID 949084, 10 pages, 2013. https://doi.org/10.1155/2013/949084

# Transmitted Waveform Design Based on Iterative Methods

Revised09 Jan 2013
Accepted24 Jan 2013
Published11 Mar 2013

#### Abstract

In intelligent radar, it is an important problem for the transmitted waveform to adapt to the environment in which radar works. In this paper, we propose mutual information model of adaptive waveform design, which can convert the problem of adaptive waveform design into the problem of optimization. We consider two situations of no clutter and clutter and use Newton method and interior point method to solve the optimization problem. Then we can draw the design criterion for the transmitted waveform in cognitive radar and get a greater mutual information from the simulation results. Finally, the whole paper is summarized.

#### 1. Introduction

In this paper, we propose mutual information model of adaptive waveform design, which can convert the problem of adaptive waveform design into the problem of optimization. We consider two situations of no clutter and clutter and use Newton method and interior point method to solve the optimization problem. Then we can draw the design criterion for the transmitted waveform in cognitive radar from the simulation results.

#### 2. Mutual Information Model of Adaptive Waveform Design

The basic parts of a radar system are illustrated in the diagram of Figure 1. The equipment is divided into several subsystems, corresponding to the usual design specialties within the radar engineering field. A radar system has a transmitter that emits radio waves called radar signals in predetermined directions. When they come into contact with an object the signals are usually reflected and/or scattered in many directions. The radar signals that are reflected back towards the transmitter are the desirable ones which make radar work.

Different from traditional radar, cognitive radar is constructed using intelligent signal processing, information feedback loop, and soft information processing. In cognitive radar, the radar continuously learns about the environment through experience gained from interactions of the receiver with the environment, the transmitter adjusts its illumination of the environment in an intelligent manner, and the whole radar system constitutes a closed-loop dynamic system. Figure 2 is the block diagram of cognitive radar .

From the block diagram of cognitive radar, we can see that constructing a waveform library is very important in cognitive radar. Through sensing the environment, cognitive radar transmits waveform suited to the working conditions. The radar returns, and environment factors can help to reconstruct the waveform library. Then the radar can select different waveforms to transmit. It forms a feedback loop, and the cycle goes on and on. Figure 3 is block diagram of waveform library.

Waveform library can store many kinds of waveforms. The design of radar waveforms has been a topic of considerable research interest for several decades. In traditional radar systems, the radar transmits single waveform. The radar is difficult to adapt to different environments. Modern radar is required to transmit different waveforms according to different environments. So a more flexible design framework is required, which should be able to synthesize waveforms that provide a smooth trade-off between competing design criteria.

Cognitive radar is the next generation radar system. Figure 4 is a basic signal-processing cycle in cognitive radar. Cognitive radar has the capability to observe and learn from the environment. It operates in closed loop, and the transmitted waveform will be adaptive. In order to achieve objectives more efficiently, the waveforms should be adapted in response to prior measurements.

Generally speaking, the waveform design is different as a result of different tasks of radar. For detection task, the optimal radar waveform should be able to put as much transmitted energy as possible into the largest mode of the target to maximize the output signal-to-noise ratio (SNR). For estimation task, the optimal radar waveform should allocate the energy between the received signal and the target signature. For other tasks, other performance is required.

Cognitive radar is an intelligent system. In different radar environments, it can synthesize different waveforms. One possible scheme is to make a trade-off among different performances. Efficient algorithms are required in the construction of cognitive radar systems. Such algorithms should provide a flexible framework that can synthesize waveforms that provide different trade-offs between a variety of performance objectives which themselves may also be adapted to the perceived nature of the environment.

We consider two situations of no clutter and clutter and set up their mutual information model, respectively.

Figure 5 is the signal model of a target in which there is no clutter.

We want to find the mutual information , that is, the mutual information between the random target impulse response and the received radar waveform. Those functions that maximize also maximize . So we maximize firstly. Assume that target is Rayleigh type and noise is Gaussian type. They are statistically independent.

Assume that represents frequency domain sampling point, is a frequency point. Let correspond to the component of with frequency components in , let correspond to the component of with frequency components in , and let correspond to the component of with frequency components in . So the overall mutual information is

Assume that the frequency interval is sufficiently small, so for , , , . is the bandwidth.

Next, we define mutual information. In probability theory and information theory, the mutual information of two random variables is expressed as the dependence of them. Mutual information can be defined from mathematics as where is joint probability distribution function, and and are marginal probability distribution function of and , respectively. Intuitively, mutual information contains the total information of and . Assume that represents the marginal entropy of , represents the conditional entropy of given . So the mutual information can also be expressed as Since , are statistically independent, so the variance of is We will now solve and , respectively. So the mutual information is given by

Referring again to the signals , , , and with frequency components confined to the interval , we have from the sampling theory that each of the signals can be represented by a sequence of samples taken at a uniform sampling rate of . Since we assume that the spectra , , and are smooth and have a constant value (at least approximately) for all , the samples of the Gaussian process sampled at a uniform rate are statistically independent.

The samples are independent, identically distributed random variables with zero mean and variance ; we note that the total energy in is

Over the time interval , this energy is evenly spread among statistically independent samples. Hence, the variance of each sample, , is

Similarly, the noise process has total energy on the interval given by

This energy is evenly distributed among the statistically independent, zero-mean samples of . Hence, the variance of each sample is

Substituting (8) and (10) into (6), we have that for each sample of and corresponding sample of , the mutual information between and is

Now these are statistically independent sample values for both and in the observation interval . Thus,

The overall mutual information is

Following we will consider the situation that there is clutter. Figure 6 is the signal model of a target ensemble in ground clutter. Assume that target is Rayleigh type, noise is Gaussian type, and clutter is Rayleigh type. They are statistically independent.

We want to find the mutual information , that is, the mutual information between the random target impulse response and the received radar waveform. Those functions that maximize also maximize . So we maximize firstly. Assume that target is Rayleigh type, noise is Gaussian type, and clutter is Rayleigh type.

The probability density distribution function of Rayleigh distribution is where is clutter amplitude, and is standard deviation of clutter. Curve of probability distribution of Rayleigh clutter is in Figure 7. Figure 8 is Rayleigh distribution clutter.

Assume that represents frequency domain sampling point, is a frequency point. Let correspond to the component of with frequency components in , let correspond to the component of with frequency components in , and let correspond to the component of with frequency components in . So the overall mutual information is

Assume that the frequency interval is sufficiently small, so for , , , , and . is the bandwidth.

Since , , and are statistically independent, so the variance of is

We will now solve and , respectively, as follows:

So the mutual information is given by

The total energy in is

Over the time interval , this energy is evenly spread among statistically independent samples. Hence, the variance of each sample, , is

Similarly, the noise process has total energy on the interval given by

This energy is evenly distributed among the statistically independent, zero-mean samples of . Hence, the variance of each sample is

Similarly,

Substituting (20), (22), and (24) into (18), we have that for each sample of and corresponding sample of , the mutual information between and is

Now these are statistically independent sample values for both and in the observation interval . Thus,

The overall mutual information is

#### 3. Optimal Waveform Design Using Newton Method and Interior Point Method

Considering the situation that there is no clutter, we can get

We can write in another form where is the length of discrete transmitted signal, and is the autocorrelation function of the transmitted signal as follows: Following we will consider some constraints.

First, in order to guarantee that radar can detect target, SNR needs to be greater than a certain threshold; that is,

Second, the energy of transmitted signal should be a fixed value; that is,

Third, the power of transmitted signal should be greater than a certain threshold; that is,

Finally, the PSD of transmitted signal should be nonnegative for all frequencies; that is,

Let the four constraints convert the constraint that contains , then the optimization problem can be expressed as where

Writing the optimization problem in a simple form, we can get where

Following we will use Newton method and interior point method to solve the optimization problem. Using a log-barrier function, the optimization can be changed into the following form:

The gradient and Hessian of can be calculated by

Using Newton method to find the optimal solution, Newton step size and Newton descent are needed to calculate. The formula of Newton step size is

The formula of Newton descent is

Then we will use interior point method to solve the optimization problem. Interior point method contains double loop. The outer loop step size is , and the inner loop (Newton loop) step size is . The algorithm can be described as follows:(1)given an initial value , execute Newton loop;(2)given an initial value , and allowable error ;(3)calculate Newton step size (4)calculate Newton descent (5)when , stop Newton loop and return optimal point ; otherwise, given initial value , renew variable ;(6)for the Newton loop returns, certify in the outer loop. If , the Newton loop returns are the global optimal point; otherwise, calculate and execute Newton loop.

For the situation there is clutter, and we can get Using the previous method, we can solve the optimization problem when there is clutter.

#### 4. Simulations

We consider a point target in the radar’s surveillance region with a known impulse response. Suppose that the frequency of the signal is normalized to be . In order to satisfy the Shannon sampling theorem, the sampling frequency is set to 2. The lengths of target impulse response and waveform vector are 63 and 63, respectively. The energy constraint is 1, and power percentage is 0.9. This frequency interval is divided into 2048 non-overlapping frequency bins. The tolerance of Newton method is 10-5. The noise variance is 0.1. The SNR threshold is -5?dB. Tolerance of Newton method and barrier method are 10-5 and 10-5, respectively. The step size increment factor for the barrier method loop is 5. Initial step size for outer loop and inner loop is 5 and 1, respectively. Parameters in backtracking line search method are 0.3 and 0.7, respectively. We use Rayleigh-type radar signature in this paper.

Figure 9 is PSD of Rayleigh target impulse. Figure 10 is PSD of transmitted waveform. Figure 11 is the superposition of PSD of Rayleigh target impulse and transmitted waveform (no clutter). The figure shows that the optimal radar waveforms will spread its energy among most of the spectral peaks of the target response. When mutual information reaches maximum, the peak of PSD of transmitted waveform changes with the peak of PSD of target impulse. So in order to maximize the mutual information in the waveform design of cognitive radar, we should transmit the waveform whose peak of PSD changes with that of target impulse.

Figure 12 is mutual information for all central points. It can be seen that when the central point can accurately approximate the optimal point, the mutual information tends to reach the maximum. The mutual information is greater than that of Bell proposed in .

Figure 13 is superposition of PSD of Rayleigh target impulse and transmitted waveform (in clutter). The figure shows that the optimal radar waveforms will spread its energy among most of the spectral peaks of the target response. When mutual information reaches maximum, the peak of PSD of transmitted waveform changes with the peak of PSD of target impulse. However, as a result of clutter influence, the trend level of the peak of PSD of transmitted waveform to that of target impulse in clutter is weak than no clutter. So the design of transmitted waveform should consider the influence of clutter. The quantitative analysis of clutter to the trend level is necessary.

Figure 14 is mutual information for all central points. It can be seen that when the central point can accurately approximate the optimal point, the mutual information tends to reach the maximum. The mutual information is also greater than that of Bell proposed in . However, the mutual information is lower than that in Figure 11 due to the influence of clutter. So we should consider the quantitative analysis of clutter to the mutual information.

#### 5. Conclusions

Cognitive radar can optimally decide or select the radar waveform for next transmission based on the observations of past radar returns. Adaptive waveform design is an important problem in cognitive radar. In this paper, we propose mutual information model of adaptive waveform design, which can convert the problem of adaptive waveform design into the problem of optimization. We consider two situations of no clutter and clutter and use Newton method and interior point method to solve the optimization problem. From the simulation results, we can see that using the IPM, the mutual information tends to reach the maximum when the central point can accurately approximate the optimal point, and the mutual information is greater than that proposed before. In the next step we should consider the quantitative analysis of clutter to the trend level and the mutual information.

#### Acknowledgments

This work was supported by the National Natural Science Foundation of China, under Grant no. 61004052 and 61104005. It was also supported by the Fundamental Research Funds for the Central Universities, under Grant no. N110323005.

1. S. Haykin, “Cognitive radar: a way of the future,” IEEE Signal Processing Magazine, vol. 23, no. 1, pp. 30–40, 2006. View at: Publisher Site | Google Scholar
2. S. Haykin, “Cognition is the key to the next generation of radar systems,” in Proceedings of the Digital Signal Processing Workshop and 5th IEEE Signal Processing Education Workshop (DSP/SPE '09), pp. 463–467, IEEE Press, 2009. View at: Google Scholar
3. I. Arasaratnam and S. Haykin, “Cubature kalman filters,” IEEE Transactions on Automatic Control, vol. 54, no. 6, pp. 1254–1269, 2009. View at: Publisher Site | Google Scholar | MathSciNet
4. S. Haykin, Y. B. Xue, and T. Davidson, Optimal WaveForm Design For Cognitive Radar, in Asilomar Conference, IEEE Press, 2008.
5. Y. Yang and R. S. Blum, “Waveform design for MIMO radar based on mutual information and minimum mean-square error estimation,” in Proceedings of the 40th Annual Conference on Information Sciences and Systems, pp. 111–116, IEEE Press, 2006. View at: Google Scholar
6. N. A. Goodman, R. P. Venkata, and A. M. Neifeld, “Adaptive waveform design and sequential hypothesis testing for target recognition with active sensors,” IEEE Journal of Selected Topics in Signal Processing, vol. 1, no. 1, pp. 105–1132, 2007. View at: Google Scholar
7. S. P. Sira, A. Papandreou-Suppappola, and D. Morrell, “Dynamic configuration of time-varying waveforms for agile sensing and tracking in clutter,” IEEE Transactions on Signal Processing, vol. 55, no. 7, part 1, pp. 3207–3217, 2007. View at: Publisher Site | Google Scholar | MathSciNet
8. Y. Yang and R. S. Blum, “MIMO radar waveform design based on mutual information and minimum mean-square error estimation,” IEEE Transactions on Aerospace and Electronic Systems, vol. 43, no. 1, pp. 330–343, 2007. View at: Google Scholar
9. A. Leshem, O. Naparstek, and A. Nehorai, “Information theoretic adaptive radar waveform design for multiple extended targets,” IEEE Journal of Selecte Topics in Signal Processing, vol. 1, no. 1, pp. 42–55, 2007. View at: Google Scholar
10. R. A. Romero and N. A. Goodman, “Waveform design in single-dependent interference and application to target recognition with multiple transmissions,” IET Radar, Sonar and Navigation, vol. 3, no. 4, pp. 328–340, 2009. View at: Google Scholar
11. K. S. Kwon, “Waveform design methods for piezo inkjet dispensers based on measured menicus motion,” Journal of Microelectromechanical Systems, vol. 18, no. 5, pp. 1118–1125, 2009. View at: Google Scholar
12. D. J. Kershaw and R. J. Evans, “Waveform selective probabilistic data association,” IEEE Transactions on Aerospace and Electronic Systems, vol. 33, no. 4, pp. 1180–1188, 1997. View at: Google Scholar
13. S. P. Sira and D. Cochran, “Adaptive waveform design for improved detection of low-RCS targets in heavy sea clutter,” IEEE Journal of Selected Topics in Signal Processing, vol. 1, no. 1, pp. 56–66, 2007. View at: Google Scholar
14. C. Rago, P. Willett, and Y. Bar-Shalom, “Detecting-tracking performance with combined waveforms,” IEEE Transactions on Aerospace and Electronic Systems, vol. 34, no. 2, pp. 612–624, 1998. View at: Google Scholar
15. M. R. Bell, “Information theory and radar waveform design,” IEEE Transactions on Information Theory, vol. 39, no. 5, pp. 1578–1597, 1993. View at: Google Scholar