?? tfarma_est_tfywu.m
字號(hào):
function AA= tfarma_est_tfywu(Ayy, MMA, N)% function AA= tfarma_est_tfywu(Ayy, MMA, N)% This file is part of the TFPM toolbox v1.0 (c)% michael.jachan@tuwien.ac.at and underlies the GPL.% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%clear;tfpm;N = 128;MAR = 2;LAR = 1;MMA = 1;LMA = LAR;re_im= 'r';mo_no= 'n';tfpm_file_gen;%-------------alpha= 1/2;beta= 1/2;Aalpha= tfarma_ambi(Aml, Bml, N, alpha);tf_show(Aalpha)Ayy= Aalpha(N/2-3*LAR+1:N/2+3*LAR+1, N/2+MMA-MAR+1:N/2+MMA+MAR+1);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% DimensionsMARmax= (size(Ayy, 2)-1)/2;LAR = (size(Ayy, 1)-1)/6;LL = 2*LAR+1;% The outputAA= zeros(LL, MARmax+1, MARmax);BB= zeros(LL, 1, MARmax);% Fill the Toeplitz/block-Toeplitz generator rPM and the RHS BMrPM= [];for mu= MARmax-1:-1:-MARmax+1%This is the right one;) rPM= [rPM; fliplr(Ayy(3*LAR+2-LL:3*LAR+LL, MARmax+mu+1).')];end;BM = [];for m= 1:MARmax BM= [BM; -Ayy(2*LAR+1:4*LAR+1, MARmax+m+1)];end;% Solve the systemAM= invert_waxkailath(rPM, BM);% Reshape AM to A and Compute Bfor k= 1:MARmax Aml_e= [[zeros(LAR, 1); 1; zeros(LAR, 1)] reshape(AM(1:LL*k, 1, k), LL, k)]; AA(:, :, k)= [Aml_e zeros(LL, MARmax-k)];% B0l= tfar_est_b0l(Ayy(:, MARmax+1-k:MARmax+1+k), Aml_e, N);% BB(:, :, k)= B0l;end;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%if(0)% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%param_norm(AA(:, :, MAR), Aml, 0)%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%end;% TEST TEST TEST TEST TEST TEST TEST TEST TEST TEST%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
?? 快捷鍵說(shuō)明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號(hào)
Ctrl + =
減小字號(hào)
Ctrl + -