24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ndadams

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 50.7
散金: 9
帖子: 106
在线: 59.3小时
虫号: 208064
注册: 2006-03-05
专业: 机械动力学

[求助] 请教：mathematica 每次的计算结果不同已有1人参与

Tee = Chop[TLee, 10*10^-8];
H = Det[Tee];
fpz = \[Omega] /. Solve[H == 0, \[Omega]]

请教各位问题：
求矩阵Tee的行列式，然后通过行列式的值为零求频率Omega。
在求行列式之前用了Chop命令 Chop[TLee, 10*10^-8];
当Chop的第二个参数（10*10^-8）取不同值时，计算结果Omega差别很大，这个是什么原因，怎么解决呢？

谢谢！

回复此楼

» 猜你喜欢

为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？已经有9人回复
版面费该交吗已经有9人回复
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低已经有13人回复
面上可以超过30页吧？已经有4人回复
“人文社科而论，许多学术研究还没有达到民国时期的水平” 已经有5人回复
什么是人一生最重要的？已经有4人回复

1楼 2015-08-24 17:20:24

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

ndadams

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 50.7
散金: 9
帖子: 106
在线: 59.3小时
虫号: 208064
注册: 2006-03-05
专业: 机械动力学

我明白了，先前的想法是不对的。原因如下：

在这里不能使用Chop指令，因为矩阵中的元素是\[Omega]\的函数，并不知道矩阵中每个元素的数值大小，Chop截断的是接近于零常数和常系数，因此很多包含\[Omega]的项因为系数较小被截断而去掉了。
导致行列式中\[Omega]的最高阶次严重降低，很多高次项被去掉了。这样，得到的多项式就和原来的大不相同，得到的解也就不同。