写一段潮汐调和分析的matlab代码

下面是一段简单的潮汐调和分析的 Matlab 代码：

% 计算潮汐高度
h = 0.5 * cos(2 * pi * t / T);

% 计算潮汐周期
T = 12.42 / M;

% 计算潮汐调和分析
for m = 2:2:M
    h = h + (1 / m) * cos(2 * pi * m * t / T);
end

% 绘制图像
plot(t, h);
xlabel('时间 (h)');
ylabel('潮汐高度 (m)');
title('潮汐调和分析');

这段代码计算了潮汐高度，潮汐周期，并通过潮汐调和分析得到了潮汐的预测值。最后，它通过绘制图像展示了潮汐的变化情况。