matlab功率谱补零,matlab如何在功率谱密度图中添加一条信度值为0.9的置信曲线

下面的数据是从1770年至1869年这100年间每年12个月所记录到的太阳黑子出现次数的平均值。

a. 用周期图法做该数据的功率谱，从该曲线您是否可以看出太阳黑子活动的周期？

b. 对该数据建立一个三阶的AR模型，分别用Levinson-Durbin法和Burg法求解AR模型的参数，并得到功率谱曲线，从该曲线来发现太阳黑子的活动周期。改变AR模型的阶次，看不同阶次下所估计出的周期的差别。

x = [101 82 66 35 31 7 20 92 154 125 85 68 38 23 10 24 83 132 131 118 90 67 60 47 41 21 16 6 4 7 14 34 45 43 48 42 28 10 8 2 0 1 5 12 14 35 46 41 30 24 16 7 4 2 8 17 36 50 62 67 71 48 28 8 13 57 122 138 103 86 63 37 24 11 15 40 62 98 124 96 66 64 54 39 21 7 4 23 55 94 96 77 59 44 47 30 16 7 37 74];

下面是matlab程序：

clear;

clc;

x = [101 82 66 35 31 7 20 92 154 125 85 68 38 23 10 24 83 132 131 118 90 67 60 47 41 21 16 6 4 7 14 34 45 43 48 42 28 10 8 2 0 1 5 12 14 35 46 41 30 24 16 7 4 2 8 17 36 50 62 67 71 48 28 8 13 57 122 138 103 86 63 37 24 11 15 40 62 98 124 96 66 64 54 39 21 7 4 23 55 94 96 77 59 44 47 30 16 7 37 74];

N=length(x)