24小时热门版块排行榜    

CyRhmU.jpeg
查看: 929  |  回复: 6

sg47

铜虫 (小有名气)

[求助] 一个关于positive definite matrix的问题

今看到一个theorem为
If matrix A (n*n) is positive semidefinite, then B'AB is also positive semidefinite independent of the rank of B (n*m).
本人尝试:
A=diag([0,0,1,1])
B=[1,0;0,1;2,0;0,2];
B'*A*B=  [4     0
                 0     4]
which is positive defnite!!!而不是个positive seimidefinite的matrix啊
求解!!!
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

hyit_lxq

木虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
soliton923: 回帖置顶 2011-12-16 22:16:42
soliton923(金币+2): 谢谢参与讨论~~~ 2011-12-16 22:16:54
sg47(金币+2): 2011-12-17 10:49:23
Proof. Denote y=Bx for any given vector x of dimensional m. The fact that A is positive semidefinite shows x'(B'AB)x=y'Ay>=0. Therefore, B'AB is positive semidefinite.
~ ~ ~
5楼2011-12-16 20:40:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

pengyehui

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
soliton923: 谢谢专家的参与 2011-12-16 22:16:31
不矛盾!
从定义看 正定  x'Ax>0  for x is not 0
   半正定 x'Ax>=0  for all x
2楼2011-12-16 18:40:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sg47

铜虫 (小有名气)

半正定 x'Ax>=0  也是x is not 0吧
3楼2011-12-16 18:59:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sg47

铜虫 (小有名气)

至少这本书上面的定义是
半正定 x'Ax>=0  for x is not 0
4楼2011-12-16 19:03:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

saladin983

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
楼上的证明足以说明问题,正定矩阵自然半正定,半正定的定义里用的是大于或者等于号,等于号不必要成立。
6楼2011-12-17 00:27:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

math2000

铁杆木虫 (职业作家)

矩阵B不要取满秩矩阵(行或列满秩)即可!得到的B‘AB就是半正定了
7楼2011-12-17 20:36:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sg47 的主题更新
信息提示
请填处理意见