[求助]多自由度矩阵特征值？

minthy · 发表于 2006-1-26 17:40

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

x

多自由度振动系统特征值可以直接用eig[(K,M)来求吗 会不会误差太大？ 还劳烦高人指教

欧阳中华 · 发表于 2006-1-26 21:09

可以的，没有问题

minthy · 发表于 2006-1-27 09:49

谢啦！ 那这个函数是用哪种方法求的 我这个是要求固有频率和振型的，怕误差太大

欧阳中华 · 发表于 2006-1-27 11:51

这个你得查一下Matlab的说明，我一般是用FORTRAN自己编的，Matlab提供的求解器也就是那几种常规的方法，一般使用时，还应该设定收敛条件，也就是迭代到什么程度停止迭代，这就直接影响到你的计算精度。你再琢磨琢磨... ... [em05][em05][em05]

mxlzhenzhu · 发表于 2013-7-12 23:52

确实有这个问题，我用eig比较和nastran的结果，发现，第一阶可以差别很大，最大甚至30Hz;

还不理解eig（）函数；不过，我打算自己安装Rao的书上算法自己搞。

ChaChing · 发表于 2013-7-13 00:09

mxlzhenzhu 发表于 2013-7-12 23:52

登录/注册后可看大图

确实有这个问题，我用eig比较和nastran的结果，发现，第一阶可以差别很大，最大甚至30Hz;

还不理解eig（ ...

个人以为不太可能!
matlab及NASTRAN都属於成熟的软件, 别再怀疑不适用
一般年轻人跟我说这些, 我会反问那麼多使用者都是呆子吗
建议先检查下建模是否正确

mxlzhenzhu · 发表于 2013-7-13 16:46

ChaChing 发表于 2013-7-13 00:09

登录/注册后可看大图

个人以为不太可能!
matlab及NASTRAN都属於成熟的软件, 别再怀疑不适用
一般年轻人跟我说这些, 我会反问 ...

建模没有问题；其实他们两个的算法是不一样的；结论：Lanczos是一种近似，我更相信matlab。

ChaChing · 发表于 2013-7-13 22:05

mxlzhenzhu 发表于 2013-7-13 16:46

登录/注册后可看大图

建模没有问题；其实他们两个的算法是不一样的；结论：Lanczos是一种近似，我更相信mat ...

...第一阶可以差别很大，最大甚至30Hz...

个人水平有限, 没真正深究过两者算法差异及优劣, 仅是图个方便的普通小工使用者而已
不过即使两者算法不一样, 个人真的还是难以相信第一阶会差到30Hz!?

mxlzhenzhu · 发表于 2013-7-14 00:30

ChaChing 发表于 2013-7-13 22:05

登录/注册后可看大图

个人水平有限, 没真正深究过两者算法差异及优劣, 仅是图个方便的普通小工使用者而已
不过即使两者算法 ...

我做过很多试验了，确实是这样，这种事情发生过两次，且只有第一阶是差别这么大，别的都正常，很少有超过4Hz的误差。

账号		自动登录	找回密码
密码			我要加入