-
Notifications
You must be signed in to change notification settings - Fork 0
/
Y_plot_all_fft_and_t_history.m
executable file
·57 lines (47 loc) · 1.27 KB
/
Y_plot_all_fft_and_t_history.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
clear all
close all
sysdir
list0=dir('g0.4a*m*r*');
n0=length(list0);
s=0;
for j=1:n0
cd (list0(j).name)
listf1=dir('0*');
n=length(listf1);
cd ('../')
str=list0(j).name;
%str(str=='r')=[];
str(str=='u')=[];
str(str=='n')=[];
str(str=='_')=[];
str(str=='z')=[];
str(str=='r')=[' '];
%str1=str;
str(str=='g')=' ';
str(str=='a')=' ';
str(str=='m')=' ';
gam=str2num(str);
G=gam(1);
A=gam(2);
m=gam(3);
for i=1:n % i=1:n
dirnamef=[list0(j).name,'/',listf1(i).name,'/'];
%v01=load([dirname,'VIV01.DAT']);
cd(dirnamef)
disp(pwd);
%[y,t]=loadv02yt(0);
fdir=pwd;
cd ../; cd ../;
mkdir(['plotting/',list0(j).name])
cd (['plotting/',list0(j).name])
ttl=['G=',num2str(G),' A_1=',num2str(A),' m=',num2str(m),' f1=',listf1(i).name,' Cylinder 2'];
pdir=pwd;
pname=[list0(j).name,'_',listf1(i).name];
f1=str2double(listf1(i).name);
ylimfft=[0 0.0005]; disp(['ylimfft=',num2str(ylimfft)])
plotfft(fdir,[],24000,pdir,[],pname,ylimfft,f1,4,[])
cd ('../')
cd ('../')
close all
end
end