24小时热门版块排行榜    

查看: 5765  |  回复: 5

aagiler

金虫 (小有名气)

[求助] 在matlab中如何提高矩阵求逆的运算精度? 已有4人参与

在matlab下,如何提高矩阵求逆的运算精度。目前用inv求逆得到的结果与理论值还有些差别,不知道有何方法可以提高运算精度。以6阶矩阵为例,要求矩阵Q的逆,该矩阵具有如下结构:
Q=[C;
     C*A;
     C*A*A;
     C*A*A*A;
     C*A*A*A*A;
     C*A*A*A*A*A]
其中
C=[1,0,1,0,1,0];
a1=2*pi*0.001;
a2=4*pi*0.001;
a3=6*pi*0.001;
A=[cos(a1) -sin(a1) 0 0 0 0;
   sin(a1)  cos(a1) 0 0 0 0;
   0 0 cos(a2) -sin(a2) 0 0;
   0 0 sin(a2)  cos(a2) 0 0;
   0 0 0 0 cos(a3) -sin(a3);
   0 0 0 0 sin(a3)  cos(a3)];
明显,A^-1的特征值均在单位圆上
D=[1;0;0;0;0;0;]
P=[Q^-1*D A*Q^-1*D A*A*Q^-1*D A*A*A*Q^-1*D A*A*A*A*Q^-1*D A*A*A*A*A*Q^-1*D]
-1表示矩阵的逆
但目前P^-1*A^-1*P的特征值与1还有些差别。
谢谢先
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

RT_Fem

铜虫 (正式写手)

【答案】应助回帖

根据矩阵相乘满足结合律,那么矩阵Q应该就是=[C;C;C;C;C;C],从而不可求逆的吧。
No sacrifice, no victory.
2楼2014-07-07 21:26:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็

【答案】应助回帖

精度要求高,而且问题比较简单,当然是用符号计算了;
符号计算mathematica更加擅长一些; matlab凑合能用;
ส็็็็็็็็็็็็็็็็็็็็
3楼2014-07-07 21:56:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ayismas

木虫 (正式写手)

【答案】应助回帖

用矩阵的分解方法,应该求的比较快,特别是对于矩阵的行列数比较大的时候非常有效。

[ 发自手机版 http://muchong.com/3g ]
4楼2014-07-07 22:36:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mathstudy

金虫 (正式写手)

【答案】应助回帖

1. 利用左除运算
2.根据你矩阵是否具有特殊的结构,尝试矩阵分解
3.从网上查一查matlab的高精度工具箱试试
5楼2014-07-08 09:51:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hy东舟

银虫 (初入文坛)

有限精度下,矩阵求逆应与矩阵条件数密切相关,如果是一个稍病态的矩阵,不能奢望能高精度计算矩阵的逆
6楼2014-07-15 11:01:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 aagiler 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化学调剂 +5 yzysaa 2026-03-21 5/250 2026-03-21 22:12 by peike
[考研] 材料求调剂 +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 22 350 本科985求调剂,求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 材料学学硕080502 337求调剂-一志愿华中科技大学 +4 顺顺顺mr 2026-03-18 5/250 2026-03-21 10:22 by luoyongfeng
[考研] 求调剂 +6 Mqqqqqq 2026-03-19 6/300 2026-03-21 08:04 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-16 3/150 2026-03-21 03:59 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 考研调剂求学校推荐 +3 伯乐29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 中南大学化学学硕337求调剂 +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[考研] 289求调剂 +6 怀瑾握瑜l 2026-03-20 6/300 2026-03-20 20:30 by 学员8dgXkO
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 0854,计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 283求调剂 +3 听风就是雨; 2026-03-16 3/150 2026-03-17 07:41 by 热情沙漠
信息提示
请填处理意见