时域同步平均算法

%%
%清除命令窗口，工作区
clc
clear all 
%%
%导入原始振动数据和方位角数据
filename1=‘A.txt‘%振动数据
data1=textread（filename1‘%n‘）;
filename2=‘B.txt‘%方位角数据
data2=textread（filename2‘%n‘）;
fs=1024;%采样频率
f=6.77;%时标所在轴频率
f_t=6.77;%感兴趣时标所在轴频率
n=length（data2）;
dic=zeros（n1）;
%%提取原时标对应振动信号的点数，存在dic矩阵中
for i=1:n
    integer=fix（（data2（i1）*fs）/（13107200））;
    dicimal=（data2（i1）*fs）/（13107200）-integer;
    if integer==（data2（i1）*fs）/（13107200）
       dic（i1）=fix（（data2（i1）*fs）/（13107200））;
    else
        dic（i1）=fix（（data2（i1）*fs）/（13107200））+1;
    end
end
%%等效时标转换
k=length（dic）;%原时标个数
dicn=fix（（f_t/f）*（k-1））+1;%等效时标个数
dic_n=zeros（dicn1）;
for j=1:dicn
    integer=fix（1+（j-1）*f/f_t）;%整数部分
    decimal=1+（j-1）*f/f_t-integer;%小数部分
    if integer~=k
       dic_n（j1）=dic（integer1）+round（decimal*（dic（integer+1）-dic（integer1）））;%round四舍五入，小数部分 
    else
    dic_n（j1）=dic（integer1）;
    end
end
dic_n=dic_n（2000:end1）;%我的数据里选择从第dic_n（2000，1）个时标开始