data_matrix.m

data_matrix.zip进行数据信号处理源程序

%形成数据矩阵
NN=6000 ;
S_NUMBER=4 ;

sample_rate=200 ;
dt=1.0/sample_rate ;

mydata=zeros(4,1) ;


acce_t=zeros(NN,1) ;
acce_ns=zeros(NN,1) ;
acce_ew=zeros(NN,1) ;
acce_8_17_ns=zeros(NN,1) ;
acce_8_17_ew=zeros(NN,1) ;

acce_matrix=zeros(NN,S_NUMBER) ;

data2=zeros(NN,1) ;
%读8月17日NS向数据
fid3=fopen('e:\calculating\dzc_8_17_ns.dat','r') ;
n=1 ;
while feof(fid3)==0
    mydata=fscanf(fid3,'%g %g %g %g\n',4) ;
    data2(n)=mydata(2) ;
    n=n+1 ;
end
for n=1:NN 
    acce_8_17_ns(n)=data2(n) ; %在数据开始前加MM个零
end
fclose(fid3) ;

data2=zeros(NN,1) ;
%读8月17日EW向数据
fid4=fopen('e:\calculating\dzc_8_17_ew.dat','r') ;
n=1 ;
while feof(fid4)==0
    mydata=fscanf(fid4,'%g %g %g %g\n',4) ;
    data2(n)=mydata(2) ;
    n=n+1 ;
end
for n=1:NN 
    acce_8_17_ew(n)=data2(n) ; %在数据开始前加MM个零
end
fclose(fid4) ;


data2=zeros(NN,1) ;
%读11月12日NS向数据
fid1=fopen('e:\calculating\dzc_11_12_ns5180.dat','r') ;
n=1 ;
while feof(fid1)==0
    mydata=fscanf(fid1,'%g %g %g %g\n',4) ;
    data2(n)=mydata(2) ;
    n=n+1 ;
end
for n=1:NN 
    acce_ns(n)=data2(n) ; %在数据开始前加MM个零
end
fclose(fid1) ;

data2=zeros(NN,1) ;
%读11月12日EW向数据
fid2=fopen('e:\calculating\dzc_11_12ew5182.dat','r') ;
n=1 ;
while feof(fid2)==0
    mydata=fscanf(fid2,'%g %g %g %g\n',4) ;
    data2(n)=mydata(2) ;
    n=n+1 ;
end
for n=1:NN 
    acce_ew(n)=data2(n) ; %在数据开始前加MM个零
end
fclose(fid2) ;

for n=1:NN
    acce_t(n)=dt*(n-1) ;
end

%形成数据矩阵
acce_matrix=[acce_8_17_ns,acce_8_17_ew,acce_ns,acce_ew] ;

fid1=fopen('e:\calculating\acce_matrix.txt','w') ;
fprintf(fid1,'%12.8g %12.8g %12.8g %12.8g\n',acce_matrix') ;
fclose(fid1) ;