理解prony 算法

刚开始阅读一些文章，在介绍prony方法时，一上来就跳到构造线性差分方程，由于好久没做动态的东西，很是想不通怎么求解。 之后翻看了一些以前的基本知识，才逐渐明白。

1. 首先prony 方法的目的是：采用 p 个具有任意幅值、相位、频率和衰减因子的指数函数的线性组合对振荡数据进行等间距采样，提取并分析振荡特征，进行模式识别 

                                                                            

其中

                                            

2. 然后是求解方法：

   这是基于一个很强的假设，信号的是线性微分模型产生的。由此导出假设，解的形式是线性差分方程：

                                                                            

  由此表达式，可以构建p维线性方程组，通过求解该方程就可以获得系数a[i]的值。

                                                        

a[i] 求出来后，就可以求解z, 即可获得alpha and f. 

通过求解如下方程，又可获得幅值和初始相位等信息。至此，p个振荡指数函数的参数计算完成，且参考文献[2]提供了matlab代码。

                                              

参考文献：

[1].竺炜, 唐颖杰, 周有庆, et al. 基于改进Prony算法的电力系统低频振荡模式识别[J]. 电网技术, 2009, 33(5):44-47.

[2].Rodríguez A F, de Santiago Rodrigo L, Guillén E L, et al. Coding Prony’s method in MATLAB and applying it to biomedical signal filtering[J]. BMC bioinformatics, 2018, 19(1): 451.   https://bmcbioinformatics.biomedcentral.com/track/pdf/10.1186/s12859-018-2473-y

[3]. https://zh.wikipedia.org/wiki/%E7%BA%BF%E6%80%A7%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B

[4]. http://course.shufe.edu.cn/ode/skja_ch8.pdf 

[5].https://blog.csdn.net/Einstellung/article/details/77199054