请问有阻尼转子系统求临界转速的编程

qianen

滑动轴承支撑的转子系统，考虑转子阻尼时求临界转速的matlab程序，通过传递矩阵求得的非线性方程为

.2662237332e-26*S^28+.2853707525e-23*S^27-.8478044388e-18*S^26-.9096790787e-15*S^25+.8994810637e-10*S^24+.9670022289e-7*S^23-.3174314383e-2*S^22-3.432281481*S^21-533.4376380*S^20+466689.3862*S^19+126737547.1*S^18-.1951799733e11*S^17-.7149478858e13*S^16+.2200636699e15*S^15+.1497645713e18*S^14+.1615371056e19*S^13-.1198032955e22*S^12-.3476596232e23*S^11+.2490798252e25*S^10+.9300526908e26*S^9-.1163039727e28*S^8-.6964934962e29*S^7-.5253088727e30*S^6+.3578229567e31*S^5+.3256779549e32*S^4-.3100845164e32*S^3-.3105418599e33*S^2+.3677934032e32*S+.3773107290e33