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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab求解正交矩阵。
111/2返回列表12下一页
查看: 4364  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

virtualzx

木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
引用回帖:
5楼: Originally posted by luzongxing at 2013-12-26 20:27:50
由于R是未知，把R中的变量当作未知数，构造新的线性方程组，就是
AX=b的模式。

现在我的问题是，我构造出AX=b的线性方程组了。A，b是一直，哪一种求解方法比较精确。

求出旋转矩阵R以后，我们发现R*R'跟单位 ...

是的，修正一个接近特殊正交的矩阵，使其变为特殊正交矩阵。
赞一下(1人) 回复此楼
7楼2013-12-27 10:03:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

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的回帖

luzongxing

银虫 (小有名气)
引用回帖:
4楼: Originally posted by virtualzx at 2013-12-27 01:45:16
嗯？奇怪呢，一般都是x.R=y求R，x和y都是3维矢量，或者A.R=B求R，A和B的所有特征值相等，没遇到过你说的9x9的情况。

既然是旋转变换，不单是正交矩阵，而且需要是特殊正交矩阵。特殊正交矩阵可以要求矩阵对数反 ...

由于R是未知，把R中的变量当作未知数，构造新的线性方程组，就是
AX=b的模式。

现在我的问题是，我构造出AX=b的线性方程组了。A，b是一直，哪一种求解方法比较精确。

求出旋转矩阵R以后，我们发现R*R'跟单位矩阵会有一个很小的偏差。怎么修正使得R*R'跟接近于单位矩阵。。
你的
Q=logm(R)
Qsym=1/2(Q-Q.')
Rsym=expm(Qsym)
是针对矩阵的微小修正方法吗？
谢谢！
赞一下 回复此楼
5楼2013-12-27 09:27:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzongxing

银虫 (小有名气)
引用回帖:
4楼: Originally posted by virtualzx at 2013-12-27 01:45:16
嗯？奇怪呢，一般都是x.R=y求R，x和y都是3维矢量，或者A.R=B求R，A和B的所有特征值相等，没遇到过你说的9x9的情况。

既然是旋转变换，不单是正交矩阵，而且需要是特殊正交矩阵。特殊正交矩阵可以要求矩阵对数反 ...

现在我的问题是，我构造出AX=b的线性方程组了。A{9*9}，b{1*9}是已知，哪一种求解方法比较精确。

上面打错字了。。
赞一下 回复此楼
6楼2013-12-27 09:29:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzongxing

银虫 (小有名气)
引用回帖:
7楼: Originally posted by virtualzx at 2013-12-27 10:03:55
是的，修正一个接近特殊正交的矩阵，使其变为特殊正交矩阵。...

就是这个。谢谢你
请问这个的数学方法的理论计算过程在哪边可以找到。
方便的话可以留个邮箱吗？

还有这个分数怎么给你。
赞一下 回复此楼
8楼2013-12-27 10:22:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
luzongxing: 金币+5, ★★★★★最佳答案 2013-12-29 21:07:45
引用回帖:
8楼: Originally posted by luzongxing at 2013-12-26 21:22:59
就是这个。谢谢你
请问这个的数学方法的理论计算过程在哪边可以找到。
方便的话可以留个邮箱吗？

还有这个分数怎么给你。...

这个是矩阵指数的一个性质，你可以查看一下矩阵指数的维基百科页面，性质那一栏里的倒数第二项（斜对称矩阵是反对称矩阵的别称）：

  • exp(X^T) = (exp X)^T，其中X^T表示X的转置。从中可以推出，如果X是对称矩阵，则e^X也是对称矩阵；如果X是斜对称矩阵，则e^X是正交矩阵。


另外如果你对群论比较熟悉的话，坐标旋转构成一个特殊正交群，而特殊正交群的李代数是反对称矩阵(就是矩阵的转置是矩阵本身的相反数)。而实矩阵的李代数是指数映射携带的。
赞一下 回复此楼
9楼2013-12-27 14:48:47
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 luzongxing 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 3/150 2026-03-15 17:32 by 小物理化学
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[基金申请] 现在如何回避去年的某一个专家，不知道名字 +3 zk200107 2026-03-12 6/300 2026-03-14 17:13 by zk200107
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 301求调剂 +3 归零lbm 2026-03-09 3/150 2026-03-14 02:20 by JourneyLucky
[考研] 材料与化工 一志愿山大 321分 求调剂 +7 每天散步 2026-03-09 8/400 2026-03-14 02:18 by JourneyLucky
[考研] 295复试调剂 +5 简木ChuFront 2026-03-09 5/250 2026-03-14 01:29 by JourneyLucky
[考研] 求调剂，一志愿江南大学环境工程085701 +3 Djdjj12 2026-03-10 4/200 2026-03-14 00:31 by JourneyLucky
[考研] b区环境工程求调剂 +4 Maps1 2026-03-10 6/300 2026-03-14 00:23 by JourneyLucky
[考研] 279求调剂 +3 抓着星星的女孩 2026-03-10 3/150 2026-03-13 23:47 by userper
[考研] 材料工程,326分，求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 0856材料与化工301求调剂 +5 奕束光 2026-03-13 5/250 2026-03-13 22:00 by 星空星月
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考研] 化工0817调剂 +8 灿若星晨 2026-03-10 8/400 2026-03-10 22:44 by 星空星月
信息提示
关闭
请填处理意见
关闭