1

作为作业的一部分，我必须推导出汽车悬架系统的运动方程。本质上，这是一个弹簧质量阻尼器问题。尚未给出汽车质量 M1、车轮质量 M2、弹簧常数 k1 和 k2 以及阻尼常数 c 的值。我已经导出了运动方程，并导出了一个将路面（输入）与最终的车身位移（输出）相关的传递函数。我必须确定这个传递函数的极点，因此我需要找到特征方程（分母）的根。问题是我没有上述变量的任何值，我试图在 MATLAB 中以符号方式分解我的四阶多项式或立即计算根。我不能假设任何值，它必须象征性地解决，

我对 MATLAB 没有太多经验，所以我不知道它的所有功能。

我试图解决的特征方程是：

(M1*M2)*s^4 + c*(M1+M2)*s^3 + ((M1*k1)+(M1*k2)+c^2+(M2*k2)-c)*s^2 + k1*c*s + ((k1*k2)-(k2^2))

先感谢您。

4

2 回答 2

2

您的方程式中有一些错误；

 c(M1+M2)*s^3 -> c*(M1+M2)*s^3
 + +k1*c*s -> + k1*c*s

但是如果你想解决多元方程，你可以这样做；

syms M1 M2 c k1 k2 s
eqn = (your equation) == 0;
roots = solve(eqn, s);

更多信息在这里：solve

于 2014-11-03T13:28:49.963 回答
0

现在您只需要按照以下步骤操作，以防您只想计算方程的根，这与前面的注释类似：

1. syms c s
2. roots=solve((M1*M2)*s^4 + c*(M1+M2)*s^3 + ((M1*k1)+(M1*k2)+c^2+(M2*k2)-c)*s^2 + k1*c*s + ((k1*k2)-(k2^2)),s)
or
roots=solve((M1*M2)*s^4 + c*(M1+M2)*s^3 + ((M1*k1)+(M1*k2)+c^2+(M2*k2)-c)*s^2 + k1*c*s +((k1*k2)-(k2^2)),c)
or 
roots=solve((M1*M2)*s^4 + c*(M1+M2)*s^3 + ((M1*k1)+(M1*k2)+c^2+(M2*k2)-c)*s^2 + k1*c*s + ((k1*k2)-(k2^2)),s,c)

取决于您想要的解决方案

于 2014-11-03T16:15:54.700 回答