% 溫度分別為 221.0645， 233.7419，250.7742，229.6129 ，（間隔為6小時）
% （1）采用線性內插、多項式內插和樣條內插方法插值到壹小時間隔，並比較不同；
% （2）對上述壹小時內插結果進行調和分析；
% （3）對上述海表溫度壹小時內插結果進行Fourier分析。
clc; clear all; close all;
t0 = [2 8 14 20];
t1 = min(t0) : 1 : max(t0);
v0 = [221.0645 233.7419 250.7742 229.6129];
p1 = polyfit(t0, v0, 1); % 線性
v11 = polyval(p1, t1);
p2 = polyfit(t0, v0, 3); % 3次多項式
v12 = polyval(p2, t1);
v13 = spline(t0, v0, t1); '樣條'
figure; hold on; box on;
plot(t0, v0, 'k*');
plot(t1, v11, 'r-');
plot(t1, v12, 'g-');
plot(t1, v13, 'b-');
legend('節點', '線性', '3次多項式', '樣條');