Perl 中的矩阵矩阵

Question

在我正在处理的 Perl 脚本中，我需要用其他几个矩阵构建一个矩阵。 我查看了 CPAN 中的几个模块（数学： ：矩阵，PDL::矩阵, Math::Cephes::Matrix< /a>），但这些似乎都不支持这一点。

在 Octave 中，这非常容易。 这是一个与我想做的事情类似的例子：

octave:1> A = [ 1, 2; 3, 4 ]
A =    
   1   2
   3   4

octave:2> B = [ 5, 6; 7, 8 ]
B =    
   5   6
   7   8

octave:3> C = [ 9, 10; 11, 12 ]
C =    
    9   10
   11   12

octave:4> D = [ 13, 14; 15, 16 ]
D =    
   13   14
   15   16

octave:5> E = [ A, B; C, D ]
E =    
    1    2    5    6
    3    4    7    8
    9   10   13   14
   11   12   15   16

似乎自己尝试这样做会很快变得混乱，这可能就是为什么这些模块不支持它的原因......还有其他人有过这样的经历吗？需要这个吗？ 你解决了吗？

Answer 1

自己滚动并不太痛苦。

use List::Util qw(max);

@A = ([1, 2], [3, 4]);
@B = ([5, 6], [7, 8]);
@C = ([9, 10], [11, 12]);
@D = ([13, 14], [15, 16]);

sub hmerge(\@\@;\@\@\@\@\@\@) {
    my @ret;
    for my $i (0 .. max map $#$_, @_) {
        push @ret, [map @{$_[$i]}, @_];
    }
    @ret;
}

@E = (hmerge(@A, @B), hmerge(@C, @D));

Answer 2

Perl 数据语言 (PDL) 版本 2.4.10 支持 pdl 构造函数，以及 append 和 glue 例程可用于将子数组粘贴在一起，如 pdl2 会话所示

pdl> $A = pdl q[ 1, 2 ; 3, 4 ];       # pdl constructor with string argument

pdl> $B = pdl q[ 5, 6 ; 7, 8 ];       # pdl constructor with string argument

pdl> $C = pdl q[ 9, 10 ; 11, 12 ];    # pdl constructor with string argument

pdl> $D = pdl q[ 13, 14 ; 15, 16];    # pdl constructor with string argument

pdl> ?vars
PDL variables in package main::

Name         Type   Dimension       Flow  State          Mem
----------------------------------------------------------------
$A           Double D [2,2]                P            0.03KB 
$B           Double D [2,2]                P            0.03KB 
$C           Double D [2,2]                P            0.03KB 
$D           Double D [2,2]                P            0.03KB 


pdl> p $A, $B, $C, $D;

[
 [1 2]
 [3 4]
]

[
 [5 6]
 [7 8]
]

[
 [ 9 10]
 [11 12]
]

[
 [13 14]
 [15 16]
]

pdl> p $AB = $A->append($B);          # concatenate horizontally (actually on dim(0))

[
 [1 2 5 6]
 [3 4 7 8]
]

pdl> p $CD = $C->append($D);          # concatenate horizontally (actually on dim(0))

[
 [ 9 10 13 14]
 [11 12 15 16]
]

pdl> p $E = $AB->glue(1,$CD);         # glue vertically (actually on dim(1))
[
 [ 1  2  5  6]
 [ 3  4  7  8]
 [ 9 10 13 14]
 [11 12 15 16]
]

： a href="http://pdl.perl.org/content/pdl-book-toc.html" rel="nofollow">PDL 书 和 PDL 邮件列表是有关 PDL 的更多信息的重要来源。

Answer 3

编辑

我误解了OP，认为他们想要迭代几个矩阵的所有可能的排列（这就是 Iterator::Array::Jagged 所做的）。

---

看看 Iterator::Array::Jagged

这是一个概要示例：

use Iterator::Array::Jagged;

# Build up a set of data:
my @data = (
  [qw/ a b /],
  [qw/ c d /],
  [qw/ e f g /]
);

# Iterator is a subref:
my $itersub = Iterator::Array::Jagged->get_iterator( @data );
while( my @set = $itersub->() )
{
  print "Next set: '" . join("&", @set) . "'\n";
}# end while()

上面代码中的示例打印以下内容：

Next set: 'a&c&e'
Next set: 'b&c&e'
Next set: 'a&d&e'
Next set: 'b&d&e'
Next set: 'a&c&f'
Next set: 'b&c&f'
Next set: 'a&d&f'
Next set: 'b&d&f'
Next set: 'a&c&g'
Next set: 'b&c&g'
Next set: 'a&d&g'
Next set: 'b&d&g'