24小时热门版块排行榜    

查看: 2237  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 4/200 2026-11-02 03:50 by xzcv0Tehd8lQ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-02 03:18 by xzcv0Tehd8lQ
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 xvJS1BYDxvBS 2026-10-31 11/550 2026-11-02 03:13 by xzcv0Tehd8lQ
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:05 by xzcv0Tehd8lQ
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +7 X5bHxOULxFqb 2026-10-30 24/1200 2026-11-02 02:39 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-02 02:34 by xzcv0Tehd8lQ
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:05 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:45 by BCdvBAEpt3zm
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +9 zuesli 2026-10-27 15/750 2026-11-01 12:59 by Equinoxhua
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 3/150 2026-11-01 12:53 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 4/200 2026-11-01 11:01 by BCdvBAEpt3zm
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 X5bHxOULxFqb 2026-10-30 26/1300 2026-11-01 10:22 by BCdvBAEpt3zm
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
信息提示
请填处理意见