24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2097  |  回复: 16
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cyq101600

铁杆木虫 (初入文坛)

[交流] 【求助】矩阵导数【已解决】已有12人参与

请问有人研究过矩阵求导问题吗?现在遇到一个矩阵导数不知如何解决:
F(X)=X-1aaT X-1,其中X是矩阵,a是列向量,(-1)表示矩阵的逆,T为矩阵转置,要F(X)对矩阵X求导。(矩阵对矩阵求导)

看了一些资料,发现矩阵求导有不同的定义方式,有常规的定义为分别对矩阵的每个元素求导平铺生成超级矩阵,如http://muchong.com/bbs/viewthread.php?tid=1450822中所讨论的。然而,根据Magnus的讨论,这种常规的定义方法没有意义,主张把矩阵先按列排成一个列向量,然后类似Jacobian对向量进行求导(详见http://en.wikipedia.org/wiki/Talk:Matrix_calculus上的讨论)。不知如何是好!

[ Last edited by cyq101600 on 2010-6-12 at 16:29 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

cyq101600

铁杆木虫 (初入文坛)

引用回帖:
Originally posted by 小雨萌萌 at 2010-06-08 08:48:51:
看你怎么定义矩阵的范数,求导数和范数有关系

谢谢参与!我想求导是和范数定义有关,但通常矩阵并没有定义什么范数下。
一个最常见的例子,F(X)=aTXb,通常给出的结果是DF(X)=abT。
但是按照Magnus主张的向量求导,结果是DF(X)=kron(aT,bT),结果都不一样。尤其在矩阵函数求导中常规的算法链式法则和乘法法则都可能不成立了,也就是Magnus指出的问题,不知道该采用哪种算法
3楼2010-06-08 09:42:39
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

cyq101600

铁杆木虫 (初入文坛)

矩阵求导很多地方有可能会遇到,给大家提供一个链接,https://ccrma.stanford.edu/~dattorro/matrixcalc.pdf 这里面的附录(Table)总结了各种求导的公式。但是要提醒的是,这里的求导都是按照平铺矩阵块这种思路的,有可能遇到乘法法则不成立的问题,这也是本贴求助的问题。另一种就是Magnus的定义,向量算法。但是不知道该用哪种定义算法。
5楼2010-06-08 11:11:11
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

cyq101600

铁杆木虫 (初入文坛)

引用回帖:
Originally posted by Pchief at 2010-06-08 12:03:02:
楼主看看我这个如何?这个是从非线性泛函分析里拿来的东东


8楼2010-06-08 13:28:16
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

cyq101600

铁杆木虫 (初入文坛)

引用回帖:
Originally posted by saladin983 at 2010-06-08 16:15:38:


这个例子是个普通的单值函数,处理起来倒不是很麻烦。我算了个2×2的例子,Kronecker积构造的导数结果是一致的。不过我觉得用这个做定义就太不直观了,单纯用于计算可能还行。

Fréchet导数拿来定义看起来 ...

是的,矩阵导数定义本身就存在很多争议,计算法则也不一。
10楼2010-06-08 17:14:20
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cyq101600 的主题更新
信息提示
请填处理意见