24小时热门版块排行榜    

查看: 2165  |  回复: 4

Medichen

铜虫 (小有名气)

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

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

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的回帖

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的回帖

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的回帖
相关版块跳转 我要订阅楼主 Medichen 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +9 wangmingqi 2026-08-28 10/500 2026-08-28 21:28 by Robin_BE
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +11 Kittylucky 2026-08-27 12/600 2026-08-28 16:26 by lyxiaomc
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 2026年叶企孙基金 +3 bud_bud 2026-08-27 6/300 2026-08-28 11:53 by bud_bud
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
请填处理意见