%%利用ARMA、AR、MA模型，以及周期圖等進(jìn)行系統(tǒng)參數(shù)估計(jì)
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）;%產(chǎn)生信號(hào)的頻域響應(yīng)
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
%采用自協(xié)方差法對(duì)AR模型參數(shù)進(jìn)行估計(jì)%
%gA1：AR模型的參數(shù)；gE1：激勵(lì)白噪聲的方差%
y1=filter（B1A1randn（1N））.*[zeros（1200）ones（1256）];????%濾波（窗函數(shù)）
[Py11F]=pcov（y145121）;%AR（4）的估計(jì)%
[Py12F]=pcov（y185121）;%AR（8）的估計(jì)%
[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

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

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

-----------?---------??----------?-----??----

?????????????????4204????????????????????1