read me.txt

用Matlab编写的功率谱估计

实验说明：
1、此程序主要用于输入数据的功率谱估计，该程序由三部分组成：burg_Kp.m，Mainfun.m，ZXG.m.其中：Mainfun.m为主函数，burg_Kp.m为Burg算法函数，ZXG.m为求输入信号自相关的函数(估计)。
2、程序中求FFT的点数可通过宏定义的H来改变，此次选取1024点。
3、程序运行的结果为：
   a.改进后的周期图法：第一点峰值频率点：89
   b.AR模型用Levinson-Durbin法，70阶时：第一点峰值频率点为：88
   c.AR模型用burg法,30阶时：第一点峰值频率点为：88
4、由于数据采样频率为Fs=1/1年，故可得到太阳黑子的活动频率设为f,周期设为   T则：
    89/1024=f/Fs;
    解得T=(1024/89)年=11.5056年
    88/1024=f/Fs;
    解得T=(1024/88)年=11.6364年
5、结论：太阳黑子的活动周期为11年左右


附：实验题目：

6. 下面的数据是从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];