在MATLAB中生成時間序列的面積圖

Question

我正在嘗試生成如下所示的數字：

圖http://www.cpc.ncep.noaa.gov/data/teledoc/nao.timeseries.gif

我在這里可以找到用於嘗試生成此圖的數據。

這是我提取的數據的子集：

d = [1950    1   0.56
1950    2   0.01
1950    3  -0.78
1950    4   0.65
1950    5  -0.50
1950    6   0.25
1950    7  -1.23
1950    8  -0.19
1950    9   0.39
1950   10   1.43
1950   11  -1.46
1950   12  -1.03
1951    1  -0.42
1951    2   0.35
1951    3  -1.47
1951    4  -0.38
1951    5  -0.50
1951    6  -1.35
1951    7   1.39
1951    8  -0.41
1951    9  -1.18
1951   10   2.54
1951   11  -0.54
1951   12   1.13
1952    1   0.57
1952    2  -1.38
1952    3  -1.97
1952    4   0.95
1952    5  -0.99
1952    6  -0.10
1952    7  -0.06
1952    8  -0.49
1952    9  -0.38
1952   10  -0.28
1952   11  -1.32
1952   12  -0.49
1953    1  -0.12
1953    2  -1.00
1953    3  -0.45
1953    4  -1.96
1953    5  -0.56
1953    6   1.41
1953    7   0.43
1953    8  -1.04
1953    9  -0.19
1953   10   1.95
1953   11   0.96
1953   12  -0.52
1954    1  -0.08
1954    2   0.40
1954    3  -1.27
1954    4   1.31
1954    5  -0.03
1954    6   0.06
1954    7  -0.57
1954    8  -2.57
1954    9  -0.28
1954   10   1.16
1954   11   0.29
1954   12   0.55
1955    1  -2.65
1955    2  -1.71
1955    3  -0.96
1955    4  -0.60
1955    5  -0.26
1955    6  -0.80
1955    7   1.78
1955    8   1.25
1955    9   0.46
1955   10  -1.09
1955   11  -1.49
1955   12   0.07
1956    1  -0.76
1956    2  -1.71
1956    3  -0.46
1956    4  -1.30
1956    5   2.10
1956    6   0.41
1956    7  -0.72
1956    8  -1.89
1956    9   0.38
1956   10   1.47
1956   11   0.40
1956   12   0.00
1957    1   0.71
1957    2  -0.32
1957    3  -1.73
1957    4   0.39
1957    5  -0.68
1957    6  -0.42
1957    7  -1.16
1957    8  -0.83
1957    9  -1.47
1957   10   1.95
1957   11   0.63
1957   12   0.02
1958    1  -1.14
1958    2  -1.64
1958    3  -2.46
1958    4   0.26
1958    5  -0.17
1958    6  -1.08
1958    7  -1.69
1958    8  -2.13
1958    9   0.08
1958   10   0.68
1958   11   1.59
1958   12  -0.74
1959    1  -1.52
1959    2   0.33
1959    3  -0.56
1959    4   0.25
1959    5   0.41
1959    6   0.71
1959    7   0.77
1959    8  -0.05
1959    9   1.00
1959   10   1.48
1959   11   0.30
1959   12   0.32
1960    1  -2.01
1960    2  -2.59
1960    3  -0.93
1960    4   1.33
1960    5   0.47
1960    6   0.10
1960    7   0.38
1960    8  -1.93
1960    9   0.53
1960   10  -1.37
1960   11  -0.67
1960   12  -0.03
1961    1  -0.03
1961    2   0.07
1961    3   0.17
1961    4  -1.83
1961    5  -0.28
1961    6   1.17
1961    7  -0.36
1961    8   1.03
1961    9   1.36
1961   10   1.07
1961   11  -0.79
1961   12  -1.46
1962    1   0.20
1962    2   0.18
1962    3  -2.99
1962    4   0.93
1962    5  -0.04
1962    6   0.47
1962    7  -2.43
1962    8   0.05
1962    9  -0.21
1962   10   0.96
1962   11  -0.38
1962   12  -1.31
1963    1  -2.98
1963    2  -1.53
1963    3  -0.85
1963    4  -1.61
1963    5   2.05
1963    6  -0.13
1963    7  -0.74
1963    8  -0.95
1963    9   1.89
1963   10   1.53
1963   11  -1.47
1963   12  -1.87
1964    1  -1.62
1964    2  -2.06
1964    3  -1.66
1964    4   0.25
1964    5   0.53
1964    6   1.61
1964    7   1.92
1964    8  -2.40
1964    9   0.34
1964   10   1.32
1964   11  -0.14
1964   12  -0.23
1965    1  -0.65
1965    2  -2.20
1965    3  -1.99
1965    4   0.64
1965    5  -0.52
1965    6   0.60
1965    7   0.35
1965    8   0.45
1965    9   0.51
1965   10   0.92
1965   11  -1.88
1965   12   1.18
1966    1  -2.54
1966    2  -2.02
1966    3   0.18
1966    4  -0.96
1966    5   0.25
1966    6   1.37
1966    7   0.35
1966    8  -2.39
1966    9  -0.29
1966   10  -0.23
1966   11  -0.18
1966   12   0.58];

我已經嘗試了下面的代碼來生成上面的圖：

time = datenum(d(:,1),d(:,2),repmat(15,size(d,1),1));
nao = d(:,3);

figure(1);
ax1 = subplot(211);
area(time(nao < 0),nao(nao < 0),'FaceColor',[86 180 233]/255,'EdgeColor','none');
hold on;
area(time(nao > 0),nao(nao > 0),'FaceColor','r','EdgeColor','none');
datetick('x','yyyy','keeplimits');

ax2 = subplot(212);
bar(time(nao < 0),nao(nao < 0),'FaceColor',[86 180 233]/255,'EdgeColor','none');
hold on;
bar(time(nao > 0),nao(nao > 0),'FaceColor','r','EdgeColor','none');
datetick('x','yyyy','keeplimits');

linkaxes([ax1 ax2], 'x');

然而，輸出與原始圖中所示的輸出不同，因為它們具有三個月的移動平均值。 理想情況下，我會看到一些看起來像我生成的每個圖的一半。 也就是說，區域圖，但在任何給定時間只有紅色或藍色。 區域功能似乎表明這些發生在同一時間，這是不正確的。 任何有關如何生成附圖的建議都將受到贊賞。

Answer 1

使用MATLAB FileExchange中的areashade ，以下內容應該可以滿足您的需求：

figure(2);
areashade(time, nao, 0, [86 180 233]/255, 'h')
areashade(time, nao, 0,'r','l')
datetick('x','yyyy','keeplimits');