arma谱估计

代码说明：

输入信号为方波信号+高斯白噪声，方波信号基 频为 1kHz，幅值为 1，高斯 白噪声方差为 0.2， 均值为 0，采样频率为 20kHz。试用周期图法和 ARMA 谱估 计法(SVD-TLS)估计此信号的功率 谱，并分析比较其结果  3. ARMA 模型 （1） 平稳 ARMA 功率谱密度：   （2） 任意 ARMA：   可写为：   其功率谱：   4. Cadzow 谱估计子： （1） Cadzow 将 ARMA 功率谱分 解成两部分   （2） ni 的计算方法   5. p 阶数确定： （1） 奇异值分解：   （2） 阶数界定，v(k)小于某值（例如，0.99995）：   6. A(i)的计算： （1） Sp 矩阵：   （2） A 的计算：