在perl中指數求一個矩陣

我正在尋找一個使用perl來求矩陣的函數。我搜查了CPAN，但似乎無法找到任何合適的東西。有誰知道一個合適的包裝，或簡單的功能？在perl中指數求一個矩陣

例如，對於下面的矩陣（$ M）：

$m = [[-1.85294,0.36417,1.34865,0.14012], 
     [0.13385,-0.70885,0.16502,0.40998], 
     [0.48588,0.16174,-0.77471,0.12709], 
     [0.13424,1.06859,0.33797,-1.54081]];

E 1米〜=

0.25438971 0.2029723 0.4580481 0.08458940 
0.07460393 0.6061030 0.1588883 0.16040370 
0.16502112 0.1557330 0.6025270 0.07671839 
0.08103971 0.4180838 0.2040153 0.29685570

來源

2012-02-23 wannymahoots

http://search.cpan.org/~ellipse/PDL-LinearAlgebra-0.06/Trans/trans.pd#mexp或http://www.davideous.com/matrixexp/matrix.pm？ – 2012-02-23 18:18:07

PDL將處理這一點，因爲@ J.D。指出。

use strict; 
use warnings; 

use PDL; 
use PDL::LinearAlgebra::Trans; 
use PDL::IO::Dumper; 

my $m = pdl [[-1.85294,0.36417,1.34865,0.14012], 
      [0.13385,-0.70885,0.16502,0.40998], 
      [0.48588,0.16174,-0.77471,0.12709], 
      [0.13424,1.06859,0.33797,-1.54081]]; 

my $new = mexp $m; 

print sdump $new;

結果：

{my($VAR1); 
my($PDL_140600590505856) = (double(0.254389714179579,0.202972289234997,0.458048094091295,0.08458940315894629,0.0746039289026906,0.606102994157082,0.158888283526803,0.160403695685552,0.165021120118487,0.155733015860466,0.602527027489139,0.07671838765811791,0.08103970821403091,0.418083774202018,0.204015347133539,0.296855703714981)->reshape(4,4)); 


$VAR1 = $PDL_140600590505856; 

}

的double()函數內的比特是結果矩陣。

來源

2012-02-23 20:14:38 friedo

這樣做。它似乎打印了三條關於被重定義的子路由器的警告。我相信它沒有什麼可擔心的，但是你知道一種防止這種情況的方法嗎？謝謝！ – wannymahoots 2012-02-23 20:59:02

奇怪。我沒有得到任何警告，在5.14.2下運行它。你使用的是什麼版本的Perl和PDL？ – friedo 2012-02-24 00:30:17

我正在運行perl v5.12.4，但可能會出現一些與我的安裝PDL（主要安裝在Ubuntu上使用apt-get，但從源安裝PDL :: LinearAlgebra :: Trans）。我會調查...... – wannymahoots 2012-02-24 12:00:13

在perl中指數求一個矩陣

回答

相關問題