原创力文档- 36
- 0
- 约4.7千字
- 约 11页
- 2017-08-05 发布于重庆
- 举报
现代信号处理大作业1217
现代信号处理大作业 姚苏洋
PAGE \* MERGEFORMAT11
《现代信号处理》大作业
姚苏洋 1130349171
1.使用matlab实现LD迭代算法
Levinson-Durbin算法
功率谱估计大致可以分为经典谱估计和现代功率谱估计，经典谱估计方法存在着以下几点缺陷：
A). 数据加窗或自相关加窗，都隐含着假定在窗外未观测到的数据或自相关系数为零，该假设不切实际。
B). 要性能好往往需要较长的数据，但实际数据长度有限
C). 窗函数容易造成谱的模糊。采用AR模型的现代谱估计方法可以克服这些不足。其中LD递推算法可以在计算机上方便实现。
LD递推算法具体计算步骤如下：
Yule-Walker方程的矩阵形式(1)所示：
(1)
系数矩阵,为Hermitian矩阵，对角线上元素相同，即为Topliez矩阵。
P-1阶Yule-Walker方程为：
(2)
其中，为误差功率。
写成联立方程：
(3)
取共轭得：
(4)
变量替换，并利用得：
(5)
表示成矩阵：
(6)
求解得：
(7)
(8)
(9)
(10)
(11)
当k=1时，即一阶递推为：
求解可得：
(12)
对于时，递推为:
, , (13)
(14)
矩阵Rx已知，可得到各阶AR模型系数为：
(15) (16)
(17)
(18)
1.2实验结果
假设p=5,使用matlab求得递推结果矩阵A如图1所示。画出每次迭代矩阵A的对应值如图2所示。
图1 LD算法得到递推结果A
图2 多次迭代输出
2.令信号由三个不同频率的复正弦信号首尾相连而成。其中，。
（1）试求x(t)的WV分布，并画出三维WV分布图
（2）指出并分析其WV分布的信号项和交叉项。
2.1 原理分析
已知WV分布公式如(19)所示：
(19)
(20)
求得的wv分布为：
(21)
其中：
(22)
2.2 实验结果
使用Matlab时频分析工具箱，可以很方便地帮助我们画出WV分布以及信号项，交叉项。将信号看作三个分段信号的叠加，分别用信号1,2,3表示，其结果分别如图3至图9所示。
图3 信号1信号项 图4 信号2信号项
图5 信号3信号项 图6 x(t) WV分布
图7 信号1 信号2 交叉项 图8 信号1 信号3 交叉项
图9 信号2 信号3 交叉项
3.一非平稳信号由两个高斯信号叠加而成：
分别求出z(t)的WV分布及模糊函数，画出二者的波形图，指出并分析其信号项和交叉项。
3.1原理分析
根据WV分布的定义式，可以得到的WV分布为;
(23)
所以的WV分布的信号项为:
(24)
交叉项为：
(25)
(26)
其中,,
所以z(t)的WV分布为：
(27)
由模糊函数的定义：
(28)
可以计算得出z1(t)的模糊函数为：
(29)
模糊函数信号项为：
(30)
交叉项为：
(31)
其中，，，，
所以可得z(t)的模糊函数为：
(32)
3.2 实验结果
图10 信号1wv分布信号项 图11 信号2 wv分布交叉项
如图12至图17所示，分别为wv分布信号项，交叉项，模糊函数，信号项，模糊函数交叉项。
图12 WV分布交叉项 图14 模糊函数
图15 模糊函数交叉项
图16 信号1 信号项 图1
文档评论（0）