× 评论会员:little_tiger 时间:2011/11/30
喜保,
伟大的工作。我只是想知道,如果你有一些这些代码的更新版本?有一个C版本?
非常感谢很多
× 评论会员:游客 时间:2011/11/30
。TharinduNishadPatikirikorala
嗨,
我用这个DLL真正有效地在一个研究项目。我想使我的工作访问,以免费的其他研究人员。但与依赖的DLL,我不知道我是否可以做不做?
有一个复制权的问题,如果我与我的DLL分发您的DLL?
谢谢你
Tharindu
× 评论会员:保罗Selormey 时间:2011/11/30
您好Tharindu,
正如你可以看到在许可部分,这是一个公共域代码,
这样你就可以在任何你喜欢它。
你不需要有一个DLL,您可以直接在项目中使用的源代码。
的源代码是很老了,你可以清洁并保持您的任何变化,这是
公共领域的东西,只要保持在源代码本身的任何牌照的通知。
最好的问候,保罗。 耶稣基督是爱!请告诉别人。
× 评论会员:游客 时间:2011/11/30
TharinduNishadPatikirikorala
感谢您的澄清
× 评论会员:muhaa222 时间:2011/11/30
喜
我用你的图书馆。这是非常好的代码{S2}
但我仍是一个问题
我想DELET DA行或一个山坳?
我没#T已经发现了这个问题soulution
CYA
muhaa
× 评论会员:andeezlemoyo 时间:2011/11/30
喜保,首先,感谢您分享您的工作
这里我要解决一个样品的线性系统:
49X - 7Y ž = 6
25X - 5年Z = 4
9X - 3Y Z = 2
我怎样才能解决与GeneralMatrix此?我将不胜感激,如果你能给我一些示例代码。
,
安迪
× 评论会员:muhaa222 时间:2011/11/30
嗨,我认为在解决是
GeneralMatrix一个
像
49 -7 1
25 -5 1
9 -3 1 GeneralMatrix B
像
6
4
2
GeneralMatrix N = A。TRANSPOSE()*一个;GeneralMatrix N = A. TRANSPOSE()* B
GeneralMatrix X = N。逆()* N
x.GetElement(0,0)是你的X
x.GetElement(1,0)是你的Ÿ
x.GetElement(2,0)是您的z
我使用的求解超定方程,但仍是工作你的问题
CYA
muhaa | LachlanGro:
× 评论会员:游客 时间:2011/11/30
您好,我有一个3x3矩阵,V:V=100101 625251900301在Matlab中,我做了QR和得到:GT,GT[Q,R]=QR(V,0〜)Q=-0.09090.87890.4683-0.56800.3405-0.7493-0.8180-0.33410.4683R=1.0E003*-1.1003-0.0396-0.001500.00730.0009000.0002使用DotNetMatrix然而,当我得到:矩阵Q=-0.09088562141317660.878886544558686-0.468292905790847-0.5680351338323540.3404887339694660.749268649265355-0.817970592718589-0.334104570207539-0.468292905790847而价值观是好的,右栏是一个负的Matlab结果。有什么建议?关于拉克兰
。保罗Selormey
× 评论会员:游客 时间:2011/11/30
LachlanGro写道:"当值都不错,右栏是所有的Matlab结果负。没错,我已经验证它。我也用在维基百科imgsrc=http://www.orcode.com/upimg/2011_11_30_20_49_53_3.gif例如矩阵]和标志不符合LachlanGro写道:"有什么建议我目前看,我不会有任何问题的根源,立即思想。最好的问候,保罗。耶稣基督是爱!请告诉别人。
LachlanGro
× 评论会员:游客 时间:2011/11/30
喜保,我贾马讨论列表张贴在一个类似的问题,显然是没有错误,只是解决的办法是"不是唯一的。"虽然我相信是这样的话,我无法找到的会是什么独特的迹象。无论哪种方式,我正是通过相同的数据到其他库(如BlueBit,MATLAB),结果是罚款。的问候,拉克兰
保罗Selormey
× 评论会员:游客 时间:2011/11/30
LachlanGro写道:我贾马讨论列表上发布了类似的问题,显然是没有错误,只是解决的办法是"不是唯一的"。信息LachlanGro写道:虽然我相信在这种情况下,我无法找到一个什么会是唯一的指示成果的基础上计算是有效的,也Q'Q=我,所以可能出现的情况。LachlanGro写道:无论哪种方式,我正是通过相同的数据到其他库(如BlueBit,MATLAB),结果是罚款目前,我不会有任何数学包,所以我只是寻找书籍和网上的例子来比较。我会考虑见惯了的算法,如果我们能够得到相同的结果如Matlab,以防止任何疑问。最好的问候,保罗。耶稣基督是爱!请告诉别人。
Kuryn
× 评论会员:游客 时间:2011/11/30
为什么使一个类的虚?这会导致性能问题。你不应该使一个类虚拟
保罗Selormey
× 评论会员:游客 时间:2011/11/30
Kuryn写道:"你不应该让一个类虚拟。哼哼,我不知道这个镇的新规则。然而,你的源代码,所以请做任何你喜欢它-它是一个公共域许可证最好的问候,保罗。耶稣基督是爱!请告诉别人。
Kuryn
× 评论会员:游客 时间:2011/11/30
你应该在需要的时候,一类虚拟。但不是数学类,因为没有什么变化。就像在一个Vector,一个dotProduct是一个dotProduct
。FueledByIgnorance
× 评论会员:游客 时间:2011/11/30
!感谢这个库,它是非常有用的我有一个问题:我在某个时刻到达矩阵是奇异的,我要执行一个逆。{C}逆方法调用解决,这将实例化一个LUDecomposition对象,并调用它的解决方法。它测试的奇点,和我得到的是奇异的"黑客帝国"的例外。我有非常有限的数学知识,所以包涵。在被认为是矩阵的奇异点,不应该由执行一个伪逆在反演方法的意见,它说:"逆(a)如果是方形的,否则伪逆"。有一个从实施缺少的一部分?感谢:保罗Selormey
您好,单一性和伪逆是不同的。要了解伪逆,看到下面的
1。
2。 {A5}
最好的问候,
保罗。
耶稣基督是爱!请告诉别人。
× 评论会员:ajwo 时间:2011/11/30
您好:
SVD分解有一个错误,从贾马继承。问题是,当n> M. V矩阵这次调整是:
if (wantv)
{
for (int k = n - 1 ; k >= 0 ; k--)
{
if ((k < nrt) & (e[k] != 0 .0 ))
{
for (int j = k + 1 ; j < n; j++) {
double t = 0 ;
这里报告错误是:
http://cio.nist.gov/esd/emaildir/lists/jama/msg00430.html
× 评论会员:保罗Selormey 时间:2011/11/30
谢谢你,我会看一看
。
最好的问候,
保罗。
耶稣基督是爱!请告诉别人
× 评论会员:iotha 时间:2011/11/30
谢谢你!
我没有数学的家伙,但我有一个问题后,另一件事(可能是一个"错误"?):{ BR}
SVD的是A [MXN] U = [MXM] ・ S [MXN] * V [NXN] ^ T - 所以应该是一个方阵。但它在这里被定义为:
M = Arg.RowDimension;
N = Arg.ColumnDimension;
NU = System.Math.Min INT(M,N);
新的双U = [M ][];{ BR}为(int i = 0;我LT,M,I)
{
U [I] =新的双[怒江]
难道不应该是这样的:
M = Arg.RowDimension;
N = Arg.ColumnDimension;
NU = System.Math.Min INT(M,N);
S =新的双[System.Math.Min(M 1,N)]
新的双U = [M ][];{ BR}为(int i = 0;我LT,M,I)
{
U [I] =新的双[M]
}
× 评论会员:celisdelorena 时间:2011/11/30
大家好,
我在C#中使用的矩阵,我想感谢和祝贺的一般矩阵作者。这是一个非常有用的和容易的类工作。
目前我正与SVD的descomposition获得的伪逆...但它不工作。在奇异值计算的代码块。没有任何人有同样的问题吗?
由于事先
× 评论会员:raz0rblade 时间:2011/11/30
您好,
我遇到的两个问题解决欠定方程组:
第一个问题是一个ArrayIndexOutOfBoundsException异常,当试图获得S矩阵。我认为这个问题是因为A [MXN] U = [MXM] ・ S [MXN] * V [NXN]"的尺寸应该是M - N,但在他们的代码N - S的大小错误, - N: / / /返回的奇异值的对角矩阵/ / /
/ / /
虚拟公共GeneralMatrix S
{
得到
{
GeneralMatrix X =新GeneralMatrix(N,N); / /这应该是(M,N)
双[] [] = X的数组
为(int i = 0;我LT,N,I)/ /这个循环应该去,直到没有N
{
(研究= 0,J LT,N,J)
{
S [I] [J] = 0.0;
}
S [I] [I] = this.s [I]
}
返回X;
}
然后有一个问题:SingularValueDecomposition构造函数的奇异值向量维
S =新的双[System.Math.Min(M 1,N)] 例如,我的方程组有六个方程和8个变量,所以M = 6和n = 8。
则S将七个大小,将会有一个奇异值并不真的存在。
除了这个,我得到不同的结果时,使用Matlab的SVD功能。的奇异值和U型矩阵是完全一样的,但在V -矩阵的最后纳米coloumns有不同的价值观,其他人一样在Matlab相同。但我没有找到一个合适的地方在代码这个问题。
请纠正我,如果我理解数学错误的东西。 {S2}
修改2008年8月29日,上午10:16
× 评论会员:mbreidt 时间:2011/11/30
不知道这个讨论板仍然活跃,但最近我发现这个库,并发现它非常有用的,所以首先感谢非常它!
不幸的是,我还发现了与奇异值计算的问题,错误的尺寸和所产生的矩阵。
如果你做一个测试的例子{A6}],这个库产生错误的结果,至于我可以告诉(FWIW,MATLAB返回相同的维基百科上的给定值)。
在运行的例子(一个4 × 5矩阵的SVD),从S矩阵的SVD对象是一个5 × 5的矩阵,而正确的结果应该是一个4 × 5的矩阵。此外,U和V矩阵中的条目似乎是在错误的顺序和/或有错误的迹象。
此外,SVD的对象返回的奇异值向量在这个例子中有五个要素,但是,根据维基百科
所以应该是最高的4个奇异值,而不是五年。
如果有人仍然是这个棋盘上,我非常喜欢听是否任何人都可以证实我的发现和/或有错误修复!对我的工作,我真的很喜欢用这个库,但只有与非方阵的作品。
谢谢!
马丁
× 评论会员:Cukaricn 时间:2011/11/30
我与特征值和特征向量之间的关系问题。我的特征向量不符合他们的特征值。相反,它们混合(A1 exaple特征值向量V3,而不是V1)。我所拥有的所有特征值都是负的,所以可以说是问题吗?
感谢
Cukaric | sergiomad:
× 评论会员:游客 时间:2011/11/30
?是为获得KarhunenLoève向量的方法我想翻译的Java{A7}这段代码,使用JAMA,但我无法找到像方法"时代()"谢谢
保罗Selormey
× 评论会员:游客 时间:2011/11/30
?sergiomad写道:可用的方法获得KarhunenLoève向量,Karhunen-Loève是没有直接关系矩阵。sergiomad写道:"但我不能找到类似的方法"次()"此端口使用。NET命名的,它应该是"乘(...)"{BR}最好的问候,保罗。耶稣基督是爱!请告诉别人