计算多自由度非经典阻尼线性系统在动力作用下的响应。 （使用纽马克的方法）
函数结果=Newmark_Linear(M,C,K,f,fs)
输入M：质量矩阵 (n*n) C:阻尼矩阵 (n*n) K：刚度矩阵 (n*n) f:外力矩阵(n,N) fs：采样频率其中n是自由度数，N是动力数据点的长度
输出： 结果：是一个结构，由Result.Displacement: 位移 (n*N) 结果.速度：速度 (n*N) Result.Acceleration: 加速度 (n*N)
参考： Chopra, Anil K.“结构动力学。理论和应用”。 地震工程（2017）。
