24小时热门版块排行榜    

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

Medichen

铜虫 (小有名气)

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

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

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 估计是周四 +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
[基金申请] 2026国自然函评费到账 +17 羊腰板 2026-08-21 18/900 2026-08-24 12:03 by iaeyuan
[基金申请] 建议基金发布提前给出明确的时间点 +12 kulium 2026-08-21 15/750 2026-08-24 12:02 by iaeyuan
[基金申请] 范进中举一文的中心思想 +6 炎黄贵胄 2026-08-22 7/350 2026-08-24 11:58 by 6543yes
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧! +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
[基金申请] 明天应该可查了!? +4 chengyan1220 2026-08-23 4/200 2026-08-24 08:47 by gloomy6159
[基金申请] 什么时候开奖? +10 CrisMessi 2026-08-18 11/550 2026-08-24 06:50 by 开心的小狮子
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +4 process2012 2026-08-23 5/250 2026-08-23 19:58 by jurkat.1640
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见