anran2000_1984 发表于 2009-3-17 16:05

提取多项式系数时出现nan,有没有解决的办法??

我用整体传递矩阵法求一双转子系统的涡动频率,用fortran编程,采用二分法搜根得到正反各16阶涡动频率,与已有结果一致。但是用matlab得到频率方程后,用solve求解则正反涡动频率各有两阶出现复根,其它实根仍然一致。频率方程是一个关于涡动频率(ω)的矩阵的行列式,实质上最后形成一个关于ω的多项式。我想验证一下哪种解法是准确的,所以想在matlab中形成频率方程之后,再用fortran求解,这样就需要提取这个多项式的系数。用sym2poly()提取后发现只有0~7次幂的系数有数值,其它57个都是"nan",有没有什么办法能够让我得到这组系数?? 请达人指点啊~~
页: [1]
查看完整版本: 提取多项式系数时出现nan,有没有解决的办法??