24小时热门版块排行榜    

查看: 3516  |  回复: 5

yaoshuiyun

新虫 (初入文坛)

[求助] 求问正交矩阵的维数以及表示 已有1人参与

弱弱的求助,一个n*n的正交矩阵的维数是n*(n-1)/2,那么怎么用n*(n-1)/2个变量表示一个n*n的矩阵呢?
现在知道n=3的情况,用(3-1)*3/2=3个数来表示,如附件所示,求n=4和n=5的情况,工科小硕万分感谢,拜托拜托。

求问正交矩阵的维数以及表示
捕获.PNG
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hank612

至尊木虫 (著名写手)

http://www.bio21.bas.bg/ibf/dpb_ ... /CD_IM_2012_AIP.pdf

http://stanford.edu/~rezab/nips2 ... ted/coordescent.pdf

http://www.wikiwand.com/en/Euler ... extrinsic_rotations

https://en.wikipedia.org/wiki/Rotation_matrix
截图从理论到实践都告诉我们如何将 SO(n)分解为n(n−1)/2 个Givens 矩阵的乘积,你自己搜搜关键词,看看对n=4,5 有没有实例操作。
求问正交矩阵的维数以及表示-1
Givens rotation compositon.png

We_must_know. We_will_know.
2楼2015-12-27 12:44:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

getengqing

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
如果:AA'=E(E为单位矩阵,A'表示“矩阵A的转置矩阵”。)或A′A=E,则n阶实矩阵A称为正交矩阵
这个矩阵中A与A'的积不是单位矩阵,显然不是正交矩阵。另外,可一眼看出A的行(或列)向量组都不是单位向量,所以A不是正交组阵。
1.方阵A正交的充要条件是A的行(列)向量组是单位正交向量组;

2.方阵A正交的充要条件是A的n个行(列)向量是n维向量空间的一组标准正交基;

3.A是正交矩阵的充要条件是:A的行向量组两两正交且都是单位向量;

4.A的列向量组也是正交单位向量组。

5.正交方阵是欧氏空间中标准正交基到标准正交基的过渡矩阵
一起交流学习/分享优秀资源
3楼2015-12-27 18:11:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yaoshuiyun

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by hank612 at 2015-12-27 12:44:05
http://www.bio21.bas.bg/ibf/dpb_files/im/mladenovpub/CD_IM_2012_AIP.pdf

http://stanford.edu/~rezab/nips2013workshop/accepted/coordescent.pdf

http://www.wikiwand.com/en/Euler_angles#/Euler_angle ...

大神,能不能帮看看n=4的结果,有没有可能显式地表示出来,还是没看懂==
4楼2015-12-27 20:05:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hank612

至尊木虫 (著名写手)

引用回帖:
4楼: Originally posted by yaoshuiyun at 2015-12-27 20:05:36
大神,能不能帮看看n=4的结果,有没有可能显式地表示出来,还是没看懂==...

http://www.cs.nthu.edu.tw/~cherung/teaching/2011anm/note07.pdf

这篇文章的第二,第三页非常清晰地给出了操作步骤。

基本上,一个n维单位向量由它与n个坐标轴的夹角唯一确定, 而作一次Givens rotation可以使向量旋转到(在指定的坐标轴上投影为零)。 所以当一个n维单位向量经过依次(n-1)个 Givens rotation, 垂直x2,...,xn 这(n-1)个坐标轴以后,自然等于(1,0,0...,0)。 经过这(n-1)个Givens矩阵作用后, O(n) 就转化成第一行为(1,000),(于是第一列也是(1,0000)),从第二行第二列开始是另一个O(n-1)矩阵。 然后如法炮制,用(n-2), (n-3),..,1 个Givens rotation 化简到单位矩阵。 这就是为什么至多  n(n-1)/2 个Givens 矩阵就一定能表示一般的SO(n)了。

你会操作n=3的情况, 其它n=4,5什么的是一模一样的操作,楼主动动手就知道了。
http://www.cs.rpi.edu/~flaherje/pdf/lin13.pdf
"Givens QR Factorization" 也可以供参考。

http://faculty.ucmerced.edu/mhya ... tures/lecture12.pdf
很经典的内容哦
We_must_know. We_will_know.
5楼2015-12-28 05:48:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiajiasnoopy

金虫 (小有名气)

引用回帖:
3楼: Originally posted by getengqing at 2015-12-27 18:11:54
如果:AA'=E(E为单位矩阵,A'表示“矩阵A的转置矩阵”。)或A′A=E,则n阶实矩阵A称为正交矩阵
这个矩阵中A与A'的积不是单位矩阵,显然不是正交矩阵。另外,可一眼看出A的行(或列)向量组都不是单位向量,所以A不 ...

$A$显然是正交矩阵,不知道你哪里看出不是的
6楼2015-12-28 14:57:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yaoshuiyun 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 14:15 by jCd0dEvKHShX
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:13 by jCd0dEvKHShX
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +5 ASdOkHsho7FD 2026-08-28 9/450 2026-08-29 14:03 by jCd0dEvKHShX
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:02 by jCd0dEvKHShX
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:51 by jCd0dEvKHShX
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +14 晴天加油 2026-08-26 15/750 2026-08-29 11:06 by biily
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见