Clear Filters
Clear Filters

how to use desoiking function?

2 views (last 30 days)
mehra
mehra on 29 Aug 2022
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

Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!