随机减量法

本程序基于随机减量法 Random Decrement Technique对白噪声激励下的结构振动响应进行处理的，得到结构自由振动响应。

%随机减量法

clear
close all
format long

fprintf（‘Samples In Test :\n‘）;
path_base = sprintf（‘test‘）;
dir_sample = dir（path_base）;
[nsampledummy] = size（dir_sample）;
nsample = nsample - 2;
for i =1:nsample
    fprintf（strcat（num2str（i）‘:     ‘））;
    fprintf（dir_sample（i+2）.name）;
    fprintf（‘\n‘）;
end
SampleNum = input（‘Please select the sample number:\n‘）;

path_sample = strcat（path_base‘\‘dir_sample（SampleNum+2）.name）; %strcat：连接字符串的函数
name_sample = dir_sample（SampleNum+2）.name;

%提取测试数据
data（:1）=load（path_sample）; 
%所选的时域样本内数据个数
ntdata = length（data（:1））; 
Cvtdata= zeros（ntdata1）;


%【1】正常信号处理:
Cvtdata = data ;
%【1】正常信号处理end;


% %【2】Narada信号转换:
% for j=1:length（data（:1））
%     data（j1） = data（j1） * 5.0 / 65535.0;  %Convert the data 转化为电压信号
% end       
% 
% %增益系数
% gain = 10;
% for j=1:ntdata      
%     Cvtdata（j1）= data（j1）/（gain*200/1000）;
% end
% Avgdata = zeros（11）;
% Avgdata = mean（ Cvtdata（:1）  ）; %求数据平均值
% %【2】Narada信号转换end;

sf=800;
np = 100;
x = Cvtdata‘;

t=0:1/sf:（np-1）/sf;    %输出时间向量
nt=length（x）;    %提取输入数据的长度
s=1.5*std（x）;    %设置截取振幅A的大小
m=0;
y=zeros（1np）;
for k=2:nt-np
    a=abs（x（k-1）-s）;
    b=abs（x（k）-s）;
    c=abs（x（k+1）-s）;
    if b        y（1:np）=y（1:np）+x（k:k+np-1）;
        m=m+1;
    end
end

y=y./m;   %对叠加结果进行平均

%绘图
plot（ty）;
xlabel（‘时间（s）‘）;
ylabel（‘幅值‘）;
grid on;
% fid=fopen（fno‘w‘）;
% for k=1:np
%    fprintf（fid‘%f\n‘y（k））;
% end
% status=fclose（fid）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1928  2015-09-16 13:05  RD法\outRDT.mat

     文件       1657  2015-10-21 13:47  RD法\RDT.m

     文件      25600  2015-10-21 12:03  RD法\test\result1_2.txt

     文件      25600  2015-10-21 12:03  RD法\test\result2_2.txt

     文件     350023  2015-09-16 13:11  RD法\x值.xlsx

     文件      80192  2007-11-18 21:58  RD法\y.mat

     目录          0  2015-10-21 13:44  RD法\test

     目录          0  2015-10-21 13:44  RD法

----------- ---------  ---------- -----  ----

               485000                    8