Least-Squares Design

  • Example data

  • Least-squares approach

  • Data and example

  • Trade-off curve

Example data

For our particular example, we have the following parameters:

  • Number of antennas: n = 16.

  • Wavelength: lambda = 8.

  • Pass-band size: Delta = pi/6. In what follows, we use the following matlab applet to define the data of our problem.

Matlab data
N = 10; % discretization parameter
n = 16; % number of antennas
lambda = 8; % wavelength
Phi = pi/6; % sidelobe parameter

Least-squares approach

We examine how least-squares can be used to strike a trade-off between the sideloble level attenuation and thermal noise power. The initial problem reads

 min_z : |z|_2 ~:~ mbox{bf Re}(D_z(0)) ge 1, ;; |D_z(phi_i)| le delta, ;; i=1,ldots,N.

The basic idea is to penalize the constraints, as follows: we choose a ‘‘trade-off parameter’’ mu>0 and solve the problem

 min_z : |z|_2^2 + mu sum_{i=1}^N |D_z(phi_i)|^2 ~:~ mbox{bf Re}(D_z(0)) ge 1.

Remember that the function D_z is linear in z:

 D_z(phi) := sum_{k=1}^n z_k exp( frac{2pi jmath k cos(phi)}{lambda} )

Hence, the above problem is a linearly constrained least-squares problem.

CVX implementation

The above problem can be solved via linear algebra (SVD) methods. Here, we do not even bother to invoke SVD — CVX will do. A CVX implementation of this problem is given below. Note how we use a new vector variable to handle the objective function. Note also that CVX understands complex variables.

CVX implementation
Angles = linspace(Phi,pi,N); % angles in the stop band
a = 2*pi*sqrt(-1)/lambda; % intermediate parameter
cvx_begin
    variable z(n,1) complex;
    variable r(N,1) complex;
    minimize( z'*z+mu*r'*r )
    subject to
        for i = 1:N,
            exp(a*cos(Angles(i))*(1:n))*z == r(i);
        end
        real( exp(a*(1:n))*z) >= 1;
cvx_end
alt text 

Antenna array design via least-squares: minimal combination of thermal noise and sidelobe level violations, computed at N=90 points, with trade-off parameter mu = 0.5. The corresponding sidelobe attenuation level is delta = -7.5386 mbox{dB}, and the thermal noise power is |z|_2 = 0.5409.

Note that the penalty approach adopted here will not enforce the desired constraints. We can only hope that, for mu large enough, all of the terms in the sum will go below the desired threshold delta. In practice, we expect that some terms might not. It seems to be hard to achieve the best trade-off this way. The following plot shows the trade-off curve that we can achieve.

alt text 

Trade-off curve achieved via via least-squares.