matlab开发-振动1的新标记。求解线性1自由度系统的纽马克方法。