rushuidetian 发表于 2008-7-9 17:33

求解矩阵行列式的方程

例如：矩阵如下：
A=; 其中x为未知量。
求解det(A)=0 ？
最好有命令，因为所处理的矩阵很大！
各位达人，路过就拔刀相助吧。:@)

[ 本帖最后由 ChaChing 于 2010-8-4 11:12 编辑 ]

sigma665 发表于 2008-7-10 09:39

syms x

y=det(A);
转为多项式，然后再求

rushuidetian 发表于 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   
要想求这个多项式的根，必须把它的系数写成向量形式，可是这样的系数写出来难度太大，请问各位大侠有何高见？

messenger 发表于 2008-7-10 19:45

回复 4楼 的帖子

没看出来有啥难度呀:@) :@)

rushuidetian 发表于 2008-7-10 20:24

问题已经基本解决

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

无水1324 发表于 2008-7-11 21:36

回复 6楼 的帖子

后面两句主要是改变现实的精度吧？

quanwang2003 发表于 2010-8-4 09:37

学习了啊啊啊啊
页: [1]
查看完整版本: 求解矩阵行列式的方程