periodic.m

arma模型的功率谱估计代码。还包括系统的极零点分析等完整软件包。

function [w, P] = Periodic(x, N0)
% 周期图法估计x(n)的功率谱,输出N0个点
N = length(x);
w = 1:1:N0;
w = (w-1)*pi/N0;
for m=1:N0
    P(m) = 0;
    for n=1:N
        P(m) = P(m)+x(n)*exp(-j*w(m)*n);
    end
    P(m) = 1/N*(abs(P(m)))^2;
end