%随机减量法模态参数识别预处理
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear
clc
close all hidden
format long
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
fni=input('随机减量法模态参数识别预处理-输入数据文件名:','s');
fid=fopen(fni,'r');
sf = fscanf(fid,'%f',1);     %采样频率
np = fscanf(fid,'%d',1);     %输出数据长度
fno = fscanf(fid,'%s',1);    %输出数据文件名
x = fscanf(fid,'%f',[1,inf]);%输入时程信号数据存成行向量
status=fclose(fid);
%建立离散输出时间向量
t=0:1/sf:(np-1)/sf;
%取输入数据长度
nt=length(x);
%设置截取振幅为输入信号标准差的1.5倍
s=1.5*std(x);
%获取输入信号的子样本函数进行叠加
m=0;
y=zeros(1,np);
for k=2:nt-np
  a=abs(x(k-1)-s);
  b=abs(x(k)-s);
  c=abs(x(k+1)-s);
  if b<a & b<c
    y(1:np)=y(1:np)+x(k:k+np-1); m=m+1;
  end
end
%对叠加结果做平均
y=y./m;
%绘制自由振动时程曲线图
plot(t,y);
xlabel('时间 (s)');  
ylabel('幅值');
grid on;                        

%打开文件输出随机减量法计算的结果               
fid=fopen(fno,'w');   
for k=1:np                           
  fprintf(fid,'%f\n',y(k));
end                           
status=fclose(fid);