wangfan1618 发表于 2007-12-30 12:23

急求:matlab中的方程求解!

在matlab中用solve求解A方程,得到B的结果,可是再把B的结果代回方程A时却不等于0
A
.221128e31-.328079e29*w+.116683e27*w^2-.245252e23*w^3+.137878e19*w^4-.277110e14*w^5+.180525e9*w^6=0
B =
[ 107.801, 187.265, 6865.04, 24601.2, 50266.2, 71474.7]
比如说把107.801代入A时,结果是A=-3.2218e+024
请高手们指点一下,如果真要求解,又该用什么命令?谢谢!

16443 发表于 2007-12-30 13:21

回复 #1 wangfan1618 的帖子

A=-3.2218e+024这么小的数了不就是零了么 。你限制一下数据有效位数不就是零了。

finalwsl 发表于 2007-12-30 15:42

回复 #2 16443 的帖子

Matlab得出的解是数值解哦.

wangfan1618 发表于 2007-12-30 15:53

回复 #2 16443 的帖子

数值不小哦,那应该表示-3.2218*10^24吧

无水1324 发表于 2007-12-30 18:46

回复 #1 wangfan1618 的帖子

这个方程比较有意思,
1、通过ezplot确实可以找到解,但是只要有微弱的差别,即使是10e-8或者更小的误差也可能使解产生很大的误差;
2、solve可能不适合于求解这种问题;
3、楼主应该将系数化简一下。

zhuofeng 发表于 2007-12-31 13:06

回复 #5 无水1324 的帖子

传说中的刚性方程?需要特殊的求解方法?

wangfan1618 发表于 2008-1-2 11:28

原帖由 无水1324 于 2007-12-30 18:46 发表 http://www.chinavib.com/forum/images/common/back.gif
这个方程比较有意思,
1、通过ezplot确实可以找到解,但是只要有微弱的差别,即使是10e-8或者更小的误差也可能使解产生很大的误差;
2、solve可能不适合于求解这种问题;
3、楼主应该将系数化简一下。
我将系数化简了,结果还是误差很大,请问用什么软件可以求到比较精确的解呢?

[ 本帖最后由 wangfan1618 于 2008-1-2 11:29 编辑 ]

messenger 发表于 2008-1-2 14:35

需要特殊的求解方法

107.80062091532734619660729972912175038721089989631, 187.26496876934531373649978739934093664093217569692, 6865.0395482175273473450129075368505200177344199174, 24601.230169777239324737935704419066288174096383437, 50266.205340974286330516204339086390930328835840085, 71474.744353423548953171936056552952050562532200369.

sigma665 发表于 2008-1-2 16:23

回复 #8 messenger 的帖子

特殊的怎么处理?
页: [1]
查看完整版本: 急求:matlab中的方程求解!