how to use desoiking function?

Hi every one,

I need to use func_despike_phasespace3d functions in order to despike my data. But because input data for this function should be with dimension (n,1), and my data had the dimension of (n,30) I got some errors, later I changed my code as following, right now am not getting any error but I am not sure if the despiking function is working properly or not. Could any one help me with this?

clear,
close all;
clc;

%% BASIC PARAMETERS/CONSTANTS
% Constant variables 

nu=1*10^-6;  % Kinematic viscosity [m^2/s]
ro=1000;
Mu=nu.*ro;
B=0.70; % Flume width(m)
g=9.807;    % Gravitational acceleration [m/s^2]

% Varying Parameters
%Profile
Q_p=89.8/3600; % m^3/h Discharge (pump monitor)
V_p=1.994; % m/s velocity (pump monitor)
h0=0.25; % Water depth(m)
Area=B*h0; % Wetted area(m^2)
P=B+2.*h0; % Wetted perimeter(m)
Rh_rad=Area/P; % Hydraulic radius(m)
V_f=Q_p/Area;

% vel_data=dir('V*.mat');
% N=length(vel_data);
% load(vel_data.name)
load('VectrinoData.235.22.Vectrino Profiler.00002.mat')
k=length(Data.Profiles_TimeStamp);
z=0.029:-0.001:0; % Starts from 
% z=0:0.001:0.029; % Starts from 
results=zeros(4,30);
ur=zeros(k,30);
vr=zeros(k,30);
w1r=zeros(k,30);
w2r=zeros(k,30);
u=zeros(k,30);
v=zeros(k,30);
w1=zeros(k,30);
w2=zeros(k,30);
w=zeros(k,30);
for np=1:30;
    ur(:,np)=Data.Profiles_VelX(:,np);
    vr(:,np)=Data.Profiles_VelY(:,np);
    w1r(:,np)=Data.Profiles_VelZ1(:,np);
    w2r(:,np)=Data.Profiles_VelZ2(:,np);
    [u(:,np), ip1] = func_despike_phasespace3d( ur(:,np),8, 2);
    [v(:,np), ip2] = func_despike_phasespace3d( vr(:,np),8, 2);
    [w1(:,np), ip3] = func_despike_phasespace3d( w1r(:,np), 8, 2);
    [w2(:,np), ip4] = func_despike_phasespace3d( w2r(:,np), 8, 2);
    w(:,np)=(w1(:,np)+w2(:,np))/2;
    results(1,np)=mean(Data.BottomCheck_BottomDistance)-0.001*np;
    results(2,np)=mean(u(:,np),1);
    results(3,np)=mean(v(:,np),1);
    results(4,np)=mean(w(:,np),1);
end

Thanks