24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1549  |  回复: 9

tang_zm

新虫 (小有名气)

[求助] 矩阵特征值计算已有1人参与

matlab有计算矩阵特征值的功能,然后数学中好像有个 LAPACK library 也可以算矩阵特征值,请教下两者有何区别,哪个更有优势些?我想操作符号矩阵,但是matlab在矩阵超过6阶后,符号矩阵就不能给出显示解了(给出一个Sizeoff函数).请问有何办法得到特征值的显示解析解???谢谢
回复此楼

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

wurongjun

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
tang_zm: 金币+10, ★★★很有帮助 2014-12-19 15:34:23
tang_zm: 回帖置顶 2014-12-19 15:48:04
引用回帖:
4楼: Originally posted by tang_zm at 2014-12-19 10:29:12
您的意思是两者均不能给出解析解吗,迄今为止,也没有能够解析搞定6阶以上符号矩阵特征值的是吧?...

应该说是:不能给出这类问题的公式解!
数值解是可以的!
3#说得对!四次以上就没有公式解了!
善恶到头终有报,人间正道是沧桑.
5楼2014-12-19 10:35:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

wurongjun

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
两个软件都可以算数值解,但是规模较大时后者更有优势!规模超大的本身是学科里面的难题!没有完全解决!
对于符号真的解析解,6阶后的一般是不能给出的!这等价于给出6次以上代数方程的解析解!
善恶到头终有报,人间正道是沧桑.
2楼2014-12-19 10:13:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

引用回帖:
2楼: Originally posted by wurongjun at 2014-12-19 10:13:30
两个软件都可以算数值解,但是规模较大时后者更有优势!规模超大的本身是学科里面的难题!没有完全解决!
对于符号真的解析解,6阶后的一般是不能给出的!这等价于给出6次以上代数方程的解析解!

好像4阶以上代数方程就没有公式解了
3楼2014-12-19 10:19:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tang_zm

新虫 (小有名气)

引用回帖:
2楼: Originally posted by wurongjun at 2014-12-19 10:13:30
两个软件都可以算数值解,但是规模较大时后者更有优势!规模超大的本身是学科里面的难题!没有完全解决!
对于符号真的解析解,6阶后的一般是不能给出的!这等价于给出6次以上代数方程的解析解!

您的意思是两者均不能给出解析解吗,迄今为止,也没有能够解析搞定6阶以上符号矩阵特征值的是吧?
4楼2014-12-19 10:29:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tang_zm

新虫 (小有名气)

引用回帖:
5楼: Originally posted by wurongjun at 2014-12-19 10:35:04
应该说是:不能给出这类问题的公式解!
数值解是可以的!
3#说得对!四次以上就没有公式解了!...

继续求助。。。,如想得到一个高阶符号矩阵的特征值特性,由于上述达人所言,无法获得特征值解析解,而只能是数值解。那么,从数值解出发,如何能很好地还原其本身的特征值特性呢?
6楼2014-12-19 15:48:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zaq123321

专家顾问 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
6楼: Originally posted by tang_zm at 2014-12-19 15:48:00
继续求助。。。,如想得到一个高阶符号矩阵的特征值特性,由于上述达人所言,无法获得特征值解析解,而只能是数值解。那么,从数值解出发,如何能很好地还原其本身的特征值特性呢?...

It had been proved by Galois theory that the general quintic equation is not solvable by radicals. So nobody dreams to solve the insolvable problem you asked since you ask for general symbolic solution of a sixth order equation unless you can prove Galois theory is wrong. Mathematicians have already wasted hundreds of years time before Galois. On numerical solution, your input is not numerical, how to get numerical solution? Consider whether your question makes sense or not in next stage.

[ 发自手机版 http://muchong.com/3g ]

» 本帖已获得的红花(最新10朵)

小木虫给我温暖,给我希望,爱就要爱小木虫。
7楼2014-12-19 16:19:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tang_zm

新虫 (小有名气)

送红花一朵
引用回帖:
7楼: Originally posted by zaq123321 at 2014-12-19 16:19:10
It had been proved by Galois theory that the general quintic equation is not solvable by radicals. So nobody dreams to solve the insolvable problem you asked since you ask for general symbolic solut ...

OK,thanks! I agree with you to some extent.
8楼2014-12-19 21:29:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tang_zm

新虫 (小有名气)

引用回帖:
7楼: Originally posted by zaq123321 at 2014-12-19 16:19:10
It had been proved by Galois theory that the general quintic equation is not solvable by radicals. So nobody dreams to solve the insolvable problem you asked since you ask for general symbolic solut ...

可能是我没表述清楚,现详述一番: 比如,一个区域内的物质特性均满足某个方程,而这个方程其中包含了多个参数。我现在将该区域离散为N个子区域,且假定在每个子区域中物质均满足原来的方程,这样我就可以得到N个线性方程组。
      这N个线性方程组的系数矩阵,理论上就该和原来方程所含的参数以及N有关,因此,其为一符号矩阵。如各位所言,对该符号矩阵求特征值的是不可能的是吧?
       但是将各个参数具体化后,就可以计算出系数矩阵的特征值了。这就是我所谓的“numerical solution”,这显然只是符号矩阵特征值的一个特例。那么,从这样的数值解出发,如何能很好地还原方程其本身的具体普遍意义的特征值特性呢(也就是符号矩阵的特征值特性)?希望继续指导。。。
9楼2014-12-19 21:49:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zaq123321

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
tang_zm: 金币+15, ★★★很有帮助, Thanks a lot! 2014-12-20 09:54:38
Simplify your problem as follows. Matrix A has a bunch of parameters in its entry and A is of 6th order. When you choose some specific values for the parameters, you can solve an eigensystem by numerical way ( hope so). Then next you expect to get some relationship  from the eigenvalues to the eigenvalue of the parametrized matrix. I think the first step you may want to go is to figure out the entries of matrix A. However, this is not possible in the current stage. The research topic is called eigensystem inverse problem. Basically it wants to reconstruct the matrix from its eigensystem. The next step is to solve parameters from matrix entries, which is usually nonlinear system.

[ 发自手机版 http://muchong.com/3g ]
小木虫给我温暖,给我希望,爱就要爱小木虫。
10楼2014-12-19 23:40:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tang_zm 的主题更新
信息提示
请填处理意见