param_phase.m

用于模拟时变非平稳的ARMA过程

function Pml= param_phase(M, L, N, beta)
% function Pml= param_phase(M, L, N, beta)
%   This file is part of the TFPM toolbox v1.0 (c)
%   [e-mail] and underlies the GPL.
% 
% Generates a parameter phase matrix Pml according to the TF shift
% parameter beta. If(beta==1/2) then the ones-matrix is obtained. 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;tfpm;
MAR  =   7;
LAR  =   5;
MMA  =   5;
LMA  =   3;
N    =  64;
re_im= 'i';
mo_no= 'i';
tfpm_file_gen;
%-------------
M= MAR;
L= LAR;
beta =-1/2;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

mm= 0:M;
ll= -L:L;
Pml= exp(j*2*pi*ll'*mm*(1/2-beta)/N);