应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab求解正交矩阵。
51/1返回列表
查看: 4734  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

luzongxing

银虫 (小有名气)

[求助] matlab求解正交矩阵。 已有1人参与

1，我们求解变换矩阵的时候得到一个AX=b的矩阵，
现在如何通过matlab得到3*3的旋转矩阵R比较精确地解，一下是我从网上查找的一些常用的解法。
%一是用伪逆法求解
% X=pinv(A)*b

% % 二是用左除法求解，
X=A\b

% %  三是用最小二乘法求解,
% X=lsqnonneg(A,b)

%四LU分解法
% [L,U]=lu (A);
% X=U\(L\b)

2.求得的矩阵R，理论上R应该是一个标准的单位正交矩阵，如果给点的A和B的对应关系不是那么精确，可能求出来的R跟正交矩阵有点偏差。
R*R'=[0.9985 -0.0085 0.0013;
           0.0045 1.0032  -0.0018;
           0.0032 0.0063  1.0024]

那么如何修正这个旋转矩阵，让R*R'的约接近于单位正交矩阵。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-12-26 13:18:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzongxing

银虫 (小有名气)
引用回帖:
9楼: Originally posted by virtualzx at 2013-12-27 14:48:47
这个是矩阵指数的一个性质，你可以查看一下矩阵指数的维基百科页面，性质那一栏里的倒数第二项（斜对称矩阵是反对称矩阵的别称）：

exp(X^T) = (exp X)^T，其中X^T表示X的转置。从中可以推出，如果X是对称矩阵， ...

我说说的的理解你看对不对。
Q=logm(R)去矩阵的对数，我还不知道有啥用。
A^T - A是斜对称矩阵。
所以Qsym=1/2*(Q-Q.')是构建斜对称矩阵，但是为什么要乘以1/2,还有那个逆前面为什么要加一个点？
如果X是斜对称矩阵，则e^X是正交矩阵。
Rsym=expm(Qsym)是吧斜对称矩阵转换为正交矩阵？

所以我到最后我还是不知道这样转换的数学依据。
就像推导公式一样。
赞一下 回复此楼
10楼2013-12-27 17:56:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

virtualzx

木虫 (著名写手)
请问你的R矩阵是从什么地方出来的，和AX=B有啥联系啊？
赞一下 回复此楼
2楼2013-12-26 14:46:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzongxing

银虫 (小有名气)
引用回帖:
2楼: Originally posted by virtualzx at 2013-12-26 14:46:14
请问你的R矩阵是从什么地方出来的，和AX=B有啥联系啊？

旋转矩阵R是点不同坐标系下的旋转变换矩阵。是3*3的矩阵、

第一问你就可以假设A是一个9*9的满秩矩阵
b是一个1*9的矩阵。来求矩阵X。
x是一个1*9的矩阵。求出来是R中的各个元素。
怎么有更高精度的算法。

第二问，由于在两个坐标系中的点对应不是那么精确，可能导致旋转矩阵算出来不是标准的单位正交矩阵，这时候需要微小的修正。
赞一下 回复此楼
3楼2013-12-26 17:28:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)
xzhdty: 2013-12-27 06:18:54
引用回帖:
3楼: Originally posted by luzongxing at 2013-12-26 04:28:23
旋转矩阵R是点不同坐标系下的旋转变换矩阵。是3*3的矩阵、

第一问你就可以假设A是一个9*9的满秩矩阵
b是一个1*9的矩阵。来求矩阵X。
x是一个1*9的矩阵。求出来是R中的各个元素。
怎么有更高精度的算法。

...

嗯？奇怪呢，一般都是x.R=y求R，x和y都是3维矢量，或者A.R=B求R，A和B的所有特征值相等，没遇到过你说的9x9的情况。

既然是旋转变换，不单是正交矩阵，而且需要是特殊正交矩阵。特殊正交矩阵可以要求矩阵对数反对称得到。就是先求矩阵自然对数，然后反对称化，再求矩阵指数就行了
Q=logm(R)
Qsym=1/2(Q-Q.')
Rsym=expm(Qsym)
赞一下 回复此楼
4楼2013-12-27 01:45:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +12 医学老男孩 2026-08-13 32/1600 2026-08-14 08:19 by 医学老男孩
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 30/1500 2026-08-14 07:55 by 医学老男孩
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +3 启萌科技 2026-08-12 11/550 2026-08-13 22:57 by 启萌科技
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +4 majunge000 2026-08-11 6/300 2026-08-13 16:50 by 四季常青藤
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +22 Lanmanbaby 2026-08-09 36/1800 2026-08-12 22:50 by sdfapple719
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +4 工厂打螺丝 2026-08-08 4/200 2026-08-12 09:18 by hansi2025
[基金申请] 是这周出结果还是下周出结果？ +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
信息提示
关闭
请填处理意见
关闭