tfarma_impr.m

来自「用于模拟时变非平稳的ARMA过程」· M 代码 · 共 44 行

function Halpha= tfarma_impr(Aml, Bml, N, alpha)
% function Halpha= tfarma_impr(Aml, Bml, N, alpha)
%   This file is part of the TFPM toolbox v0.5 (c)
%   [e-mail] and underlies the GPL.
% 
% Computes the impulse response (underspread approximation) of a
% TFARMA model. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;tfpm;
N    =  64;
MAR  =   0;
LAR  =   0;
MMA  =   5;
LMA  =   3;
re_im= 'r';
mo_no= 'm';
tfpm_file_gen;
%-------------
%Bml= Bml= [1 1]
alpha= 1/2;

Anm= param_expand(Aml.*param_phase(MAR, LAR, N, alpha), N);
Bnm= param_expand(Bml.*param_phase(MMA, LMA, N, alpha), N);
H= ker_to_lag(tvarma_impr(Anm, Bnm), -1, alpha);
figure(97);mesh(abs(H))
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Lalpha= tfarma_weyl(Aml, Bml, N, alpha);
Halpha= nk_to_nm(Lalpha);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(98);mesh(abs(Halpha))
figure(99);mesh(abs(H-Halpha))
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%