谐波恢复的ARMA建模

利用MATLAB软件，在原有的ARMA模型的基础上进行改进，实现谐波恢复的ARMA建模

%%利用ARMA、AR、MA模型，以及周期图等进行系统参数估计
N=456;
B1=[1 0.3544 0.3508 0.1736 0.2401]; 
A1=[1 -1.3817 1.5632 -0.8843 0.4096];
w=linspace（0pi512）;
H1=freqz（B1A1w）;%产生信号的频域响应
Ps1=abs（H1）.^2;
SPy11=0;%20次AR（4）
SPy12=0;%20次AR（8）
SPy13=0;%20次AR周期图
SPy14=0;%20次ARMA（44）
SPy15=0;%20次ARMA（88）
VSPy11=0;%20次AR（4）
VSPy12=0;%20次AR（8）
VSPy13=0;%20次AR周期图
VSPy14=0;%20次ARMA（44）
VSPy15=0;%20次ARMA（88）
for k=1:20
%采用自协方差法对AR模型参数进行估计%
%gA1：AR模型的参数；gE1：激励白噪声的方差%
y1=filter（B1A1randn（1N））.*[zeros（1200）ones（1256）];    %滤波（窗函数）
[Py11F]=pcov（y145121）;%AR（4）的估计%
[Py12F]=pcov（y185121）;%AR（8）的估计%
[Py13F]=periodogram（y1[]5121）;     %求功率谱
SPy11=SPy11+Py11;
SPy12=SPy12+Py12;
SPy13=SPy13+Py13;
VSPy11=VSPy11+abs（Py11）.^2;
VSPy12=VSPy12+abs（Py12）.^2;
VSPy13=VSPy13+abs（Py13）.^2;
fig

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       4204  2017-12-01 20:24  ARMA.m

----------- ---------  ---------- -----  ----

                 4204                    1