采用Wilson-θ法计算其位移响应，并作出前10s内每个质点的时间位移曲线。m1=m2=1，c1=c2=1，k1=1,k2=5,f1=0,f2=sin(w2*t)，w1=w21,x1=x1'=x2=x2'=0,matlab代码

好的，根据您的要求，我会为您解答这个问题。首先，我们需要根据题目给出的参数，建立系统的动力学方程：

m1 * x1'' + c1 * (x1' - x2') + k1 * (x1 - x2) = f1
m2 * x2'' + c2 * (x2' - x1') + k2 * (x2 - x1) = f2

其中，x1、x1'、x1''分别表示质点1的位移、速度和加速度，x2、x2'、x2''分别表示质点2的位移、速度和加速度。f1和f2分别表示两个质点受到的外力。

根据Wilson-θ法，我们可以将动力学方程离散化，得到如下的求解公式：

M * (X - Xn) / (θ * Δt)^2 + C * (X - Xn) / (θ * Δt) + K * X = F

其中，M、C、K和F分别表示质量矩阵、阻尼矩阵、刚度矩阵和外力矩阵。X和Xn分别表示当前时刻和上一时刻的位移矩阵。θ和Δt分别表示Wilson-θ法中的参数和时间步长。在此问题中，我们可以取θ=0.5，Δt=0.01s。

根据上述公式，我们可以编写如下的Matlab代码：