声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1110|回复: 5

[综合讨论] 求解det([h-I]2*2)=0

[复制链接]
发表于 2010-4-12 17:23 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
已知:
h=[28.61031789247    -28.89049270091    3599.14530104906;
   -14.77850691994     14.92324472504   -1859.13590304261;
   -0.02616072627      -0.02758480103     56.52067419820];
I=[f 0 0;0 f 0;0 0 f];
并且:det([h-I]2*2)=0 其中[h-I]2*2表示矩阵[h-I]的任意一个2*2阶子矩阵
问:怎么求解f(matlab代码)

[ 本帖最后由 ChaChing 于 2010-4-19 14:38 编辑 ]
回复
分享到:

使用道具 举报

发表于 2010-4-13 16:56 | 显示全部楼层
LZ想要的结果是一个f值么?还是每2*2矩阵算出一个f值?每2*2矩阵得到一个f值如下所示:
clear;clc
syms f
h=[28.61031789247    -28.89049270091    3599.14530104906;
   -14.77850691994     14.92324472504   -1859.13590304261;
   -0.02616072627      -0.02758480103     56.52067419820];

l=[f 0 0;0 f 0;0 0 f]; a=h-l;
a1=a(1:2,1:2); a2=a(1:2,2:3);
a3=a(2:3,1:2); a4=a(2:3,2:3);
f1=solve(det(a1)); f2=solve(det(a2)); f3=solve(det(a3)); f4=solve(det(a4));

不知是不是LZ想要的。。

[ 本帖最后由 ChaChing 于 2010-4-13 20:58 编辑 ]

评分

1

查看全部评分

发表于 2010-4-13 20:57 | 显示全部楼层
个人水平专业有限, 真的看不清楼主说的
若如2F所说, 尚缺a(1:2,[1,3]),a(2:3,[1,3]),a([1,3],1:2),a([1,3],2:3),a([1,3],[1,3])
但明显无一f值可以满足所有任意一个2*2阶子矩阵det=0!
 楼主| 发表于 2010-4-14 16:44 | 显示全部楼层

回复 沙发 maigicku 的帖子

你的方法求解的f是单个的,每次都不相等。这个问题我已经解决,大家一块分享:
思路:det([h-I]2*2)=0 其中[h-I]2*2表示矩阵[h-I]的任意一个2*2阶子矩阵,其中[h-I]2*2有9个2*2阶子矩阵;6个为f的一次函数,3个为f的二次函数(这里不用):
A=[h(1,3); h(2,3); h(3,1); h(2,1); h(3,2); h(1,2)];
B=[h(1,3)*h(2,2)-h(1,2)*h(2,3);  h(1,1)*h(2,3)-h(1,3)*h(2,1);
     h(2,2)*h(3,1)-h(2,1)*h(3,2);   h(2,1)*h(3,3)-h(2,3)*h(3,1);
     h(1,1)*h(3,2)-h(1,2)*h(3,1);   h(1,2)*h(3,3)-h(3,2)*h(1,3)];
f=A\B;

[ 本帖最后由 ChaChing 于 2010-4-19 13:10 编辑 ]

评分

1

查看全部评分

发表于 2010-4-15 00:16 | 显示全部楼层
看完LZ的解法, 有点了解LZ要什么了!
看来我跟maigicku都误解LZ在1F, "det([h-I]2*2)=0 其中[h-I]2*2表示矩阵[h-I]的任意一个2*2阶子矩阵"的含义了!?
 楼主| 发表于 2010-4-19 09:44 | 显示全部楼层

回复 5楼 ChaChing 的帖子

不好意思,可能是我的问题,帖子写的很草,很让人费解,见谅!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 13:29 , Processed in 0.069882 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表