24小时热门版块排行榜    

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

Medichen

铜虫 (小有名气)

[求助] |A+a*B|对a的导数怎么求? 已有2人参与

如题,A和B均为已知的对称矩阵,a为一标量,||表示行列式。
那么|A+a*B|对a的导数怎么求呢?
如果特殊一点的话,A为一对角阵,那么根据特征分解,其整个行列式对a的导数很容易通过特征值表示出来。但对于更一般的对称阵A貌似通过特征分解不好求,请问大家有其他解法吗? 谢谢大家了
(ps. 这不是数学作业题,而是我在做算法时遇到的一个问题)
回复此楼
笨鸟先飞!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zaq123321

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
Medichen: 金币+5, ★★★★★最佳答案 2015-10-13 21:09:54
Jacobi's formula see https://en.m.wikipedia.org/wiki/Jacobi%27s_formula

[ 发自手机版 http://muchong.com/3g ]
小木虫给我温暖,给我希望,爱就要爱小木虫。
3楼2015-10-13 18:12:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

0404600213

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
我的思路是这样的
根据导数的定义|A+aB|的导数应该是(|A+(a+da)B|-|A+aB|)/da当da趋于0时的极限

可以把|A+(a+da)B|=|(A+aB)+daB|用行列式加法展开成若干个行列式的和,而且每个行列式中任意一列要么是|A+aB|当中的对应列要么是|daB|当中的对应列
易知,当上述行列式中有两个或以上|daB|中的对应列时,这个行列式就是da的高阶无穷小;当上述行列式中没有|daB|中的对应列时,就是|A+aB|

所以|A+aB|的导数就是n个行列式的和,期中每个行列式中都有一列是|B|当中的对应列,其他列都是|A+aB|当中的对应列

因为上面的过程中并没有用到对称这个条件,所以这个结果应该还可以化简。



PS:另外一个思路是|A+aB|一定是一个关于a的n次多项式,可以考虑它的n个根……但是我没有想到比较好的方法。
2楼2015-10-13 17:47:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Medichen

铜虫 (小有名气)

引用回帖:
3楼: Originally posted by zaq123321 at 2015-10-13 18:12:00
Jacobi's formula see https://en.m.wikipedia.org/wiki/Jacobi%27s_formula

非常感谢!要得就是这个答案!
笨鸟先飞!
4楼2015-10-13 21:09:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Medichen

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by 0404600213 at 2015-10-13 17:47:22
我的思路是这样的
根据导数的定义|A+aB|的导数应该是(|A+(a+da)B|-|A+aB|)/da当da趋于0时的极限

可以把|A+(a+da)B|=|(A+aB)+daB|用行列式加法展开成若干个行列式的和,而且每个行列式中任意一列要么是|A+aB|当中 ...

谢谢你!楼下给的链接有直接的公式可以用。
笨鸟先飞!
5楼2015-10-13 21:26:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +5 scalable 2026-08-24 5/250 2026-08-24 17:03 by spiderzhuzhi
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 19/950 2026-08-24 16:52 by iaeyuan
[基金申请] 能否退出参与的面上项目解除限项 +19 koalala 2026-08-24 21/1050 2026-08-24 16:48 by Wei_ren
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 我面上完蛋了 +8 且听虎啸 2026-08-20 9/450 2026-08-24 13:52 by zzuzxg
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见