I need to put all vectors in the same window with a single origin
i have this code and i want to put all vectors in the same window with the same origin to draw the mean flow of fluid :
clear all
close all
% Initiate time counter
tic
%L’origine est à préciser
x01 = 7.75; %
y01 = 132;%
titre = ‘250rpm_bas_pos1’; % 05/02/2024
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Specify working directory and get total n images
% Right-click the document tab, and then select Copy Full Path to Clipboard
path = ‘D:PIV_Donnée250rpm_bas_pos1’; % % chemin pour lire les données
files=dir(fullfile(path,’*.csv’));
nb_images = size(files,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n_figure=0;
%Nom du fichier à lire
file = (‘3 (‘);
% Première lecture d’un champ de vecteurs vitesses pour définir les tailles
num=1;
% concatenation pour le nom de fichier
% name=strcat(file,num2str(num),’).csv’);
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
[m,n]=size(DATA);
L = 127; %127;
C=m/L; %127; % nombre de lignes
x_DATA=DATA(:,1);
x_mat=reshape(x_DATA,L,C);
x_vec=x_mat(:,1);
clear x_DATA
y_DATA=DATA(:,2);
y_mat=reshape(y_DATA,L,C);
y_vec=y_mat(1,:)’;
clear y_DATA
u=zeros(L,C,nb_images);
v=zeros(L,C,nb_images);
u_DATA=DATA(:,3);
u_mat=reshape(u_DATA,L,C);
u(:,:,1)=u_mat(:,:);
clear u_DATA
v_DATA=DATA(:,4);
v_mat=reshape(v_DATA,L,C);
v(:,:,1)=v_mat(:,:);
clear v_DATA
% Lecture des vecteurs vitesses
for num=2:nb_images
num=num-1;
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
u_DATA=DATA(:,3);
u_mat=reshape(u_DATA,L,C);
u(:,:,num)=u_mat(:,:);
clear u_DATA
v_DATA=DATA(:,4);
v_mat=reshape(v_DATA,L,C);
v(:,:,num)=v_mat(:,:);
clear v_DATA
end
clear u_mat v_mat
%%%%%%
% Initiate time counter
tic
%L’origine est à préciser
x02 = 7.75; %
y02 = 150;%
titre = ‘250rpm_middle_pos1’; % 06/02/2024
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Specify working directory and get total n images
% Right-click the document tab, and then select Copy Full Path to Clipboard
path = ‘D:PIV_Donnée250rpm_middle_pos1’; % % chemin pour lire les données
files=dir(fullfile(path,’*.csv’));
nb_images = size(files,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n_figure=0;
%Nom du fichier à lire
file = (‘7 (‘);
% Première lecture d’un champ de vecteurs vitesses pour définir les tailles
num=1;
% concatenation pour le nom de fichier
% name=strcat(file,num2str(num),’).csv’);
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
[m,n]=size(DATA);
L = 127; %127;
C=m/L; %127; % nombre de lignes
x1_DATA=DATA(:,1);
x1_mat=reshape(x1_DATA,L,C);
x1_vec=x1_mat(:,1);
clear x1_DATA
y1_DATA=DATA(:,2);
y1_mat=reshape(y1_DATA,L,C);
y1_vec=y1_mat(1,:)’;
clear y1_DATA
u1=zeros(L,C,nb_images);
v1=zeros(L,C,nb_images);
u1_DATA=DATA(:,3);
u1_mat=reshape(u1_DATA,L,C);
u1(:,:,1)=u1_mat(:,:);
clear u1_DATA
v1_DATA=DATA(:,4);
v1_mat=reshape(v1_DATA,L,C);
v1(:,:,1)=v1_mat(:,:);
clear v1_DATA
% Lecture des vecteurs vitesses
for num=2:nb_images
num=num-1;
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
u1_DATA=DATA(:,3);
u1_mat=reshape(u1_DATA,L,C);
u1(:,:,num)=u1_mat(:,:);
clear u1_DATA
v1_DATA=DATA(:,4);
v1_mat=reshape(v1_DATA,L,C);
v1(:,:,num)=v1_mat(:,:);
clear v1_DATA
end
clear u1_mat v1_mat
%%%%
% Initiate time counter
tic
%L’origine est à préciser
x03 = 7.75; %
y03 =419.5;%
titre = ‘250rpm_haut_pos1’; % 06/02/2024
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Specify working directory and get total n images
% Right-click the document tab, and then select Copy Full Path to Clipboard
path = ‘D:PIV_Donnée250rpm_haut_pos1’; % % chemin pour lire les données
files=dir(fullfile(path,’*.csv’));
nb_images = size(files,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n_figure=0;
%Nom du fichier à lire
file = (‘5 (‘);
% Première lecture d’un champ de vecteurs vitesses pour définir les tailles
num=1;
% concatenation pour le nom de fichier
% name=strcat(file,num2str(num),’).csv’);
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
[m,n]=size(DATA);
L = 127; %127;
C=m/L; %127; % nombre de lignes
x2_DATA=DATA(:,1);
x2_mat=reshape(x2_DATA,L,C);
x2_vec=x2_mat(:,1);
clear x2_DATA
y2_DATA=DATA(:,2);
y2_mat=reshape(y2_DATA,L,C);
y2_vec=y2_mat(1,:)’;
clear y2_DATA
u2=zeros(L,C,nb_images);
v2=zeros(L,C,nb_images);
u2_DATA=DATA(:,3);
u2_mat=reshape(u2_DATA,L,C);
u2(:,:,1)=u2_mat(:,:);
clear u2_DATA
v2_DATA=DATA(:,4);
v2_mat=reshape(v2_DATA,L,C);
v2(:,:,1)=v2_mat(:,:);
clear v2_DATA
% Lecture des vecteurs vitesses
for num=2:nb_images
num=num-1;
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
u2_DATA=DATA(:,3);
u2_mat=reshape(u2_DATA,L,C);
u2(:,:,num)=u2_mat(:,:);
clear u2_DATA
v2_DATA=DATA(:,4);
v2_mat=reshape(v2_DATA,L,C);
v2(:,:,num)=v2_mat(:,:);
clear v2_DATA
end
clear u2_mat v2_mat
%%%%
% Initiate time counter
tic
%L’origine est à préciser
x04 = 7.75; %
y04 = 450;%
titre = ‘250Rpm_haut_surfacelibre_pos1’; % 06/02/2024
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Specify working directory and get total n images
% Right-click the document tab, and then select Copy Full Path to Clipboard
path = ‘D:PIV_Donnée250Rpm_haut_surfacelibre_pos1’; % % chemin pour lire les données
files=dir(fullfile(path,’*.csv’));
nb_images = size(files,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n_figure=0;
%Nom du fichier à lire
file = (‘8 (‘);
% Première lecture d’un champ de vecteurs vitesses pour définir les tailles
num=1;
% concatenation pour le nom de fichier
% name=strcat(file,num2str(num),’).csv’);
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
[m,n]=size(DATA);
L = 127; %127;
C=m/L; %127; % nombre de lignes
x3_DATA=DATA(:,1);
x3_mat=reshape(x3_DATA,L,C);
x3_vec=x3_mat(:,1);
clear x3_DATA
y3_DATA=DATA(:,2);
y3_mat=reshape(y3_DATA,L,C);
y3_vec=y3_mat(1,:)’;
clear y3_DATA
u3=zeros(L,C,nb_images);
v3=zeros(L,C,nb_images);
u3_DATA=DATA(:,3);
u3_mat=reshape(u3_DATA,L,C);
u3(:,:,1)=u3_mat(:,:);
clear u3_DATA
v3_DATA=DATA(:,4);
v3_mat=reshape(v3_DATA,L,C);
v3(:,:,1)=v3_mat(:,:);
clear v3_DATA
% Lecture des vecteurs vitesses
for num=2:nb_images
num=num-1;
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
u3_DATA=DATA(:,3);
u3_mat=reshape(u3_DATA,L,C);
u3(:,:,num)=u3_mat(:,:);
clear u3_DATA
v3_DATA=DATA(:,4);
v3_mat=reshape(v3_DATA,L,C);
v3(:,:,num)=v3_mat(:,:);
clear v3_DATA
end
clear u3_mat v3_mati have this code and i want to put all vectors in the same window with the same origin to draw the mean flow of fluid :
clear all
close all
% Initiate time counter
tic
%L’origine est à préciser
x01 = 7.75; %
y01 = 132;%
titre = ‘250rpm_bas_pos1’; % 05/02/2024
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Specify working directory and get total n images
% Right-click the document tab, and then select Copy Full Path to Clipboard
path = ‘D:PIV_Donnée250rpm_bas_pos1’; % % chemin pour lire les données
files=dir(fullfile(path,’*.csv’));
nb_images = size(files,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n_figure=0;
%Nom du fichier à lire
file = (‘3 (‘);
% Première lecture d’un champ de vecteurs vitesses pour définir les tailles
num=1;
% concatenation pour le nom de fichier
% name=strcat(file,num2str(num),’).csv’);
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
[m,n]=size(DATA);
L = 127; %127;
C=m/L; %127; % nombre de lignes
x_DATA=DATA(:,1);
x_mat=reshape(x_DATA,L,C);
x_vec=x_mat(:,1);
clear x_DATA
y_DATA=DATA(:,2);
y_mat=reshape(y_DATA,L,C);
y_vec=y_mat(1,:)’;
clear y_DATA
u=zeros(L,C,nb_images);
v=zeros(L,C,nb_images);
u_DATA=DATA(:,3);
u_mat=reshape(u_DATA,L,C);
u(:,:,1)=u_mat(:,:);
clear u_DATA
v_DATA=DATA(:,4);
v_mat=reshape(v_DATA,L,C);
v(:,:,1)=v_mat(:,:);
clear v_DATA
% Lecture des vecteurs vitesses
for num=2:nb_images
num=num-1;
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
u_DATA=DATA(:,3);
u_mat=reshape(u_DATA,L,C);
u(:,:,num)=u_mat(:,:);
clear u_DATA
v_DATA=DATA(:,4);
v_mat=reshape(v_DATA,L,C);
v(:,:,num)=v_mat(:,:);
clear v_DATA
end
clear u_mat v_mat
%%%%%%
% Initiate time counter
tic
%L’origine est à préciser
x02 = 7.75; %
y02 = 150;%
titre = ‘250rpm_middle_pos1’; % 06/02/2024
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Specify working directory and get total n images
% Right-click the document tab, and then select Copy Full Path to Clipboard
path = ‘D:PIV_Donnée250rpm_middle_pos1’; % % chemin pour lire les données
files=dir(fullfile(path,’*.csv’));
nb_images = size(files,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n_figure=0;
%Nom du fichier à lire
file = (‘7 (‘);
% Première lecture d’un champ de vecteurs vitesses pour définir les tailles
num=1;
% concatenation pour le nom de fichier
% name=strcat(file,num2str(num),’).csv’);
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
[m,n]=size(DATA);
L = 127; %127;
C=m/L; %127; % nombre de lignes
x1_DATA=DATA(:,1);
x1_mat=reshape(x1_DATA,L,C);
x1_vec=x1_mat(:,1);
clear x1_DATA
y1_DATA=DATA(:,2);
y1_mat=reshape(y1_DATA,L,C);
y1_vec=y1_mat(1,:)’;
clear y1_DATA
u1=zeros(L,C,nb_images);
v1=zeros(L,C,nb_images);
u1_DATA=DATA(:,3);
u1_mat=reshape(u1_DATA,L,C);
u1(:,:,1)=u1_mat(:,:);
clear u1_DATA
v1_DATA=DATA(:,4);
v1_mat=reshape(v1_DATA,L,C);
v1(:,:,1)=v1_mat(:,:);
clear v1_DATA
% Lecture des vecteurs vitesses
for num=2:nb_images
num=num-1;
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
u1_DATA=DATA(:,3);
u1_mat=reshape(u1_DATA,L,C);
u1(:,:,num)=u1_mat(:,:);
clear u1_DATA
v1_DATA=DATA(:,4);
v1_mat=reshape(v1_DATA,L,C);
v1(:,:,num)=v1_mat(:,:);
clear v1_DATA
end
clear u1_mat v1_mat
%%%%
% Initiate time counter
tic
%L’origine est à préciser
x03 = 7.75; %
y03 =419.5;%
titre = ‘250rpm_haut_pos1’; % 06/02/2024
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Specify working directory and get total n images
% Right-click the document tab, and then select Copy Full Path to Clipboard
path = ‘D:PIV_Donnée250rpm_haut_pos1’; % % chemin pour lire les données
files=dir(fullfile(path,’*.csv’));
nb_images = size(files,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n_figure=0;
%Nom du fichier à lire
file = (‘5 (‘);
% Première lecture d’un champ de vecteurs vitesses pour définir les tailles
num=1;
% concatenation pour le nom de fichier
% name=strcat(file,num2str(num),’).csv’);
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
[m,n]=size(DATA);
L = 127; %127;
C=m/L; %127; % nombre de lignes
x2_DATA=DATA(:,1);
x2_mat=reshape(x2_DATA,L,C);
x2_vec=x2_mat(:,1);
clear x2_DATA
y2_DATA=DATA(:,2);
y2_mat=reshape(y2_DATA,L,C);
y2_vec=y2_mat(1,:)’;
clear y2_DATA
u2=zeros(L,C,nb_images);
v2=zeros(L,C,nb_images);
u2_DATA=DATA(:,3);
u2_mat=reshape(u2_DATA,L,C);
u2(:,:,1)=u2_mat(:,:);
clear u2_DATA
v2_DATA=DATA(:,4);
v2_mat=reshape(v2_DATA,L,C);
v2(:,:,1)=v2_mat(:,:);
clear v2_DATA
% Lecture des vecteurs vitesses
for num=2:nb_images
num=num-1;
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
u2_DATA=DATA(:,3);
u2_mat=reshape(u2_DATA,L,C);
u2(:,:,num)=u2_mat(:,:);
clear u2_DATA
v2_DATA=DATA(:,4);
v2_mat=reshape(v2_DATA,L,C);
v2(:,:,num)=v2_mat(:,:);
clear v2_DATA
end
clear u2_mat v2_mat
%%%%
% Initiate time counter
tic
%L’origine est à préciser
x04 = 7.75; %
y04 = 450;%
titre = ‘250Rpm_haut_surfacelibre_pos1’; % 06/02/2024
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Specify working directory and get total n images
% Right-click the document tab, and then select Copy Full Path to Clipboard
path = ‘D:PIV_Donnée250Rpm_haut_surfacelibre_pos1’; % % chemin pour lire les données
files=dir(fullfile(path,’*.csv’));
nb_images = size(files,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n_figure=0;
%Nom du fichier à lire
file = (‘8 (‘);
% Première lecture d’un champ de vecteurs vitesses pour définir les tailles
num=1;
% concatenation pour le nom de fichier
% name=strcat(file,num2str(num),’).csv’);
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
[m,n]=size(DATA);
L = 127; %127;
C=m/L; %127; % nombre de lignes
x3_DATA=DATA(:,1);
x3_mat=reshape(x3_DATA,L,C);
x3_vec=x3_mat(:,1);
clear x3_DATA
y3_DATA=DATA(:,2);
y3_mat=reshape(y3_DATA,L,C);
y3_vec=y3_mat(1,:)’;
clear y3_DATA
u3=zeros(L,C,nb_images);
v3=zeros(L,C,nb_images);
u3_DATA=DATA(:,3);
u3_mat=reshape(u3_DATA,L,C);
u3(:,:,1)=u3_mat(:,:);
clear u3_DATA
v3_DATA=DATA(:,4);
v3_mat=reshape(v3_DATA,L,C);
v3(:,:,1)=v3_mat(:,:);
clear v3_DATA
% Lecture des vecteurs vitesses
for num=2:nb_images
num=num-1;
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
u3_DATA=DATA(:,3);
u3_mat=reshape(u3_DATA,L,C);
u3(:,:,num)=u3_mat(:,:);
clear u3_DATA
v3_DATA=DATA(:,4);
v3_mat=reshape(v3_DATA,L,C);
v3(:,:,num)=v3_mat(:,:);
clear v3_DATA
end
clear u3_mat v3_mat i have this code and i want to put all vectors in the same window with the same origin to draw the mean flow of fluid :
clear all
close all
% Initiate time counter
tic
%L’origine est à préciser
x01 = 7.75; %
y01 = 132;%
titre = ‘250rpm_bas_pos1’; % 05/02/2024
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Specify working directory and get total n images
% Right-click the document tab, and then select Copy Full Path to Clipboard
path = ‘D:PIV_Donnée250rpm_bas_pos1’; % % chemin pour lire les données
files=dir(fullfile(path,’*.csv’));
nb_images = size(files,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n_figure=0;
%Nom du fichier à lire
file = (‘3 (‘);
% Première lecture d’un champ de vecteurs vitesses pour définir les tailles
num=1;
% concatenation pour le nom de fichier
% name=strcat(file,num2str(num),’).csv’);
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
[m,n]=size(DATA);
L = 127; %127;
C=m/L; %127; % nombre de lignes
x_DATA=DATA(:,1);
x_mat=reshape(x_DATA,L,C);
x_vec=x_mat(:,1);
clear x_DATA
y_DATA=DATA(:,2);
y_mat=reshape(y_DATA,L,C);
y_vec=y_mat(1,:)’;
clear y_DATA
u=zeros(L,C,nb_images);
v=zeros(L,C,nb_images);
u_DATA=DATA(:,3);
u_mat=reshape(u_DATA,L,C);
u(:,:,1)=u_mat(:,:);
clear u_DATA
v_DATA=DATA(:,4);
v_mat=reshape(v_DATA,L,C);
v(:,:,1)=v_mat(:,:);
clear v_DATA
% Lecture des vecteurs vitesses
for num=2:nb_images
num=num-1;
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
u_DATA=DATA(:,3);
u_mat=reshape(u_DATA,L,C);
u(:,:,num)=u_mat(:,:);
clear u_DATA
v_DATA=DATA(:,4);
v_mat=reshape(v_DATA,L,C);
v(:,:,num)=v_mat(:,:);
clear v_DATA
end
clear u_mat v_mat
%%%%%%
% Initiate time counter
tic
%L’origine est à préciser
x02 = 7.75; %
y02 = 150;%
titre = ‘250rpm_middle_pos1’; % 06/02/2024
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Specify working directory and get total n images
% Right-click the document tab, and then select Copy Full Path to Clipboard
path = ‘D:PIV_Donnée250rpm_middle_pos1’; % % chemin pour lire les données
files=dir(fullfile(path,’*.csv’));
nb_images = size(files,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n_figure=0;
%Nom du fichier à lire
file = (‘7 (‘);
% Première lecture d’un champ de vecteurs vitesses pour définir les tailles
num=1;
% concatenation pour le nom de fichier
% name=strcat(file,num2str(num),’).csv’);
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
[m,n]=size(DATA);
L = 127; %127;
C=m/L; %127; % nombre de lignes
x1_DATA=DATA(:,1);
x1_mat=reshape(x1_DATA,L,C);
x1_vec=x1_mat(:,1);
clear x1_DATA
y1_DATA=DATA(:,2);
y1_mat=reshape(y1_DATA,L,C);
y1_vec=y1_mat(1,:)’;
clear y1_DATA
u1=zeros(L,C,nb_images);
v1=zeros(L,C,nb_images);
u1_DATA=DATA(:,3);
u1_mat=reshape(u1_DATA,L,C);
u1(:,:,1)=u1_mat(:,:);
clear u1_DATA
v1_DATA=DATA(:,4);
v1_mat=reshape(v1_DATA,L,C);
v1(:,:,1)=v1_mat(:,:);
clear v1_DATA
% Lecture des vecteurs vitesses
for num=2:nb_images
num=num-1;
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
u1_DATA=DATA(:,3);
u1_mat=reshape(u1_DATA,L,C);
u1(:,:,num)=u1_mat(:,:);
clear u1_DATA
v1_DATA=DATA(:,4);
v1_mat=reshape(v1_DATA,L,C);
v1(:,:,num)=v1_mat(:,:);
clear v1_DATA
end
clear u1_mat v1_mat
%%%%
% Initiate time counter
tic
%L’origine est à préciser
x03 = 7.75; %
y03 =419.5;%
titre = ‘250rpm_haut_pos1’; % 06/02/2024
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Specify working directory and get total n images
% Right-click the document tab, and then select Copy Full Path to Clipboard
path = ‘D:PIV_Donnée250rpm_haut_pos1’; % % chemin pour lire les données
files=dir(fullfile(path,’*.csv’));
nb_images = size(files,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n_figure=0;
%Nom du fichier à lire
file = (‘5 (‘);
% Première lecture d’un champ de vecteurs vitesses pour définir les tailles
num=1;
% concatenation pour le nom de fichier
% name=strcat(file,num2str(num),’).csv’);
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
[m,n]=size(DATA);
L = 127; %127;
C=m/L; %127; % nombre de lignes
x2_DATA=DATA(:,1);
x2_mat=reshape(x2_DATA,L,C);
x2_vec=x2_mat(:,1);
clear x2_DATA
y2_DATA=DATA(:,2);
y2_mat=reshape(y2_DATA,L,C);
y2_vec=y2_mat(1,:)’;
clear y2_DATA
u2=zeros(L,C,nb_images);
v2=zeros(L,C,nb_images);
u2_DATA=DATA(:,3);
u2_mat=reshape(u2_DATA,L,C);
u2(:,:,1)=u2_mat(:,:);
clear u2_DATA
v2_DATA=DATA(:,4);
v2_mat=reshape(v2_DATA,L,C);
v2(:,:,1)=v2_mat(:,:);
clear v2_DATA
% Lecture des vecteurs vitesses
for num=2:nb_images
num=num-1;
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
u2_DATA=DATA(:,3);
u2_mat=reshape(u2_DATA,L,C);
u2(:,:,num)=u2_mat(:,:);
clear u2_DATA
v2_DATA=DATA(:,4);
v2_mat=reshape(v2_DATA,L,C);
v2(:,:,num)=v2_mat(:,:);
clear v2_DATA
end
clear u2_mat v2_mat
%%%%
% Initiate time counter
tic
%L’origine est à préciser
x04 = 7.75; %
y04 = 450;%
titre = ‘250Rpm_haut_surfacelibre_pos1’; % 06/02/2024
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Specify working directory and get total n images
% Right-click the document tab, and then select Copy Full Path to Clipboard
path = ‘D:PIV_Donnée250Rpm_haut_surfacelibre_pos1’; % % chemin pour lire les données
files=dir(fullfile(path,’*.csv’));
nb_images = size(files,1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n_figure=0;
%Nom du fichier à lire
file = (‘8 (‘);
% Première lecture d’un champ de vecteurs vitesses pour définir les tailles
num=1;
% concatenation pour le nom de fichier
% name=strcat(file,num2str(num),’).csv’);
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
[m,n]=size(DATA);
L = 127; %127;
C=m/L; %127; % nombre de lignes
x3_DATA=DATA(:,1);
x3_mat=reshape(x3_DATA,L,C);
x3_vec=x3_mat(:,1);
clear x3_DATA
y3_DATA=DATA(:,2);
y3_mat=reshape(y3_DATA,L,C);
y3_vec=y3_mat(1,:)’;
clear y3_DATA
u3=zeros(L,C,nb_images);
v3=zeros(L,C,nb_images);
u3_DATA=DATA(:,3);
u3_mat=reshape(u3_DATA,L,C);
u3(:,:,1)=u3_mat(:,:);
clear u3_DATA
v3_DATA=DATA(:,4);
v3_mat=reshape(v3_DATA,L,C);
v3(:,:,1)=v3_mat(:,:);
clear v3_DATA
% Lecture des vecteurs vitesses
for num=2:nb_images
num=num-1;
name=fullfile(path,strcat(file,num2str(num),’).csv’));
DATA = csvread(name,9);
u3_DATA=DATA(:,3);
u3_mat=reshape(u3_DATA,L,C);
u3(:,:,num)=u3_mat(:,:);
clear u3_DATA
v3_DATA=DATA(:,4);
v3_mat=reshape(v3_DATA,L,C);
v3(:,:,num)=v3_mat(:,:);
clear v3_DATA
end
clear u3_mat v3_mat vector, superposition MATLAB Answers — New Questions
