cal.asv

kalman滤波

clear all;

% Calculate responses of a 3-story shear building subject to White Noise Excitation on Floor 
n=3;

%Input data

% Integrates ordinary differential equations y' = f(t,y) by Runge-Kutta, medium order method.
y0=zeros(2*n,1);   % initial conditions
OPTIONS = [];
[t,y]=ode45(@linear3dof,te,y0,OPTIONS,n,duration);

% Acceleration response
ll=length(y);
for it=1:ll 
   xx(:,it)=linear3dof(t(it),y(it,:)',n,duration);
end

xt=xx';
acc=xt(:,n+1:2*n);

save response.mat acc
save xstate.mat y