function [d,v,a]=linear(para,np,per,xg,t,dt)
st=length(xg);
m=1;%先假定為單位質量；
k1=m*(2*pi/per)^2;%初始彈性剛度，計算時頻反應譜是一個系列；
%此為恢復力模型輸入參數；
eta=0.05;
txg=zeros(1,2);
x1=zeros(1,3);
x2=zeros(1,3);
x1(3)=-xg(1);
for nt=1:st    
	if (nt==1)%初始時刻位移和速度為零，加速度為負初始時刻地震加速度；
        d(1)=0;
        v(1)=0;
        a(1)=-xg(1);
        ff(1)=0;
	else
		txg(1)=xg(nt-1);
		txg(2)=xg(nt);
        [x2]=newmark_single(dt,k1,m,eta,txg,x1);
        f2=k1*x2(1);
		d(nt)=x2(1);
		v(nt)=x2(2);
		a(nt)=x2(3);
        x1=x2;
        ff(nt)=f2;
    end  
end