clear;clc
for k=1:1:15 %% 此处的15记得修改,15个开边界点
filename=sprintf("data_%d.mat",k);
load(filename);
Tidal(:,k)=TimeSeries(1,:);
end
%%输出
fid=fopen("Elevation.dat","wt");
[p,q] = size(Tidal)
for i=1:1:p
for j=1:1:q
if j==q
fprintf(fid,"%10.6f\n",Tidal(i,j));
else
fprintf(fid,"%10.6f\t",Tidal(i,j));
end
end
end
fclose(fid);