二自由度振动计算matlab