MA模型的建立

clc;
clear;
close all;

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;
figure（1）
plot（w./（2*pi）Ps1FPy11）;
legend（‘真实功率谱‘‘20次AR（4）估计图‘）;
hold on;
figure（2）
plot（w./（2*pi）Ps1FPy12）;
legend（‘真实功率谱‘‘20次AR（8）估计图‘）;
hold on;
figure（3）
plot（w./（2*pi）Ps1FPy13）;
legend（‘真实功率谱‘‘20次周期图法