现代法谱估计(4)Music算法MATLAB及Python实现
原理原博客是用word编辑成的,插入了很多Mathtype公式,而CSDN不支持此格式,转换格式后出现了很多麻烦和错误,故此博客直接放了截图,请见谅。程序和结果MATLAB程序:结果:Python程序:结果:分析由上图可见,我给程序输入的N为256,信号中f1=0.1,f2=0.13,为了看得明显,我把f1幅值调大了,在图中我们可以看到对应的位置出现了峰值,寻峰得到频率估计均为0.1和0.3左右。程序中我取的 的范围是0到π弧度,所以也只是画出了0.1到0.5