1월부터 11월까지의 조위 자료가 담긴 파일이 11개가 있습니다. 현재 코드에서는 모든 파일을 하나로 합쳐 조화분석하였는데, 저는 월별로 조화분석을 하여 월별 분조의 진폭 변화를 plot 하고 싶습니다. 하지만 M2, S2, K1, O1 분조만 추출하여 plot 하는 방법과 월별로 plot 하는 방법이 감이 잘 잡히지 않습니다. 이 두가지 방법에 대해서 도움을 주시면 감사할 것 같습니다.
latitude = 34.7797222222222;
filename = sprintf("%d.txt", month);
f = readtable(filename, 'Delimiter', '\t', 'VariableNamingRule', 'preserve');
time = datetime(f{:, 1}, 'InputFormat', 'yyyy/MM/dd HH:mm:ss');
combined_amplitude = [combined_amplitude; amplitude];
combined_time = [combined_time; time];
start_time = combined_time(1);
[NAME, FREQ, TIDECON, ~] = t_tide(combined_amplitude, 'interval', time_interval, 'start time', start_time, 'latitude', latitude);
코드는 여기까지 작성하였으며 자료는 아래와 같이 출력됩니다. *Q1 아래로도 여러 분조들이 있습니다.
tide freq amp amp_err pha pha_err snr
SSA 0.0002282 1.7574 3.367 357.82 153.18 0.27
MSM 0.0013098 1.4068 3.492 157.39 164.77 0.16
MM 0.0015122 0.2377 2.915 42.53 219.09 0.0066
MSF 0.0028219 1.9428 3.460 74.66 130.33 0.32
MF 0.0030501 0.8132 3.360 131.62 174.29 0.059
ALP1 0.0343966 0.3861 1.115 335.87 164.88 0.12
2Q1 0.0357064 0.9802 1.606 247.44 104.05 0.37
SIG1 0.0359087 0.2855 1.068 307.22 229.19 0.071
*Q1 0.0372185 3.5933 1.573 161.17 28.45 5.2
