Kossel_studienarbeit.pdf 1.4 MB - Institut für Strömungsmechanik ...

eadspec1d.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
% lese 1D-Spektrum aus spec1d.plt (1D-Rechnung)
% bzw. Spectrum1D_in.plt (2D-Rechnung)
% TK
function [data,header] = readspec1d(pathname);
if (nargin < 1) | (isempty(pathname))
% Menu zur Wahl des Directorys
pathname = uigetdir('','Verzeichnis wählen');
if pathname == 0
data = 0; header = 0;
return
end
end
fid = fopen(strcat(pathname, '\spec1d.plt')); % 1D
D = 3;
NHeaderLines = 5;
if fid == -1
fid = fopen(strcat( pathname,'\Spectrum1D_in.plt'));
D = 2;
NHeaderLines = 4;
if fid == -1
data = 0; header = 0;
return
end
end
for i = 1 : NHeaderLines
header{i,1} = fgetl(fid);
end
datat = fscanf(fid,'%f',[D inf])';
data = datat(:,1:2);
fclose(fid);
% 2D
readsubstrif.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
% lese sea_substr.plt bzw. sea2d_substr.plt
% berechne integrale Seegangslasten
% TK
function [IF,header] = readsubstrif(pathname);
if (nargin < 1) | (isempty(pathname))
% Menu zur Wahl des Directorys
pathname = uigetdir('','Verzeichnis wählen');
if pathname == 0
data = 0; header = 0;
return
end
end
fid = fopen(strcat(pathname,'\sea_substr.plt'));
if fid == -1
fid = fopen(strcat(pathname,'\sea2d_substr.plt'));
if fid == -1
data = 0; header = 0;
return
end
end
NHeaderLines = 19;
for i = 1 : NHeaderLines
header{i,1} = fgetl(fid);
end
i = 1;
j = 1;
line = 0;
while line ~= -1
line = fgetl(fid);
if line(1) == 'Z' | line == -1
29