重采样，大侠帮帮忙 看不懂~~

xmuszq

重采样，大侠帮帮忙 看不懂~~
这是一个用转速a1来对振动信号a2进行重采样的matlab程序，
重采样部分没看懂是啥意思
大侠帮帮指点指点




load F:\F\ss2.txt;
a1=ss2(:,1);
a2=ss2(:,2);
%[b,a]=butter(8,1000/16384);
%a2=filter(b,a,a21);
k=zeros(1,40);
j=1;
for i=1:1:65535
   if a1(i)>-3 & a1(i+1)<-3
      k(j)=i;
      j=j+1;
   end   
end   
tim=zeros(1,j-2);
tim1=zeros(1,j-2);
% 获取每段脉冲长度
for i=1:1:j-2
   tim(i)=k(i+1)-k(i);
end   
% 重采样前滤波
%[b,a]=butter(8,1500/16384);
%a21=filter(b,a,a2);
a21=a2;
%a21=a2;
kk=1024;
speed=1564.5/4;  % 分析轴转速
data12=zeros(1,kk);


% 重采样
for mm=1:1:25
  tim(mm)=round(tim(mm)*4009.3/4/speed);% 脉冲变换后每段长度
  data1=a21(k(mm):(k(mm)+tim(mm)-1));  %截取的原数据
  k(mm+1)=k(mm)+tim(mm);
  delta=1/32768*(tim(mm)/kk); %重采样后时间间隔
  nn=tim(mm);

  for i=1:1:kk
    shike=round(i*tim(mm)/kk);
    data11=0;
    for j=-5:1:5
        %yy=j;
  %      canshu=pi*((i-1)*delta-(shike-1+j)/32768)/(1/32768);
        canshu=pi*((i)*delta-(shike+j)/32768)/(1/32768);
        
       if ((canshu==0))
          data11 =data11+a21(k(mm)+shike+j)*1;  
       else
         %hh=sin(canshu)/canshu;
         data11 =data11+a21(k(mm)+shike+j)*sin(canshu)/canshu;
       end

    end
    data12(i)=data12(i)+data11;
  end
end
data13=data12/mm; %重采样后平均