生成潮位矩阵数据

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);