[现代信号处理.docxVIP

一、有两个ARMA过程，其中信号1是宽带信号，信号2是窄带信号，分别用AR谱估计算法、ARMA谱估计算法和周期图算法估计其功率谱。产生信号1的系统函数为激励白噪声的方差为1。产生信号2的系统函数为激励白噪声的方差为1。每次实验使用的数据长度为256。对信号1，分别使用AR(4)，AR(8)，ARMA(4,4)和ARMA(8,8)模型进行谱估计，对AR方法采用自协方差算法，对ARMA算法采用改进的Yule-Walker方程算法，也用周期图法作谱估计。做20次独立实验，将20次实验结果画在一张图上，观察谱估计的随机分布性质，另将20次的平均值和真实谱画在一张图上进行比较。对信号2，分别使用AR(4)，AR(8)，AR(12)，AR(16)，ARMA(4,2)，ARMA(8,4)和ARMA(12,6)模型进行谱估计，对AR方法采用自协方差算法，对ARMA算法采用改进的Yule-Walker方程算法，也用周期图法作谱估计。做20次独立实验，将20次实验结果画在一张图上，观察谱估计的随机分布性质，另将20次的平均值和真实谱画在一张图上进行比较。对各种算法的性能进行比较分析。解:(1 )很多随机过程可以由或近似由均值为零、方差为 的白噪声序列u(n)经过具有有理传输函数H(z)的ARMA线性系统来得到。称该随机过程为ARMA过程。任意平稳ARMA过程，其功率谱密度有如下形式： (1)则x(n)可用ARMA(p,q)模型描述，即 (2)则可以根据给出的信号1的系统函数来进行计算。仿真图形结果下所示：信号1的20次谱估计总图如图一所示： A B C D 图一 信号1的20次谱估计信号1的20次谱估计平均如下图二所示： a b C d 图二信号1的20次谱估计平均上图中蓝色曲线表示理论值，红色曲线表示估计值。信号1的20次周期图法图谱估计如图三所示： a b图三 信号1的20次周期图法图 (2)仿真结果如下所示： 信号2的20次谱估计总图如图四所示 a b C d e f G图四 信号2的20次谱估计信号2的20次谱估计平均如图五所示：图六信号2的20次谱估计平均信号2的周期法谱估计如图七所示：3 各个算法分析，对于信号1，有图形可以看出，AR模型，估计所得值