sunspot.m

现代信号谱估计算法的matlab代码 调用函数代码

1869    60.9  59.9  52.7  41.0 103.9 108.4  59.2  79.6  80.6  59.3  78.1  104.3
1870    77.3 114.9 157.6 160.0 176.0 135.6 132.4 153.8 136.0 146.4 147.5 130.0
1871    88.3 125.3 143.2 162.4 145.5  91.7 103.0 110.1  80.3  89.0 105.4  90.4
1872    79.5 120.1  88.4 102.1 107.6 109.9 105.5  92.9 114.6 102.6 112.0  83.9
1873    86.7 107.0  98.3  76.2  47.9  44.8  66.9  68.2  47.1  47.1  55.4  49.2
1874    60.8  64.2  46.4  32.0  44.6  38.2  67.8  61.3  28.0  34.3  28.9  29.3
1875    14.6  21.5  33.8  29.1  11.5  23.9  12.5  14.6   2.4  12.7  17.7   9.9
1876    14.3  15.0  30.6   2.3   5.1   1.6  15.2   8.8   9.9  14.3   9.9   8.2
1877    24.4   8.7  11.9  15.8  21.6  14.2   6.0   6.3  16.9   6.7  14.2   2.2
1878     3.3   6.6   7.8   0.1   5.9   6.4   0.1   0.0   5.3   1.1   4.1   0.5
1879     1.0   0.6   0.0   6.2   2.4   4.8   7.5  10.7   6.1  12.3  13.1   7.3
1880    24.0  27.2  19.3  19.5  23.5  34.1  21.9  48.1  66.0  43.0  30.7  29.6
1881    36.4  53.2  51.5  51.6  43.5  60.5  76.9  58.4  53.2  64.4  54.8  47.3
1882    45.0  69.5  66.8  95.8  64.1  45.2  45.4  40.4  57.7  59.2  84.4  41.8
1883    60.6  46.9  42.8  82.1  31.5  76.3  80.6  46.0  52.6  83.8  84.5  75.9
1884    91.5  86.9  87.5  76.1  66.5  51.2  53.1  55.8  61.9  47.8  36.6  47.2
1885    42.8  71.8  49.8  55.0  73.0  83.7  66.5  50.0  39.6  38.7  30.9  21.7
1886    29.9  25.9  57.3  43.7  30.7  27.1  30.3  16.9  21.4   8.6   0.3  13.0
1887    10.3  13.2   4.2   6.9  20.0  15.7  23.3  21.4   7.4   6.6   6.9  20.7
1888    12.7   7.1   7.8   5.1   7.0   7.1   3.1   2.8   8.8   2.1  10.7   6.7
1889     0.8   8.5   6.7   4.3   2.4   6.4   9.4  20.6   6.5   2.1   0.2   6.7
1890     5.3   0.6   5.1   1.6   4.8   1.3  11.6   8.5  17.2  11.2   9.6   7.8
1891    13.5  22.2  10.4  20.5  41.1  48.3  58.8  33.0  53.8  51.5  41.9  32.5
1892    69.1  75.6  49.9  69.6  79.6  76.3  76.5 101.4  62.8  70.5  65.4  78.6
1893    75.0  73.0  65.7  88.1  84.7  89.9  88.6 129.2  77.9  80.0  75.1  93.8
1894    83.2  84.6  52.3  81.6 101.2  98.9 106.0  70.3  65.9  75.5  56.6  60.0
1895    63.3  67.2  61.0  76.9  67.5  71.5  47.8  68.9  57.7  67.9  47.2  70.7
1896    29.0  57.4  52.0  43.8  27.7  49.0  45.0  27.2  61.3  28.7  38.0  42.6
1897    40.6  29.4  29.1  31.0  20.0  11.3  27.6  21.8  48.1  14.3   8.4  33.3
1898    30.2  36.4  38.3  14.5  25.8  22.3   9.0  31.4  34.8  34.4  30.9  12.6
1899    19.5   9.2  18.1  14.2   7.7  20.5  13.5   2.9   8.4  13.0   7.8  10.5
1900     9.4  13.6   8.6  16.0  15.2  12.1   8.3   4.3   8.3  12.9   4.5   0.3
1901     0.2   2.4   4.5   0.0  10.2   5.8   0.7   1.0   0.6   3.7   3.8   0.0
1902     5.5   0.0  12.4   0.0   2.8   1.4   0.9   2.3   7.6  16.3  10.3   1.1
1903     8.3  17.0  13.5  26.1  14.6  16.3  27.9  28.8  11.1  38.9  44.5  45.6
1904    31.6  24.5  37.2  43.0  39.5  41.9  50.6  58.2  30.1  54.2  38.0  54.6
1905    54.8  85.8  56.5  39.3  48.0  49.0  73.0  58.8  55.0  78.7 107.2  55.5
1906    45.5  31.3  64.5  55.3  57.7  63.2 103.6  47.7  56.1  17.8  38.9  64.7
1907    76.4 108.2  60.7  52.6  42.9  40.4  49.7  54.3  85.0  65.4  61.5  47.3
1908    39.2  33.9  28.7  57.6  40.8  48.1  39.5  90.5  86.9  32.3  45.5  39.5
1909    56.7  46.6  66.3  32.3  36.0  22.6  35.8  23.1  38.8  58.4  55.8  54.2
1910    26.4  31.5  21.4   8.4  22.2  12.3  14.1  11.5  26.2  38.3   4.9   5.8
1911     3.4   9.0   7.8  16.5   9.0   2.2   3.5   4.0   4.0   2.6   4.2   2.2
1912     0.3   0.0   4.9   4.5   4.4   4.1   3.0   0.3   9.5   4.6   1.1   6.4
1913     2.3   2.9   0.5   0.9   0.0   0.0   1.7   0.2   1.2   3.1   0.7   3.8
1914     2.8   2.6   3.1  17.3   5.2  11.4   5.4   7.7  12.7   8.2  16.4  22.3
1915    23.0  42.3  38.8  41.3  33.0  68.8  71.6  69.6  49.5  53.5  42.5  34.5
1916    45.3  55.4  67.0  71.8  74.5  67.7  53.5  35.2  45.1  50.7  65.6  53.0
1917    74.7  71.9  94.8  74.7 114.1 114.9 119.8 154.5 129.4  72.2  96.4 129.3
1918    96.0  65.3  72.2  80.5  76.7  59.4 107.6 101.7  79.9  85.0  83.4  59.2
1919    48.1  79.5  66.5  51.8  88.1 111.2  64.7  69.0  54.7  52.8  42.0  34.9
1920    51.1  53.9  70.2  14.8  33.3  38.7  27.5  19.2  36.3  49.6  27.2  29.9
1921    31.5  28.3  26.7  32.4  22.2  33.7  41.9  22.8  17.8  18.2  17.8  20.3
1922    11.8  26.4  54.7  11.0   8.0   5.8  10.9   6.5   4.7   6.2   7.4  17.5
1923     4.5   1.5   3.3   6.1   3.2   9.1   3.5   0.5  13.2  11.6  10.0   2.8
1924     0.5   5.1   1.8  11.3  20.8  24.0  28.1  19.3  25.1  25.6  22.5  16.5
1925     5.5  23.2  18.0  31.7  42.8  47.5  38.5  37.9  60.2  69.2  58.6  98.6
1926    71.8  69.9  62.5  38.5  64.3  73.5  52.3  61.6  60.8  71.5  60.5  79.4
1927    81.6  93.0  69.6  93.5  79.1  59.1  54.9  53.8  68.4  63.1   67.2  45.2
1928    83.5  73.5  85.4  80.6  77.0  91.4  98.0  83.8  89.7  61.4  50.3  59.0
1929    68.9  62.8  50.2  52.8  58.2  71.9  70.2  65.8  34.4  54.0  81.1 108.0
1930    65.3  49.9  35.0  38.2  36.8  28.8  21.9  24.9  32.1  34.4  35.6  25.8
1931    14.6  43.1  30.0  31.2  24.6  15.3  17.4  13.0  19.0  10.0  18.7  17.8
1932    12.1  10.6  11.2  11.2  17.9  22.2   9.6   6.8   4.0   8.9   8.2  11.0
1933    12.3  22.2  10.1   2.9   3.2   5.2   2.8   0.2   5.1   3.0   0.6   0.3
1934     3.4   7.8   4.3  11.3  19.7   6.7   9.3   8.3   4.0   5.7   8.7  15.4
1935    18.6  20.5  23.1  12.2  27.3  45.7  33.9  30.1  42.1  53.2  64.2  61.5
1936    62.8  74.3  77.1  74.9  54.6  70.0  52.3  87.0  76.0  89.0 115.4 123.4
1937   132.5 128.5  83.9 109.3 116.7 130.3 145.1 137.7 100.7 124.9  74.4  88.8
1938    98.4 119.2  86.5 101.0 127.4  97.5 165.3 115.7  89.6  99.1 122.2  92.7
1939    80.3  77.4  64.6 109.1 118.3 101.0  97.6 105.8 112.6  88.1  68.1  42.1
1940    50.5  59.4  83.3  60.7  54.4  83.9  67.5 105.5  66.5  55.0  58.4  68.3
1941    45.6  44.5  46.4  32.8  29.5  59.8  66.9  60.0  65.9  46.3  38.4  33.7
1942    35.6  52.8  54.2  60.7  25.0  11.4  17.7  20.2  17.2  19.2  30.7  22.5
1943    12.4  28.9  27.4  26.1  14.1   7.6  13.2  19.4  10.0   7.8  10.2  18.8
1944     3.7   0.5  11.0   0.3   2.5   5.0   5.0  16.7  14.3  16.9  10.8  28.4
1945    18.5  12.7  21.5  32.0  30.6  36.2  42.6  25.9  34.9  68.8  46.0  27.4
1946    47.6  86.2  76.6  75.7  84.9  73.5 116.2 107.2  94.4 102.3 123.8 121.7
1947   115.7 133.4 129.8 149.8 201.3 163.9 157.9 188.8 169.4 163.6 128.0 116.5
1948   108.5  86.1  94.8 189.7 174.0 167.8 142.2 157.9 143.3 136.3  95.8 138.0
1949   119.1 182.3 157.5 147.0 106.2 121.7 125.8 123.8 145.3 131.6 143.5 117.6
1950   101.6  94.8 109.7 113.4 106.2  83.6  91.0  85.2  51.3  61.4  54.8  54.1
1951    59.9  59.9  55.9  92.9 108.5 100.6  61.5  61.0  83.1  51.6  52.4  45.8
1952    40.7  22.7  22.0  29.1  23.4  36.4  39.3  54.9  28.2  23.8  22.1  34.3
1953    26.5   3.9  10.0  27.8  12.5  21.8   8.6  23.5  19.3   8.2   1.6   2.5
1954     0.2   0.5  10.9   1.8   0.8   0.2   4.8   8.4   1.5   7.0   9.2   7.6
1955    23.1  20.8   4.9  11.3  28.9  31.7  26.7  40.7  42.7  58.5  89.2  76.9
1956    73.6 124.0 118.4 110.7 136.6 116.6 129.1 169.6 173.2 155.3 201.3 192.1
1957   165.0 130.2 157.4 175.2 164.6 200.7 187.2 158.0 235.8 253.8 210.9 239.4
1958   202.5 164.9 190.7 196.0 175.3 171.5 191.4 200.2 201.2 181.5 152.3 187.6
1959   217.4 143.1 185.7 163.3 172.0 168.7 149.6 199.6 145.2 111.4 124.0 125.0
1960   146.3 106.0 102.2 122.0 119.6 110.2 121.7 134.1 127.2  82.8  89.6  85.6
1961    57.9  46.1  53.0  61.4  51.0  77.4  70.2  55.8  63.6  37.7  32.6  39.9
1962    38.7  50.3  45.6  46.4  43.7  42.0  21.8  21.8  51.3  39.5  26.9  23.2
1963    19.8  24.4  17.1  29.3  43.0  35.9  19.6  33.2  38.8  35.3  23.4  14.9
1964    15.3  17.7  16.5   8.6   9.5   9.1   3.1   9.3   4.7   6.1   7.4  15.1
1965    17.5  14.2  11.7   6.8  24.1  15.9  11.9   8.9  16.8  20.1  15.8  17.0
1966    28.2  24.4  25.3  48.7  45.3  47.7  56.7  51.2  50.2  57.2  57.2  70.4
1967   110.9  93.6 111.8  69.5  86.5  67.3  91.5 107.2  76.8  88.2  94.3 126.4
1968   121.8 111.9  92.2  81.2 127.2 110.3  96.1 109.3 117.2 107.7  86.0 109.8
1969   104.4 120.5 135.8 106.8 120.0 106.0  96.8  98.0  91.3  95.7  93.5  97.9
1970   111.5 127.8 102.9 109.5 127.5 106.8 112.5  93.0  99.5  86.6  95.2  83.5
1971    91.3  79.0  60.7  71.8  57.5  49.8  81.0  61.4  50.2  51.7  63.2  82.2
1972    61.5  88.4  80.1  63.2  80.5  88.0  76.5  76.8  64.0  61.3  41.6  45.3
1973    43.4  42.9  46.0  57.7  42.4  39.5  23.1  25.6  59.3  30.7  23.9  23.3
1974    27.6  26.0  21.3  40.3  39.5  36.0  55.8  33.6  40.2  47.1  25.0  20.5
1975    18.9  11.5  11.5   5.1   9.0  11.4  28.2  39.7  13.9   9.1  19.4   7.8
1976     8.1   4.3  21.9  18.8  12.4  12.2   1.9  16.4  13.5  20.6   5.2  15.3
1977    16.4  23.1   8.7  12.9  18.6  38.5  21.4  30.1  44.0  43.8  29.1  43.2
1978    51.9  93.6  76.5  99.7  82.7  95.1  70.4  58.1 138.2 125.1  97.9 122.7
1979   166.6 137.5 138.0 101.5 134.4 149.5 159.4 142.2 188.4 186.2 183.3 176.3
1980   159.6 155.0 126.2 164.1 179.9 157.3 136.3 135.4 155.0 164.7 147.9 174.4
1981   114.0 141.3 135.5 156.4 127.5  90.9 143.8 158.7 167.3 162.4 137.5 150.1
1982   111.2 163.6 153.8 122.0  82.2 110.4 106.1 107.6 118.8  94.7  98.1 127.0
1983    84.3  51.0  66.5  80.7  99.2  91.1  82.2  71.8  50.3  55.8  33.3  33.4
1984    57.0  85.4  83.5  69.7  76.4  46.1  37.4  25.5  15.7  12.0  22.8  18.7
1985    16.5  15.9  17.2  16.2  27.5  24.2  30.7  11.1   3.9  18.6  16.2  17.3
1986     2.5  23.2  15.1  18.5  13.7   1.1  18.1   7.4   3.8  35.4  15.2   6.8
1987    10.4   2.4  14.7  39.6  33.0  17.4  33.0  38.7  33.9  60.6  39.9  27.1
1988    59.0  40.0  76.2  88.0  60.1 101.8 113.8 111.6 120.1 125.1 125.1 179.2 
1989   161.3 165.1 131.4 130.6 138.5 196.2 126.9 168.9 176.7 159.4 173.0 165.5];
%
%1990   177.3 130.5 140.3 139.8 132.0 105.2 147.0

% we generate sunspot numbers as the annual sunspot activity, found by summing
% the monthly numbers
year=x(:,1);
sspot=x(:,2:13);
sunspot=sum(sspot')';