数学建模社区-数学中国

标题: 【求助】怎么让mathematica真正算出结果？ [打印本页]

作者: stalco 时间: 2007-11-17 12:32
标题: 【求助】怎么让mathematica真正算出结果？

想用函数MatrixConditionNumber[]求矩阵条件数
不知道为什么结果不返回一个数而仍是表达式？
等待高人解答
bow

例子：
In[30]:=
MatrixConditionNumber[{{1,0},{0,1}},2]

Out[30]=
MatrixConditionNumber[{{1,0},{0,1}},2]

作者: 00704050 时间: 2007-11-17 18:38
用N[[表达式],位数]

作者: madio 时间: 2007-11-18 01:33

主要是因为这个函数所在的包并不是被直接载入Mathematica的核心中，需要你先运行包文件！另外我运行你的例子确实不给计算，你可以试试别的例子！下面是包文件，解压缩到Mathematica的AddOns\Packages目录下，运行这个包文件就可以使用你提到的函数了！[attach]3676[/attach]
你可以试试这个例子：

In[108]:= BandedMatrix[n_] := Module[
   {A, spar, diag, udiag, udiag2, ldiag, ldiag2},
   diag = Table[{i, i} -> 14., {i, 1, n}];
   udiag = Table[{i, i + 1} -> 3., {i, 1, n - 1}];
   udiag2 = Table[{i, i + 8} -> -1., {i, 1, n - 8}];
   ldiag = Table[{i + 1, i} -> 5., {i, 1, n - 1}];
   ldiag2 = Table[{i + 500, i} -> -4., {i, 1, n - 500}];
   spar = Flatten[{diag, udiag, udiag2, ldiag, ldiag2}];
   A = SparseArray[spar, {n, n}]
   ];

In[109]:= MatrixConditionNumber[BandedMatrix[100]]

Out[109]= 4.6

matrixmanipulation.rar.rar

6.42 KB, 下载次数: 8, 下载积分: 体力 -2 点

matrixmanipulation.rar

欢迎光临数学建模社区-数学中国 (http://www.madio.net/)