GetPC_GaussianFiltering.m 508 B

123456789101112131415
  1. function z = GetPC_GaussianFiltering(x, window, gaussFilter)
  2. for i=1:(window-1)/2
  3. temp = x((1:i+(window-1)/2),:);
  4. z(i,:) = (temp'*gaussFilter((window+1)/2+1-i:window,i))';
  5. end
  6. for i=(window+1)/2:size(x,1)-(window-1)/2
  7. temp = x((i-(window-1)/2:i+(window-(window+1)/2)),:);
  8. z(i,:) = (temp'*gaussFilter(:,(window+1)/2))';
  9. end
  10. j=1;
  11. for i=size(x,1)-(window-1)/2+1:size(x,1)
  12. temp=x((end-window+1+j:end),:);
  13. z(i,:) = (temp'*gaussFilter(1:window-j,(window+1)/2-j))'; j=j+1;
  14. end