Homec4science
Diffusion PP routine (Snowy-2.0) Browse

R11173/read_v1.m
f14a1b0eee55Snowy-2.0

Last Change
Raw File

read_v1.m
View Options

This document is not UTF8. It was detected as Shift JIS and converted to UTF8 for display.
function [nb,dat_cp,freq,time,MP,n11,Q11,sigma,nm,hm,Hp,alpha,Em,rho,Pp,Qp,coeff_n,coeff_q,channel,folder,folder_results, folder_save, F_pump, F_pump2, nbsensor, length_sig,date1]= read_v1(name,coeff,projectname)
%% IMPORTANT : MUST BE UPDATED IF THE NAME OF THE ROOTS IS MODIFIED !!!
% the function mfilename doesnt work in the standalone app, if the code is
% run by the standalone app ( if isdeployed), il faut donner directement le
% nom du dossier ou sont les codes matlab.
%% 3.1 - Definition of the folder variable. The intermediate files will be stored somewhere in the folder variable
if isdeployed
folder = '\\lmhnas2\gem\07_Software\02_Matlab\PP modeles\Matlab';
else
% On ajoute au matlab path le chemin d'acc鑚 du code actuel.
folder = fileparts(which(mfilename));
addpath(genpath(folder));
end
% creates a folder where the pdf and the .fig will be saved.At the end of
% the code, all these files will be copied to the working folder
folder_save=fullfile(folder, 'Results','Saves\');
%Where Results will be saved
slash=find(name=='\');
l=max(slash);
folder_results=name(1:(l-1));
%% 3.2 - Creates the internal save folder (if they don't already exist)
if ~exist([folder '\Results'],'dir')
mkdir(folder,'Results')
end
if ~exist([folder '\Results\Saves'],'dir')
mkdir(folder,'Results\Saves')
end
if ~exist([folder '\Results\images'],'dir')
mkdir(folder,'Results\images')
end
if ~exist([folder '\Results\images\' projectname],'dir')
mkdir(folder,['Results\images\',projectname])
end
% Creates the folder where all the final files will be copied.
%THis will be at the same location than the tdms files but in others
%folders.
if ~exist([folder_results '\PDF'],'dir')
mkdir(folder_results,'PDF')
end
if ~exist([folder_results '\EXCEL'],'dir')
mkdir(folder_results,'EXCEL')
end
if ~exist([folder_results '\EXCEL\signal_post'],'dir')
mkdir([folder_results,'\EXCEL'],'signal_post')
end
if ~exist([folder_results '\FIG'],'dir')
mkdir(folder_results,'FIG')
end
Files=dir(name); %Where all the tdms are.
%% 3.3 get the date of measurement. This date will then be written on the top
% right of the pdf files.
[Y, M, D, ~, ~, ~]= datevec(Files(end).date); % Gets the date of measurement of the last tdms file
Y_s=sprintf('%.20f',Y);
idx=regexp(Y_s,'0+','start');
Y_s(idx(end):end)=[];
M_s=sprintf('%.20f',M);
idx=regexp(M_s,'0+','start');
M_s(idx(end):end)=[];
if M<10
M_s = ['0' M_s]; % adds a zero to the date to have format xx.xx.xxxx instead of x.x.xxxx
end
D_s=sprintf('%.20f',D);
idx=regexp(D_s,'0+','start');
D_s(idx(end):end)=[];
if D<10
D_s = ['0' D_s]; % adds a zero to the date to have format xx.xx.xxxx instead of x.x.xxxx
end
date2 = [D_s,M_s,Y_s];
date1=date2(1:end-1);
L1=length(Files);
%% 3.4 - Delete index files if necessary
compt=1;
for k=1:L1
FileNames=Files(k).name;
d=strfind(FileNames,'index');
if isempty(d) ==0
ind_del(compt) = k;
compt=compt+1;
end
end
A= exist('ind_del')
if A == 1
Files(ind_del)=[];
end
%% 3.5 - Read and store the data from the tdms files :
L1 =length(Files);
for k=3:L1 % We start with 3 to avoid "." and ".." always present from function dir
FileNames=Files(k).name;
totalname=fullfile(name,FileNames);
MyData(k-2) = TDMS_readTDMSFile(totalname); % Mydata contains all the tdms
end
nb = L1-2; %nb of measurements points
length_sig = length(cell2mat(MyData(1).data(19)));
nbsensor= length(MyData(1).chanNames{1,1})-16;
%% Order My data : ReadTDMS doesn't read the file in proper order we want to reorder them according to their MP files
for i=1:nb
ind_order(i)=cell2mat(MyData(i).data(5));
end
MyData2=MyData;
for i =1:nb
[~, ind_min]=min(ind_order);
MyData2(i)=MyData(ind_min);
ind_order(ind_min)=99999;
end
MyData=MyData2;
%% 3.6 - Convert the cell containig all data in a matrix
%Find where the data begins in the tdms document :
start = find(strcmp(MyData(1).chanNames{1,1},'Channel 1'));
disp(start)
dat= ones(nb,nbsensor,length_sig);
disp(length_sig)
for i=1:nb % nb being the number of measurements points
for j=1:nbsensor
dat(i,j,:) = cell2mat(MyData(i).data(start+1+j));
end
end
% dat contains now all the pressure fluctuation data.
%% Center and permute the signal
moy= mean(dat,3);
dat_c= dat-moy; % data centered
dat_cp=permute(dat_c,[3,1,2]); %data centered and permuted
%% 3.7 - get all interesting data :
freq = cell2mat(MyData(1).data(3)); % frequency of acquisition
time= cell2mat(MyData(1).data(4)); % time of aquisition
%preallocating
MP=zeros(1,nb);
n11=zeros(1,nb);
Q11=zeros(1,nb);
sigma=zeros(1,nb);
nm=zeros(1,nb);
Hp=zeros(1,nb);
alpha=zeros(1,nb);
Em=zeros(1,nb);
rho=zeros(1,nb);
Pp=zeros(1,nb);
Qp=zeros(1,nb);
for j = 1 : nb
MP(j)= round(cell2mat(MyData(j).data(5))); % MP
n11(j)= round(cell2mat(MyData(j).data(6)),3); % n11
Q11(j)= round(cell2mat(MyData(j).data(7)),3); %q11
sigma(j)= round(cell2mat(MyData(j).data(8)),3); % sigma
nm(j)= round(cell2mat(MyData(j).data(9)),1); % rotation speed
Hp(j)= round(cell2mat(MyData(j).data(10)),2); % prototype head
alpha(j)= round(cell2mat(MyData(j).data(11)),1); % angle
Em(j)=cell2mat(MyData(j).data(10));1; % Energy
rho(j)= (cell2mat(MyData(j).data(13)));
Hp(j)= (cell2mat(MyData(j).data(14)));
Pp(j)= (cell2mat(MyData(j).data(15)));
Qp(j)= (cell2mat(MyData(j).data(16)));
F_pump(j)=(cell2mat(MyData(j).data(17)))/60;
F_pump2(j)=(cell2mat(MyData(j).data(18)))/60;
end
hm=round(Em./9.8063,2);
% %% 02.11.2020 - modification Zhen'an
% % Read the results Excelfiles
% B = readcell(fullfile(folder_results,'724-05_ZhenAn_MAT.xlsx'));
% B(1,1)={'MP'}; % Renane MP in results file
%
% %% Check for blank lines, suppress them if any :
% row_missing=cellfun(@(x) all(ismissing(x)), B(:,1));
% B=B(~row_missing,:);
%
% %% Convert to double to transpose the reynold and the pressure to desired units
% B(cellfun(@(x) any(ismissing(x)), B)) = {1000000}; % On convertit on num駻ique les cases qui ont du texte pour garder le bon format.
% B2=B(3:end,:);
% idx = find(cellfun(@isempty,B2));
% idx2 = find(cellfun(@ischar,B2));
% B2(idx)={1}; % They might be some empty elements preventing from the cell2mat function to work proprely
% B2(idx2)={1}; % They might be some char elements preventing from the cell2mat function to work proprely
% B3=cell2mat(B2);
% [line,~]=find((B3(:,1))==MP);
%
% Em= B3(line,5)'; % Energy
% hm= round(B3(line,41),3)';
% n11= round(B3(line,15),3)';
% Q11=round(B3(line,19),4)';
% Qp=B3(line,32)';
%% 3.8 - Write proprely the chosen dimensionless coefficient in a string variable
%Relevant speed and discharge coefficients :
if coeff == 'n11 Q11'
coeff_n = 'n_{11}';
coeff_q='Q_{11}';
else if coeff == 'ned Qed'
coeff_n = 'n_{\text{ED}}';
coeff_q='Q_{\text{ED}}';
else if coeff == 'End Qnd'
coeff_n = 'E_{\text{nD}}';
coeff_q='Q_{\text{nD}}';
else
coeff_n = '\psi';
coeff_q='\varphi';
end
end
end
%% find where the name of the sensor is in myfile :
add = find(strcmp(MyData(1).chanNames{1,2},'Sensor Name'));
channel= (MyData(1).data(add+start+nbsensor+2)); % accces each by channel{1}{i}
%% replace _ by - :
Index_under = strfind(channel{:},'_');
Index2=find(~cellfun(@isempty,Index_under));
under=cell2mat(Index_under(~cellfun('isempty',Index_under)));
for i=1:length(Index2)
channel{1,1}{Index2(i)}(under(i))='-'
end
save([folder_save 'read']) % for debugging purposes
end
c4science · Help