Any pointer for analyzing this matrix?

  • Context: Graduate 
  • Thread starter Thread starter divB
  • Start date Start date
  • Tags Tags
    Matrix
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
3 replies · 2K views
divB
Messages
85
Reaction score
0
Row subset of matrix with lowest condition number

Hi,

I have a matrix with [itex]N[/itex] rows and [itex]K\cdot M \ll N[/itex] columns.
Is there a generic way to choose [itex]L[/itex] rows s.t. the condition number of the matrix is minimized?
I would appreciate any pointer or any hints how I could approach the problem.
Currently I am brute forcing all possible subsets and seek for the set with the smallest condition number but this is nearly undeasible.

I do not think it helps but if it does, this is the structure of the matrix: Is is a pointwise multiplication of 2 matrices. The first summand consists of [itex]K[/itex] vertical "ribbons" of size [itex]N \times M[/itex]. Each consists of [itex]M[/itex] equal columns where the columns are the DFT of a power of a discrete random(like) signal (I do not know if this matters ...)
The second summand consists of [itex]K[/itex] equal Vandermonde matrices (to be precicse, the first [itex]M[/itex] columns of the DFT matrix).
Formally:

[tex] \left[\begin{matrix} \mathbf{A}_1 & \mathbf{A}_2 & \cdots & \mathbf{A}_K \end{matrix}\right] * \left[\underbrace{\begin{matrix}\mathbf{V} & \mathbf{V} & \cdots & \mathbf{V} \end{matrix}}_{K-\mathrm{times}}\right][/tex]

where [itex]*[/itex] denotes point-wise multiplication.

[tex] \mathbf{A}_k = \left[\underbrace{\begin{matrix}<br /> \mathrm{DFT}(x_n^{k-1})^T & \cdots & \mathrm{DFT}(x_n^{k-1})^T<br /> \end{matrix}}_{M-\mathrm{times}}\right][/tex]

[tex] \mathbf{V} = \left[\begin{matrix}<br /> 1 & 1 & 1 & \cdots & 1 \\<br /> 1 & u^1 & u^2 & \cdots & u^M \\<br /> 1 & u^2 & u^4 & \cdots & u^{2M} \\<br /> \vdots & \ddots & \ddots & \ddots & \vdots \\<br /> 1 & u^N & u^{2N} & \cdots & u^{NM}<br /> \end{matrix}\right][/tex]

Thanks,
divB
 
Last edited:
Physics news on Phys.org
divB said:
I have a matrix with [itex]N[/itex] rows and [itex]K\cdot M \ll N[/itex] columns.
Is there a generic way to choose [itex]L[/itex] rows s.t. the condition number of the matrix is minimized?

If you choose [itex]L \ne K\cdot M[/itex], I'm curious how you define the "condition number" of the resulting non-square matrix.
 
Hi Stephen,

Thanks. I have no formal requirement for this. The system should be solved via LS and stability and accuracy increases as the condition number decreases. Which definition is exactly used is not so relevant to me; the ratio between largest and smallest singular value works pretty well (cond in MATLAB) and of course works for any, non-square matrix.

divB