Finding Jordan Canonical Form (V and J) of a big square matrix

Question

Rahul Singh 2015-6-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/224759-finding-jordan-canonical-form-v-and-j-of-a-big-square-matrix

编辑： Rahul Singh 2015-6-25

While using [V J] = jordan(A); for a 33x33 matrix A, I am getting the following error.

Code:

if true
  % code
end
N = 33;
W = zeros(N,N);
for i=1:N 
  for j = 1:N
      if(i ~= j)
          if(round(rand - 0.1))
              W(i,j)  = round(5*rand);
          end
      end
  end
end
[V J] = jordan(W);

Command Window:

if true
  % code
end
Warning: Found roots of the minimal polynomial that cannot be determined
in terms of radicals. [linalg::jordanForm] 
Error using mupadmex
Error in MuPAD command: The operand is invalid. [_index]
Evaluating: symobj::jordan
Error in sym/mupadmexnout (line 879)
      out = mupadmex(fcn,args{:});
Error in sym/jordan (line 34)
      [Vsym,Jsym] = mupadmexnout('symobj::jordan',A,'All');
Error in jordan (line 25)
 [V,J] = jordan(sym(A));

When I increase the size of the matrix A, I get the following error

if true
  % code
end
Error using mupadmex
Error in MuPAD command: Similarity matrix too large.
Error in sym/mupadmexnout (line 879)
      out = mupadmex(fcn,args{:});
Error in sym/jordan (line 34)
      [Vsym,Jsym] = mupadmexnout('symobj::jordan',A,'All');
Error in jordan (line 25)
 [V,J] = jordan(sym(A));

How can I get the jordan decomposition of bigger matrices say a maximum of 200x200 matrix??

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Rahul Singh 2015-6-22

Please help!

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Mukul Rao 2015-6-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/224759-finding-jordan-canonical-form-v-and-j-of-a-big-square-matrix#answer_183731

在 MATLAB Online 中打开

The Jordan function has an imposed size limit to help prevent exceedingly long calculations. In order to get around the error, execute this function instead:

>> feval(symengine, 'linalg::jordanForm', A, 'All')

Where "A" is the matrix which you are analyzing. Note that while this will bypass the size limitation, the calculation can take a very long time. If you are just interested in the eigenvalues and eigenvectors of the matrix, consider using the "eig" function instead.