param_stabilize.m

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

function [Dml, lambdamax, mm, nrFDIR]= param_stabilize(Cml, N, lambda, alg)
% function [Dml, lambdamax, mm, nrFDIR]= param_stabilize(Cml, N, lambda, alg)
%   This file is part of the TFPM toolbox v1.0 (c)
%   [e-mail] and underlies the GPL.
% 
% alg==1: shrink
% alg==2: saturate

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;tfpm;
MAR  =   3;
LAR  =   2;
MMA  =   3;
LMA  =   2;
N    = 256;
re_im= 'i';
mo_no= 'n';
tfpm_file_gen;
%-------------
alpha= 1/2;
lambda= .98;
alg= 1;
%Cml= (1.1*Aml-.1*Bml);
Cml= Bml;Cml(:, 2:end)= Cml(:, 2:end)+.8*randn(2*LAR+1, MAR);
%[Dml, lambdamax, mm, nrFDIR]= param_stabilize(Cml, N, lambda)
%figure(21);mesh(fliplr(abs(tfarma_impr(Cml, 1, N, 1/2))))

Cml= [
   0.0878 + 0.0785i   0.5858 - 0.4004i   0.0862 + 0.0260i   0.2836 - 0.0528i
   0.0044 - 0.1536i   0.1840 + 0.1776i   0.1798 + 0.0571i   0.7029 + 0.1465i
   1.7227             0.3453 - 0.5094i   0.0713 - 0.1039i  -0.2028 + 0.3373i
   0.0044 + 0.1536i   0.3803 + 0.0471i   0.9098 + 0.2566i   0.4255 - 0.1673i
   0.0878 - 0.0785i   0.1162 + 0.1160i   0.4867 - 0.1522i   0.3405 + 0.1902i
];
%Cml= param_hermite(Cml);
TDIR= param_tdir(Cml, N);
FDIR= param_fdir(Cml, N);
%figure(21);clf;plot(abs(TDIR)  );axis([1 N .5 1.5])
%figure(22);clf;plot(abs(FDIR).');axis([1 N .5 1.5])
max(max(abs(TDIR)))
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
alg= 1;%!!!!!!!!!!!!!!!!!!!!
Cml0= Cml;
[M, L]= param_dim(Cml);
if(~M)
   TDIRopt= 0;
else
   TDIRopt= 1;
end;

if(~L)
   FDIRopt= 0;
else
   FDIRopt= 1;
end;
   FDIRopt= 0;

[TDIRnm, TDIRn0, Cnm]= param_tdir(Cml, N);
[FDIRkl, FDIRkL, Ckl]= param_fdir(Cml, N);
abst= abs(TDIRnm);
absf= abs(FDIRkl);
exceedt1= abst<lambda;
exceedt2= (lambda<=abst) & (abst<=1);
exceedt3= (1<abst) & (abst<=1/lambda);
exceedt4= 1/lambda<abst;
exceedf1= absf<lambda;
exceedf2= (lambda<=absf) & (absf<=1);
exceedf3= (1<absf) & (absf<=1/lambda);
exceedf4= 1/lambda<absf;
lambdamax= max(max(abst));
mm= 0;
%while( (lambdamax>=lambda) | sum(sum(exceedf2+exceedf3)) )
while( (lambdamax>=lambda) )
   mm= mm+1;
   [ mm lambdamax ];
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(1);clf;
%subplot(3, 2, 1);cla;
plot(abs(Cnm), 'k', 'Linewidth', 2);
axis([1 N 0 3]);
grid;
set(gca, 'XTick', [1 N/4 N/2 3*N/4 N])
set(gca, 'XTickLabel', [0 N/4-1 N/2-1 3*N/4-1 N-1])
set(gca, 'XTickLabel', ['V'; 'W'; 'X'; 'Y'; 'Z'])
set(gca, 'YTick', 0:4)
set(gca, 'YTickLabel', ['A'; 'B'; 'C'; 'D'; 'E'])
%text(225, 2.6, 'a')
%ylabel('cmna')
title('a')
print -deps2 -r300 /users/mjachan/tex/prj/TF/figures/TDIR1.eps

%subplot(3, 2, 2);cla;
clf;hold on
plot(abs(TDIRnm), 'k.');
plot(ones(N, 1), 'k')
axis([1 N .4 1.2]);
grid;box
set(gca, 'XTick', [1 N/4 N/2 3*N/4 N])
set(gca, 'XTickLabel', [0 N/4-1 N/2-1 3*N/4-1 N-1])
set(gca, 'XTickLabel', ['V'; 'W'; 'X'; 'Y'; 'Z'])
set(gca, 'YTick', 0.4:.2:1.2)
set(gca, 'YTickLabel', ['55'; '66'; '77'; '88'; '99'])
%text(225, 1.1, 'b')
%ylabel('rmnb')
print -deps2 -r300 /users/mjachan/tex/prj/TF/figures/TDIR2.eps
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
sum(sum(exceedt1))
sum(sum(exceedt2))
sum(sum(exceedt3))
sum(sum(exceedt4))
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Change
if(TDIRopt)
   if(alg==1)
      TDIR= lambda   *(exceedt1.*TDIRnm) + ...
	    lambda^2 *(exceedt2.*TDIRnm) + ...
            lambda^2./conj(TDIRnm).*(exceedt3) + ...
            lambda  ./conj(TDIRnm).*(exceedt4);
   else
      TDIR= (exceedt1.*TDIRnm) + ...
	    lambda*(exceedt2.*TDIRnm./abs(TDIRnm)) + ...
            lambda*(exceedt3.*TDIRnm./abs(TDIRnm)) + ...
            1./conj(TDIRnm).*(exceedt4);
   end;
else
   TDIR= TDIRnm;
end;
   Dnm= param_invtdir(TDIR, TDIRn0);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clf;%subplot(3, 2, 3);cla;
plot(abs(Dnm), 'k', 'Linewidth', 2);
axis([1 N 0 3]);
grid;
set(gca, 'XTick', [1 N/4 N/2 3*N/4 N])
set(gca, 'XTickLabel', [0 N/4-1 N/2-1 3*N/4-1 N-1])
set(gca, 'XTickLabel', ['V'; 'W'; 'X'; 'Y'; 'Z'])
set(gca, 'YTick', 0:4)
set(gca, 'YTickLabel', ['A'; 'B'; 'C'; 'D'; 'E'])
%text(225, 2.6, 'c')
%ylabel('cmnc')
title('b')
print -deps2 -r300 /users/mjachan/tex/prj/TF/figures/TDIR3.eps

clf;hold on;%subplot(3, 2, 4);cla;hold on
plot(abs(TDIR), 'k.');
plot(ones(N, 1), 'k')
axis([1 N .4 1.2]);
grid;box
set(gca, 'XTick', [1 N/4 N/2 3*N/4 N])
set(gca, 'XTickLabel', [0 N/4-1 N/2-1 3*N/4-1 N-1])
set(gca, 'XTickLabel', ['V'; 'W'; 'X'; 'Y'; 'Z'])
set(gca, 'YTick', 0.4:.2:1.2)
set(gca, 'YTickLabel', ['55'; '66'; '77'; '88'; '99'])
%text(225, 1.1, 'd')
%ylabel('rmnd')
print -deps2 -r300 /users/mjachan/tex/prj/TF/figures/TDIR4.eps
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

   Dml= param_compress_nm(Dnm, L);
   [TDIRnm, TDIRn0, Dnm]= param_tdir(Dml, N);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clf;%subplot(3, 2, 5);cla;
plot(abs(Dnm), 'k', 'Linewidth', 2);
axis([1 N 0 3]);
grid;
set(gca, 'XTick', [1 N/4 N/2 3*N/4 N])
set(gca, 'XTickLabel', [0 N/4-1 N/2-1 3*N/4-1 N-1])
set(gca, 'XTickLabel', ['Q'; 'R'; 'S'; 'T'; 'U'])
set(gca, 'YTick', 0:4)
set(gca, 'YTickLabel', ['A'; 'B'; 'C'; 'D'; 'E'])
%text(225, 2.6, 'e')
%ylabel('cmne')
title('c')
print -deps2 -r300 /users/mjachan/tex/prj/TF/figures/TDIR5.eps

clf;hold on;%subplot(3, 2, 6);cla;hold on
plot(abs(TDIRnm), 'k.');
plot(ones(N, 1), 'k')
axis([1 N .4 1.2]);
grid;box
set(gca, 'XTick', [1 N/4 N/2 3*N/4 N])
set(gca, 'XTickLabel', [0 N/4-1 N/2-1 3*N/4-1 N-1])
set(gca, 'XTickLabel', ['Q'; 'R'; 'S'; 'T'; 'U'])
set(gca, 'YTick', 0.4:.2:1.2)
set(gca, 'YTickLabel', ['55'; '66'; '77'; '88'; '99'])
%text(225, 1.1, 'f')
%ylabel('rmnf')
print -deps2 -r300 /users/mjachan/tex/prj/TF/figures/TDIR6.eps

[TDIRnm, TDIRn0, Cnm]= param_tdir(Dml, N);
[FDIRkl, FDIRkL, Ckl]= param_fdir(Dml, N);
abst= abs(TDIRnm);
absf= abs(FDIRkl);
exceedf1= absf<lambda;
exceedf2= (lambda<=absf) & (absf<=1);
exceedf3= (1<absf) & (absf<=1/lambda);
exceedf4= 1/lambda<absf;
max(max(abst))
sum(sum(exceedf2+exceedf3))

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%













   [FDIRkl, FDIRkL, Dkl]= param_fdir(Dml, N);
   absf= abs(FDIRkl);
   exceedf1= absf<lambda;
   exceedf2= (lambda<=absf) & (absf<=1);
   exceedf3= (1<absf) & (absf<=1/lambda);
   exceedf4= 1/lambda<absf;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
figure(mm);clf;
subplot(3, 2, 1);cla;
plot(abs(Dkl).', 'k.');
axis([1 N 0 800]);
grid;
set(gca, 'XTick', [1 N/4 N/2 3*N/4 N])
set(gca, 'XTickLabel', [0 N/4-1 N/2-1 3*N/4-1 N-1])
set(gca, 'YTick', 0:200:800)
set(gca, 'YTickLabel', ['A'; 'B'; 'C'; 'D'; 'E'])
text(225, 700, 'a')
ylabel('ckla')

subplot(3, 2, 2);cla;
plot(abs(FDIRkl).', 'k.');
axis([1 N .8 1.2]);
grid;
set(gca, 'XTick', [1 N/4 N/2 3*N/4 N])
set(gca, 'XTickLabel', [0 N/4-1 N/2-1 3*N/4-1 N-1])
set(gca, 'YTick', 0.8:.1:1.2)
set(gca, 'YTickLabel', ['55'; '66'; '77'; '88'; '99'])
text(225, 1.15, 'b')
ylabel('sklb')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
sum(sum(exceedf1))
sum(sum(exceedf2))
sum(sum(exceedf3))
sum(sum(exceedf4))
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Change
if(FDIRopt)
   FDIR= lambda     .*(exceedf1.*FDIRkl)              + ...
	 lambda^2   .*(exceedf2.*FDIRkl./abs(FDIRkl)) + ...
	 lambda^(-2).*(exceedf3.*FDIRkl./abs(FDIRkl)) + ...
         lambda^(-1).*(exceedf4.*FDIRkl);
else
   FDIR= FDIRkl;
end;
   Ekl= param_invfdir(FDIR, FDIRkL);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(3, 2, 3);cla;
plot(abs(Ekl).', 'k.');
axis([1 N 0 800]);
grid;
set(gca, 'XTick', [1 N/4 N/2 3*N/4 N])
set(gca, 'XTickLabel', [0 N/4-1 N/2-1 3*N/4-1 N-1])
set(gca, 'YTick', 0:200:800)
set(gca, 'YTickLabel', ['A'; 'B'; 'C'; 'D'; 'E'])
text(225, 700, 'c')
ylabel('cklc')

subplot(3, 2, 4);cla;
plot(abs(FDIR).', 'k.');
axis([1 N .8 1.2]);
grid;
set(gca, 'XTick', [1 N/4 N/2 3*N/4 N])
set(gca, 'XTickLabel', [0 N/4-1 N/2-1 3*N/4-1 N-1])
set(gca, 'YTick', 0.8:.1:1.2)
set(gca, 'YTickLabel', ['55'; '66'; '77'; '88'; '99'])
text(225, 1.15, 'd')
ylabel('skld')
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

   Eml= param_compress_kl(Ekl, M);
   Eml(:, 1)= (Eml(:, 1)+flipud(conj(Eml(:, 1))))/2;
   [FDIRkl, FDIRkL, Ekl]= param_fdir(Eml, N);
   [TDIRnm, TDIRn0, Enm]= param_tdir(Eml, N);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(3, 2, 5);cla;
plot(abs(Ekl).', 'k.');
axis([1 N 0 800]);
grid;
set(gca, 'XTick', [1 N/4 N/2 3*N/4 N])
set(gca, 'XTickLabel', [0 N/4-1 N/2-1 3*N/4-1 N-1])
set(gca, 'YTick', 0:200:800)
set(gca, 'YTickLabel', ['A'; 'B'; 'C'; 'D'; 'E'])
text(225, 700, 'e')
ylabel('ckle')

subplot(3, 2, 6);cla;
plot(abs(FDIRkl).', 'k.');
axis([1 N .8 1.2]);
grid;
set(gca, 'XTick', [1 N/4 N/2 3*N/4 N])
set(gca, 'XTickLabel', [0 N/4-1 N/2-1 3*N/4-1 N-1])
set(gca, 'YTick', 0.8:.1:1.2)
set(gca, 'YTickLabel', ['55'; '66'; '77'; '88'; '99'])
text(225, 1.15, 'f')
ylabel('sklf')
%print -deps2 -tiff -r300 /users/mjachan/tex/prj/TF/figures/FDIR.eps

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

   abst= abs(TDIRnm);
   absf= abs(FDIRkl);
   exceedt1= abst<lambda;
   exceedt2= (lambda<=abst) & (abst<=1);
   exceedt3= (1<abst) & (abst<=1/lambda);
   exceedt4= 1/lambda<abst;
   exceedf1= absf<lambda;
   exceedf2= (lambda<=absf) & (absf<=1);
   exceedf3= (1<absf) & (absf<=1/lambda);
   exceedf4= 1/lambda<absf;
   lambdamax= max(max(abst));

   Cml= Eml;Cnm= Enm;
end
Dml= Cml;

[FDIRkl, FDIRkL, Ckl]= param_fdir(Dml, N);
abst= abs(TDIRnm);
absf= abs(FDIRkl);
exceedf1= absf<lambda;
exceedf2= (lambda<=absf) & (absf<=1);
exceedf3= (1<absf) & (absf<=1/lambda);
exceedf4= 1/lambda<absf;
nrFDIR= sum(sum(exceedf2+exceedf3));

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
mm
TDIR= param_tdir(Dml, N);
FDIR= param_fdir(Dml, N);
figure(23);clf;plot(abs(TDIR)  );axis([1 N .5 1.5])
figure(24);clf;plot(abs(FDIR).');axis([1 N .5 1.5])
figure(25);mesh(abs(rot90(tfarma_impr(Dml, 1, N, 1/2))))

TDIR0= param_tdir(Cml0, N);
FDIR0= param_fdir(Cml0, N);

[max(max(abs(TDIR0))) max(max(abs(TDIR)))]
[sum(sum((abs(FDIR0)>lambda)&(abs(FDIR0)<1/lambda)))  sum(sum((abs(FDIR)>lambda)&(abs(FDIR)<1/lambda)))]
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(4, 2, 7);cla;
plot(abs(Dnm), 'k.');
axis([1 N 0 4]);
grid;
set(gca, 'XTick', [1 N/4 N/2 3*N/4 N])
set(gca, 'XTickLabel', [0 N/4-1 N/2-1 3*N/4-1 N-1])
set(gca, 'YTick', 0:4)
text(225, 3.5, 'g')


subplot(4, 2, 8);cla;
plot(abs(TDIRnm), 'k.');
axis([1 N .4 1.2]);
grid;
set(gca, 'XTick', [1 N/4 N/2 3*N/4 N])
set(gca, 'XTickLabel', [0 N/4-1 N/2-1 3*N/4-1 N-1])
set(gca, 'YTick', 0:.2:1.2)
text(225, 1.1, 'h')
%print -deps2 -tiff -r300 /users/mjachan/tex/prj/TF/figures/TDIR.eps

max(max(abs(param_tdir(Dml, N))))
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear;tfpm;
MAR  =   3;
LAR  =   2;
MMA  =   3;
LMA  =   2;
N    = 256;
re_im= 'i';
mo_no= 'n';
tfpm_file_gen;
%-------------
alpha= 1/2;
lambda= .98;
MM= 100;
for mm= 1:MM
   Cml= Bml;
   Cml(:, 2:end)= Cml(:, 2:end)+.8*randn(2*LAR+1, MAR);
   [Dml1, lambdamax1, mm1, nrFDIR1]= param_stabilize(Cml, N, lambda, 1);
   if(nrFDIR1)
      nrFDIR1
   end;
%   [Dml2, lambdamax2, mm2, nrFDIR2]= param_stabilize(Cml, N, lambda, 2);
end;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%