作者celestialgod (天)
看板MATLAB
标题Re: [讨论] 取出 Jordan blocks
时间Mon Aug 24 16:16:55 2015
※ 引述《math99 (新世界)》之铭言：
: 我有一个 Jordan form , 我想要取出各别的 Jordan blocks
: 请教 matlab 有什麽方便的做法吗？
: 目前我做法只能用 for loop 一个一个取, 希望能够不要用到 loop
: 例
: Input:
: A = [2 0 0 0;
: 0 1 1 0;
: 0 0 1 0;
: 0 0 0 2]
: Output:
: J1 = [2], J2=[1,1;0,1], J3=[2];
: 感谢
indices = mat2cell(1:size(A, 1), 1, [1, 2, 1]);
out = cellfun(@(x) A(x,x), indices, 'UniformOutput', false);
out会有三个element, [2], [1,1;0,1], [2]
至於[1, 2, 1]，目前没有想到更好的方法去抓出来
後来想到用下面的方法来做，目前测试还未有BUG出现
有其他问题请麻烦通知我
% matrix generation
bksize = randi(3, 4, 1);
As = arrayfun(@(x) rand(x), bksize, 'UniformOutput', false);
A = blkdiag(As{:});
% extract blocks
[I, J] = ind2sub(size(A), find(A ~= 0));
I = sort(I);
repeatLength = diff([0; find(I(2:end) ~= I(1:end-1)); length(I)]);
repeatLoc = [find(repeatLength(2:end) ~= repeatLength(1:end-1)); ...
length(repeatLength)];
bksize = repeatLength(repeatLoc);
if sum(bksize) ~= size(A,1)
tmp = diff([0; repeatLoc]) ./ bksize;
bksize = cell2mat(arrayfun(@(x, y) x * ones(y, 1), bksize, tmp, ...
'UniformOutput', false));
end
indeices = mat2cell(1:size(A, 1), 1, bksize);
out = cellfun(@(x) A(x,x), indeices, 'UniformOutput', false);
isequal(As, out') % true
我稍微改写了一下 放到matlab的file exchange去了~~
网址：
http://tinyurl.com/oev2p98
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 123.205.27.107
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/MATLAB/M.1440404218.A.9A1.html
1F：推 sunev: jordan form只要判断diag(A,1)即可 08/24 18:43
我已经忘记jordan form是啥了XDDD 四年前学的东西 没用到过~"~
写成general form就都可以用了~~~
2F：推 sunev: A=[1 0 1;0 1 0;0 0 1]; ?08/25 12:47
这个我没办法QQ
3F：→ math99: 感谢 好漂亮用法08/25 17:58
4F：→ math99: s大举得例我应该不会遇到08/25 17:59
5F：推 math99: 在请教一下 cellfun 参数 'UniformOutput' 作用是什麽?08/25 18:03
如果每一个输出不是长度为一的向量，就要把这个选项设定false
※ 编辑: celestialgod (123.205.27.107), 08/25/2015 18:13:54
6F：推 math99: ok, thanks 08/25 18:16