clear;tfpm;N    = 256;MAR  =   3;LAR  =   2;MMA  =   0;LMA  = LAR;re_im= 'r';mo_no= 'n';F    =   2;filename= sprintf('vtfar%02d%02d%02d%04d', MAR, LAR, F, N)load(filename);%-------------alpha= 1/2;beta = 1/2;AML0= AML;B0L0= B0L;E0= randn(N, F);X= real(vtfarma_gen(E0, AML, B0L, F, alpha));x= X(:, 1);y= X(:, 2);[Aml_x, B0l_x]= tfar(x, 5, 3);[Aml_y, B0l_y]= tfar(y, 5, 3);Lx= tfarma_weyl(Aml_x, B0l_x, N, alpha);Ly= tfarma_weyl(Aml_y, B0l_y, N, alpha);Ax= ambi_est(x, x, -1, alpha);%[Psi, Mask, v2]= tf_multiwin(N, 10, 10, 0, 2, 1);RxW= ml_to_nk(Ax.*conj(Psi));Rx= ml_to_nk(Ax);IDij= real(sum(1./sum(Lx).*conj(sum(Rx)))/sum(1./sum(Ly).* conj(sum(Rx))))TFIDij= real(sum(sum(1./Lx.*conj(Rx)))/sum(sum(1./Ly.*conj(Rx))))[AML_X, B0L_X]= vtfar(X, 5, 3);[Mx, Lx]= param_dim(Aml_x)[My, Ly]= param_dim(Aml_y)[MX, LX]= param_dim(AML_X)b0x= real(param_expand(B0l_x, N)).^2;b0y= real(param_expand(B0l_y, N)).^2;bji= real(param_expand(reshape(B0L_X(2, 1, :, :), 2*LX+1, 1), N)).^2;bjj= real(param_expand(reshape(B0L_X(2, 2, :, :), 2*LX+1, 1), N)).^2;GCij= log(b0y./(bji + bjj))figure(1);clf;subplot(4, 1, 1);plot(X(:, 1))subplot(4, 1, 2);plot(X(:, 2))subplot(4, 1, 3);plot(X)subplot(4, 1, 4);plot(GCij)