声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3181|回复: 6

[编程技巧] 求解矩阵行列式的方程

[复制链接]
发表于 2008-7-9 17:33 | 显示全部楼层 |阅读模式

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

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

x
例如:矩阵如下:
A=[1-x  2  3; 4  5-x  6;7  8  9-x]; 其中x为未知量。
求解det(A)=0 ?
最好有命令,因为所处理的矩阵很大!
各位达人,路过就拔刀相助吧。:@)

[ 本帖最后由 ChaChing 于 2010-8-4 11:12 编辑 ]
回复
分享到:

使用道具 举报

发表于 2008-7-10 09:39 | 显示全部楼层
syms x

y=det(A);
转为多项式,然后再求
 楼主| 发表于 2008-7-10 17:53 | 显示全部楼层

疑惑

按3楼主任的方法做了之后所得的多项式如下:
x^9-66673120054560363/70368744177664*x^8+478235594689933356941139765688113/1237940039285380274899124224*x^7-7732838775237209051331383658906252999057465836775/87112285931760246646623899502532662132736*x^6+154841367116801208993346117297884763304560427032637411291003654703/12259964326927110866866776217202473468949912977468817408*x^5-991370917947471657696594091089401583212241133477063514104443509640105905113318845/862718293348820473429344482784628181556388621521298319395315527974912*x^4+2018220659752477454755431848782449482456285745422725432899954704980952483679023948410686140801517/30354201441027016733116592294117482916287606860189680019559568902170379456331382784*x^3-4999916680097234771047347051000043575998204769269400667849414416414869450798775578249266358281255799707220373651/2135987035920910082395021706169552114602704522356652769947041607822219725780640550022962086936576*x^2+27035406820626855544158110585561249944892437204337245033252956767912206945201984942841013446820015534257017928143443481740398121/601226901190101306339707032778070279008174732520529886901066488712245510429339761526706943586500787976175353856*x-14883951798095742968031210595923755991325385553888278458807761464439893464583530132918466118066273183748589034654065872762128907000183245734195/42307582002575910332922579714097346549017899709713998034217522897561970639123926132812109468141778230245837569601494931472384     
要想求这个多项式的根,必须把它的系数写成向量形式,可是这样的系数写出来难度太大,请问各位大侠有何高见?
发表于 2008-7-10 19:45 | 显示全部楼层

回复 4楼 的帖子

没看出来有啥难度呀:@) :@)
 楼主| 发表于 2008-7-10 20:24 | 显示全部楼层

问题已经基本解决

代码如下:
syms x
y=det(A);
[x]=solve(y,x);
vpa(x)
vpa(x,5)

评分

1

查看全部评分

发表于 2008-7-11 21:36 | 显示全部楼层

回复 6楼 的帖子

后面两句主要是改变现实的精度吧?
发表于 2010-8-4 09:37 | 显示全部楼层
学习了啊啊啊啊
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-3 01:33 , Processed in 0.083281 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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